Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,626 G→C 4.4% intergenic (+98/+96) purD → / ← sdaA phosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,566,6260G→C 4.4% 9.4 ‑0.4 ‑0.1158intergenic (+98/+96)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (4/3):  ref base (99/52):  total (103/55)
Fisher's exact test strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.61e-01

Only 100 of 170 total aligned reads displayed.

ACGAAAC.GGGGGCGCTGCCCCC.CTCCTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTCAGC  >  NC_014034/2566530‑2566726
                                                                                                    |                                                                                                    
 cGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                                                                                                    <  1:3404108/101‑1 (MQ=21)
  gAAACGGGGGGCGCTGCTCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                                                                                   <  2:3806296/101‑1 (MQ=21)
  gAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                                                                                   <  2:1593263/101‑1 (MQ=21)
    aacggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata                                                                                                 <  1:3581109/97‑1 (MQ=25)
      cggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc                                                                                               <  2:1120311/99‑1 (MQ=25)
      cagggggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggg                                                                                         >  1:3951563/3‑101 (MQ=0)
       ggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc                                                                                              <  2:1571825/100‑1 (MQ=25)
       agggggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg                                                                                        >  1:4076720/2‑101 (MQ=0)
         ggggCGCTGCCCCCTCTTGGC..CTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggggga                                                                                      <  1:3959563/101‑3 (MQ=2)
          gggCGCTGCCCCCTCTTGGC..CTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag                                                                                     >  1:601326/1‑98 (MQ=2)
          gggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag                                                                                     >  1:1254032/1‑98 (MQ=2)
           ggCGCTGCCCCCTCTTGGC..CTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggagc                                                                                    <  1:1048067/101‑5 (MQ=2)
           ggCGCTGCCCCCTCTTGGC..CTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggagc                                                                                    <  1:2334609/101‑5 (MQ=2)
             cgcTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCCAAATATCCCggggg                                                                                        >  2:3022218/1‑101 (MQ=25)
                      ccctcttggccTTCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt                                                                                <  1:1052427/91‑1 (MQ=38)
                      ccctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt                                                                                <  1:93659/91‑1 (MQ=255)
                           ttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg                                                                           <  1:3070202/96‑1 (MQ=255)
                             ggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggc                                                                         >  1:1210820/4‑101 (MQ=255)
                              gccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggag                                                                        >  2:3850891/3‑99 (MQ=255)
                               ccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggcag                                                                       >  2:1717474/2‑99 (MQ=255)
                                   cGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGc                                                                   <  2:2144179/101‑1 (MQ=255)
                                    ggCCAATTCACCCCCCGGGATATTTTGGGCGAGATAAAGAAGGGGTTTTTCTTTTCTCTTTCCCGAAAAATCACGGGGGGGGCCGCCTCGTGGCggggggg                                                                  >  2:756706/1‑99 (MQ=16)
                                     gCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAg                                                                 <  2:1958095/101‑1 (MQ=255)
                                       cAAGTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAgcg                                                               <  1:1876199/101‑1 (MQ=255)
                                        aaTTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAgcgc                                                              >  2:1742483/1‑101 (MQ=255)
                                         aTTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc                                                             >  1:2704649/1‑101 (MQ=255)
                                          ttCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccc                                                            >  1:652542/1‑101 (MQ=255)
                                            cACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccccc                                                          <  1:3077981/101‑1 (MQ=255)
                                            cACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccccc                                                          <  1:2029001/101‑1 (MQ=255)
                                             aCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCg                                                         >  1:584338/1‑101 (MQ=255)
                                              ccccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCCGCGCCCCCGa                                                        >  2:2440648/1‑101 (MQ=255)
                                              ccccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGa                                                        <  1:3964337/101‑1 (MQ=255)
                                               cccccGAGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGaa                                                       <  2:2040485/101‑1 (MQ=255)
                                                 cccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgc                                                     <  1:1742602/101‑1 (MQ=255)
                                                  ccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcg                                                    >  2:1190061/1‑101 (MQ=255)
                                                   cGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcgc                                                   >  1:1676698/1‑101 (MQ=255)
                                                   cGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcgc                                                   >  1:946245/1‑101 (MQ=255)
                                                    gggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGcc                                                  <  1:3779274/101‑1 (MQ=255)
                                                     ggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCCGCGCCCCCGAAGCGccc                                                 >  2:2508330/1‑101 (MQ=255)
                                                     ggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGccc                                                 >  1:3557518/1‑101 (MQ=255)
                                                     ggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGAAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGccc                                                 <  2:680417/101‑1 (MQ=255)
                                                      gATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCg                                                <  1:2636550/101‑1 (MQ=255)
                                                       atatTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGa                                               >  1:3851558/1‑101 (MQ=255)
                                                       atatTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGa                                               >  2:2199887/1‑101 (MQ=255)
                                                         atTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                                             >  1:2516173/1‑101 (MQ=255)
                                                          tttCGGACAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgc                                            <  2:436788/101‑1 (MQ=255)
                                                           ttCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcg                                           <  2:2617717/101‑1 (MQ=255)
                                                           ttCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcg                                           <  2:3428342/101‑1 (MQ=255)
                                                             cGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcgcg                                         >  1:3450870/1‑101 (MQ=255)
                                                              gggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcgg                                        >  2:3547478/1‑101 (MQ=255)
                                                              gggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcgg                                        >  2:848799/1‑101 (MQ=255)
                                                               ggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggc                                       >  2:185544/1‑101 (MQ=255)
                                                                  aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg                                    >  2:2700590/1‑101 (MQ=255)
                                                                  aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg                                    >  1:117707/1‑101 (MQ=255)
                                                                     aTAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtt                                 <  2:1457617/101‑1 (MQ=255)
                                                                      tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGGttc                                >  1:3181312/1‑100 (MQ=255)
                                                                       aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc                               >  1:794897/1‑101 (MQ=255)
                                                                       aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc                               >  2:3148243/1‑101 (MQ=255)
                                                                       aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc                               >  1:1983891/1‑101 (MQ=255)
                                                                       aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc                               >  1:2248743/1‑101 (MQ=255)
                                                                        aagaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTccc                              >  1:3903167/1‑101 (MQ=255)
                                                                          gaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccg                            >  2:2955844/1‑101 (MQ=255)
                                                                          gaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccg                            <  2:1222457/101‑1 (MQ=255)
                                                                           aagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTtcccgc                           >  2:4124807/1‑101 (MQ=255)
                                                                           aagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgc                           >  1:3795794/1‑101 (MQ=255)
                                                                             ggTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgccc                         <  1:3823933/101‑1 (MQ=255)
                                                                               tCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCaa                       <  1:1717589/101‑1 (MQ=255)
                                                                                 ttttCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAagg                     >  2:48972/1‑101 (MQ=255)
                                                                                 ttttCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAagg                     <  2:1249088/101‑1 (MQ=255)
                                                                                 ttttCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAagg                     <  1:2272398/101‑1 (MQ=255)
                                                                                  tttCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggg                    >  1:3251960/1‑100 (MQ=255)
                                                                                  tttCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAagga                    >  1:2588182/1‑101 (MQ=255)
                                                                                  tttCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAagga                    >  1:3121277/1‑101 (MQ=255)
                                                                                   ttCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggag                   >  2:1206568/1‑101 (MQ=255)
                                                                                   ttCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggag                   >  2:3637307/1‑101 (MQ=255)
                                                                                    tCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggagg                  >  2:3367572/1‑101 (MQ=255)
                                                                                    tCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggagg                  >  1:338741/1‑101 (MQ=255)
                                                                                     ccTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAggg                 >  1:3230237/1‑101 (MQ=255)
                                                                                     ccTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAggg                 >  2:2052951/1‑101 (MQ=255)
                                                                                     ccTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAggg                 <  1:4224027/101‑1 (MQ=255)
                                                                                      cTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggaagat                >  2:2132525/1‑97 (MQ=255)
                                                                                      cTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgc                <  1:3297777/101‑1 (MQ=255)
                                                                                      cTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgc                <  2:2559851/101‑1 (MQ=255)
                                                                                      cTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgc                >  1:2963084/1‑101 (MQ=255)
                                                                                       tttCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgcg               >  2:1659022/1‑101 (MQ=255)
                                                                                       tttCATCTTGCCCAAAATACCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgcg               <  2:3979986/101‑1 (MQ=255)
                                                                                        ttCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgcgc              >  1:4099437/1‑101 (MQ=255)
                                                                                         tCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCa             >  1:2518967/1‑101 (MQ=255)
                                                                                         tCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCa             >  2:3624076/1‑101 (MQ=255)
                                                                                         tCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCa             <  2:3091450/101‑1 (MQ=255)
                                                                                         tCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCa             >  2:3088180/1‑101 (MQ=255)
                                                                                          cATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCACGGAGGGCGCaa            >  2:575171/1‑101 (MQ=255)
                                                                                          cATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCaa            >  1:2456653/1‑101 (MQ=255)
                                                                                          cATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCaa            >  1:1837436/1‑101 (MQ=255)
                                                                                           aTCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAAc           >  2:2391960/1‑101 (MQ=255)
                                                                                            tCTTGCCGGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAAcc          <  2:1676583/101‑1 (MQ=255)
                                                                                            tCTTGCCCGAAATACCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAAcc          <  1:2700796/101‑1 (MQ=255)
                                                                                              ttGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGc        >  2:1657159/1‑101 (MQ=255)
                                                                                               tGCACGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGcc       >  2:3979110/1‑101 (MQ=255)
                                                                                                    gAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTcagc  <  2:639166/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
ACGAAAC.GGGGGCGCTGCCCCC.CTCCTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTCAGC  >  NC_014034/2566530‑2566726

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG