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,5230AG96.5% 270.3 / ‑5.6 87intergenic (‑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/1);  major base G (33/50);  minor base T (3/0);  total (36/51)
Fisher's exact test for biased strand distribution p-value = 6.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCG  >  NC_000913/2662376‑2662668
                                                                                                                                                   |                                                                                                                                                 
tACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCgg                                                                                                                                                >  1:668620/1‑151 (MQ=255)
tACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCCTATCCGACCTGCgg                                                                                                                                                >  1:70257/1‑151 (MQ=255)
   gcatCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTTTCCGACCCTCGGtt                                                                                                                                              >  2:1468315/4‑150 (MQ=255)
    cgTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCt                                                                                                                                            <  1:338586/151‑1 (MQ=255)
       cGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTTTCCGCCCTGCGGTTCt                                                                                                                                            >  2:1115983/1‑148 (MQ=255)
        gCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGt                                                                                                                                          >  2:685375/1‑149 (MQ=255)
          aTCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCAc                                                                                                                                      <  2:1369842/151‑1 (MQ=255)
            cAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACt                                                                                                                                     <  1:1186724/150‑1 (MQ=255)
            cAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTg                                                                                                                                    <  1:662688/151‑1 (MQ=255)
              ggCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGt                                                                                                                                   <  2:658921/150‑1 (MQ=255)
                   cGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGcc                                                                                                                              >  1:1324948/1‑150 (MQ=255)
                   cGGCTGTCGGATGCGGCGTGAACGCCTTACCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGGTGCGGCGTGGACGCCTTTTTCGACCTTCGGTTCTGTTCACTGTTGGc                                                                                                                               >  2:1197996/1‑149 (MQ=14)
                     gCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGcc                                                                                                                              >  1:1258057/1‑148 (MQ=255)
                     gCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGa                                                                                                                           <  1:1154873/151‑1 (MQ=255)
                        gTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGAt                                                                                                                          <  2:729406/149‑1 (MQ=255)
                        gTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGAt                                                                                                                          <  2:511683/149‑1 (MQ=255)
                        gTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGAt                                                                                                                          <  2:860521/149‑1 (MQ=255)
                          cGGATGCGGCGTGAATGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATaa                                                                                                                        >  1:162307/1‑149 (MQ=255)
                          cGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAg                                                                                                                       <  1:882530/150‑1 (MQ=255)
                          cGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGa                                                                                                                      <  2:894985/151‑1 (MQ=255)
                           ggATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGa                                                                                                                      >  1:740268/1‑150 (MQ=255)
                           ggATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGGGAACCCCTTATCCGCCCTTCGGTTTTGTTCACTGTAGGCCCGATAAGAc                                                                                                                     >  2:1360503/1‑151 (MQ=14)
                                cggcgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCAtt                                                                                                                >  1:392902/1‑151 (MQ=255)
                                 ggcgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCCtt                                                                                                                >  2:1390715/1‑150 (MQ=255)
                                   cgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAcg                                                                                                             <  1:1077836/151‑1 (MQ=255)
                                       aaCGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGta                                                                                                         <  2:634815/151‑2 (MQ=255)
                                       aaCGCCTGATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTc                                                                                                         <  2:1279902/151‑1 (MQ=255)
                                         cGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGTACCCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACTCATTACGCGt                                                                                                          >  2:720294/1‑148 (MQ=255)
                                           ccTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCAt                                                                                                     <  1:108046/151‑1 (MQ=255)
                                               aTCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                                                                  <  2:462778/150‑1 (MQ=38)
                                               aTCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                                                                  <  2:1069828/150‑1 (MQ=38)
                                                 ccGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGc                                                                                                <  2:70248/150‑1 (MQ=38)
                                                  cGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCa                                                                                               >  1:583765/1‑150 (MQ=38)
                                                    aCCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACg                                                                                            <  1:1026538/151‑1 (MQ=25)
                                                     ccTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACgg                                                                                           <  1:223923/151‑1 (MQ=25)
                                                       tACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCCGCGGTTCTGTTCACTGTACGCCTGAGAAGACGCATTACGCGTCGCATCAGGCAACCGc                                                                                          >  2:171537/1‑150 (MQ=255)
                                                        aCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCt                                                                                         <  2:197734/150‑1 (MQ=18)
                                                        aCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCt                                                                                         <  2:331787/150‑1 (MQ=18)
                                                        aCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTg                                                                                        <  1:1086356/151‑1 (MQ=18)
                                                          ggTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGtt                                                                                      >  1:644771/1‑151 (MQ=17)
                                                          ggTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGtt                                                                                      <  2:630718/151‑1 (MQ=17)
                                                           gTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGGCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTggt                                                                                      >  2:723162/1‑148 (MQ=9)
                                                                gTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTAtt                                                                                   <  1:1455376/148‑1 (MQ=21)
                                                                gTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTAtt                                                                                   <  1:1417449/148‑1 (MQ=21)
                                                                    aCTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGTACGCCTTATCCGACCTGCGGTTCTGTTCCCTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGcc                                                                              >  2:1080200/1‑149 (MQ=255)
                                                                    aCTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTTATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCtt                                                                            >  1:647197/1‑151 (MQ=255)
                                                                    aCTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCtt                                                                            >  2:1026716/1‑151 (MQ=37)
                                                                    aCTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCtt                                                                            <  1:685451/151‑1 (MQ=37)
                                                                       gTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTg                                                                           <  1:124164/149‑1 (MQ=38)
                                                                          ggCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTAcc                                                                      <  1:1274533/151‑1 (MQ=255)
                                                                               gATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTAt                                                                 <  1:1471021/151‑1 (MQ=255)
                                                                                aTAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTAtt                                                                <  2:1387319/151‑1 (MQ=255)
                                                                                aTAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTAtt                                                                >  1:401137/1‑151 (MQ=255)
                                                                                  aaGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTTGTACCTTTTTTct                                                              >  2:1082972/1‑151 (MQ=255)
                                                                                     aCGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctgct                                                           <  2:81506/151‑1 (MQ=255)
                                                                                      cGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctgct                                                           <  1:793392/150‑1 (MQ=255)
                                                                                         aTTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTTTTCACTGTAGGCCTGATAAGACGCATTACTCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTTACTTTATTctgctcccc                                                       >  1:1377455/1‑147 (MQ=14)
                                                                                             cgcgTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCTTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATa                                                   <  2:1103324/151‑1 (MQ=255)
                                                                                              gcgTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAg                                                  <  2:49350/151‑1 (MQ=255)
                                                                                                gTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGa                                                 <  1:1233931/150‑1 (MQ=255)
                                                                                                 tCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGa                                                 >  1:227188/1‑149 (MQ=255)
                                                                                                  cGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAg                                               <  1:395148/150‑1 (MQ=255)
                                                                                                  cGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAg                                               <  1:401476/150‑1 (MQ=255)
                                                                                                  cGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGc                                              <  1:1162801/151‑1 (MQ=255)
                                                                                                       cAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGaa                                         <  1:536670/151‑1 (MQ=255)
                                                                                                       cAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGaa                                         <  1:1190168/151‑1 (MQ=255)
                                                                                                       cAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGa                                          >  2:249622/1‑150 (MQ=255)
                                                                                                            aaCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATccc                                     >  2:1402124/1‑150 (MQ=255)
                                                                                                            aaCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTAACTTTATTTTTCTGCTCAATAGAAGCCAGAATccc                                     >  1:1250144/1‑150 (MQ=255)
                                                                                                              cGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCcgcg                                  <  2:9993/151‑1 (MQ=255)
                                                                                                                gCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCgcgc                                 >  1:545208/1‑150 (MQ=255)
                                                                                                                  tGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAgg                              >  2:975386/1‑151 (MQ=255)
                                                                                                                   gTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAg                               <  2:406874/149‑1 (MQ=255)
                                                                                                                     cGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGa                             >  2:538498/1‑149 (MQ=255)
                                                                                                                       gATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATAtt                         <  1:33571/151‑1 (MQ=255)
                                                                                                                          gcggcgTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCa                       >  1:834981/1‑150 (MQ=255)
                                                                                                                           cggcgTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAAc                     >  2:1057976/1‑151 (MQ=255)
                                                                                                                             gcgTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTc                   <  2:286539/151‑1 (MQ=255)
                                                                                                                             gcgTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTc                   <  1:1468423/151‑1 (MQ=255)
                                                                                                                                 gAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTg                <  2:663233/150‑1 (MQ=255)
                                                                                                                                   aCGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGTTAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCt             >  1:958629/1‑151 (MQ=255)
                                                                                                                                     gCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCttt           <  2:923582/151‑1 (MQ=255)
                                                                                                                                        ttATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTcc         >  2:712301/1‑150 (MQ=255)
                                                                                                                                         tATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTcc         <  1:631590/149‑1 (MQ=255)
                                                                                                                                         tATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTcc         <  1:1390819/149‑1 (MQ=255)
                                                                                                                                            ccGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgcg    <  2:28094/151‑1 (MQ=255)
                                                                                                                                              gACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgcgcg  >  2:1389222/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCG  >  NC_000913/2662376‑2662668

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