Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,668,951 T→G R214R (CGA→CGC)  ynfL ← LysR family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,668,9510TG98.1% 545.1 / ‑3.7 156R214R (CGA→CGC) ynfLLysR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base G (95/57);  total (98/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 156 total aligned reads displayed.

CAGGCAAAATTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGTCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTCAGCCAGCGTTACATT  >  NC_000913/1668854‑1669046
                                                                                                 |                                                                                               
cAGGCAAAATTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCg                                                                                               <  1:3014391/100‑1 (MQ=255)
  ggCAAAATTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCa                                                                                             >  1:677561/1‑100 (MQ=255)
  ggCAAAATTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCa                                                                                             >  1:5179408/1‑100 (MQ=255)
   gCAAAATTGAAACACCCAGATCGGCGGCAACCAGACCGATGATGGTCAGTGCCTGGCACACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCAt                                                                                            >  2:5480392/1‑100 (MQ=255)
   gCAAAATTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCAt                                                                                            <  1:138507/100‑1 (MQ=255)
    cAAAATTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATc                                                                                           >  2:3763606/1‑100 (MQ=255)
    cAAAATTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATc                                                                                           >  1:4693449/1‑100 (MQ=255)
    cAAAATTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATc                                                                                           >  2:1443649/1‑100 (MQ=255)
     aaaaTTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCa                                                                                          <  1:3158962/100‑1 (MQ=255)
     aaaaTTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCa                                                                                          >  1:3016154/1‑100 (MQ=255)
        aTTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGcc                                                                                       >  2:5453845/1‑100 (MQ=255)
        aTTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGcc                                                                                       <  1:2067601/100‑1 (MQ=255)
        aTTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGcc                                                                                       >  1:899187/1‑100 (MQ=255)
          tGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCg                                                                                     <  1:2175062/100‑1 (MQ=255)
            aaaCACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCgag                                                                                   <  1:4282921/100‑1 (MQ=255)
            aaaCACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCgag                                                                                   <  1:3798467/100‑1 (MQ=255)
             aaCACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGTCGCATCAGCCCgaga                                                                                  >  1:886454/1‑100 (MQ=255)
             aaCACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCgaga                                                                                  >  1:144673/1‑100 (MQ=255)
              acacCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGaa                                                                                 <  1:2678151/100‑1 (MQ=255)
              acacCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGaa                                                                                 <  2:1727931/100‑1 (MQ=255)
                acCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAata                                                                               >  1:1318517/1‑100 (MQ=255)
                acCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAata                                                                               >  2:785453/1‑100 (MQ=255)
                 cccAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAatat                                                                              >  2:2925541/1‑100 (MQ=255)
                 cccAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAatat                                                                              <  2:3369448/100‑1 (MQ=255)
                 cccAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAatat                                                                              >  2:2007583/1‑100 (MQ=255)
                  ccAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATc                                                                             >  2:797667/1‑100 (MQ=255)
                     gACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTc                                                                          <  1:5003231/100‑1 (MQ=255)
                      aCCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCa                                                                         <  1:3962554/100‑1 (MQ=255)
                       ccGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCAt                                                                        >  1:3817964/1‑100 (MQ=255)
                         ggcggAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATAc                                                                      <  1:1179006/100‑1 (MQ=255)
                         ggcggAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATAc                                                                      >  1:192444/1‑100 (MQ=255)
                            ggAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTACCGGCGCCTCAGCCCGAGAATATCGTCATACAGc                                                                   >  2:3706597/1‑100 (MQ=255)
                              aaaCCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGccc                                                                 >  2:4014959/1‑100 (MQ=255)
                                 ccAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGt                                                              >  1:1797800/1‑100 (MQ=255)
                                 ccAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGt                                                              >  2:3180282/1‑100 (MQ=255)
                                  cAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTc                                                             >  2:2164617/1‑100 (MQ=255)
                                   aGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTcc                                                            >  2:415492/1‑100 (MQ=255)
                                    gACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTccc                                                           >  1:3617502/1‑100 (MQ=255)
                                      ccGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGa                                                         >  1:228221/1‑100 (MQ=255)
                                          tgatgGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCGGTCCCGACGTg                                                     >  2:3594371/1‑100 (MQ=255)
                                             tgGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCgg                                                  >  2:1396572/1‑100 (MQ=255)
                                              ggTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGa                                                 >  2:3276300/1‑100 (MQ=255)
                                               gTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGAt                                                >  1:278777/1‑100 (MQ=255)
                                               gTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATCGTAACGACGCAGCAGCCCGAGAATATATTAATACATCCATGTCCCGACGTGCGGAt                                                >  2:3105425/1‑100 (MQ=255)
                                                tCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATc                                               >  1:4236244/1‑100 (MQ=255)
                                                tCATTGCCTCGCCCACCTCCTGAGCGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATc                                               <  1:1633871/100‑1 (MQ=255)
                                                 cATTGCCTCGCCCACCTCCTGAGTGATGGCGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCa                                              >  2:4441981/1‑100 (MQ=255)
                                                 cATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCa                                              >  1:5257244/1‑100 (MQ=255)
                                                  nttGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCaa                                             >  1:4427540/2‑100 (MQ=255)
                                                  aTTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCaa                                             <  2:1059931/100‑1 (MQ=255)
                                                    tGCCTCGCTCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGGGCGGATCaaaa                                           >  2:5276644/1‑100 (MQ=255)
                                                    tGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCaaac                                           >  2:2890163/1‑99 (MQ=255)
                                                    tGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCaaaa                                           <  2:3986088/100‑1 (MQ=255)
                                                     gCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCaaaaa                                          >  1:3276551/1‑100 (MQ=255)
                                                     gCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCaaaaa                                          >  1:608669/1‑100 (MQ=255)
                                                       cTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAg                                        <  1:4019984/100‑1 (MQ=255)
                                                       cTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAg                                        >  2:4445986/1‑100 (MQ=255)
                                                       cTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAACATCGTCATACAGCCCTGTCCCGACCTGGGGATCAACAAAg                                        >  2:1146299/1‑100 (MQ=255)
                                                         cGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGAc                                      >  1:2469158/1‑100 (MQ=255)
                                                           cccACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACaa                                    >  2:1584135/1‑100 (MQ=255)
                                                             cACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAg                                  >  2:4206442/1‑100 (MQ=255)
                                                              aCCTCCTGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAgg                                 >  2:2356786/1‑100 (MQ=255)
                                                                ctcctGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGt                               <  1:452721/100‑1 (MQ=255)
                                                                ctcctGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGt                               >  2:1703422/1‑100 (MQ=255)
                                                                 tcctGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGtt                              >  1:2845236/1‑100 (MQ=255)
                                                                  cctGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTc                             <  1:2878587/100‑1 (MQ=255)
                                                                   ctGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGtaga                            >  1:5345555/1‑97 (MQ=255)
                                                                   ctGAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCt                            >  1:1401600/1‑100 (MQ=255)
                                                                     gAGTGATGACGGGGGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCttt                          <  2:4729349/100‑1 (MQ=255)
                                                                     gAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCttt                          >  2:2121826/1‑100 (MQ=255)
                                                                     gAGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCttt                          >  2:2643573/1‑100 (MQ=255)
                                                                      aGTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTc                         <  2:887312/100‑1 (MQ=255)
                                                                       gTGATGACGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGATTATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGGTCTTTCg                        >  1:5370479/1‑100 (MQ=255)
                                                                         gatgaCGGGCGTCAAATGGTAACGTCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGcc                      >  2:4360607/1‑100 (MQ=255)
                                                                         gatgaCGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGcc                      <  1:2447043/100‑1 (MQ=255)
                                                                         gatgaCGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGcc                      >  2:3638785/1‑100 (MQ=255)
                                                                          atgaCGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCa                     >  2:5041601/1‑100 (MQ=255)
                                                                            gaCGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGt                   >  2:2533045/1‑100 (MQ=255)
                                                                             aCGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGtt                  <  2:2945362/100‑1 (MQ=255)
                                                                             aCGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGACCTGTGCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGtt                  <  2:1087747/100‑1 (MQ=255)
                                                                              cGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTc                 <  1:1036968/100‑1 (MQ=255)
                                                                              cGGGCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTc                 <  1:1038244/100‑1 (MQ=255)
                                                                                ggCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTcag               <  2:63991/100‑1 (MQ=255)
                                                                                 gCGTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTcagc              >  2:4726141/1‑100 (MQ=255)
                                                                                  cGTCAAACGGTAACGGCTCATCAGCCCGAGAACATCGTCATACAACCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTcagcc             <  1:5109541/100‑1 (MQ=255)
                                                                                  cGTCAAACGGTAACGGCCCATCAGCCCGAGAACATCGTCATACATCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTcagcc             <  1:5112153/100‑1 (MQ=255)
                                                                                   gTCAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTcagcca            >  1:1680755/1‑100 (MQ=255)
                                                                                     cAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTcagccagc          <  1:5527487/100‑1 (MQ=255)
                                                                                     cAAATGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAAGGTTCTTTCGCCAGTTcagccagc          >  1:1501836/1‑100 (MQ=255)
                                                                                       aaTGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTCAGCCAGCGt        >  1:335531/1‑100 (MQ=255)
                                                                                       aaTGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTCAGCCAGCGt        <  1:3280311/100‑1 (MQ=255)
                                                                                        aTGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTCAGCCAGCGtt       >  1:4046920/1‑100 (MQ=255)
                                                                                         tGGTACCGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATTAAAAAAGACAAAGGGTTCTTTCGCCAGTTCAGCCAGCGTTa      <  1:5277150/100‑1 (MQ=255)
                                                                                         tGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTCAGCCAGCGTTa      <  2:2351252/100‑1 (MQ=255)
                                                                                         tGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATAAAAAAAGACAAAGGGTTCTTTCGCCAGTTCAGCCAGCGTTa      <  2:2818153/100‑1 (MQ=255)
                                                                                         tGGTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATAAAAAAAGACAAAGGGTTCTTTCGCCAGTTCAGCCAGCGTTa      <  2:2003064/100‑1 (MQ=255)
                                                                                           gTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTCAGCCAGCGTTACa    >  2:2404213/1‑100 (MQ=255)
                                                                                           gTAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTCAGCCAGCGTTACa    >  2:1009889/1‑100 (MQ=255)
                                                                                            tAACGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTCAGCCAGCGTTACAt   >  1:1947852/1‑100 (MQ=255)
                                                                                             aaCGGCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGCGTTCTTTCGCCAGTTCAGCCAGCGCTACAtt  >  2:308886/1‑100 (MQ=255)
                                                                                                 |                                                                                               
CAGGCAAAATTGAAACACCCAGACCGGCGGAAACCAGACCGATGATGGTCATTGCCTCGCCCACCTCCTGAGTGATGACGGGCGTCAAATGGTAACGTCGCATCAGCCCGAGAATATCGTCATACAGCCCTGTCCCGACGTGCGGATCAAAAAAGACAAAGGGTTCTTTCGCCAGTTCAGCCAGCGTTACATT  >  NC_000913/1668854‑1669046

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