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‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑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‑‑ACa                                                                                            >  2:2135555/1‑100 (MQ=255)
              cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACaa                                                                                           >  2:1580979/1‑100 (MQ=255)
                  gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATc                                                                                       >  1:2114859/1‑100 (MQ=255)
                    aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGt                                                                                     <  1:1201430/100‑1 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACACAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGc                                                                                   <  1:2682228/100‑1 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGc                                                                                   <  1:1628590/100‑1 (MQ=255)
                       cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCa                                                                                  >  2:1863234/1‑100 (MQ=255)
                       cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAAGACTTATCCGGCCT‑‑ACAAAATCGTGc                                                                                   <  1:266293/99‑1 (MQ=255)
                            gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTc                                                                             <  1:1164818/100‑1 (MQ=255)
                            gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTc                                                                             <  2:1442715/100‑1 (MQ=255)
                              acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCaa                                                                           <  1:748110/100‑1 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAAt                                                                          >  2:1036276/1‑100 (MQ=255)
                                      aacGGGTATACAGATTGTTGACAAAGTGCGCATTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGc                                                                   >  2:1960069/3‑100 (MQ=255)
                                      aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGc                                                                   <  2:2662823/100‑1 (MQ=255)
                                       cggggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACAtt                                                                     <  1:319334/95‑1 (MQ=255)
                                              taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAAtt‑‑‑tt                                                      >  1:2285169/1‑100 (MQ=255)
                                              taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAAtt‑‑‑tt                                                      >  2:1833524/1‑100 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑g                                                  >  2:1486724/1‑100 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑g                                                  <  1:1070657/100‑1 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑g                                                  <  1:767779/100‑1 (MQ=255)
                                                 aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑Gt                                                 <  1:1056153/100‑1 (MQ=255)
                                                 aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTCAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑Gt                                                 <  1:1863226/100‑1 (MQ=255)
                                                  gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTa                                                <  1:2253249/100‑1 (MQ=255)
                                                   aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAg                                               >  2:2234618/1‑100 (MQ=255)
                                                    ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAgg                                              <  1:2278694/100‑1 (MQ=255)
                                                       ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGcc                                            <  1:536930/99‑1 (MQ=255)
                                                          aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGAt                                        <  1:1023825/100‑1 (MQ=255)
                                                          aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGAt                                        >  2:1373498/1‑100 (MQ=255)
                                                           cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATa                                       >  1:2423367/1‑100 (MQ=255)
                                                           cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATa                                       <  2:1553936/100‑1 (MQ=255)
                                                              aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGc                                    <  1:1719736/100‑1 (MQ=255)
                                                               gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCg                                   <  2:2066772/100‑1 (MQ=255)
                                                               gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCg                                   >  1:8102/1‑100 (MQ=255)
                                                                 gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                                 >  1:1134315/1‑100 (MQ=255)
                                                                   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                               >  1:1899801/1‑100 (MQ=255)
                                                                   gcTTTGTTTATGCCGGATGCGGCCTTAACCCCATAGCCGGCCT‑‑ACAAAGTCGTGCAAACTCCGTACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATCAGCGGAgc                               <  1:211091/100‑1 (MQ=18)
                                                                     attttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATTGTGTAAATTCAATAAATTGATCGATC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAgcgc                             >  2:796438/10‑100 (MQ=18)
                                                                     attttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcgc                             >  2:1839977/10‑100 (MQ=18)
                                                                      ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCa                            >  2:582712/1‑100 (MQ=255)
                                                                       ttcgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCGT‑‑ACGAAATCGTGCTAATTCAATATATTGTGGGCTC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCAt                           >  2:2311783/8‑100 (MQ=18)
                                                                       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                           >  1:287684/1‑100 (MQ=255)
                                                                        ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACCTTGCTAATTCAATATATTGCAGGGAC‑‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATc                          >  1:2351461/7‑100 (MQ=18)
                                                                        ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                         >  2:2138085/7‑100 (MQ=21)
                                                                        gatctggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACTTTGCAAATTCAATATATTGCAATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATc                          <  2:1402921/94‑1 (MQ=18)
                                                                        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                          >  2:501784/1‑100 (MQ=255)
                                                                        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                          <  1:2447633/100‑1 (MQ=255)
                                                                         tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATc                          >  1:1920233/6‑100 (MQ=21)
                                                                          gtaggCCGGATAAGGCGTTCACTCCGCATCCGGCCT‑‑ACAAAAGCACGCAAATTCAATATATTGCAGAGAT‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                        >  1:373862/5‑100 (MQ=14)
                                                                           taggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                       >  2:621294/4‑100 (MQ=17)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                       >  2:739538/4‑100 (MQ=25)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                       >  2:2291973/4‑100 (MQ=21)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                       >  2:2519674/4‑100 (MQ=21)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                       <  1:1396944/97‑1 (MQ=21)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                       >  2:1634854/4‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                      >  2:2882801/3‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                      >  2:889574/3‑100 (MQ=17)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAGCGTGCAAATTCAATACATTGCATGGGCC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                      >  1:2129657/3‑100 (MQ=17)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                      >  1:1492369/3‑100 (MQ=17)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAACCTGCAAATTCAATATATTGCAGGATA‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                      >  1:150656/3‑100 (MQ=17)
                                                                            aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                      <  2:199529/100‑1 (MQ=255)
                                                                             ggctggataaggcgTTCACGCCGCATCCAGCCT‑‑ACAAAATTGTGCAAAATCAATGGATTGCACAGCT‑‑‑AGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                     >  1:1699792/11‑100 (MQ=25)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                      >  2:2212556/2‑100 (MQ=21)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                     >  2:2741719/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACTTTGCAAATTCAATATATTGCAATCTC‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                     >  2:2588977/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACTTTGCAAATTCAATATATTGCAATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAgatc                     >  1:1402919/2‑97 (MQ=18)
                                                                              gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                    >  1:2211422/1‑100 (MQ=255)
                                                                              gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                    >  2:1317900/1‑100 (MQ=255)
                                                                              gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATGCAATACATTGCAGAATT‑‑‑TTC‑‑CTAGGCCTGATAAGCGGTGCGCATCAGGCAt                    >  2:1143506/1‑100 (MQ=255)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACGAAATCGTGCTAATTCAATATATTGTGGGCTC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                    >  2:715891/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                    >  2:2906945/1‑99 (MQ=21)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATTGTGTAAATTCAATAAATTGATCGATC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                    >  2:2513525/1‑99 (MQ=21)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                    >  1:1597626/1‑99 (MQ=21)
                                                                              gCCGGATAAGGCGTTAACGCCGCATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                    >  2:39523/1‑100 (MQ=18)
                                                                               ccggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                   >  2:191299/9‑98 (MQ=17)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                   >  1:1658508/9‑98 (MQ=21)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                   >  2:2369487/9‑98 (MQ=21)
                                                                               ccggataaggcgTTAACGCCGCATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                   >  2:278676/9‑100 (MQ=18)
                                                                               ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGCGCGCATCAGGCAt                    >  2:1286012/1‑99 (MQ=255)
                                                                               ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                   >  1:1448217/1‑100 (MQ=255)
                                                                               ccGGATAAGGCGTTAACGCCGCATCCGGCATAAACAAAATCGTGCAAAATCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                     >  1:2065620/1‑100 (MQ=17)
                                                                                cggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAATTCTGCAAATTCAACAAATTGTGAATC‑‑‑CCTT‑‑GTAGGCCTGATAAGCGTAACGCATCAGGCAttt                  >  1:2752226/8‑99 (MQ=18)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                 >  2:300180/7‑100 (MQ=18)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                 >  1:1188413/7‑100 (MQ=17)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                 >  2:2682618/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAGATCGTGCAAATTCAATATATTGCATGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                 >  1:1979177/7‑100 (MQ=17)
                                                                                 ggataaggcgTTCACGCCGCATCCGACCC‑‑ACAAAGTCTTGCAAAATCAATAAATTGCAGGAACG‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                 >  2:88122/7‑96 (MQ=25)
                                                                                    taaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAGCACGCAAACTCAATATATTGCAGAGAT‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttagc              >  1:609151/4‑97 (MQ=17)
                                                                                    tGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttttt              >  1:1390426/1‑98 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg         <  2:2413054/100‑1 (MQ=255)
                                                                                           ttAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc       <  2:1209265/100‑1 (MQ=255)
                                                                                           ttAACGCCTTACCCGGCCT‑‑ACAAAATCGTGCATATGCAAGACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGTGCGCATCAGGCATTtttgcctttgcc       <  2:174846/100‑1 (MQ=21)
                                                                                                gCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg  <  1:2565758/100‑1 (MQ=255)
                                                                                                gCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg  >  2:1394979/1‑100 (MQ=255)
                                                                                                   |                                                                                                          
TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑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: