Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 39.5% intergenic (‑195/+113) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9850TC39.5% 218.6 / 240.0 205intergenic (‑195/+113)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (57/67);  new base C (71/10);  total (128/77)
Fisher's exact test for biased strand distribution p-value = 6.46e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-03

Only 100 of 205 total aligned reads displayed.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAG  >  NC_000913/3705887‑3706080
                                                                                                  |                                                                                                         
 ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                         >  1:2711336/1‑100 (MQ=255)
  tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTCGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                         >  1:2530029/1‑99 (MQ=255)
    ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                      >  2:1056155/1‑100 (MQ=255)
    ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                      <  2:110610/100‑1 (MQ=255)
     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCgg                                                                                                     <  2:1273241/100‑1 (MQ=37)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                   <  2:2627334/100‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                   <  2:1920883/100‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                   <  2:353427/100‑2 (MQ=37)
        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                  >  2:1290454/1‑100 (MQ=255)
        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                  <  2:670031/100‑3 (MQ=37)
         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCAT‑‑a                                                                                               <  1:2392408/100‑1 (MQ=37)
           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑a                                                                                            >  2:2135555/1‑100 (MQ=255)
              cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑aaaa                                                                                         >  2:114740/1‑100 (MQ=255)
                gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATc                                                                                       >  1:2114859/1‑100 (MQ=255)
                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGt                                                                                     <  1:1201430/100‑1 (MQ=255)
                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGc                                                                                   <  1:1628590/100‑1 (MQ=255)
                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTc                                                                             <  1:1164818/100‑1 (MQ=255)
                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTc                                                                             <  2:1442715/100‑1 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCa                                                                            >  2:111010/1‑100 (MQ=255)
                            acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCaa                                                                           <  1:748110/100‑1 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAAt                                                                          >  2:1036276/1‑100 (MQ=255)
                                   cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTg                                                                    >  1:1714925/1‑100 (MQ=255)
                                    aacGGGTATACAGATTGTTGACAAAGTGCGCATTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGc                                                                   >  2:1960069/3‑100 (MQ=255)
                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGc                                                                   <  2:2680935/100‑1 (MQ=255)
                                     cggggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACAtt                                                                     <  1:319334/95‑1 (MQ=255)
                                     ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCa                                                                  <  1:2901766/100‑1 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAt‑t‑‑‑tt                                                      >  1:2285169/1‑100 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑g                                                  <  1:767779/100‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑g                                                  <  1:1070657/100‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑g                                                  >  2:1486724/1‑100 (MQ=255)
                                               aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑Gt                                                 <  1:1056153/100‑1 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTa                                                <  1:2253249/100‑1 (MQ=255)
                                                 aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAg                                               >  2:2234618/1‑100 (MQ=255)
                                                  ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAgg                                              <  1:2278694/100‑1 (MQ=255)
                                                     ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGcc                                            <  1:536930/99‑1 (MQ=255)
                                                        aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGAt                                        <  1:1023825/100‑1 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATa                                       >  1:2423367/1‑100 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATa                                       <  2:1553936/100‑1 (MQ=255)
                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCg                                   >  1:8102/1‑100 (MQ=255)
                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCg                                   <  2:2066772/100‑1 (MQ=255)
                                                               gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                                 >  1:1134315/1‑100 (MQ=255)
                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                <  1:2875506/100‑1 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCCTTAACCCCATAGCCGGCCT‑‑AC‑AAAGTCGTGCAAACTCCGTACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATCAGCGGAgc                               <  1:211091/100‑1 (MQ=18)
                                                                   attttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATTGTGTAAATTCAATAAATTGAT‑CGAT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAgcgc                             >  2:796438/10‑100 (MQ=18)
                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCa                            >  2:582712/1‑100 (MQ=255)
                                                                     tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                           >  2:1362652/8‑100 (MQ=18)
                                                                     ttcgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCGT‑‑AC‑GAAATCGTGCTAATTCAATATATTGTG‑GGCT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCAt                           >  2:2311783/8‑100 (MQ=18)
                                                                      ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACCTTGCTAATTCAATATATTGCA‑GGGA‑C‑‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATc                          >  1:2351461/7‑100 (MQ=18)
                                                                      ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                         >  2:2138085/7‑100 (MQ=21)
                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                          <  1:2447633/100‑1 (MQ=255)
                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                          >  2:501784/1‑100 (MQ=255)
                                                                       tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑AAT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATc                          >  1:1920233/6‑100 (MQ=21)
                                                                       tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCa                         <  2:1706794/95‑1 (MQ=18)
                                                                        gtaggCCGGATAAGGCGTTCACTCCGCATCCGGCCT‑‑AC‑AAAAGCACGCAAATTCAATATATTGCA‑GAGA‑T‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                        >  1:373862/5‑100 (MQ=14)
                                                                         taggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                       >  2:621294/4‑100 (MQ=17)
                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                       >  2:2291973/4‑100 (MQ=21)
                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                       >  2:739538/4‑100 (MQ=25)
                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                       >  2:2519674/4‑100 (MQ=21)
                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                       >  2:1396946/4‑100 (MQ=21)
                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                       <  1:1396944/97‑1 (MQ=21)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                      >  2:2882801/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                      >  2:889574/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCGTGCAAATTCAATACATTGCATGGGC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                      >  1:2129657/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCGTGCAAATTCAATACATTGCATGGGC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                      >  2:1880716/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                      >  1:1492369/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAAACCTGCAAATTCAATATATTGCA‑GGAT‑A‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                      >  1:150656/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑‑AAACCGCGCAAATTCAATATATTGCG‑GAGA‑A‑‑‑AAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                     >  1:2602641/3‑99 (MQ=21)
                                                                          aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                      <  2:199529/100‑1 (MQ=255)
                                                                           ggctggataaggcgTTCACGCCGCATCCAGCCT‑‑AC‑AAAATTGTGCAAAATCAATGGATTGCA‑CAGC‑T‑‑‑AGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                     >  1:1699792/11‑100 (MQ=25)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑AAT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                      >  2:2212556/2‑100 (MQ=21)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                     >  2:2741719/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑AAT‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAgatc                     >  1:1402919/2‑97 (MQ=18)
                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                    >  1:2211422/1‑100 (MQ=255)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                   >  1:1882686/1‑100 (MQ=18)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑GAAATCGTGCTAATTCAATATATTGTG‑GGCT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                    >  2:715891/1‑99 (MQ=18)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                    >  2:2906945/1‑99 (MQ=21)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                    >  2:1764218/1‑99 (MQ=25)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATTGTGTAAATTCAATAAATTGAT‑CGAT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                    >  2:2513525/1‑99 (MQ=21)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                    >  1:1597626/1‑99 (MQ=21)
                                                                            gCCGGATAAGGCGTTAACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                    >  2:39523/1‑100 (MQ=18)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                   >  1:1658508/9‑98 (MQ=21)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                   >  2:2369487/9‑98 (MQ=21)
                                                                             ccggataaggcgTTAACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                   >  2:278676/9‑100 (MQ=18)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                   >  1:1448217/1‑100 (MQ=255)
                                                                             ccGGATAAGGCGTTAACGCCGCATCCGGCATAAAC‑AAAATCGTGCAAAATCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                     >  1:2065620/1‑100 (MQ=17)
                                                                              cggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttat                 >  2:185518/8‑98 (MQ=18)
                                                                              cggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtt                  >  2:2710258/8‑97 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                 >  2:300180/7‑100 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                 >  2:1379629/7‑100 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                 >  1:1188413/7‑100 (MQ=17)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                 >  1:2675926/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAGATCGTGCAAATTCAATATATTGCA‑TGAG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                 >  1:1979177/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCAT‑‑AC‑AAAAGCATGAAAATTCAATATATTGCA‑GGAG‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                 >  2:1835503/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGACCC‑‑AC‑AAAGTCTTGCAAAATCAATAAATTGCAGGAAC‑G‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                 >  2:88122/7‑96 (MQ=25)
                                                                                gataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAACCTTGCTAATTCAATATATTGCA‑GGGC‑C‑‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagcttt                >  1:1840469/6‑95 (MQ=18)
                                                                                gATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                >  2:1616832/1‑100 (MQ=255)
                                                                                  taaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGA‑T‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttagc              >  1:609151/4‑97 (MQ=17)
                                                                                  tGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttttt              >  1:1390426/1‑98 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg         <  2:2413054/100‑1 (MQ=255)
                                                                                            aCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc    <  2:2587543/100‑1 (MQ=255)
                                                                                                  |                                                                                                         
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAG  >  NC_000913/3705887‑3706080

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑