Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,662,468 A→G 7.6% intergenic (‑337/+115) 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,4680AG 7.6% 351.9 / 16.1 132intergenic (‑337/+115)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 (67/53);  new base G (5/5);  total (72/60)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 135 total aligned reads displayed.

GCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCT  >  NC_000913/2662319‑2662614
                                                                                                                                                     |                                                                                                                                                  
gcgTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGTACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTa                                                                                                                                                    >  2:219695/1‑150 (MQ=255)
gcgTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTa                                                                                                                                                    >  1:1192762/1‑150 (MQ=255)
gcgTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTa                                                                                                                                                    >  1:1196656/1‑150 (MQ=255)
        gCCATTTCCGCCCTGGGGTTCTGTACAATGTAGGCCTGATACGACGCATGACGGGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATTAGACGCATTACGCGTCGc                                                                                                                                            <  2:2077924/150‑1 (MQ=14)
        gCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGc                                                                                                                                            >  1:1027655/1‑150 (MQ=255)
        gCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGc                                                                                                                                            >  1:2279267/1‑150 (MQ=255)
        gCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCAGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGc                                                                                                                                            >  1:2279175/1‑150 (MQ=255)
        gCCATATCCGCCCCGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGc                                                                                                                                            <  2:2283732/150‑1 (MQ=255)
            tctCCGCCGTGGGGTTTTTTACTCTGTAGTCCTGTTTAGACGCATTAGGGGTCGAATGTGGCAACGGCTGTCGGTTGCGGCGTTAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCa                                                                                                                                        <  1:937010/148‑1 (MQ=11)
            tatCCGCCCTGTTTTTCTGTACAATTTAGCCCTGATATGACTATTTACGCGTCGCATAAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCa                                                                                                                                        <  2:1433617/150‑1 (MQ=11)
            tatCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCa                                                                                                                                        <  1:876076/150‑1 (MQ=39)
            tatCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCa                                                                                                                                        <  2:575051/150‑1 (MQ=39)
             atCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                                                                                                       <  2:1377257/150‑1 (MQ=39)
              tCCGCCCTGGGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAgg                                                                                                                                      <  1:191504/150‑1 (MQ=21)
              tCCGCCCTCGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAgg                                                                                                                                      <  1:687707/150‑1 (MQ=255)
               ccGCCCTGGGGTTGTTTACACTGTTGGCCTGATAAGCGGCATTACGGGTGGCATCCGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATGCGACGCGTAATGCGTCTTATCAGGc                                                                                                                                     <  2:1215930/150‑1 (MQ=14)
               ccGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGc                                                                                                                                     >  1:45774/1‑150 (MQ=39)
                 gCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGAATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCaa                                                                                                                                   <  1:2639449/150‑1 (MQ=38)
                 gCCCGGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCaa                                                                                                                                   <  1:869798/150‑1 (MQ=25)
                         tttgtGTACACTGTAGGCCTGATAAGACGCTTTTCGCGTCGCAGCTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTc                                                                                                                           <  2:348346/146‑1 (MQ=11)
                         gTTCTGTCCCCGGTCGGCTTTATAGGACGATTTACGGGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCTTTACGCGTCGCATCTGGCAACGGCTGTc                                                                                                                           <  2:1288958/150‑1 (MQ=11)
                         gTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTc                                                                                                                           <  2:2405141/150‑1 (MQ=32)
                         gTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTc                                                                                                                           <  2:348534/150‑1 (MQ=32)
                          ttCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGttg                                                                                                                          >  1:2069084/1‑148 (MQ=14)
                          ttCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCg                                                                                                                          >  1:156629/1‑150 (MQ=33)
                               tACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGGTgc                                                                                                                     >  1:2090316/1‑150 (MQ=21)
                               tACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATgc                                                                                                                     >  1:2148596/1‑150 (MQ=37)
                               tACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATgc                                                                                                                     >  1:2150178/1‑150 (MQ=37)
                               tACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATgc                                                                                                                     >  1:1211648/1‑150 (MQ=37)
                               tACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTAATCAGACCTACGGTTCTGTACACTGTAGGCCTGATAAGACGCAATACGCGTCGCATCAGGCAACGGCTGtgatacgc                                                                                                                     >  1:2148597/1‑143 (MQ=11)
                                 cacTGTCGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTCATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATACGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATgcgg                                                                                                                   <  2:96983/150‑1 (MQ=255)
                                 cacTGTAGGCCTGATATGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATgcgg                                                                                                                   <  2:928895/150‑1 (MQ=37)
                                 cacTGTAGGCCTGATAAGACGCTTTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATgcgg                                                                                                                   <  2:1619879/150‑1 (MQ=37)
                                    tGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGATTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGt                                                                                                                >  1:563416/1‑150 (MQ=255)
                                         gCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACg                                                                                                           >  1:1699802/1‑150 (MQ=255)
                                          ccTGATGCGACGCGTAATGCGTCTTATCAGGCAACGGCTTTCGGATGCGGCGTGAACGCCATATGCGACCTACGGTTCTGTTCACTGTAGGCCAGATAACACGCATTACGCGTCGCAACACGCAACCGCTGTCCGATGCGGCGGGAAccc                                                                                                          >  2:2138347/1‑148 (MQ=11)
                                          ccTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGc                                                                                                          >  2:1456170/1‑150 (MQ=255)
                                              aTACGACGATTTACGCGTCGCATCTGGAAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTa                                                                                                      <  2:740935/150‑1 (MQ=14)
                                              aTAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTa                                                                                                      <  2:740872/150‑1 (MQ=255)
                                                aacactaaTTACCCGTCGAATAAGGAAACGGATGTCGGATGCGGCTTCAACGCCTTATCCGACCTACGTTTCTGTACACTGTAGGCCTGATAAGACGAATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATc                                                                                                    <  2:2368174/143‑1 (MQ=14)
                                                aaGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATc                                                                                                    <  2:2165619/150‑1 (MQ=255)
                                                aaGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATc                                                                                                    <  2:201077/150‑1 (MQ=255)
                                                 aGACGCAATACGCGTCGCATCAGGCAACGGCAGTCGGATGCGGCGTGAACGACTTAACCGACCTACGGTTCTGATCACAGAAGGCCTGATAAGACGCAATACGCGACCCAGCAGGCAACCGCTGTCGGAAGCGGCGTGAACGccagaacc                                                                                                   >  2:2212094/1‑144 (MQ=14)
                                                   aCGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGa                                                                                                 >  1:1754284/1‑150 (MQ=255)
                                                   aCGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCAGa                                                                                                 >  1:1526986/1‑150 (MQ=255)
                                                       gtaatgcgTCTTATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTGCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTa                                                                                             >  1:2637046/6‑150 (MQ=255)
                                                       gtaatgcgTCTTATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTGCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTa                                                                                             >  1:1940985/6‑150 (MQ=255)
                                                       gtaatgcgTCTTATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTGCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTa                                                                                             >  1:1757204/6‑150 (MQ=255)
                                                       gtaatgcgTCTTATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTGCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTa                                                                                             >  1:522161/6‑150 (MQ=255)
                                                            gcgTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGtt                                                                                        >  1:2479060/1‑150 (MQ=255)
                                                               tCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTg                                                                                     <  2:617715/150‑1 (MQ=255)
                                                                  cATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTc                                                                                  >  1:1076091/1‑150 (MQ=255)
                                                                  cATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTc                                                                                  >  1:818707/1‑150 (MQ=255)
                                                                  cATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTAGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTc                                                                                  >  1:1822210/1‑150 (MQ=255)
                                                                  cATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGGTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTc                                                                                  >  1:2624567/1‑150 (MQ=255)
                                                                      aGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGAGGCGGCGGGTACGCCTGATCCGACCGGCGGTTCTgggcagtg                                                                              >  2:1204291/1‑143 (MQ=14)
                                                                        gtaaCTGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTCCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTa                                                                            <  2:2472388/148‑1 (MQ=255)
                                                                        gCAAGGGGTGTCGGATGCGGCGTCAACGCATTATCCGACCTCGGGTTCTGTTCACTGTAGGCCTTATAAGACGCTTTACGCGTCGCACCTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTa                                                                            <  1:2370992/150‑1 (MQ=11)
                                                                        gCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCAGTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTGCTGTTCACTGtc                                                                            >  2:2637244/1‑149 (MQ=255)
                                                                        gCAACGGCTGTCGGATGCGGCGTGAACGCCATATCCGACCTACGGTTATGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTAGGATGCGGCGTGAACGCCTTATCCGACCAACGGTTCTGTTCACTGTa                                                                            >  1:470621/1‑150 (MQ=255)
                                                                         caccGGGTGTCGGATGCGGCGTGAACGCCTTATCCGACCTCCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGGGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAg                                                                           <  2:2696435/147‑1 (MQ=255)
                                                                           ccGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTCGGTTTCTGTTCACTGTAGGCCTGATAAGACGCTTTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGc                                                                         <  2:1187321/149‑1 (MQ=255)
                                                                                    ggctGCGGCGTGAACGCCTTTTCCGCCCTCAGTTTCTGTTCACTGTACGCCTGATATGACGCATTACGCGTCGCATCAGGGAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGa                                                                <  2:427204/147‑1 (MQ=11)
                                                                                    ggATGCGGCGTGAACGCCTTATCCGACCTTCGGTTTGGTTCACTGTTGGCCTGATACGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGa                                                                <  1:1860370/150‑1 (MQ=14)
                                                                                     gATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGAc                                                               >  2:2120555/1‑150 (MQ=255)
                                                                                          ggcgTGAACGCCTTTTCCGACCTGCGGTTCTGTGCACTGTAGGCCTGATATGACGAATTACGCGTCGGATCATGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCAtt                                                          <  2:1070029/150‑1 (MQ=11)
                                                                                            cgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAc                                                        >  1:2225905/1‑150 (MQ=255)
                                                                                                   gCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGc                                                 >  1:1393324/1‑150 (MQ=255)
                                                                                                   gCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGc                                                 >  2:2237839/1‑150 (MQ=255)
                                                                                                     cTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCAt                                               >  1:877806/1‑150 (MQ=40)
                                                                                                     cTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCAt                                               >  1:1333939/1‑150 (MQ=40)
                                                                                                     cTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCAt                                               >  1:374929/1‑150 (MQ=40)
                                                                                                     cTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCAt                                               >  1:374959/1‑150 (MQ=40)
                                                                                                      ttATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATc                                              >  1:2449884/1‑150 (MQ=40)
                                                                                                      ttATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATc                                              >  1:1119930/1‑150 (MQ=40)
                                                                                                      ttATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTAATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATc                                              >  1:2449853/1‑150 (MQ=255)
                                                                                                       tATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTGCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCa                                             <  2:522161/150‑1 (MQ=255)
                                                                                                       tATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTGCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCa                                             <  2:1940985/150‑1 (MQ=255)
                                                                                                         tCCGACCTACGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAgg                                           <  1:2320105/150‑1 (MQ=21)
                                                                                                            gACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCaa                                        >  2:2722614/1‑150 (MQ=38)
                                                                                                            gACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTAATCCGACCTACGGTTCTGTTCACAGTAGGCCTGATAAGACGCATAACGCGTCGCAGCAGGCac                                        >  2:226857/1‑149 (MQ=255)
                                                                                                             aCCTACGATTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAAc                                       <  2:563416/150‑1 (MQ=25)
                                                                                                                tACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGc                                    >  1:2168415/1‑150 (MQ=34)
                                                                                                                tACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGc                                    >  1:2168430/1‑150 (MQ=34)
                                                                                                                tACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGc                                    >  1:2168507/1‑150 (MQ=34)
                                                                                                                        tGTTCACTTTTGGCCTGATAATACGCATTACGCTTCGAATCTGTCAAGGGCTGTCGGATGCGGCGTTAACGCCTTATCCGACCTACGTTTCTTTTCACTTTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTc                            <  2:45774/150‑1 (MQ=11)
                                                                                                                          ttCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCgg                          >  2:261423/1‑150 (MQ=37)
                                                                                                                          ttCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCgg                          >  2:1742977/1‑150 (MQ=37)
                                                                                                                              cTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCtt                      >  1:698599/1‑150 (MQ=38)
                                                                                                                               tGTTGGCCTGATAAGACGAATTCCGCGTCGCATCTGGCACCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTg                     <  2:811057/150‑1 (MQ=14)
                                                                                                                                gttggCCTGAAAAGACGATTTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTCAACGCCTTATCCGACCTAGGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGt                    <  2:2315626/147‑1 (MQ=14)
                                                                                                                                gTAGGCCTGATAAGCCGCATTACGCGTCGCCTCAGGCAACGGCGGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGt                    <  2:406718/150‑1 (MQ=255)
                                                                                                                                   ggCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTAc                 <  2:1962458/150‑1 (MQ=255)
                                                                                                                                   ggCCGGATAAGACGCTCTCCGCGCCCCCTCTGGCAACGGCGGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTACACGGTAGGCCTGATAAGACGCATGACGCGTCGCATCCGGCAACGGCTGTTATTCGGCCTTGTTAc                 <  2:2090316/150‑1 (MQ=14)
                                                                                                                                    gCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTAcc                >  2:2443509/1‑150 (MQ=255)
                                                                                                                                    gCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGAGAAGACGCATTACGCGTCGCATCAGGCCACGGCTGTTCTTCGGCCTTGGTAcc                >  2:2017365/1‑150 (MQ=255)
                                                                                                                                          tAAGACGCATTACGCGGCGGTTCTGGCACGGGCTGTCGGAGGCGGCGTGAACGCCTTTTCCGAGCTACGGTTGTGTTCACTGTAGGCCTGATAAGACGCTTTAGGCGTCGGATCTGGCAAGGGCTGTTATTCGGCCTTGTTACCTTTAtt          <  2:1654696/150‑1 (MQ=17)
                                                                                                                                              aCGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctgc      >  1:1483984/1‑150 (MQ=255)
                                                                                                                                              aCGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctgc      >  2:1405312/1‑150 (MQ=255)
                                                                                                                                                  aTTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTgctgct  >  1:2319339/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCT  >  NC_000913/2662319‑2662614

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: