Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,978 A→C 26.8% intergenic (‑188/+120) 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,9780AC26.8% 558.2 / 212.9 285intergenic (‑188/+120)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base A (81/112);  new base C (62/2);  total (171/114)
Fisher's exact test for biased strand distribution p-value = 6.83e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-03

Only 100 of 289 total aligned reads displayed.

GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705879‑3706077
                                                                                                   |                                                                                                            
   gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACGAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                           <  1:2789006/100‑1 (MQ=255)
   gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                           <  1:1300910/100‑1 (MQ=255)
        gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                      >  2:3754103/1‑100 (MQ=255)
        gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                      <  1:532975/100‑1 (MQ=255)
        gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                      >  1:3018742/1‑100 (MQ=255)
        gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                      <  1:1774437/100‑1 (MQ=255)
         ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                     <  2:1280650/100‑1 (MQ=255)
          tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                    >  1:150553/1‑100 (MQ=255)
           gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                   >  1:830343/1‑100 (MQ=255)
           gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                   >  2:1896680/1‑100 (MQ=255)
           gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                   <  1:2589120/100‑1 (MQ=255)
             cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                 >  2:3629323/1‑100 (MQ=255)
              ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                <  1:1981352/100‑1 (MQ=37)
              ggATGGCGCTGCACGCAGACGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTGGTTTATGCCGTTTGCGGCTTTAACGCCtttttcg                                                                                                  >  2:1585974/1‑93 (MQ=25)
                aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                               <  2:3391859/99‑1 (MQ=255)
                aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                              <  2:1868115/100‑1 (MQ=255)
                aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                              <  2:2541420/100‑1 (MQ=255)
                 tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑a                                                                                           <  1:697369/100‑1 (MQ=255)
                   tcgcTGCACGCAGAGGCACAGCGCCACCGGGGATACAGATGGTTGACAAAGTGCGCGTTGGGTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACa                                                                                        >  1:3145757/2‑100 (MQ=21)
                     gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACaaa                                                                                      >  1:1588550/1‑100 (MQ=255)
                      cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACaaaa                                                                                     >  2:2187834/1‑100 (MQ=255)
                      cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATACCGGATCCGGCGTTAACGCCTTATCCGCCCT‑‑ACaaa                                                                                      >  2:2724628/1‑99 (MQ=255)
                       tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAAt                                                                                    <  1:3629351/100‑1 (MQ=255)
                       tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAAt                                                                                    >  1:3256633/1‑100 (MQ=255)
                        gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATc                                                                                   >  1:1916829/1‑100 (MQ=255)
                         gaCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACACAGTACGCTATGTTTATGCCGAATGCGTCTTTAACGCCTTATCCGTCCT‑‑Acgca                                                                                       >  2:1485222/2‑93 (MQ=25)
                         cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCg                                                                                  <  2:2286037/100‑1 (MQ=255)
                         cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCg                                                                                  <  1:2869382/100‑1 (MQ=255)
                          aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGt                                                                                 <  2:834424/100‑1 (MQ=255)
                           cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGt                                                                                 >  2:1686935/1‑99 (MQ=255)
                           cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTg                                                                                <  1:3787153/100‑1 (MQ=255)
                            gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGc                                                                               >  1:3551049/1‑100 (MQ=255)
                                  gCACAGCGCTACCGGGTATTCAGATGGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTc                                                                         <  2:1924443/100‑1 (MQ=255)
                                  gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTc                                                                         <  1:1653585/100‑1 (MQ=255)
                                   cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCa                                                                        >  1:2314827/1‑100 (MQ=255)
                                    acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCaa                                                                       <  1:1145316/100‑1 (MQ=255)
                                           cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTg                                                                <  2:1773465/100‑1 (MQ=255)
                                              cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAg                                                             <  2:1039507/100‑1 (MQ=255)
                                               gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGa                                                            >  1:727708/1‑100 (MQ=255)
                                                 gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAATCGTGCAAATTCAATACATTGCAGAAtt                                                        >  1:1377222/1‑100 (MQ=255)
                                                   ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAACTTCAATACATTGCAGAAtt‑‑‑t                                                    <  2:2411721/100‑1 (MQ=255)
                                                    taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAAtt‑‑‑tt                                                   >  1:929310/1‑100 (MQ=255)
                                                    taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAAtt‑‑‑tt                                                   >  2:1259000/1‑100 (MQ=255)
                                                    taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAAtt‑‑‑tt                                                   >  1:1104764/1‑100 (MQ=255)
                                                     aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTc                                                  <  2:296087/100‑1 (MQ=255)
                                                      cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑g                                               <  2:1370813/100‑1 (MQ=255)
                                                           cataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGc                                          >  2:3227649/6‑100 (MQ=18)
                                                              tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTg                                       >  1:2901318/1‑100 (MQ=255)
                                                                   aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAg                                  <  1:1216033/100‑1 (MQ=255)
                                                                      tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                               <  2:1329633/100‑1 (MQ=255)
                                                                       gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                              <  1:3423430/100‑1 (MQ=255)
                                                                         gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                            <  1:2997335/100‑1 (MQ=255)
                                                                          cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                           >  2:2019414/1‑100 (MQ=255)
                                                                          cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                           >  2:1503135/1‑100 (MQ=255)
                                                                            ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCa                         <  2:1123426/100‑1 (MQ=255)
                                                                              ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATc                       >  1:3864568/7‑100 (MQ=21)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                       >  1:875368/1‑100 (MQ=255)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                       <  2:3528445/100‑1 (MQ=255)
                                                                               tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                      <  1:3694331/100‑1 (MQ=255)
                                                                               tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                      >  1:307054/1‑100 (MQ=255)
                                                                                gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAATCGTGCAAATTCAACATATTGCAATTCT‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                     >  1:2104335/5‑100 (MQ=18)
                                                                                gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                     >  2:2324274/5‑100 (MQ=18)
                                                                                 taggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  1:1637322/4‑100 (MQ=18)
                                                                                 taggCCGGATAAGGCGTTCACTCCGCATCCGGCCT‑‑ACAAAAGCACGCAAATTCAATATATTGCAGAGAT‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                    >  2:813508/4‑100 (MQ=14)
                                                                                 taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                    >  1:2696846/4‑100 (MQ=18)
                                                                                 tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                    <  2:3018724/100‑1 (MQ=255)
                                                                                  aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  2:210147/3‑100 (MQ=17)
                                                                                  aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  1:70706/3‑100 (MQ=17)
                                                                                  aggCCGGATAAGGCGTTAACGCCGCATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                   >  1:1418089/3‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                 >  1:2442379/2‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:1631156/2‑100 (MQ=17)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACTTTGCAAATTCAATATATTGC‑‑AATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  1:1823785/2‑100 (MQ=18)
                                                                                    gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                 >  1:2328186/1‑99 (MQ=18)
                                                                                    gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAATCGTGCAAATTCAACATATTGCAATTCT‑‑‑CTT‑‑GTAGGTCTGATAAGCGTAGCGCATCAGGctg                 >  1:563498/1‑98 (MQ=18)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                 >  2:3093156/1‑99 (MQ=18)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                 >  2:592406/1‑99 (MQ=21)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAACCTGCAAATTCAATATATTGCAGGATA‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                 >  1:3472357/1‑99 (MQ=17)
                                                                                     ctggataaggcgTTCACGCCGCATCCAGCCT‑‑ACAAAATTGTGCAAAATCAATGGATTGCACAGCT‑‑‑AGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaga                >  1:2728497/9‑98 (MQ=25)
                                                                                     ccggataaggcgTTTACGCCTTATCCGGCCT‑‑ACAAAACAGTGCAAATTCAATATATTGCAATCTT‑‑‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaag                >  1:3911787/9‑98 (MQ=18)
                                                                                     ccggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                >  1:1746472/9‑98 (MQ=17)
                                                                                     ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAAACATGCAAATTCAATAGATTGCAGAGAT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  1:201161/1‑100 (MQ=17)
                                                                                      cggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAATTGTGTAAATTCAATAACTTGATCGATC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt               >  2:670402/8‑100 (MQ=21)
                                                                                       ggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAAAATGTGCAAATTCAATAAATTGC‑‑AATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt               >  2:3871751/7‑97 (MQ=21)
                                                                                       ggataaggcgTTCACGCCGCATCCTGCCT‑‑ACAAAAGCACGCAAATTCAATATATTGCAGAGAT‑‑‑GGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttt              >  1:2198275/7‑100 (MQ=14)
                                                                                       ggataaggcgTTCACGCCGCATCCTGCCT‑‑ACAAAAGCACGCAAATTCAATATATTGCAGAGAT‑‑‑GGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttt              >  1:2191830/7‑100 (MQ=14)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt              >  2:1516479/7‑100 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt              >  2:3946351/7‑100 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct              >  1:3203297/7‑96 (MQ=21)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt              >  1:212303/7‑100 (MQ=18)
                                                                                        gataaggcgTTTACGCCTTATCCGGCCT‑‑ACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt             >  2:1845336/6‑100 (MQ=18)
                                                                                        gataaggcgTTAACGCCGCATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt             >  2:2615176/6‑100 (MQ=18)
                                                                                         aTGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg            <  1:3804108/100‑1 (MQ=255)
                                                                                          taaggcgTTTACGCCGCATCCGGCCT‑‑ACAAAAGCGTGCAAATTCAATACATTGCATGGGC‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattttac           >  2:1429123/4‑98 (MQ=18)
                                                                                          tGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAACTCGTGTAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc           <  2:2668498/100‑1 (MQ=255)
                                                                                          tGCGGCGTTAACGCCTTATCCGGCCT‑‑AAAAAATCGGGCAAAGTCAATACATTGCAGCATC‑‑‑GTC‑‑GTAGGCCTGATAAGCGGCCCGCAGCAGGCATTcttgc           >  2:3577533/1‑100 (MQ=17)
                                                                                            cggcgTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttttttct         >  2:1180819/1‑96 (MQ=255)
                                                                                             ggcgTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt        >  2:1973394/1‑100 (MQ=255)
                                                                                              gcgTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt       >  2:2926462/1‑100 (MQ=255)
                                                                                                gTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc     <  1:905809/100‑1 (MQ=255)
                                                                                                  tAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa   <  1:1116322/100‑1 (MQ=255)
                                                                                                   |                                                                                                            
GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705879‑3706077

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: