Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,048 G→T 22.8% intergenic (‑258/+50) 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,706,0480GT22.8% 523.8 / 128.4 215intergenic (‑258/+50)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (79/82);  new base T (47/2);  total (131/84)
Fisher's exact test for biased strand distribution p-value = 2.88e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 220 total aligned reads displayed.

GCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705950‑3706146
                                                                                                        |                                                                                                  
gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGa                                                                                                   >  1:1695818/1‑100 (MQ=255)
  gcTTTGTTTATGCCGTATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCATGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGa                                                                                                   <  1:845356/98‑1 (MQ=255)
  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                 <  2:1180077/100‑1 (MQ=255)
   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                <  2:1907433/100‑1 (MQ=255)
   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                 <  1:1219329/99‑1 (MQ=255)
      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                                                                                             <  1:3753949/100‑1 (MQ=255)
      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                                                                                             <  2:2250984/100‑1 (MQ=255)
       ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATc                                                                                            >  1:1135662/7‑100 (MQ=18)
         ttATGCCGTATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                                                          <  2:615806/100‑1 (MQ=255)
         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                          >  1:1949991/5‑100 (MQ=18)
         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                         >  2:3120682/5‑100 (MQ=21)
          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                                         >  1:1550646/1‑100 (MQ=255)
           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                        >  1:3201060/3‑100 (MQ=17)
           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                        >  2:3365943/3‑100 (MQ=17)
           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTCCGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                        >  1:722526/3‑100 (MQ=18)
           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGc                                                                                        >  2:1029086/3‑100 (MQ=18)
           aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                                        >  2:2151747/3‑100 (MQ=18)
            ggcaggataaggcgTTCACGCCGCATCCTGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATGGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:1258959/11‑100 (MQ=14)
            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  1:471172/2‑100 (MQ=17)
            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:1784483/2‑100 (MQ=21)
             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                      >  1:3347711/1‑99 (MQ=21)
             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                      >  2:3302709/1‑99 (MQ=17)
              ccggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                      >  2:939175/9‑98 (MQ=18)
              ccGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                     >  1:929404/1‑98 (MQ=18)
               cggataaggcgTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                                                    >  2:335103/8‑97 (MQ=18)
               cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                                                                                    >  2:2844637/1‑100 (MQ=255)
                ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  1:2524882/7‑100 (MQ=17)
                ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCGGGCaattc                                                                                   >  2:635156/7‑96 (MQ=21)
                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                                                   >  2:3181850/7‑100 (MQ=17)
                ggataaggcgTTCACGCCGCATCCAGCCTACAAAATTGTGCAAAATCAATGGATTGCACAGCTAGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagatt                                                                                   >  1:2434254/7‑96 (MQ=25)
                ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                                                                                   >  1:1090170/1‑100 (MQ=255)
                 taggcaaacatcacttcgaaaacgtctttatcgaataattaaaagccatacaagactgatggcaaaggcaaaaatGCCTGATAAGCGGAGCGCATCAGGc                                                                                        >  1:1390709‑M1/76‑100 (MQ=255)
                 gataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgcg                                                                                 >  1:3655984/6‑97 (MQ=21)
                 gataaggcgTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGCAATCTTCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaagctt                                                                                  >  2:3568092/6‑95 (MQ=21)
                  aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  1:1032340/100‑1 (MQ=255)
                  aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  2:70289/100‑1 (MQ=255)
                   taaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgc                                                                                >  2:2685643/4‑100 (MQ=17)
                     gctcgTTAACGGCGTATGCGGCCTACAAAATCTTGCACAGTCAATACATTGCAGAGTTTTC‑‑GTATGCCTGATAAGCGGAGCGCATCAGGCATGtttgc                                                                                <  1:3087845/95‑1 (MQ=17)
                     cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                              >  2:151445/1‑100 (MQ=255)
                       gcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt                                                                            <  2:1597117/100‑1 (MQ=255)
                       aacatcacttcgaaaacgtctctatcgaataattaaaagccatacaagactgatggcaaaggcaaaaatGCCTGATAAGCGGAGCGCATCAGGCAttttt                                                                                  >  1:179055‑M1/70‑100 (MQ=255)
                        cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                           >  2:3285306/1‑100 (MQ=255)
                        acatcacttcgaaaacgtctctatcgaataattaaaagccatacaagactgatggcaaaggcaaaaatGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 >  1:3021037‑M1/69‑100 (MQ=255)
                        acatcacttcgaaaacgtctctatcgaataattaaaagccatacaagactgatggcaaaggcaaaaatGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 >  1:2595677‑M1/69‑100 (MQ=255)
                          ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                                         <  2:2016149/100‑1 (MQ=255)
                             aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                      <  2:2523203/100‑1 (MQ=255)
                             aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                      <  2:3495213/100‑1 (MQ=255)
                              cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                      >  1:3679804/1‑99 (MQ=255)
                                ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATCTTC‑‑GTGGCCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                                                   >  1:164981/1‑100 (MQ=255)
                                   tATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt                                                                >  2:194285/1‑100 (MQ=255)
                                   tATCCGGACTACAAAATCGTGCAAATTCAATCCATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                                                 <  1:535698/99‑1 (MQ=255)
                                        ggCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                                           >  1:396405/1‑100 (MQ=255)
                                        ggCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                                           <  1:2333133/100‑1 (MQ=255)
                                        ggCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                                           >  2:1337342/1‑100 (MQ=255)
                                          ccTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                         <  2:2366167/100‑1 (MQ=255)
                                            tACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                       <  2:1478849/100‑1 (MQ=255)
                                             acacAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCTTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                       <  1:2758435/99‑1 (MQ=255)
                                             aCAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                                      >  2:984993/1‑100 (MQ=255)
                                                 aaaTCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                                                    >  1:3335552/1‑99 (MQ=255)
                                                    tCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatta                                                <  2:2307533/100‑1 (MQ=255)
                                                       tGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  2:1067822/100‑1 (MQ=255)
                                                       tGCAAATGCAATACATTGCAGAATTTTA‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTACCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  1:2252360/100‑1 (MQ=255)
                                                        gCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                            <  1:3550230/100‑1 (MQ=255)
                                                            aTTCAATACATTGTAGAATTTTC‑‑GTAGCTCTGATAAGCGGATCGTATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGTCTTTTAATTATTc                                             <  1:2350371/95‑1 (MQ=37)
                                                             ttCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga                                       >  2:1874643/1‑100 (MQ=255)
                                                             ttCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga                                       <  2:882728/100‑1 (MQ=255)
                                                             ttCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga                                       <  1:1297889/100‑1 (MQ=255)
                                                                    cATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt                                >  1:2555635/1‑100 (MQ=255)
                                                                      ttGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc                              <  1:265547/100‑1 (MQ=255)
                                                                          aGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                          <  1:1335982/100‑1 (MQ=255)
                                                                            gAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                         >  2:881959/1‑100 (MQ=255)
                                                                            gAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                         <  1:198666/100‑1 (MQ=255)
                                                                               aTTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa                       <  1:66391/100‑1 (MQ=255)
                                                                                    tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGtt                   <  1:2844620/100‑1 (MQ=255)
                                                                                    tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGtt                   >  2:3051554/1‑100 (MQ=255)
                                                                                     c‑‑GTAGGCTTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGttt                  <  1:555104/100‑1 (MQ=255)
                                                                                        gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                 <  1:3483059/100‑1 (MQ=255)
                                                                                          aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  1:1400031/100‑1 (MQ=255)
                                                                                           ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt              >  2:969394/1‑100 (MQ=255)
                                                                                           ggCCTGATAAGCGGAGCGCATCAGGCATTTATGACTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt              >  1:2501592/1‑100 (MQ=255)
                                                                                             ccTGATAAGCGGAGCGCATCAGGCATTTTTTTTTTTTTTTTCAGTCTTGTATGGCTTTTAATTATTCTTTAGAGAAGTTTTCGAAGTGATGTTTGCCt              >  1:770004/1‑98 (MQ=25)
                                                                                             CCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTaa            <  2:2873811‑M1/100‑1 (MQ=255)
                                                                                              CTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt           <  2:1638078‑M1/100‑1 (MQ=255)
                                                                                              CTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt           >  2:3732197‑M1/1‑100 (MQ=255)
                                                                                              CTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt           >  2:608271‑M1/1‑100 (MQ=255)
                                                                                               TGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTGGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg          <  1:1866607‑M1/100‑1 (MQ=255)
                                                                                               TGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg          >  1:3452776‑M1/1‑100 (MQ=255)
                                                                                                gATAAGCGGAGCGCATCAGGCATTTTTTTTTTTTTTATCAGTCTTGTATGGCTTTTAATTATTTTATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt           >  1:1230176/1‑98 (MQ=255)
                                                                                                GATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATgg         >  2:1401882‑M1/1‑100 (MQ=255)
                                                                                                GATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATgg         >  1:2198281‑M1/1‑100 (MQ=255)
                                                                                                GATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATgg         <  2:772033‑M1/100‑1 (MQ=255)
                                                                                                GATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATgg         >  2:2816142‑M1/1‑100 (MQ=255)
                                                                                                  TAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt       <  2:1611590‑M1/100‑1 (MQ=255)
                                                                                                   aaGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATc      <  2:1695825/100‑1 (MQ=255)
                                                                                                    aGCGGAGCGCATCAGGCAGTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGGTAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa     >  1:749820/1‑100 (MQ=255)
                                                                                                      cGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAc   <  1:2167308/100‑1 (MQ=255)
                                                                                                      cGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAc   <  2:396407/100‑1 (MQ=255)
                                                                                                      cGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAc   >  2:3526130/1‑100 (MQ=255)
                                                                                                       ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc  <  2:1163944/100‑1 (MQ=255)
                                                                                                       ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc  <  1:2823039/100‑1 (MQ=255)
                                                                                                        |                                                                                                  
GCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705950‑3706146

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