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

GTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGG  >  NC_000913/2662381‑2662664
                                                                                                                                               |                                                                                                                                             
gTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCt                                                                                                                                        <  2:653056/150‑1 (MQ=255)
 tCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTg                                                                                                                                       <  1:287490/150‑1 (MQ=255)
  cGCATCAGGCACCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGCCCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCt                                                                                                                                        <  1:1090418/148‑1 (MQ=255)
  cGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGtt                                                                                                                                     <  2:149776/151‑1 (MQ=255)
       aaGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGTTGCGGCGTGGACGCCTTATCCGA‑CCTGCGGTTCTTGTCACt                                                                                                                                 >  1:1073951/2‑150 (MQ=17)
           cAACGGCTGTCGGATTCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCAGCTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAgg                                                                                                                            <  1:957201/151‑1 (MQ=255)
           cAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAgg                                                                                                                            <  2:824521/151‑1 (MQ=255)
            aaCGGCTGTCGGATGCGGCGTGCACGCCTTATCCCCCCTCCCGTTCTGTTCACTGTAGGCCTGGTAAGCCGCATTACGCGTCGCATCAGGCAACCGCTGTTGGTCCCCGCGTTTACCCCTTTTCCGC‑CCCGCCGTTCTGTTccccgtccg                                                                                                                            >  1:1094754/1‑142 (MQ=2)
             aCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGcc                                                                                                                          <  2:36038/151‑1 (MQ=255)
             aCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGCTTGCGGCGTGCACCCCTTTTCCGATCCTCCGG‑TCCGTTCAc                                                                                                                                  >  1:683730/1‑143 (MQ=14)
                   gTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATa                                                                                                                     <  2:496242/150‑1 (MQ=255)
                     cGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGa                                                                                                                  <  2:10287/151‑1 (MQ=255)
                       gATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGGTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACg                                                                                                                <  2:985976/151‑1 (MQ=255)
                             gcgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAc                                                                                                          <  1:440523/151‑1 (MQ=255)
                               gTGAACGCCTTATCCGCCCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAcgc                                                                                                        <  1:504720/151‑1 (MQ=25)
                                 gAACGCCTTATCGGGCCTGCTGTTCTATTCACTGTAGGCCTGAGTAGACGCATTACGCTTCGCGTCAGTCAACGTCTGTCTGATGCGGCGTGAACGCCTTATCCGA‑CGTGCAGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAcgcg                                                                                                       <  1:881817/150‑1 (MQ=11)
                                 gAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAcgcg                                                                                                       <  2:364937/150‑1 (MQ=255)
                                 gAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAcgcg                                                                                                       <  2:1126386/150‑1 (MQ=255)
                                 gAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGt                                                                                                      <  2:582455/151‑1 (MQ=255)
                                  aaCGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTc                                                                                                     <  2:125704/151‑1 (MQ=255)
                                      ccTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGGTCTGTTCACAGGAGGCCCGATAAGACGCATTACGCGTcccct                                                                                                 >  1:1195533/1‑147 (MQ=14)
                                         tATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCa                                                                                               <  2:709695/150‑1 (MQ=38)
                                           tCCGAGCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAgg                                                                                             <  1:890380/150‑1 (MQ=255)
                                           tCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGCCGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGCTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                                                              <  2:233981/149‑1 (MQ=255)
                                             cGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCaa                                                                                          <  2:545116/151‑1 (MQ=37)
                                                ccTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTTTGTTCACTGTAGGGCTGATAAGACGCATTACGCGTCTCATCAGGCAACgg                                                                                       >  1:1108452/1‑151 (MQ=255)
                                                ccTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCCGCAGTTGTGGTCACAGTAGGCCTGATAAGACGCATTACCCGTCGCATCAGGc                                                                                            >  1:367431/1‑146 (MQ=11)
                                                 cTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGc                                                                                      >  2:1035849/1‑151 (MQ=21)
                                                   aCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCt                                                                                     <  1:917674/150‑1 (MQ=18)
                                                   aCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTg                                                                                    <  2:149405/151‑1 (MQ=18)
                                                      gTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTa                                                                                 <  2:218576/151‑1 (MQ=18)
                                                      gTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTa                                                                                 <  2:1067478/151‑1 (MQ=18)
                                                      gTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTa                                                                                 <  1:1171248/151‑1 (MQ=18)
                                                          tGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACCGTACGCCCGAGAAGATGCATTACCCCGCGCGTCAGGGaagcgcgtgtattc                                                                              >  1:595565/1‑138 (MQ=11)
                                                           gTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTAtt                                                                               <  2:1146023/148‑1 (MQ=21)
                                                               acggtcggCCTGATAAGACGCATTCCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAATGGCTGTTATTCGGCCtt                                                                        <  1:129220/145‑1 (MQ=255)
                                                               aCTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCtt                                                                        >  2:994076/1‑151 (MQ=37)
                                                                cTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCtt                                                                        <  1:555243/150‑1 (MQ=38)
                                                                cTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGAGAAGACGCATTACGCGTCGCCTCAGGCAACGGCTGTTATTCGGCCtt                                                                        >  1:318448/1‑150 (MQ=255)
                                                                  gTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGtt                                                                     <  2:1171600/151‑1 (MQ=255)
                                                                  gTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGt                                                                      <  2:232986/150‑1 (MQ=255)
                                                                  gTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGt                                                                      <  2:838793/150‑1 (MQ=255)
                                                                  gTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGt                                                                      <  2:412356/150‑1 (MQ=255)
                                                                          gATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTAt                                                             <  1:671754/151‑1 (MQ=255)
                                                                            tAAGACGCATTGCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGG‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTAtt                                                            <  1:22449/150‑1 (MQ=255)
                                                                            tAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTAtt                                                            >  2:644542/1‑150 (MQ=255)
                                                                                      tACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCa                                                  <  2:212817/150‑1 (MQ=255)
                                                                                          cgTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAAt                                                >  2:430681/1‑148 (MQ=255)
                                                                                           gTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGaa                                            <  2:1182901/151‑1 (MQ=255)
                                                                                             cGCATCAGGCAACGTCTGTCGGATTCGGCGTGAACGCCTTATCCTG‑CCTGCGGTTCTGTTCACTGTTGGCCTGATAATACGCATTTCGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGa                                             <  1:35531/148‑1 (MQ=18)
                                                                                             cGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGc                                          <  2:1144080/151‑1 (MQ=255)
                                                                                                   aGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGGAGGCCGGATAAGGCGCATTTCGCGTCGCATCAGGGAACGGCGGTTATTCGGTCTTGTTACCTTTTTTTTGCTGCGCGATAGAAGCCAGaa                                     >  1:1293511/1‑150 (MQ=17)
                                                                                                    ggCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGGTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGaa                                     >  2:343822/1‑149 (MQ=255)
                                                                                                     gCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATcc                                  <  1:755756/151‑1 (MQ=255)
                                                                                                       aaCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATccc                                 >  1:332317/1‑150 (MQ=255)
                                                                                                       aaCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCcg                                <  1:650984/151‑1 (MQ=255)
                                                                                                        aCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCCATAGAAGCCAGAATCCcg                                >  2:753628/1‑150 (MQ=255)
                                                                                                         cGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCcgcg                              <  1:577556/151‑1 (MQ=255)
                                                                                                         cGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCcgcg                              <  1:1037817/151‑1 (MQ=255)
                                                                                                         cGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCcgc                               <  1:320481/150‑1 (MQ=255)
                                                                                                          ggCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCgcgc                             >  2:1082709/1‑151 (MQ=255)
                                                                                                                  gATGCGGCGTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATAtt                     <  1:865649/151‑1 (MQ=255)
                                                                                                                     gcggcgTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCaa                  <  2:142228/151‑1 (MQ=255)
                                                                                                                     gcggcgTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCa                   >  2:247097/1‑150 (MQ=255)
                                                                                                                     gcggcgTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCa                   <  2:750806/150‑1 (MQ=255)
                                                                                                                      cggcgTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCaa                  >  1:828910/1‑150 (MQ=255)
                                                                                                                      cggcgTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAAc                 >  1:868968/1‑151 (MQ=255)
                                                                                                                         cgTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTcc              >  2:1091120/1‑151 (MQ=255)
                                                                                                                          gTGAACGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCt             >  2:1199385/1‑151 (MQ=255)
                                                                                                                               cGCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCtt        >  2:201631/1‑151 (MQ=255)
                                                                                                                                gCCTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCttt       <  1:288391/151‑1 (MQ=255)
                                                                                                                                  cTTATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTcc     >  2:772592/1‑151 (MQ=255)
                                                                                                                                    tATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTcc     <  2:696266/149‑1 (MQ=255)
                                                                                                                                    tATCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCg    <  1:592623/150‑1 (MQ=255)
                                                                                                                                     aTCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCggg  >  1:239078/1‑151 (MQ=255)
                                                                                                                                     aTCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCggg  >  1:646684/1‑151 (MQ=255)
                                                                                                                                     aTCCGA‑CCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCgg   <  2:713416/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                             
GTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGA‑CCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGG  >  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: ‑