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,5230AG99.2% 446.0 / NA 133intergenic (‑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 (57/75);  minor base T (1/0);  total (58/75)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 133 total aligned reads displayed.

TACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCA  >  NC_000913/2662376‑2662670
                                                                                                                                                   |                                                                                                                                                   
tACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCgg                                                                                                                                                  >  2:63094/1‑151 (MQ=255)
 aCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGt                                                                                                                                                 <  2:281365/151‑1 (MQ=255)
  cgcgTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGtt                                                                                                                                                <  1:1289333/151‑1 (MQ=255)
     gTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGGGGCGGGAACGCCTTATTCGGCCTGCGGtt                                                                                                                                                >  2:311524/1‑148 (MQ=255)
         cATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCa                                                                                                                                         <  1:737743/151‑1 (MQ=255)
            ctggtaaCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTAGGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCa                                                                                                                                         <  4:396836/143‑1 (MQ=255)
            cAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACt                                                                                                                                       <  2:477050/150‑1 (MQ=255)
               gCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGt                                                                                                                                     >  2:20818/1‑149 (MQ=255)
                    ggCTGTCGGATGCGGCGTGAACGCCTTATCCGCCCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTg                                                                                                                                      <  1:218543/143‑1 (MQ=255)
                        gTCGGATGCGGCGTGAACGCCTTATCCGCCCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATa                                                                                                                           <  4:307673/150‑1 (MQ=25)
                         tCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTTACTGTAGGCCTGAt                                                                                                                            >  2:1641302/1‑148 (MQ=255)
                         tCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCGCTGTAGGCCTGATa                                                                                                                           >  3:77076/1‑149 (MQ=255)
                           ggATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGa                                                                                                                        >  2:142907/1‑150 (MQ=255)
                           ggATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGAc                                                                                                                       >  3:440460/1‑151 (MQ=255)
                             aTGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGc                                                                                                                     <  3:328626/151‑1 (MQ=255)
                               gcggcgTGAACGCCTTATCCGCCCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCAt                                                                                                                   <  1:431722/151‑1 (MQ=25)
                               gcggcgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCa                                                                                                                    <  2:774934/150‑1 (MQ=255)
                                   cgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAcg                                                                                                               <  2:759779/151‑1 (MQ=255)
                                   cgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAcg                                                                                                               <  4:382889/151‑1 (MQ=255)
                                    gTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAcgc                                                                                                              <  2:376193/151‑1 (MQ=255)
                                        aCGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCg                                                                                                          <  2:373796/151‑1 (MQ=255)
                                          gCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCa                                                                                                        <  4:358767/151‑1 (MQ=255)
                                           ccTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCAt                                                                                                       <  1:1316063/151‑1 (MQ=255)
                                              tATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGCCCTGCGGTTCTGTTCACTGTAGGCCGGATAAGACGCATTACCCGTCGCATc                                                                                                      >  4:73078/1‑149 (MQ=255)
                                                 ccGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTTCGCGTCGCATCAGGc                                                                                                  >  1:534512/1‑150 (MQ=255)
                                                 ccGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCa                                                                                                 <  2:995634/151‑1 (MQ=38)
                                                  cGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCaa                                                                                                <  2:1087995/151‑1 (MQ=37)
                                                      cTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTACGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACCGc                                                                                            >  4:397189/1‑151 (MQ=19)
                                                        aCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTTCGCGTCGCATCAGGCAACGGCTg                                                                                          >  1:450941/1‑151 (MQ=18)
                                                        aCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGc                                                                                            <  1:175214/149‑1 (MQ=18)
                                                        aCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTg                                                                                          >  2:1544317/1‑151 (MQ=18)
                                                         cGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGt                                                                                         >  2:1150229/1‑151 (MQ=17)
                                                          ggTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGGCTGAGAAGACGCATTACCCGTCGCATCAGGCAACGGCTggt                                                                                        >  1:1385776/1‑149 (MQ=255)
                                                          ggTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTggt                                                                                        >  3:132270/1‑149 (MQ=14)
                                                          ggTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGt                                                                                         <  1:1003028/150‑1 (MQ=14)
                                                           gTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTa                                                                                       <  2:75439/151‑1 (MQ=18)
                                                            ttCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTAt                                                                                      >  2:759059/1‑151 (MQ=18)
                                                              cTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGGGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGACTGATAAGACGCATTACGCGGCGCATCAGGCAAAGGCTGTTAtt                                                                                     >  2:1540803/1‑150 (MQ=255)
                                                               tGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTc                                                                                    >  2:1681097/1‑150 (MQ=25)
                                                               tGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCCGCGTGAACGCCTTATCCGCTCCTCGGTTTCATTCACCGTAGGCCTCATAAGACGCATTACGCCACCCCTTCGGCAACCGCTGGTATTc                                                                                    >  1:1185847/1‑150 (MQ=11)
                                                                gTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCgg                                                                                  >  2:1199904/1‑151 (MQ=37)
                                                                gTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCgg                                                                                  <  2:403615/151‑1 (MQ=37)
                                                                 ttCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGc                                                                                 >  2:1408694/1‑151 (MQ=37)
                                                                  tCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGc                                                                                 <  2:1211381/150‑1 (MQ=25)
                                                                  tCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGc                                                                                 <  1:1001551/150‑1 (MQ=25)
                                                                   cACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCt                                                                               >  2:135931/1‑151 (MQ=37)
                                                                   cACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCt                                                                               <  2:1515908/151‑1 (MQ=37)
                                                                   cACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCt                                                                               <  1:701707/151‑1 (MQ=37)
                                                                    aCTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGcc                                                                                >  2:1003882/1‑149 (MQ=25)
                                                                       gTAGGCCTGATAAGACGCATTACGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGt                                                                            <  1:16411/150‑1 (MQ=255)
                                                                       gTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACTTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGtt                                                                           <  1:1075338/151‑1 (MQ=255)
                                                                       gTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTg                                                                             <  1:1349108/149‑1 (MQ=38)
                                                                       gTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGt                                                                            <  1:1000489/150‑1 (MQ=255)
                                                                       gTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTCTTATTCGGCCTTGt                                                                            >  2:344866/1‑150 (MQ=255)
                                                                         aGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTAc                                                                         <  3:35737/151‑1 (MQ=255)
                                                                          ggCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTAcc                                                                        <  1:1676409/151‑1 (MQ=255)
                                                                           gCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTAcc                                                                        >  2:307631/1‑150 (MQ=255)
                                                                           gCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTAcc                                                                        >  2:697508/1‑150 (MQ=255)
                                                                           gCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCt                                                                       <  4:153110/151‑1 (MQ=255)
                                                                            ccTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCtt                                                                      >  2:318121/1‑151 (MQ=255)
                                                                             cTGATAAGCCGCATTACGCGTCGCATCCGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCt                                                                       <  1:350986/149‑1 (MQ=255)
                                                                               gATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTa                                                                    <  2:63063/150‑1 (MQ=255)
                                                                               gATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTAt                                                                   <  2:425909/151‑1 (MQ=255)
                                                                                   aGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctg                                                               <  2:1650091/151‑1 (MQ=255)
                                                                                    gACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctg                                                               >  2:1690359/1‑150 (MQ=255)
                                                                                      cGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTTGTACCTTTATTcttc                                                              >  4:560701/1‑147 (MQ=255)
                                                                                      cGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctgctg                                                            <  1:444365/151‑1 (MQ=255)
                                                                                        cATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTgctgct                                                          <  1:137310/151‑1 (MQ=255)
                                                                                         aTTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTc                                                         >  3:371303/1‑151 (MQ=255)
                                                                                           tACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCa                                                        <  1:1165220/150‑1 (MQ=255)
                                                                                           tACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCa                                                        <  1:227677/150‑1 (MQ=255)
                                                                                            aCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAAt                                                      <  4:519188/151‑1 (MQ=255)
                                                                                            aCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAAt                                                      <  1:675676/151‑1 (MQ=255)
                                                                                               cgTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGa                                                   <  1:1574015/151‑1 (MQ=255)
                                                                                               cgTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGa                                                   <  4:566145/151‑1 (MQ=255)
                                                                                               cgTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGa                                                   <  2:1705900/151‑1 (MQ=255)
                                                                                                gTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGa                                                   >  3:533744/1‑150 (MQ=255)
                                                                                                gTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGGTATTCGGCCTTTTTAACTTTATTCTGCTTCTCCATAg                                                    >  4:172197/1‑149 (MQ=17)
                                                                                                   gCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGcc                                               <  2:1072710/151‑1 (MQ=255)
                                                                                                        aGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGaa                                           >  2:652321/1‑150 (MQ=255)
                                                                                                          gCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATcc                                        <  1:142880/151‑1 (MQ=255)
                                                                                                              cGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCcgcg                                    <  2:1198326/151‑1 (MQ=255)
                                                                                                              cGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCcgcg                                    >  2:854792/1‑151 (MQ=255)
                                                                                                                 cTGTCGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCgcgc                                   >  4:588239/1‑149 (MQ=255)
                                                                                                                     cGGATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGat                              <  2:1451988/150‑1 (MQ=255)
                                                                                                                      ggATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGatat                            >  2:568027/1‑151 (MQ=255)
                                                                                                                      ggATGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGatat                            >  1:748961/1‑151 (MQ=255)
                                                                                                                        aTGCGGCGTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTc                          <  2:2015/151‑1 (MQ=255)
                                                                                                                          gcggcgTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCa                         <  1:715624/150‑1 (MQ=255)
                                                                                                                           cggcgTGAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCa                         >  1:537528/1‑149 (MQ=255)
                                                                                                                                 gAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTgg                 >  1:533174/1‑151 (MQ=255)
                                                                                                                                 gAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTgg                 >  2:638685/1‑151 (MQ=255)
                                                                                                                                 gAACGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTg                  <  1:1663122/150‑1 (MQ=255)
                                                                                                                                    cGCCTTATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTTCTGGCtt              <  1:1232737/151‑1 (MQ=255)
                                                                                                                                        ttATCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTcc           >  1:1586674/1‑150 (MQ=255)
                                                                                                                                          aTCCGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCggg        >  2:765680/1‑151 (MQ=255)
                                                                                                                                             cGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgcgc     <  1:813765/151‑1 (MQ=255)
                                                                                                                                             cGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgcgc     >  2:1204389/1‑151 (MQ=255)
                                                                                                                                             cGACCTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGgcgc     >  2:1591627/1‑151 (MQ=255)
                                                                                                                                                ccTGCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCa  <  4:5010/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCTGGCTTTCCGGGCGCGCA  >  NC_000913/2662376‑2662670

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