Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 32.6% 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,9850TC32.6% 379.7 / 244.5 261intergenic (‑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 (101/75);  new base C (77/8);  total (178/83)
Fisher's exact test for biased strand distribution p-value = 1.66e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.46e-04

Only 100 of 265 total aligned reads displayed.

TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705886‑3706084
                                                                                                   |                                                                                                             
tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                               <  1:3238809/100‑1 (MQ=255)
     ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                          <  2:622861/100‑1 (MQ=255)
       ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                        >  1:548012/1‑100 (MQ=255)
       ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGGATACGGATTGTTGACAAATTGCGCGTTGTTTATGCCGGATGCGGCGTTAACGCTTTATCCgg                                                                                                         >  1:2311048/1‑99 (MQ=25)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                       <  1:2948775/100‑1 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                       <  2:3165294/100‑1 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                       <  2:3063968/100‑2 (MQ=37)
         aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                      <  1:1033321/100‑1 (MQ=255)
          tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                                     <  2:1193575/100‑1 (MQ=255)
          tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGCTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCTTa                                                                                                     >  1:1599127/1‑100 (MQ=255)
            gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                                  <  1:1800330/100‑1 (MQ=255)
             cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                                                 >  1:1218937/1‑100 (MQ=255)
             cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaac                                                                                                  <  2:1732678/100‑7 (MQ=25)
                 gCACGCAGAGGCGCAGCGCCACCGGGTATACAGATTGTGGGCAACGTGCGCTCCGTTTATGCCGGATGCGGCGTTAACGCCTTCTCCGGCCTAC‑AAAATc                                                                                             >  1:1184230/1‑100 (MQ=21)
                 gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaacaaag                                                                                              <  2:415097/100‑11 (MQ=21)
                   aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                                           >  1:2882100/1‑100 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                          <  1:242753/100‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                          <  1:3354362/100‑1 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                         <  2:185679/100‑1 (MQ=255)
                      cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                                        >  2:1268479/1‑100 (MQ=255)
                          ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                                    >  1:3179000/1‑100 (MQ=255)
                           gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                                   <  2:3398772/100‑1 (MQ=255)
                            cacaGCGCTTCCGGGTATACGGGTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGAGGCGGTAACGCCTTATACGGCCTAC‑ACAATCTTGCAAATTCa                                                                                  >  1:3848114/1‑100 (MQ=17)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                                  >  2:947064/1‑100 (MQ=255)
                               aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATa                                                                               >  2:887035/1‑100 (MQ=255)
                                gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTACAAATTCAATa                                                                               <  1:2250118/99‑1 (MQ=255)
                                  gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                                            >  2:1978827/1‑100 (MQ=255)
                                  gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                                            >  2:693932/1‑100 (MQ=255)
                                  gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                                            >  2:238435/1‑100 (MQ=255)
                                       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑g                                                                    <  1:3087587/100‑1 (MQ=255)
                                        gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑Ga                                                                   >  2:3861235/1‑100 (MQ=255)
                                        gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑Ga                                                                   >  2:1027286/1‑100 (MQ=255)
                                              aCAGATTGTTGACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAACTCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTc                                                         >  1:2583587/1‑100 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑g                                                      >  2:2095549/1‑100 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑g                                                      <  1:1987050/100‑1 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑g                                                      >  2:3746948/1‑100 (MQ=255)
                                                  aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAg                                                   >  1:2609808/1‑100 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGcc                                                >  2:2157340/1‑100 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGcc                                                <  2:2797448/100‑1 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGcc                                                <  2:3608908/100‑1 (MQ=255)
                                                      ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCt                                               >  2:3259422/1‑100 (MQ=255)
                                                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGa                                             >  1:2017776/1‑100 (MQ=255)
                                                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGa                                             >  1:2098217/1‑100 (MQ=255)
                                                         aCAAAGTGCGCGTTGTTTATGCCGGATGCGCCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTCCCAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑TTCGGCCTGAt                                            >  1:839425/1‑100 (MQ=17)
                                                         aCAAAGTGCGATTTGTTTATGCCGGATGCGGCGATAACGCCTTATCCGCCCTAC‑AAAATCGTGCAAATTCAATCCATTGC‑‑A‑GAAT‑T‑‑‑CGC‑‑GTAGGCCTGAt                                            >  1:2736989/1‑100 (MQ=17)
                                                           aaaGTGCGCTTTGTTTATGCCGTATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GATT‑T‑‑‑TTC‑‑GTAGGCCTGATa                                           <  1:2920906/99‑1 (MQ=255)
                                                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGc                                        <  2:785460/100‑1 (MQ=255)
                                                                gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                                     >  1:3663284/1‑100 (MQ=255)
                                                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                    >  1:3030765/1‑100 (MQ=255)
                                                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                    >  2:4818/1‑100 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                   <  2:1758935/100‑1 (MQ=255)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                  >  1:939811/1‑100 (MQ=255)
                                                                    tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcgc                                 >  1:1466270/1‑100 (MQ=255)
                                                                     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                >  2:23519/1‑100 (MQ=255)
                                                                       ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GAAA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATc                              >  1:650746/7‑100 (MQ=18)
                                                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                            >  2:743044/1‑100 (MQ=255)
                                                                          taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑‑‑AAT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAg                            >  2:3685799/4‑100 (MQ=21)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGC‑‑A‑ATA‑‑T‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:2084815/4‑100 (MQ=18)
                                                                          caTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAC‑AAGATCGTGCAAATACAATATATTGC‑‑A‑TGAG‑C‑‑‑CAC‑‑GTAGGACTGGTAAGCGGacagcaaca                             >  1:3426646/2‑90 (MQ=17)
                                                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGC‑‑ATGGGC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:156068/3‑100 (MQ=17)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑AGGGAC‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc                          >  2:2768106/3‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAACCGCGCAAATTCAATATATTGCGGA‑GAAA‑A‑‑‑T‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                         >  1:1755912/3‑99 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GGAC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:1100467/3‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                          >  2:2069224/3‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GAAT‑CC‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                        >  1:434921/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GGAC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:244738/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGC‑‑A‑‑ATA‑T‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcga                        >  2:2219880/2‑98 (MQ=21)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGA‑‑T‑CGAT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:2770704/2‑100 (MQ=21)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑AGGGAC‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                         >  2:1209031/2‑100 (MQ=18)
                                                                             gtcggataaggcgTTCACGCCGCATCCGGCCCAC‑AAAGTCTTGCAAAATCAATAAATTGC‑‑A‑GGAA‑C‑‑‑GAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                        >  2:3486112/10‑99 (MQ=18)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTCTTCCAAATTCAATATATTGC‑‑A‑GAAA‑T‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                        >  2:2890111/1‑99 (MQ=18)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                       >  2:2177810/1‑100 (MQ=18)
                                                                             gCCGGATAAGGCGTTCACTCCGCATCCGGCCTAC‑AAAAGCACGCAAATTCAATATATTGC‑‑A‑GAGA‑T‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                           >  1:880512/1‑97 (MQ=11)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TGAG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  1:3000439/1‑99 (MQ=17)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑A‑ATCT‑C‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                        >  2:2519857/1‑99 (MQ=18)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GGAG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                       >  1:3099349/1‑100 (MQ=21)
                                                                             gCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        >  1:3835333/1‑100 (MQ=17)
                                                                             gCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        >  2:1564141/1‑100 (MQ=18)
                                                                              ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGC‑‑ATGGGC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                       >  1:2456973/9‑98 (MQ=17)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑AAT‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                       >  2:2548885/9‑98 (MQ=18)
                                                                               cggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GAATCC‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttg                     >  1:3559356/8‑99 (MQ=18)
                                                                               cGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                      <  2:3037115/100‑1 (MQ=255)
                                                                               cGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                      >  1:2409372/1‑100 (MQ=255)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GGAC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:1094732/7‑100 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CCAT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattc                     >  2:2635096/7‑96 (MQ=21)
                                                                                 gataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GAAA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                    >  1:2218479/6‑100 (MQ=17)
                                                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  1:2886082/100‑1 (MQ=255)
                                                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  1:1268471/100‑1 (MQ=255)
                                                                                   tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                  >  2:2250515/1‑100 (MQ=255)
                                                                                     cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                >  2:2692477/1‑100 (MQ=255)
                                                                                       gcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt              >  2:2657045/1‑100 (MQ=255)
                                                                                        cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttn             <  1:3774941/100‑2 (MQ=255)
                                                                                        cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             <  2:3315362/100‑1 (MQ=255)
                                                                                         gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            >  2:1423063/1‑100 (MQ=255)
                                                                                         gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             >  1:86674/1‑99 (MQ=255)
                                                                                           tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa          >  2:2298617/1‑100 (MQ=255)
                                                                                             aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  2:293755/100‑1 (MQ=255)
                                                                                             aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  2:884729/100‑1 (MQ=255)
                                                                                               gCCTTCTCCGGCCTCC‑AAAATCGTGCAAATTCCATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTCGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg      >  1:1652789/1‑100 (MQ=21)
                                                                                               gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg      <  2:1511309/100‑1 (MQ=255)
                                                                                                   |                                                                                                             
TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705886‑3706084

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: ‑