Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,662,533 T→A 5.4% intergenic (‑402/+50) 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,5330TA 5.4% 1021.2 / 5.7 326intergenic (‑402/+50)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 T (180/118);  new base A (4/15);  total (184/142)
Fisher's exact test for biased strand distribution p-value = 1.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 334 total aligned reads displayed.

CATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGA  >  NC_000913/2662385‑2662682
                                                                                                                                                    |                                                                                                                                                     
 aTCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCAGACCTACGGTTCTGTTCa                                                                                                                                                     >  1:4873786/1‑150 (MQ=255)
  tCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAATACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCAc                                                                                                                                                    >  1:1169228/1‑150 (MQ=255)
    cggCACCGGGGGTCGGATGCGGCGTCACCGCCTTAGCCGACCTCCGGTTGTGCTCACTGTCGGCCGGATTCGACGCATTCGGCGGCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTg                                                                                                                                                  <  1:1784200/149‑1 (MQ=14)
    aGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTg                                                                                                                                                  >  2:2935720/1‑150 (MQ=255)
     ggCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGt                                                                                                                                                 >  2:2741382/1‑150 (MQ=255)
      gCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATACGACCAACGGTTCTGTTCACTGTa                                                                                                                                                >  1:1922239/1‑150 (MQ=255)
        accGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAgg                                                                                                                                              <  2:5175726/148‑1 (MQ=255)
        aaCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCAACGGTTCTGTTCACTGTAgg                                                                                                                                              >  2:5424176/1‑150 (MQ=255)
           ggCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTCCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCt                                                                                                                                           <  2:2450362/150‑1 (MQ=255)
           ggCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCt                                                                                                                                           <  2:4485496/150‑1 (MQ=255)
           ggCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTATGTTCACTGTATGCCAGATAAGAAGCAATACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATACGACCTACGGTTCTGTTCACTGTAGGcca                                                                                                                                           >  1:5637961/1‑149 (MQ=14)
            gCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTCCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTg                                                                                                                                          <  2:5642023/150‑1 (MQ=255)
                tCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATaa                                                                                                                                      >  1:4022924/1‑150 (MQ=255)
                 cGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATAAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAg                                                                                                                                     <  2:875378/150‑1 (MQ=255)
                  ggATGCGGCGTGAACGCCTTATCCGACCTCCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGa                                                                                                                                    <  2:3917032/150‑1 (MQ=255)
                  ggATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGa                                                                                                                                    >  1:616566/1‑150 (MQ=255)
                   gATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGAc                                                                                                                                   >  2:5515815/1‑150 (MQ=255)
                   gATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGAc                                                                                                                                   >  1:4398237/1‑150 (MQ=255)
                        ggcgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCAtt                                                                                                                              >  1:2654996/1‑150 (MQ=255)
                              aaCGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGt                                                                                                                        >  1:2029732/1‑150 (MQ=255)
                              aaCGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGt                                                                                                                        <  2:4244281/150‑1 (MQ=255)
                               aCGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAATGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCAGTACGCGTc                                                                                                                       >  1:1637256/1‑150 (MQ=255)
                                  ccTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCa                                                                                                                    >  2:3640919/1‑150 (MQ=255)
                                  ccTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCa                                                                                                                    >  1:3694557/1‑150 (MQ=255)
                                  ccTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCa                                                                                                                    >  2:608804/1‑150 (MQ=255)
                                   cTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCAt                                                                                                                   >  2:1121693/1‑150 (MQ=40)
                                   cTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCAt                                                                                                                   >  1:4619164/1‑150 (MQ=40)
                                      aTCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                                                                                >  1:2904474/1‑150 (MQ=39)
                                      aTCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                                                                                >  2:2571941/1‑150 (MQ=39)
                                      aTCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                                                                                >  1:4302354/1‑150 (MQ=39)
                                      aTCCGACCTACGGTTCTGTACACTGTAGGCCTGCTAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                                                                                <  2:3033207/150‑1 (MQ=21)
                                       tCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAgg                                                                                                               >  2:4190044/1‑150 (MQ=39)
                                            acttcgtttCTGTTCACTGTAGGCCTGATAAGACTCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACg                                                                                                          <  2:5922748/143‑1 (MQ=12)
                                             cTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACgg                                                                                                         >  1:2856970/1‑150 (MQ=35)
                                              tACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGc                                                                                                        >  2:1384926/1‑150 (MQ=34)
                                                 tgTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAATGGCTGt                                                                                                     <  2:5691121/149‑1 (MQ=17)
                                                   ttCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTa                                                                                                   >  2:5477812/1‑150 (MQ=33)
                                                   ttCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTa                                                                                                   >  1:5911619/1‑150 (MQ=33)
                                                      tGTTCACTGTAGGCTTGATAAGAAGATTTTCGCGTGGCATTTGGCACCGGCTGTCGGATGCGGCGTGAACGCATTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTc                                                                                                <  2:2809402/150‑1 (MQ=11)
                                                        ttCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCACCGGCTGTTATTCgg                                                                                              >  1:1985134/1‑150 (MQ=37)
                                                        ttCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCgg                                                                                              >  2:1019391/1‑150 (MQ=37)
                                                         acACGGTAGGCCGGAGACGACGCATTACGCGTCGCCTCTGGCACCTGCTGTCGGATGCGGCGTGAACGCCTTAGCCGACCTACGGTTCTGTCCACTGTAGGCCTGATAAGACGCTGTAGGCGTCGCATCAGGCAACGGCTGTTATTCGGc                                                                                             <  2:5228923/149‑1 (MQ=11)
                                                          cACTTTATGCATTATAATACGAATTACGCGTCGCAACAGGCACCGGCTGTGTTATGCGGCGTGAACGCCTTATCCGACCTCCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGcc                                                                                            <  2:829574/150‑1 (MQ=11)
                                                            cTGTAGGGGTGATTTGAGGCTTTTCGGGTCGCTTCTGGGTACGGCTGTCGGTTGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATACGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCtt                                                                                          <  2:3872160/150‑1 (MQ=11)
                                                             tGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTg                                                                                         >  2:2257396/1‑150 (MQ=39)
                                                             tGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTg                                                                                         >  2:2075170/1‑150 (MQ=39)
                                                              gTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGt                                                                                        >  1:1449545/1‑150 (MQ=40)
                                                                  gCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTAcc                                                                                    >  1:4487423/1‑150 (MQ=255)
                                                                  gCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTAcc                                                                                    >  1:3388703/1‑150 (MQ=255)
                                                                  gCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTAcc                                                                                    >  1:2788505/1‑150 (MQ=255)
                                                                  gCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTAcc                                                                                    >  2:150492/1‑150 (MQ=255)
                                                                  gCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTAcc                                                                                    >  1:1739664/1‑150 (MQ=255)
                                                                   ccTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCt                                                                                   >  2:505618/1‑150 (MQ=255)
                                                                    cTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCtt                                                                                  >  2:4905719/1‑150 (MQ=255)
                                                                     tGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCTCTGTAGGCCTGTTAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCttt                                                                                 >  1:5513502/1‑150 (MQ=255)
                                                                       aTAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTAt                                                                               >  1:2520155/1‑150 (MQ=255)
                                                                           gACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctg                                                                           >  1:1159138/1‑150 (MQ=255)
                                                                             cGCTTTACGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctgct                                                                         <  2:4516920/150‑1 (MQ=255)
                                                                             cGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctgct                                                                         <  2:456341/150‑1 (MQ=255)
                                                                              gCATTAGGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctgctg                                                                        <  2:4324377/150‑1 (MQ=255)
                                                                              gCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctgctg                                                                        >  2:1990146/1‑150 (MQ=255)
                                                                               cATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtgctgc                                                                       >  2:4612389/1‑150 (MQ=255)
                                                                                aTTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTgctgct                                                                      >  1:4952805/1‑150 (MQ=255)
                                                                                        tCGCATCTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGaa                                                              <  2:1178003/150‑1 (MQ=255)
                                                                                        tCGCATCCGGCACCGGGTGTGGGATGGGGCGTTAACGCCTTATCCGACATAGGGTTCTTTACTTTGTAGGCCTTATAACCCGAATTAGGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGaa                                                              <  1:2456969/150‑1 (MQ=18)
                                                                                        tCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGaa                                                              >  2:4750567/1‑150 (MQ=255)
                                                                                            ctCAGGCACCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGCCCTTGGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCa                                                          <  2:5342578/149‑1 (MQ=255)
                                                                                            aTCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCa                                                          >  1:420116/1‑150 (MQ=255)
                                                                                               aGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATACGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAAGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGaa                                                       >  1:535078/1‑150 (MQ=255)
                                                                                                 gcacgggCTGTCGTTTGGGGCGTGACCGCCTTATCCGACCTCGGTTTCTTTTCAGTTTTGGCCTGATAAGACGCATTCCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATc                                                     <  2:2654996/145‑1 (MQ=21)
                                                                                                 gCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATc                                                     >  1:3164653/1‑150 (MQ=255)
                                                                                                   accGGGTGTCGGATGCGGCGTGAACGCCTTATCCGACCTCGGGTTCTGTTCACTGTTGGCCTGATAAGACGCTTTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATccc                                                   <  1:5189728/148‑1 (MQ=21)
                                                                                                        cTGTCGGATGCGGCGTGAACGCCTTATCCGACCTCCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCa                                              <  2:2536173/150‑1 (MQ=255)
                                                                                                         tggcGGATGCGGCGTGAACGCCTTACCCGACCTCCGGTTCTGTCCACTGTCGGCCTTATATCACGCATTCCGCGCCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAg                                             <  1:21434/147‑1 (MQ=21)
                                                                                                         tGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAg                                             >  2:3333002/1‑150 (MQ=255)
                                                                                                           tCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGa                                           >  2:1687842/1‑150 (MQ=255)
                                                                                                           tCGGATAAGGCGTTCACGCCGCATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGa                                           >  1:3754608/1‑150 (MQ=25)
                                                                                                           tCGGAGGCGGCGTCAACGCCTTCTCCGACCTACGGTTCTGTTCACTGTCGGCCTGATAAGACGCATTCCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGa                                           <  1:4192601/150‑1 (MQ=25)
                                                                                                            cGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGat                                          >  1:2418003/1‑150 (MQ=255)
                                                                                                              gATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGatat                                        >  2:177452/1‑150 (MQ=255)
                                                                                                                 gcggcgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTTGGGCTGTTAAGACGCATTACGCGTCGCTTCAGGCAACGGCTGTTATTCGGCGTTGTTAGGTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTct                                     >  2:3800267/1‑149 (MQ=21)
                                                                                                                     cttGAACGCTTTATCCGACCTAGGTTTCTTTTCAGTGTTGGCCTTTTATGACGATTTACGCGTCGCATCTGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTc                                 <  1:2296883/148‑1 (MQ=255)
                                                                                                                     cgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTc                                 >  1:5959745/1‑150 (MQ=255)
                                                                                                                       tGATCGGCTTATCCGACCTCCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCt                               <  2:4407156/150‑1 (MQ=255)
                                                                                                                       tGAACGCCTTAGCCGCCCTAGGGTTCTGTACACTGTTGGCCTTATAAGACGCATTAGGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCt                               <  1:3565647/150‑1 (MQ=25)
                                                                                                                       ggaccGCATTACCCGACCTTCGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAGGCCAGAATCCCGCGCAGGATATTCAACTCCt                               <  1:4122972/146‑1 (MQ=25)
                                                                                                                         accGCCTTACCCGCCCTGCGGTTCTGTACACTGTCGCCCTTATAAGACGCAGTCCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTgg                             <  1:4685825/148‑1 (MQ=37)
                                                                                                                           cGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTATGCTGCTCAATAGAATCCAGAATCCCGCGCAGGATATTCAACTCCTGGca                           >  2:1334081/1‑149 (MQ=255)
                                                                                                                            gCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCtt                          >  1:5895898/1‑150 (MQ=255)
                                                                                                                             ccTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCGGTTATGCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCCGAATCCCGCGCAGGATATTCAACTCCTGGCgtt                         >  2:4467545/1‑150 (MQ=255)
                                                                                                                                  tCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCggg                    >  2:3474739/1‑150 (MQ=255)
                                                                                                                                      acgtcgggTTCTTTTCACTGTAGGCTTGATACGACGCATTACGCGTCGCATCAGGAAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgcgcg                <  1:3880540/144‑1 (MQ=38)
                                                                                                                                        cTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCa              >  2:4055360/1‑150 (MQ=255)
                                                                                                                                         tACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCACTAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCAc             >  1:1024474/1‑150 (MQ=255)
                                                                                                                                           cGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACgg           >  1:832729/1‑150 (MQ=255)
                                                                                                                                              ttCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGtt                                                                                       <  2:726694/71‑1 (MQ=255)
                                                                                                                                                 tGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTCTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAAc     >  1:1330247/1‑150 (MQ=255)
                                                                                                                                                  gTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACa    <  1:1384926/150‑1 (MQ=255)
                                                                                                                                                   ttCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAg   >  2:5141952/1‑150 (MQ=255)
                                                                                                                                                   ttCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAg   >  2:3209196/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCACGGGTAAACAGA  >  NC_000913/2662385‑2662682

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