Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,229,177 A→C 8.8% intergenic (+112/+261) yohP → / ← dusC uncharacterized protein/tRNA‑dihydrouridine synthase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,229,1770AC 8.8% 599.4 / 6.6 287intergenic (+112/+261)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base A (51/167);  new base C (22/9);  total (94/193)
Fisher's exact test for biased strand distribution p-value = 3.17e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 291 total aligned reads displayed.

ACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229029‑2229326
                                                                                                                                                      |                                                                                                                                                        
  tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAggg                                                                                                                                                      <  2:6399559/150‑1 (MQ=255)
    gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCgggggg                                                                                                                                                    >  1:5266612/1‑150 (MQ=255)
       gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcgg                                                                                                                                                <  1:7908418/150‑1 (MQ=255)
       gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcgg                                                                                                                                                <  2:6908876/150‑1 (MQ=255)
       gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcgg                                                                                                                                                <  1:2863628/150‑1 (MQ=255)
         gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcgggc                                                                                                                                              >  2:1120895/1‑148 (MQ=255)
          aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCa                                                                                                                                             >  2:1514918/1‑150 (MQ=255)
            tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGccca                                                                                                                                            >  2:7807960/1‑148 (MQ=255)
            tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGcccaa                                                                                                                                           >  2:4477426/1‑147 (MQ=255)
             ggagTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATc                                                                                                                                          <  1:4151959/147‑1 (MQ=255)
             ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCCGCCTGTTGCGCTCCTGTCATGGGCATTGCTGCATATACTCACTGGCGCAAGGAGCGCGCGgggggggggggccc                                                                                                                                          >  1:2624352/1‑144 (MQ=255)
             ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATc                                                                                                                                          <  2:7571445/150‑1 (MQ=255)
             ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATc                                                                                                                                          <  2:5027420/150‑1 (MQ=255)
             ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATc                                                                                                                                          <  1:1325096/150‑1 (MQ=255)
              gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcg                                                                                                                                         >  1:4246991/1‑150 (MQ=255)
                 tATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccg                                                                                                                                      <  2:3684106/150‑1 (MQ=255)
                   tttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCACTagccgcc                                                                                                                                    >  2:750504/1‑150 (MQ=255)
                      aaGATGATATTTTAAAATAAATTAATGTCATCAGGTCCGAAAATAACTAGAATATTTAAGTCTCTCATCATGTTGCGCTCCTGTCATGTGCAATGCTTCAGATACTCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAAACgccgccccc                                                                                                                                 >  1:6629440/1‑150 (MQ=255)
                       aGATGATATTTTAAAATTAATTCATGGCATCAGGTCCGAAAATAACGAGAATATTACAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGccc                                                                                                                                <  1:9102740/150‑1 (MQ=255)
                        gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                               <  1:9223225/150‑1 (MQ=255)
                          tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCACCTGTCATGTGCATTGCTTCATAAAATCACTGGCGCAAGGAGCGCGGAGCGGGCGGGCCAACCACcccccccccc                                                                                                                              >  2:4885707/1‑150 (MQ=255)
                           tatatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTg                                                                                                                            <  1:594500/149‑1 (MQ=255)
                            atatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCCACGCccccccccccccc                                                                                                                           >  1:2124092/1‑147 (MQ=255)
                            atatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCCATCCCCCCCCCCCCCCGc                                                                                                                           >  2:9128263/1‑150 (MQ=38)
                            atatTTTAAAATTAATGAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGc                                                                                                                           <  2:6953545/150‑1 (MQ=255)
                             tatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCa                                                                                                                          <  1:8387184/150‑1 (MQ=255)
                                tttAAAATTAATTAATGTCATCAGGTCCGATAATAACGAGAATATTTCAGTCTCTCTTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAccc                                                                                                                       <  2:3852495/150‑1 (MQ=255)
                                tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCTTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCACCCCCTGCAccc                                                                                                                       >  1:10087770/1‑150 (MQ=255)
                                tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGGCAAAAACCCCCCCCCCCCCCAccc                                                                                                                       >  1:333236/1‑150 (MQ=255)
                                tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCACCCCCAGCAccc                                                                                                                       >  1:302191/1‑150 (MQ=255)
                                     aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCggg                                                                                                                  <  2:223538/150‑1 (MQ=255)
                                       ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGGGGCCCCTAGCCGCCCCCCCCCCCCCCCCCCGGcg                                                                                                                >  2:432490/1‑149 (MQ=255)
                                       ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGct                                                                                                                <  2:4092700/150‑1 (MQ=255)
                                       ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGct                                                                                                                <  2:7944853/150‑1 (MQ=255)
                                       ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCAGCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAATAGccgcccccccccgcccccccccgcg                                                                                                                >  1:6851257/1‑145 (MQ=255)
                                         aattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGGGGGCCAAACCACCCCCCCCCCCCCAcccccccccca                                                                                                               >  2:2904461/1‑144 (MQ=255)
                                          attattGTCATCTGGTCCTAAAATAACTAGAATTTTTCAGTCTCTCTTCCTGTTGCGCTCCTTTCATGTGAATTGATTAATATAATCACTGGCGCAAGGAGGGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTg                                                                                                             <  2:2610371/150‑1 (MQ=21)
                                           ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGAGGGGGCGGCCCAAAGCCCCCGCCCCCACCCCCCCCCGGcacagg                                                                                                            >  1:1712356/1‑145 (MQ=255)
                                            taatGTCATCTGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGc                                                                                                           <  2:9662687/150‑1 (MQ=255)
                                              atGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAg                                                                                                                                                        >  1:3685400/1‑104 (MQ=255)
                                                     caGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCGGGGGGCGGCCCAAAGCCGCCGCCCCCCCCACCCCCCGGCACAGGCGCACaaaaa                                                                                                 >  2:6352831/1‑149 (MQ=25)
                                                     caGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCAACGCCGCCGCCCCCAGCACCCCCCGGCGCGGGCGGACaaaaa                                                                                                 >  2:6352938/1‑149 (MQ=25)
                                                        gTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGGCCCGCCCCCCCCCCCCCCCCACCCCCCCCCCCGGGCGAGCaaaaaacc                                                                                              >  1:7982163/1‑146 (MQ=255)
                                                          ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGGCCAAAGCCCCCCCCCCCCCCCCCCCCCCCCCCAGGCGACCaaaaacaccc                                                                                            >  1:7226858/1‑144 (MQ=255)
                                                          ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgccg                                                                                            <  2:1653011/150‑1 (MQ=255)
                                                            gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCATCGCCGCCCCCCCCAGCCCCCCCCCGCGCTGGCGCAACaaaacgccccg                                                                                           >  2:10124606/1‑143 (MQ=25)
                                                            gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCCCCCCCGGGCTCTGGCGAACAAAAACGCCGCt                                                                                          >  2:8478222/1‑150 (MQ=255)
                                                                   aCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGt                                                                                   <  2:6695654/150‑1 (MQ=255)
                                                                     gagaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCAGCAACCCCCGGCTATGGCGAACAAAAAAGCCGCGGCGCGGTGc                                                                                 >  2:3436666/1‑150 (MQ=21)
                                                                         ctatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTc                                                                             <  2:6720933/149‑1 (MQ=255)
                                                                         atatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCAGTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTc                                                                             <  2:5581928/150‑1 (MQ=255)
                                                                           atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTCCATATAATCACTGGCGCACGGAGCGCGCCGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:8796637/150‑1 (MQ=255)
                                                                           atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATATGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:6005747/150‑1 (MQ=255)
                                                                               cAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTa                                                                       <  1:6308374/150‑1 (MQ=255)
                                                                               cAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTa                                                                       <  2:4219207/150‑1 (MQ=255)
                                                                               cAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTa                                                                       <  2:5242443/150‑1 (MQ=255)
                                                                                 gTCTCTCATCGTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAAACACTGGCGAATGGAGCGCGCAGGGGGCGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATc                                                                     <  1:10225172/150‑1 (MQ=18)
                                                                                 gTCTCTCAGCCGGGTGCGCTCCGGTCATGTGCATTGCTTCAGCTAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCCGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATc                                                                     <  2:10024429/150‑1 (MQ=17)
                                                                                 gTCTCTAATCCTGTTGCGCCCCTCTCATGTGAATTGCTTAATATAAACACTGGCCCAAGGAGCGCGCTGGGGGCGCCAAATCGCCGCCCCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATc                                                                     <  1:4765437/150‑1 (MQ=17)
                                                                                     ctcATCCTGTGGCGCTCCTGTCTTGTGAATTGCTTCATATAATCACTGGCGCACGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCtt                                                                 <  1:9561980/150‑1 (MQ=18)
                                                                                      tcATCCTGCTGCGCTCCTGTCATGTGCATTGCTTCATATAATCAGTGGCGCAAGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  1:11463/150‑1 (MQ=255)
                                                                                        aTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                              <  2:8344443/150‑1 (MQ=255)
                                                                                        aTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                              <  2:7380154/150‑1 (MQ=255)
                                                                                        aTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                              <  2:1721545/150‑1 (MQ=255)
                                                                                        aTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                              <  2:6890879/150‑1 (MQ=255)
                                                                                        aTCCTGTTGCGCTCCTGTCATGTGCATTGCTACATATAAGCAGTGGCGCACGGCGCGCGCTGGGGGCGCCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                              <  1:4808619/150‑1 (MQ=17)
                                                                                            tGTTGCGCTCCTGTCATGTGCATTGCTTAATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTa                                                          <  2:1962384/150‑1 (MQ=255)
                                                                                             gTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAc                                                         <  2:7129948/150‑1 (MQ=255)
                                                                                               tGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCg                                                       <  2:1307880/150‑1 (MQ=255)
                                                                                                 cgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAAACGCCCCCCCCCCCAGCACCCCCCGGCGCAGGCGCACAAAAACGCCCCGGCGCGGGGCCCCCGGCGTAACCCCTACCGCTAccccg                                                     >  2:6511822/1‑147 (MQ=17)
                                                                                                   cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                   <  2:5768213/150‑1 (MQ=255)
                                                                                                     ccTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCgg                                                 <  1:1514918/150‑1 (MQ=255)
                                                                                                      cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCggg                                                <  2:1823515/150‑1 (MQ=255)
                                                                                                       tGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGgc                                               <  1:8814541/150‑1 (MQ=255)
                                                                                                        gTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGgcg                                              <  1:2142951/150‑1 (MQ=255)
                                                                                                           cggtgCAGTGCGTCAGAGACTCACTGGCGCCCGGAGCGCGCCGGGGGCGGCCACTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGACCAAAATCGCCGCTGCGCGGGGCCCTCGGCTTCTCCCTTACGGCGACCGGGTCGGGCGCGa                                           <  2:8907746/148‑1 (MQ=11)
                                                                                                           aTGTGCATTGCTTCATATACTCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  2:9449583/150‑1 (MQ=255)
                                                                                                           aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCATGGGGCGGCCAAACGCCGCCGCCCCCAGCCCCCCCCGGCTATGGCGAACAAAAACCCCGCGGCGCGGTGGCCACGGCATAACCCATAACGCTACCCGGTCGGGcgcgc                                           >  2:2504346/1‑149 (MQ=255)
                                                                                                           aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGGCCAACGCCGCCCCCCCCCGCACCCCCCGGCGCAGGGGGCCAAAAAAGCCGCGGCGCGGGGCCCACGGGGGATCCCCGAAGGCGACCCGGGGGGGcgcgc                                           >  2:9862290/1‑149 (MQ=12)
                                                                                                           aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  2:972547/150‑1 (MQ=255)
                                                                                                            tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCAAACGCCGCCCCCCCCAGCACCCCCGGGCTCTGGCGCACAAAATCCCCGCTGCGCCGTGCCCTCGGCTTATCCCATAACGCGACCCGGTCGGGcgcgcg                                          >  2:3289100/1‑148 (MQ=17)
                                                                                                            tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGACCAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                          <  2:9051928/150‑1 (MQ=255)
                                                                                                            tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                          <  2:191608/150‑1 (MQ=255)
                                                                                                              tgCATTGCTTCATATAATCAGTGGCGAAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGt                                        <  1:4728962/150‑1 (MQ=255)
                                                                                                              tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGt                                        <  2:6795639/150‑1 (MQ=255)
                                                                                                               gCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                       <  2:5479274/150‑1 (MQ=255)
                                                                                                                cATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                      <  2:6634453/150‑1 (MQ=255)
                                                                                                                  ttGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACGCCCCCCCCCCCAGCACCCCCCGGCTCAGGCGAAAAAAAAAGCCCCGGCGCGGTGCCCACGGCGTGTACCGTAACGCAAACCGGTCGGGCGCGCGGGAACa                                    >  2:10101488/1‑150 (MQ=9)
                                                                                                                     cTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                 <  1:5672492/150‑1 (MQ=40)
                                                                                                                      ttCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCCCCCCCAGCACCCCCCGGCTCTGGCGAACAAAAAACCCGCGGCCCGGTGCCCACGGCGTATCCCATAAGGCTAACCGGTCGGGCGCGAGGTAACATacc                                >  2:10289735/1‑150 (MQ=14)
                                                                                                                       tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  2:7248947/150‑1 (MQ=39)
                                                                                                                       tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:1855658/150‑1 (MQ=39)
                                                                                                                                    ggCGCCTGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc                  <  1:8626612/150‑1 (MQ=12)
                                                                                                                                      cgcAAGGAGCGCGCAGGGGGGCGCCCAAACCCCCCCCCCCCACCACCCCCCGGCGCAGGCGAACAAAAAAGCCCCGGCGCGGGGGCCCCCGCGGAACACAAAAGGCGACCCGGGCGGGCGCGAGGGAACAAACCAGTAAAACCCGccccc                >  2:3232206/1‑148 (MQ=2)
                                                                                                                                          aGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  2:4663737/150‑1 (MQ=32)
                                                                                                                                           ggAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCACGGGCTCTGGCGAACAAAATCGCCGCTTCGCTGTGCCCTCGGCTTATCCCTTACTGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc           <  1:4954180/150‑1 (MQ=255)
                                                                                                                                             aGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcac         <  1:974349/150‑1 (MQ=32)
                                                                                                                                                 cgcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc     <  1:5454848/150‑1 (MQ=32)
                                                                                                                                                      aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:4850231/1‑150 (MQ=32)
                                                                                                                                                      aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:52904/1‑150 (MQ=32)
                                                                                                                                                      |                                                                                                                                                        
ACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229029‑2229326

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