Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 39.8% intergenic (‑194/+114) 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,9840TG39.8% 247.6 / 280.7 205intergenic (‑194/+114)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/66);  new base G (71/10);  total (128/77)
Fisher's exact test for biased strand distribution p-value = 7.49e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.91e-02

Only 100 of 206 total aligned reads displayed.

TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAG  >  NC_000913/3705885‑3706080
                                                                                                   |                                                                                                          
     gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCATCGGGTATACAGATTGTTGTCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                       <  1:2903641/100‑1 (MQ=255)
      ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                      <  2:110610/100‑1 (MQ=255)
      ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                      >  2:1056155/1‑100 (MQ=255)
         gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                   <  2:353427/100‑2 (MQ=37)
          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                  <  2:2556677/100‑3 (MQ=37)
          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                  <  2:670031/100‑3 (MQ=37)
           tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCAT‑‑a                                                                                               <  1:2392408/100‑1 (MQ=37)
             gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑a                                                                                            >  2:2135555/1‑100 (MQ=255)
              cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑aa                                                                                           >  2:1580979/1‑100 (MQ=255)
                  gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATc                                                                                       >  1:2114859/1‑100 (MQ=255)
                    aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGt                                                                                     <  1:1201430/100‑1 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACACAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGc                                                                                   <  1:2682228/100‑1 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGc                                                                                   <  1:1628590/100‑1 (MQ=255)
                       cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCa                                                                                  >  2:1863234/1‑100 (MQ=255)
                       cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAAGACTTATCCGGCCT‑‑AC‑AAAATCGTGc                                                                                   <  1:266293/99‑1 (MQ=255)
                            gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTc                                                                             <  1:1164818/100‑1 (MQ=255)
                            gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTc                                                                             <  2:1442715/100‑1 (MQ=255)
                              acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCaa                                                                           <  1:748110/100‑1 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAAt                                                                          >  2:1036276/1‑100 (MQ=255)
                                      aacGGGTATACAGATTGTTGACAAAGTGCGCATTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGc                                                                   >  2:1960069/3‑100 (MQ=255)
                                      aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGc                                                                   <  2:2662823/100‑1 (MQ=255)
                                       cggggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACAtt                                                                     <  1:319334/95‑1 (MQ=255)
                                              taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAtt‑‑‑t‑t                                                      >  1:2285169/1‑100 (MQ=255)
                                              taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAtt‑‑‑t‑t                                                      >  2:1833524/1‑100 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑g                                                  >  2:1486724/1‑100 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑g                                                  <  1:1070657/100‑1 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑g                                                  <  1:767779/100‑1 (MQ=255)
                                                 aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑Gt                                                 <  1:1056153/100‑1 (MQ=255)
                                                 aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTCAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑Gt                                                 <  1:1863226/100‑1 (MQ=255)
                                                  gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTa                                                <  1:2253249/100‑1 (MQ=255)
                                                   aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAg                                               >  2:2234618/1‑100 (MQ=255)
                                                    ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAgg                                              <  1:2278694/100‑1 (MQ=255)
                                                       ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGcc                                            <  1:536930/99‑1 (MQ=255)
                                                          aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGAt                                        <  1:1023825/100‑1 (MQ=255)
                                                          aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGAt                                        >  2:1373498/1‑100 (MQ=255)
                                                           cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATa                                       >  1:2423367/1‑100 (MQ=255)
                                                           cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATa                                       <  2:1553936/100‑1 (MQ=255)
                                                              aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGc                                    <  1:1719736/100‑1 (MQ=255)
                                                               gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCg                                   <  2:2066772/100‑1 (MQ=255)
                                                               gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCg                                   >  1:8102/1‑100 (MQ=255)
                                                                 gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGa                                 >  1:1134315/1‑100 (MQ=255)
                                                                   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAgc                               >  1:1899801/1‑100 (MQ=255)
                                                                   gcTTTGTTTATGCCGGATGCGGCCTTAACCCCATAGCCGGCCT‑‑AC‑AAAGTCGTGCAAACTCCGTACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATCAGCGGAgc                               <  1:211091/100‑1 (MQ=18)
                                                                     attttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATTGTGTAAATTCAATAAATTGAT‑CGATC‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAgcgc                             >  2:796438/10‑100 (MQ=18)
                                                                     attttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAgcgc                             >  2:1839977/10‑100 (MQ=18)
                                                                      ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCa                            >  2:582712/1‑100 (MQ=255)
                                                                       ttcgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCGT‑‑AC‑GAAATCGTGCTAATTCAATATATTGTG‑GGCTC‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCAt                           >  2:2311783/8‑100 (MQ=18)
                                                                       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCAt                           >  1:287684/1‑100 (MQ=255)
                                                                        ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACCTTGCTAATTCAATATATTGCA‑GGGAC‑‑‑T‑AT‑‑GTAGGCCTGATAAGCATAGCGCATc                          >  1:2351461/7‑100 (MQ=18)
                                                                        ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAGT‑‑‑T‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                         >  2:2138085/7‑100 (MQ=21)
                                                                        gatctggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑AATC‑‑‑T‑CCGTGTAGGCCTGATAAGCGTAGCGCATc                          <  2:1402921/94‑1 (MQ=18)
                                                                        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATc                          >  2:501784/1‑100 (MQ=255)
                                                                        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATc                          <  1:2447633/100‑1 (MQ=255)
                                                                         tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑AATTCAAC‑TT‑‑GTAGGCCTGATAAGCGCAGCGCATc                          >  1:1920233/6‑100 (MQ=21)
                                                                          gtaggCCGGATAAGGCGTTCACTCCGCATCCGGCCT‑‑AC‑AAAAGCACGCAAATTCAATATATTGCA‑GAGAT‑‑‑T‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                        >  1:373862/5‑100 (MQ=14)
                                                                           taggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGTT‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                       >  2:621294/4‑100 (MQ=17)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTTA‑‑‑T‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                       >  2:739538/4‑100 (MQ=25)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTTA‑‑‑T‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                       >  2:2291973/4‑100 (MQ=21)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CCATC‑‑‑C‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                       >  2:2519674/4‑100 (MQ=21)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CCATC‑‑‑C‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                       <  1:1396944/97‑1 (MQ=21)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGAAAATTCAATATATTGCA‑GGAGC‑‑‑T‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                       >  2:1634854/4‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAAC‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                      >  2:2882801/3‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAAC‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                      >  2:889574/3‑100 (MQ=17)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCGTGCAAATTCAATACATTGCATGGGCC‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                      >  1:2129657/3‑100 (MQ=17)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGTT‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                      >  1:1492369/3‑100 (MQ=17)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAAACCTGCAAATTCAATATATTGCA‑GGATA‑‑‑C‑GT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                      >  1:150656/3‑100 (MQ=17)
                                                                            aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                      <  2:199529/100‑1 (MQ=255)
                                                                             ggctggataaggcgTTCACGCCGCATCCAGCCT‑‑AC‑AAAATTGTGCAAAATCAATGGATTGCA‑CAGCT‑‑‑A‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                     >  1:1699792/11‑100 (MQ=25)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑AATTCAAC‑TT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                      >  2:2212556/2‑100 (MQ=21)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GGACC‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                     >  2:2741719/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGCA‑ATCTC‑‑‑C‑GT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                     >  2:2588977/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑AATC‑‑‑T‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAgatc                     >  1:1402919/2‑97 (MQ=18)
                                                                              gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                    >  1:2211422/1‑100 (MQ=255)
                                                                              gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                    >  2:1317900/1‑100 (MQ=255)
                                                                              gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATGCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑CTAGGCCTGATAAGCGGTGCGCATCAGGCAt                    >  2:1143506/1‑100 (MQ=255)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑GAAATCGTGCTAATTCAATATATTGTG‑GGCTC‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                    >  2:715891/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTTA‑‑‑T‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                    >  2:2906945/1‑99 (MQ=21)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATTGTGTAAATTCAATAAATTGAT‑CGATC‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                    >  2:2513525/1‑99 (MQ=21)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CCATC‑‑‑C‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                    >  1:1597626/1‑99 (MQ=21)
                                                                              gCCGGATAAGGCGTTAACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                    >  2:39523/1‑100 (MQ=18)
                                                                               ccggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGTT‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                   >  2:191299/9‑98 (MQ=17)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CCATC‑‑‑C‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                   >  1:1658508/9‑98 (MQ=21)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CCATC‑‑‑C‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                   >  2:2369487/9‑98 (MQ=21)
                                                                               ccggataaggcgTTAACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                   >  2:278676/9‑100 (MQ=18)
                                                                               ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGCGCGCATCAGGCAt                    >  2:1286012/1‑99 (MQ=255)
                                                                               ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                   >  1:1448217/1‑100 (MQ=255)
                                                                               ccGGATAAGGCGTTAACGCCGCATCCGGCATAAAC‑AAAATCGTGCAAAATCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                     >  1:2065620/1‑100 (MQ=17)
                                                                                cggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTG‑T‑GAATC‑‑‑CCTT‑‑GTAGGCCTGATAAGCGTAACGCATCAGGCAttt                  >  1:2752226/8‑99 (MQ=18)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAAC‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                 >  2:300180/7‑100 (MQ=18)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGTT‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                 >  1:1188413/7‑100 (MQ=17)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GGACC‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                 >  2:2682618/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAGATCGTGCAAATTCAATATATTGCA‑TGAGC‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                 >  1:1979177/7‑100 (MQ=17)
                                                                                 ggataaggcgTTCACGCCGCATCCGACCC‑‑AC‑AAAGTCTTGCAAAATCAATAAATTGCAGGAACG‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                 >  2:88122/7‑96 (MQ=25)
                                                                                    taaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGAT‑‑‑C‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttagc              >  1:609151/4‑97 (MQ=17)
                                                                                    tGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttttt              >  1:1390426/1‑98 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg         <  2:2413054/100‑1 (MQ=255)
                                                                                           ttAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc       <  2:1209265/100‑1 (MQ=255)
                                                                                           ttAACGCCTTACCCGGCCT‑‑AC‑AAAATCGTGCATATGCAAGACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGTGCGCATCAGGCATTtttgcctttgcc       <  2:174846/100‑1 (MQ=21)
                                                                                                gCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg  <  1:2565758/100‑1 (MQ=255)
                                                                                                gCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg  >  2:1394979/1‑100 (MQ=255)
                                                                                                   |                                                                                                          
TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAG  >  NC_000913/3705885‑3706080

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