Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,662,480 A→T 7.1% intergenic (‑349/+103) iscR ← / ← trmJ isc operon transcriptional repressor; suf operon transcriptional activator; oxidative stress‑ and iron starvation‑inducible; autorepressor/tRNA mC32,mU32 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,662,4800AT 7.1% 679.1 / 3.9 241intergenic (‑349/+103)iscR/trmJisc operon transcriptional repressor; suf operon transcriptional activator; oxidative stress‑ and iron starvation‑inducible; autorepressor/tRNA mC32,mU32 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (139/71);  new base T (2/23);  total (142/99)
Fisher's exact test for biased strand distribution p-value = 1.57e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 245 total aligned reads displayed.

TATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAA  >  NC_000913/2662331‑2662629
                                                                                                                                                     |                                                                                                                                                     
tatCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATACGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCa                                                                                                                                                       <  1:1520577/150‑1 (MQ=255)
 ctCCGCCCTGGGGTTCTGTACACTGTCGGCCTGATAAGCCGCATTACGCGTCGCATCAGCCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                                                                                                                      <  1:1074618/149‑1 (MQ=255)
 atCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACCGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                                                                                                                      <  1:3997341/150‑1 (MQ=255)
  tCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAgg                                                                                                                                                     >  2:743899/1‑150 (MQ=39)
   ccGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGc                                                                                                                                                    >  1:966121/1‑150 (MQ=39)
   ccGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGc                                                                                                                                                    >  2:2415987/1‑150 (MQ=39)
    cGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCa                                                                                                                                                   >  2:1427527/1‑150 (MQ=38)
    cGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCa                                                                                                                                                   >  1:542906/1‑150 (MQ=38)
     gCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCaa                                                                                                                                                  <  2:791259/150‑1 (MQ=38)
      cccTGGGGTTCTGTACACTGTCGGCCTGATAAGACGAATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAAc                                                                                                                                                 <  1:815321/150‑1 (MQ=255)
      cccTGGGGTTCTGTACACTGTAGGCCTGATAAGACTCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAAc                                                                                                                                                 <  2:2773848/150‑1 (MQ=255)
       ccTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACg                                                                                                                                                >  1:775367/1‑150 (MQ=35)
       ccTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACg                                                                                                                                                >  1:72291/1‑150 (MQ=35)
             gTTCTGTACACTGTAGGCCTGATATGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTc                                                                                                                                          <  2:2158067/150‑1 (MQ=17)
                cTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCCGAAGCGCCGGGAACGCCTTATCCGACCAACGGTTCTGGCCACAGTAGGCCTCATAATAAGCATGACGAGTCGCATCTGGCAACGGCTGTccta                                                                                                                                       >  1:1986598/1‑147 (MQ=11)
                  gTACACTGTTGGCTTGATATGACGCTTTTGGGGTCGCATCTGGCAACGGGTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTTTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCTGGCAACGGCTGTCGGATg                                                                                                                                     <  2:1840218/150‑1 (MQ=11)
                   tACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATgc                                                                                                                                    >  1:1077299/1‑150 (MQ=37)
                   tACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTAGGATgc                                                                                                                                    >  1:3145816/1‑150 (MQ=21)
                    acacTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATgcg                                                                                                                                   >  2:1917146/1‑150 (MQ=37)
                    acacTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGGCGGAGgcg                                                                                                                                   >  2:3959676/1‑150 (MQ=17)
                        tGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGAGGCCGCGt                                                                                                                               >  2:2899527/1‑150 (MQ=21)
                             gCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACg                                                                                                                          >  2:3398131/1‑150 (MQ=255)
                               cTGATACGACGCATTCCGCGGCGCCGCCGGCACCGGCTGTCGGCTGCGGCGTGCCCGCCTGATCCGCCCTCCGGTTCTGTTCCCTGTAGGCCTGATACGACGCATTACGCGTCGCCTCCGGCAACGGCTGTCGGATGCGGCGTGAACGcc                                                                                                                        <  2:3964531/150‑1 (MQ=14)
                                       aCGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGGTCACTGTAGGCCTGATAAGACGCATTACGCGTCTCATCAGGAAACGGCAGTAGGATGCGGAGGGGACGCCTTATCCGa                                                                                                                >  2:125649/1‑150 (MQ=14)
                                            ttACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTAc                                                                                                           >  1:2220244/1‑150 (MQ=255)
                                                gggTCGGATCATGCACCGGCTGCCGGATGCGGCGTGAACGCCTTATCCGACCTCCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGtt                                                                                                       <  2:281712/148‑1 (MQ=14)
                                                gcgTCGAATCTGGCAAGGGCTGTCGGTTGGGGCGTGAACGCCTTACCCGCCCTAGGGTTCTGTTCACTGTAGCCCTGATAAGACGAATTACGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGtt                                                                                                       <  1:3925905/150‑1 (MQ=14)
                                                     gCATCTTGCAAGGGCTGTCGGTTGCGGGGTGACCGCATTATCCGACCTTGGTTTCTGTTCACTGTAGGCCTTATACGACTCTTTACGGGCCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGtt                                                                                                  <  1:2405000/150‑1 (MQ=14)
                                                     gCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGtt                                                                                                  <  2:2280982/150‑1 (MQ=255)
                                                     gCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGtt                                                                                                  <  2:2974123/150‑1 (MQ=255)
                                                        tCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGCCGTGAACGCCTTATCCGACCTACGGTTCTGTTCAc                                                                                               >  2:3849367/1‑150 (MQ=255)
                                                           ggCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGt                                                                                            >  1:118709/1‑150 (MQ=255)
                                                            gCAACGGGTGGCGGATGCGGCGTGAACGCCTTATCCGACCTCCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTa                                                                                           <  2:1217244/150‑1 (MQ=255)
                                                            gCAACGGCTGTCGGATGCGGCGTGAACGCATTATCCTACCTACGGTTCTGTTCAATTTAGGCCTTATAAGACGCTTTACGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTa                                                                                           <  1:3239669/150‑1 (MQ=17)
                                                            gCAACGGCGGTCGGATGCGGCGTGAACGCATTACCCGACCTAGGTTTCTGTTCACTGTAGGCCTGTTAAGACGCATTACGCGTCGCATATGGAAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTa                                                                                           <  1:604194/150‑1 (MQ=14)
                                                             cAACGGCTGTCGGTTGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAg                                                                                          <  1:2709666/150‑1 (MQ=255)
                                                                cGGCTGTCGGATGCGGCGTGAACGCCTTCTCCGACCTCGGGTTCTGTACACTGTAGGCCTGATACGCCGCATTACGCGTCGCAGCTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGcc                                                                                       <  2:3194778/150‑1 (MQ=11)
                                                                cGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGcc                                                                                       <  2:459351/150‑1 (MQ=255)
                                                                 ggCTGTCGGATGCGGCGTGAACGCCATATACGACCAACCGTTATGTTCACTGTAGGCCAGATAAGACGCAATAAGCGTCGCATCAGGCAACGGCAGTCGGATGCGGCGTGAAAGCCTTAAACGACCAACGGTTCTGTTCACTGTAGGccg                                                                                      >  1:3467356/1‑149 (MQ=14)
                                                                   cTGTGGTATGCGGCGTGAACGCCTTATCCGACCTAGGGTTCTGTTCACTGTTGGACTTATATGACGCATTAGGCGTCGAATCAGGAAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGa                                                                                    <  1:2071213/150‑1 (MQ=11)
                                                                    tGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGAt                                                                                   >  2:2219504/1‑150 (MQ=255)
                                                                    tGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGAt                                                                                   >  1:632008/1‑150 (MQ=255)
                                                                    tGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCGACGGTTCTGTTCACTGTAGGCCTGAt                                                                                   >  2:1656133/1‑150 (MQ=255)
                                                                      tCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTAGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATaa                                                                                 >  2:3639162/1‑150 (MQ=255)
                                                                        ggATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGa                                                                               >  1:3833543/1‑150 (MQ=255)
                                                                              ggcgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCAtt                                                                         >  2:3912851/1‑150 (MQ=255)
                                                                              ggcgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCAtt                                                                         >  2:1868217/1‑150 (MQ=255)
                                                                              ggcgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCAtt                                                                         >  2:3378891/1‑150 (MQ=255)
                                                                               gcgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTa                                                                        >  2:2046003/1‑150 (MQ=255)
                                                                                  tGAACGCCTTATCCGACCTCCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAcgc                                                                     <  1:2753589/150‑1 (MQ=255)
                                                                                  tGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAcgc                                                                     >  1:2462302/1‑150 (MQ=255)
                                                                                  tGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGGTCACTGTAGGCCTGATAAGACGCATTAcgc                                                                     >  2:2266002/1‑150 (MQ=255)
                                                                                     aCGCCTAATCCGACCTCCGGGTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTc                                                                  <  2:2423593/150‑1 (MQ=255)
                                                                                      cGCCTTATCCGACCTTCGGTTCTGTTCACTTTAGGCATGATTAGACGCATTACGCGTCGCATAAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCg                                                                 <  1:1330608/150‑1 (MQ=14)
                                                                                      cGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCg                                                                 >  1:2158234/1‑150 (MQ=255)
                                                                                      cGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCg                                                                 >  1:1470209/1‑150 (MQ=255)
                                                                                        ccTTATCCGACCTCGGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATGCGACGCGTAATGCGTctta                                                               <  2:407362/150‑4 (MQ=11)
                                                                                         cTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGGTTGTACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAATACTCATTTCGCGTCGctt                                                              >  2:2611983/1‑148 (MQ=14)
                                                                                          ttATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATACGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATc                                                             <  1:2973060/150‑1 (MQ=255)
                                                                                           tATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCa                                                            <  2:3537038/150‑1 (MQ=39)
                                                                                            aTCCGACCTTCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                           <  1:3701239/150‑1 (MQ=37)
                                                                                            aTCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                           >  1:3518788/1‑150 (MQ=39)
                                                                                              ccGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGc                                                         >  1:2806703/1‑150 (MQ=39)
                                                                                                 aCCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAAc                                                      >  1:3960274/1‑150 (MQ=37)
                                                                                                   cTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACgg                                                    >  1:3699054/1‑150 (MQ=35)
                                                                                                       ggTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGt                                                >  1:3410789/1‑150 (MQ=32)
                                                                                                           cTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTAtt                                            >  1:2863314/1‑150 (MQ=35)
                                                                                                            tggccACTGTCGGCCTTATAACACGAATAACGCGTCGCAACTGGCAACGGCTGTCGGATGCGGCGTTAACGCCTTATCCGCCCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTCCGCGTCGCATCAGCCAACGGCTGTTATTc                                           <  2:1103080/146‑1 (MQ=11)
                                                                                                             gtacACTGTCGGCCTGCCAAGACGCATTACGGGTCGCATCCGGCAACGGCTGTCGGAGGCGGCGCGAACGCCTTATCCGACCTCCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCg                                          <  2:1379230/147‑1 (MQ=11)
                                                                                                              ttCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGACCGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCgg                                         >  1:82321/1‑150 (MQ=37)
                                                                                                              ttCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCgg                                         >  1:1680594/1‑150 (MQ=37)
                                                                                                              ttCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCgg                                         >  1:2147615/1‑150 (MQ=37)
                                                                                                              ttCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCgg                                         >  1:2696724/1‑150 (MQ=37)
                                                                                                                         ccTGATACGACGATTTACGGGTGGCCTCTTGAAACGGGTGTCGGATGCGGCGTGAACGCCTTAACCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCt                              <  2:855636/150‑1 (MQ=14)
                                                                                                                          cTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCtt                             >  2:2889884/1‑150 (MQ=255)
                                                                                                                          cTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCtt                             <  2:542906/150‑1 (MQ=255)
                                                                                                                          cTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCtt                             >  1:3500996/1‑150 (MQ=255)
                                                                                                                          cTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCtt                             >  1:1721770/1‑150 (MQ=255)
                                                                                                                           tGATACGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTAGGGTTCTGTTCACTGTAGGCCTTATAAGACGAATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCttt                            <  2:143716/150‑1 (MQ=255)
                                                                                                                            gATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTa                           <  2:3086419/150‑1 (MQ=255)
                                                                                                                            gATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTa                           >  1:3109471/1‑150 (MQ=255)
                                                                                                                            gATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTa                           >  1:2021901/1‑150 (MQ=255)
                                                                                                                            gATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTa                           >  2:3769234/1‑150 (MQ=255)
                                                                                                                                tcaCGCTTTACGCGTCGCATTTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTTGGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTct                       <  2:1158681/148‑1 (MQ=17)
                                                                                                                                  aCGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctgc                     >  1:1274707/1‑150 (MQ=255)
                                                                                                                                     cagtACGCGTCGCATCAGGCAACGGCTGTGGGATGCGGCGTGACCGCCTTCTCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtgctgc                  <  2:3031766/147‑1 (MQ=255)
                                                                                                                                      aTTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTgctgct                 >  2:3263671/1‑150 (MQ=255)
                                                                                                                                       ttACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTc                >  1:2566816/1‑150 (MQ=255)
                                                                                                                                       ttACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTc                >  2:679639/1‑150 (MQ=255)
                                                                                                                                       ttACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTc                >  1:566966/1‑150 (MQ=255)
                                                                                                                                           gcttCGAATAAGGAAAGGGCTGTGGTAGGCGGGGTGAACCCCTTATCCTACCTCCGTTTAGGTTCACTGCCGGCCTGATAAGCCTCTTTCCGCGCCGAATAAGCCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATa            <  1:3780258/147‑1 (MQ=255)
                                                                                                                                           gcgTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATa            <  2:632008/150‑1 (MQ=255)
                                                                                                                                             gTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGa          >  1:220791/1‑150 (MQ=255)
                                                                                                                                              tCGCATCTGGCATCGGTTGTGGGATGCGGCGTTAACGCATTATCCGACCTCGGGTTCTTTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGaa         <  1:3030207/150‑1 (MQ=18)
                                                                                                                                              tCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGaa         >  1:3076239/1‑150 (MQ=255)
                                                                                                                                              tCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGaa         >  1:870811/1‑150 (MQ=255)
                                                                                                                                                  aTCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTCCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCa     <  2:3373364/150‑1 (MQ=255)
                                                                                                                                                  aTCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCa     >  1:1897155/1‑150 (MQ=255)
                                                                                                                                                   tCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAg    >  1:2269714/1‑150 (MQ=255)
                                                                                                                                                     aGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTCCGGTTCTGTCCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGaa  <  2:3674050/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAA  >  NC_000913/2662331‑2662629

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