Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,028 T→C 18.6% intergenic (‑238/+70) 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,0280TC18.6% 579.3 / 122.0 276intergenic (‑238/+70)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (106/107);  new base C (47/2);  total (167/109)
Fisher's exact test for biased strand distribution p-value = 1.25e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 290 total aligned reads displayed.

GTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATG  >  NC_000913/3705928‑3706127
                                                                                                     |                                                                                                     
  ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAACTTCAATACATTGCAGAAttt                                                                                                      <  2:2411721/100‑1 (MQ=255)
   taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAtttt                                                                                                     >  1:163557/1‑100 (MQ=255)
     cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                                                                                 <  2:3715525/100‑1 (MQ=255)
     cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                                                                                 <  2:2290613/100‑1 (MQ=255)
     cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                                                                                 <  2:1370813/100‑1 (MQ=255)
       gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTa                                                                                               <  2:1916820/100‑1 (MQ=255)
          tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGc                                                                                            <  2:2698205/100‑1 (MQ=255)
          cataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGc                                                                                            >  2:3227649/6‑100 (MQ=18)
           gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGcc                                                                                           <  1:317377/100‑1 (MQ=255)
           gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGcc                                                                                           <  2:2211694/100‑1 (MQ=255)
            ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCt                                                                                          <  2:2931999/100‑1 (MQ=255)
              gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGa                                                                                        >  2:2556248/1‑100 (MQ=255)
              gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGa                                                                                        <  1:3650780/100‑1 (MQ=255)
                  aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAg                                                                                    <  1:1216033/100‑1 (MQ=255)
                    gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCg                                                                                  <  1:2122398/100‑1 (MQ=255)
                     tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCgg                                                                                 <  1:239360/100‑1 (MQ=255)
                     tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCgg                                                                                 <  2:1329633/100‑1 (MQ=255)
                        gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgc                                                                              >  2:905804/1‑100 (MQ=255)
                        gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgc                                                                              <  1:2997335/100‑1 (MQ=255)
                            tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGATCGATCCAC‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                                          >  2:2006015/8‑100 (MQ=21)
                             ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATc                                                                         >  1:3864568/7‑100 (MQ=21)
                             gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                                         >  1:875368/1‑100 (MQ=255)
                                taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                      >  1:2696846/4‑100 (MQ=18)
                                taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                      >  1:890024/4‑100 (MQ=18)
                                taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGATCGATCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                      <  1:2006025/97‑1 (MQ=21)
                                tATGCCGGCTGCGGCGTTAACGCCTTCTCCGGCCTAC‑TAAATCGTGCAAATTCAATACATTGCATAATTTTC‑‑GTATTGCTGATAAGCGGAGCGCATCAgg                                                                      <  2:2712831/100‑1 (MQ=17)
                                tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                      <  2:3018724/100‑1 (MQ=255)
                                tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                      <  2:150553/100‑1 (MQ=255)
                                 aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                     >  1:1742172/3‑100 (MQ=18)
                                 aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGCA‑ATATGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcg                                                                    >  2:1141131/3‑99 (MQ=21)
                                 aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                     >  2:793094/3‑100 (MQ=17)
                                 aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                     >  1:2884127/3‑100 (MQ=21)
                                 aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                     >  1:1418089/3‑100 (MQ=18)
                                  ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                    >  1:599739/2‑100 (MQ=18)
                                  ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                    >  2:20689/2‑100 (MQ=17)
                                  ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑AATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                    >  1:2194927/2‑100 (MQ=18)
                                  ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑AATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                    >  1:1823785/2‑100 (MQ=18)
                                   gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                   >  1:2328186/1‑99 (MQ=18)
                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                   >  2:3093156/1‑99 (MQ=18)
                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                                   >  1:1109645/1‑99 (MQ=25)
                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                   >  2:592406/1‑99 (MQ=21)
                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCAGGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                   >  1:3892116/1‑99 (MQ=17)
                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCAGGATACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                   >  1:3472357/1‑99 (MQ=17)
                                    ctggataaggcgTTCACGCCGCATCCAGCCTAC‑AAAATTGTGCAAAATCAATGGATTGCACAGCTAGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaga                                                                  >  1:2728497/9‑98 (MQ=25)
                                    ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                                  >  1:915809/9‑98 (MQ=21)
                                    ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCAGGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                  >  1:3074817/9‑98 (MQ=18)
                                    ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGccta                                                                  <  2:1588541/100‑4 (MQ=255)
                                     cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAACTTGATCGATCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                                                                 >  2:670402/8‑100 (MQ=21)
                                      ggataaggcgTTCACGCCGCATCCTGCCTAC‑AAAAGCACGCAAATTCAATATATTGCAGAGATGGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttt                                                                >  1:2198275/7‑100 (MQ=14)
                                      ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                >  2:3946351/7‑100 (MQ=18)
                                       gataaggcgTTTACGCCTTATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                               >  2:1845336/6‑100 (MQ=18)
                                       gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattct                                                               >  2:1274928/6‑98 (MQ=21)
                                       gataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                                                               >  2:2615176/6‑100 (MQ=18)
                                        aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                              <  1:1586966/100‑1 (MQ=255)
                                        aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                              <  2:1883779/100‑1 (MQ=255)
                                         taaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATGGGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattttac                                                             >  2:1429123/4‑98 (MQ=18)
                                          gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                            <  1:985615/100‑1 (MQ=255)
                                           cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                           >  1:3242183/1‑100 (MQ=255)
                                           cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                           >  2:35405/1‑100 (MQ=255)
                                           cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                           >  2:1493828/1‑100 (MQ=255)
                                            ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                                          >  2:1973394/1‑100 (MQ=255)
                                                ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                      <  1:700237/100‑1 (MQ=255)
                                                 tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                     <  1:1128293/100‑1 (MQ=255)
                                                  aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                                    >  1:2086203/1‑100 (MQ=255)
                                                    cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                                  <  1:3749156/100‑1 (MQ=255)
                                                    cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                                  <  1:3306700/100‑1 (MQ=255)
                                                     gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                                 <  2:1244451/100‑1 (MQ=255)
                                                      ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                                <  1:818842/100‑1 (MQ=255)
                                                       cTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc                                               >  1:2925998/1‑100 (MQ=255)
                                                         tATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTTCTTTTTTATCAGTCtt                                             >  2:3352744/1‑100 (MQ=18)
                                                         tATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt                                             >  1:128132/1‑100 (MQ=255)
                                                           tccttccTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGt                                           <  2:3884239/95‑1 (MQ=255)
                                                           tCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGt                                           >  1:1961991/1‑100 (MQ=255)
                                                             cGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                                         <  1:2476189/100‑1 (MQ=255)
                                                               tccTAT‑AAAACCGTTCAAATTCAATACATTGCAGAAGTTTC‑‑GTAGGCCTGATAAGCGGATCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                        <  1:3691515/98‑1 (MQ=25)
                                                                 cTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                                     <  1:1564799/100‑1 (MQ=255)
                                                                 cTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                                     <  1:1342409/100‑1 (MQ=255)
                                                                  tAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                    >  1:175621/1‑100 (MQ=255)
                                                                   aC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                   >  1:1037090/1‑100 (MQ=255)
                                                                    c‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                                  >  1:2366901/1‑100 (MQ=255)
                                                                      aaaaTCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                                 <  1:1657849/100‑1 (MQ=255)
                                                                       aaaTCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTaa                                <  1:4210/100‑1 (MQ=255)
                                                                        aaTCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                               <  1:526428/100‑1 (MQ=255)
                                                                          tCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatta                             <  1:2187846/100‑1 (MQ=255)
                                                                             tGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                          <  2:3824882/100‑1 (MQ=255)
                                                                             tGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                          >  1:1639256/1‑100 (MQ=255)
                                                                              gCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                         >  2:1628334/1‑100 (MQ=255)
                                                                               cAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                         <  2:3489588/99‑1 (MQ=255)
                                                                               cAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa                        >  2:3193994/1‑100 (MQ=255)
                                                                               cAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa                        >  1:2316532/1‑100 (MQ=255)
                                                                                 aaTTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa                      <  2:2113225/100‑1 (MQ=255)
                                                                                           aTTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGtttt            >  1:3035376/1‑100 (MQ=255)
                                                                                            ttGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc           >  2:3787056/1‑100 (MQ=255)
                                                                                              gCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa         >  2:1654540/1‑100 (MQ=255)
                                                                                              gCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa         >  2:1723217/1‑100 (MQ=255)
                                                                                               cAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa        >  1:3116603/1‑100 (MQ=255)
                                                                                               cAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa        >  2:2324992/1‑100 (MQ=255)
                                                                                                aGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg       <  2:455336/100‑1 (MQ=255)
                                                                                                    ttttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGAt   >  1:3590499/1‑100 (MQ=255)
                                                                                                    ttttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGAt   <  2:1104756/100‑1 (MQ=255)
                                                                                                     |                                                                                                     
GTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATG  >  NC_000913/3705928‑3706127

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