Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,662,523 A→G intergenic (‑392/+60) 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,5230AG98.7% 250.7 / NA 77intergenic (‑392/+60)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 (0/0);  major base G (25/51);  minor base C (1/0);  total (26/51)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGG  >  NC_000913/2662381‑2662664
                                                                                                                                               |                                                                                                                                             
gTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCt                                                                                                                                        <  2:653056/150‑1 (MQ=255)
 tCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTg                                                                                                                                       <  1:287490/150‑1 (MQ=255)
  cGCATCAGGCACCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGCCCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCt                                                                                                                                        <  1:1090418/148‑1 (MQ=255)
  cGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGtt                                                                                                                                     <  2:149776/151‑1 (MQ=255)
       aaGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGTTGCGGCGTGGACGCCTTATCCGACCTGCGGTTCTTGTCACt                                                                                                                                 >  1:1073951/2‑150 (MQ=17)
           cAACGGCTGTCGGATTCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCAGCTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAgg                                                                                                                            <  1:957201/151‑1 (MQ=255)
           cAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAgg                                                                                                                            <  2:824521/151‑1 (MQ=255)
            aaCGGCTGTCGGATGCGGCGTGCACGCCTTATCCCCCCTCCCGTTCTGTTCACTGTAGGCCTGGTAAGCCGCATTACGCGTCGCATCAGGCAACCGCTGTTGGTCCCCGCGTTTACCCCTTTTCCGCCCCGCCGTTCTGTTccccgtccg                                                                                                                            >  1:1094754/1‑142 (MQ=2)
             aCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGcc                                                                                                                          <  2:36038/151‑1 (MQ=255)
             aCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGCTTGCGGCGTGCACCCCTTTTCCGATCCTCCGGTCCGTTCAc                                                                                                                                  >  1:683730/1‑143 (MQ=14)
                   gTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATa                                                                                                                     <  2:496242/150‑1 (MQ=255)
                     cGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGa                                                                                                                  <  2:10287/151‑1 (MQ=255)
                       gATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGGTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACg                                                                                                                <  2:985976/151‑1 (MQ=255)
                             gcgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAc                                                                                                          <  1:440523/151‑1 (MQ=255)
                               gTGAACGCCTTATCCGCCCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAcgc                                                                                                        <  1:504720/151‑1 (MQ=25)
                                 gAACGCCTTATCGGGCCTGCTGTTCTATTCACTGTAGGCCTGAGTAGACGCATTACGCTTCGCGTCAGTCAACGTCTGTCTGATGCGGCGTGAACGCCTTATCCGACGTGCAGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAcgcg                                                                                                       <  1:881817/150‑1 (MQ=11)
                                 gAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAcgcg                                                                                                       <  2:364937/150‑1 (MQ=255)
                                 gAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAcgcg                                                                                                       <  2:1126386/150‑1 (MQ=255)
                                 gAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGt                                                                                                      <  2:582455/151‑1 (MQ=255)
                                  aaCGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTc                                                                                                     <  2:125704/151‑1 (MQ=255)
                                      ccTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGGTCTGTTCACAGGAGGCCCGATAAGACGCATTACGCGTcccct                                                                                                 >  1:1195533/1‑147 (MQ=14)
                                         tATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCa                                                                                               <  2:709695/150‑1 (MQ=38)
                                           tCCGAGCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAgg                                                                                             <  1:890380/150‑1 (MQ=255)
                                           tCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGCCGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGCTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                                                              <  2:233981/149‑1 (MQ=255)
                                             cGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCaa                                                                                          <  2:545116/151‑1 (MQ=37)
                                                ccTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTTTGTTCACTGTAGGGCTGATAAGACGCATTACGCGTCTCATCAGGCAACgg                                                                                       >  1:1108452/1‑151 (MQ=255)
                                                ccTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCCGCAGTTGTGGTCACAGTAGGCCTGATAAGACGCATTACCCGTCGCATCAGGc                                                                                            >  1:367431/1‑146 (MQ=11)
                                                 cTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGc                                                                                      >  2:1035849/1‑151 (MQ=21)
                                                   aCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCt                                                                                     <  1:917674/150‑1 (MQ=18)
                                                   aCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTg                                                                                    <  2:149405/151‑1 (MQ=18)
                                                      gTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTa                                                                                 <  2:218576/151‑1 (MQ=18)
                                                      gTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTa                                                                                 <  2:1067478/151‑1 (MQ=18)
                                                      gTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTa                                                                                 <  1:1171248/151‑1 (MQ=18)
                                                          tGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACCGTACGCCCGAGAAGATGCATTACCCCGCGCGTCAGGGaagcgcgtgtattc                                                                              >  1:595565/1‑138 (MQ=11)
                                                           gTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTAtt                                                                               <  2:1146023/148‑1 (MQ=21)
                                                               acggtcggCCTGATAAGACGCATTCCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAATGGCTGTTATTCGGCCtt                                                                        <  1:129220/145‑1 (MQ=255)
                                                               aCTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCtt                                                                        >  2:994076/1‑151 (MQ=37)
                                                                cTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCtt                                                                        <  1:555243/150‑1 (MQ=38)
                                                                cTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGAGAAGACGCATTACGCGTCGCCTCAGGCAACGGCTGTTATTCGGCCtt                                                                        >  1:318448/1‑150 (MQ=255)
                                                                  gTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGtt                                                                     <  2:1171600/151‑1 (MQ=255)
                                                                  gTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGt                                                                      <  2:232986/150‑1 (MQ=255)
                                                                  gTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGt                                                                      <  2:838793/150‑1 (MQ=255)
                                                                  gTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGt                                                                      <  2:412356/150‑1 (MQ=255)
                                                                          gATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTAt                                                             <  1:671754/151‑1 (MQ=255)
                                                                            tAAGACGCATTGCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGGCCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTAtt                                                            <  1:22449/150‑1 (MQ=255)
                                                                            tAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTAtt                                                            >  2:644542/1‑150 (MQ=255)
                                                                                      tACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCa                                                  <  2:212817/150‑1 (MQ=255)
                                                                                          cgTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAAt                                                >  2:430681/1‑148 (MQ=255)
                                                                                           gTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGaa                                            <  2:1182901/151‑1 (MQ=255)
                                                                                             cGCATCAGGCAACGTCTGTCGGATTCGGCGTGAACGCCTTATCCTGCCTGCGGTTCTGTTCACTGTTGGCCTGATAATACGCATTTCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGa                                             <  1:35531/148‑1 (MQ=18)
                                                                                             cGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGc                                          <  2:1144080/151‑1 (MQ=255)
                                                                                                   aGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGGAGGCCGGATAAGGCGCATTTCGCGTCGCATCAGGGAACGGCGGTTATTCGGTCTTGTTACCTTTTTTTTGCTGCGCGATAGAAGCCAGaa                                     >  1:1293511/1‑150 (MQ=17)
                                                                                                    ggCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGGTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGaa                                     >  2:343822/1‑149 (MQ=255)
                                                                                                     gCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATcc                                  <  1:755756/151‑1 (MQ=255)
                                                                                                       aaCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATccc                                 >  1:332317/1‑150 (MQ=255)
                                                                                                       aaCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCcg                                <  1:650984/151‑1 (MQ=255)
                                                                                                        aCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCCATAGAAGCCAGAATCCcg                                >  2:753628/1‑150 (MQ=255)
                                                                                                         cGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCcgcg                              <  1:577556/151‑1 (MQ=255)
                                                                                                         cGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCcgcg                              <  1:1037817/151‑1 (MQ=255)
                                                                                                         cGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCcgc                               <  1:320481/150‑1 (MQ=255)
                                                                                                          ggCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCgcgc                             >  2:1082709/1‑151 (MQ=255)
                                                                                                                  gATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATAtt                     <  1:865649/151‑1 (MQ=255)
                                                                                                                     gcggcgTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCaa                  <  2:142228/151‑1 (MQ=255)
                                                                                                                     gcggcgTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCa                   >  2:247097/1‑150 (MQ=255)
                                                                                                                     gcggcgTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCa                   <  2:750806/150‑1 (MQ=255)
                                                                                                                      cggcgTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCaa                  >  1:828910/1‑150 (MQ=255)
                                                                                                                      cggcgTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAAc                 >  1:868968/1‑151 (MQ=255)
                                                                                                                         cgTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTcc              >  2:1091120/1‑151 (MQ=255)
                                                                                                                          gTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCt             >  2:1199385/1‑151 (MQ=255)
                                                                                                                               cGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCtt        >  2:201631/1‑151 (MQ=255)
                                                                                                                                gCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCttt       <  1:288391/151‑1 (MQ=255)
                                                                                                                                  cTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTcc     >  2:772592/1‑151 (MQ=255)
                                                                                                                                    tATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTcc     <  2:696266/149‑1 (MQ=255)
                                                                                                                                    tATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCg    <  1:592623/150‑1 (MQ=255)
                                                                                                                                     aTCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCggg  >  1:239078/1‑151 (MQ=255)
                                                                                                                                     aTCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCggg  >  1:646684/1‑151 (MQ=255)
                                                                                                                                     aTCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCgg   <  2:713416/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                             
GTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGG  >  NC_000913/2662381‑2662664

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: