Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,224,483 C→T S219F (TCC→TTC)  metH → cobalamin‑dependent methionine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,224,4830CT97.1% 220.8 / NA 70S219F (TCC→TTC) metHcobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (37/31);  minor base G (0/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGC  >  NC_000913/4224237‑4224725
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
aCCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTGCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCggg                                                                                                                                                                                                                                                <  2:310682/251‑1 (MQ=255)
    cACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGgcgcc                                                                                                                                                                                                                                            >  1:303902/1‑250 (MQ=255)
              aTTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTcc                                                                                                                                                                                                                                    <  2:127788/249‑1 (MQ=255)
                ttCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGGTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTTAGCTGCCGATTATGATCTCCGGCACCCTCACCGACGCCTTCGGgcgccccaccttcgggc                                                                                                                                                                                                                                >  2:246055/1‑238 (MQ=255)
                  ctcCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGc                                                                                                                                                                                                                                >  2:148062/1‑249 (MQ=255)
                      ggACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGacca                                                                                                                                                                                                                          <  1:397235/251‑1 (MQ=255)
                            cAACGATCCGGCATTTCGTACTATCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGacg                                                                                                                                                                                                                    >  2:90347/1‑249 (MQ=255)
                            cAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGACGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGCTTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAACACTGGGCGTTGAGCTGCCGGTTCTGATCTCCGGCACCCTAACCGCCGCCTTCGGGCGCACGCGCTCAGGGCCGaccacaga                                                                                                                                                                                                                      >  2:389767/1‑246 (MQ=255)
                             aaCGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGc                                                                                                                                                                                                                   >  1:457719/1‑251 (MQ=255)
                                     ggCATTTCGTAATATCACTTTTGACGGGCTTGTTGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCTGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTTAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACa                                                                                                                                                                                                           <  1:422521/251‑1 (MQ=255)
                                      gCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACa                                                                                                                                                                                                           >  1:182410/1‑250 (MQ=255)
                                      gCATTGCGTGATATCCCTGTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGACGGTGGGGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACTCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTGCaa                                                                                                                                                                                                          <  2:256733/251‑1 (MQ=255)
                                            cGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCAtt                                                                                                                                                                                                    >  1:126327/1‑251 (MQ=255)
                                            cGTAATATCACTGTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATGATGAGCTCCGTCACCATCACCGACGCCTGCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCAtt                                                                                                                                                                                                    <  1:439450/251‑1 (MQ=255)
                                               aaTATCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTgc                                                                                                                                                                                                  <  2:143908/250‑1 (MQ=255)
                                                atatCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAACGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGCTTGAAACCGTTTTCGCCACCCTTAACGCCAAAGCGGCGGTGTTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCCCCATCCCCGACGCCTTCGGGCGCACGCCCTCCGGGCCGACCCCCGAAGCCTTTTCCAACTCCTTgcgc                                                                                                                                                                                                >  2:224238/1‑251 (MQ=255)
                                                      cTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGcc                                                                                                                                                                                          >  1:167551/1‑251 (MQ=255)
                                                          ggACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAg                                                                                                                                                                                      <  1:150133/250‑1 (MQ=255)
                                                           gACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTTAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGc                                                                                                                                                                                     <  2:408444/251‑1 (MQ=255)
                                                              gggCTGGTGGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGctct                                                                                                                                                                                  <  1:148062/251‑1 (MQ=255)
                                                                  tggtggCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGAcc                                                                                                                                                                              >  1:376461/1‑251 (MQ=255)
                                                                   ggtggCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCCCTGACCt                                                                                                                                                                             >  2:348256/1‑251 (MQ=255)
                                                                      ggcggcTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCttt                                                                                                                                                                           <  1:133748/250‑1 (MQ=255)
                                                                         ggcTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACGCATTGCGCCACGCCGAAGCTCTGACCGTTGGcc                                                                                                                                                                       <  2:86443/251‑1 (MQ=255)
                                                                             tATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCCGTATTTGCGGTGAAAACGGAGTTTGAAGAGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCCTCACCGACGCCTTCGGGCGCACGCTATCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGGAGCTCTGACCTTTGGCCTGa                                                                                                                                                                    >  1:443277/1‑250 (MQ=255)
                                                                                cGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACtg                                                                                                                                                                >  1:257266/1‑251 (MQ=255)
                                                                                  agagGCCAGCAAAGCGCTGGTGGAAGGTGGCGCGGATCGGATCCTGATTGAAACCGTTTTCGACCCCCTTAACGCCAAATCGGCGGTATTTGCGGTGAAAACGGAGGTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCGCCGGCACCATCACCGACTCCTTCGGGCGCACGCTCTCCGGGCAGGCCACCGATGCATTTTACAACTCATGGCGCCACGCCGAAGCTCTGACCTTTGGCCGGAACtgtg                                                                                                                                                              <  2:179827/251‑1 (MQ=255)
                                                                                              aGCGCTGGTGGAAGGAGGCGCGGATATGATCCTGATTGAACCCGTTTTCGACACCCTCAACGACAACGCGGCGGCATTTGCGGTGCAAACGGAGTTTGAACCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGGTCTCCGGGCAGACCACCGGAGCATTTTACAACTCATTGCGCCCCGCCGGCGCTCTGACCTTTGGCCTGAACCGTGCCCTGGGGCCCg                                                                                                                                                  >  2:205380/1‑251 (MQ=255)
                                                                                                       ggAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTgc                                                                                                                                         >  1:18534/1‑251 (MQ=255)
                                                                                                            gTGGCGCGGATCTGATCCTGATTGACACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTAATTGCCCCACGCCGAAGCTCTGACCTTTGGCCTGAACTCTGCGCTGGGGCCCGATGAACTGCGcccg                                                                                                                                    >  2:87307/1‑249 (MQ=255)
                                                                                                               gcgcGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTAAGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTAc                                                                                                                                 <  1:45529/251‑1 (MQ=255)
                                                                                                                   ggATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGc                                                                                                                             <  2:18534/251‑1 (MQ=255)
                                                                                                                      tCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAgg                                                                                                                          >  1:60322/1‑251 (MQ=255)
                                                                                                                         gATCCTGATTGAAACCGTTTTCGACAACCTTACCGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGc                                                                                                                       >  2:326305/1‑251 (MQ=255)
                                                                                                                            ccTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTg                                                                                                                     <  1:396692/250‑1 (MQ=255)
                                                                                                                                tttGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCAGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGCAGCTCTGACCTTTGGCCTTAACTGTGCGCTGTGGCCCGATGAACTTCGCCAGTTCGTGCCTGAGCTGTCACg                                                                                                                >  2:334289/2‑251 (MQ=255)
                                                                                                                                  tGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGCAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGCAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGGTGAACGGCGCCAGTACGGGCAGGAGCTGTCACggg                                                                                                              >  2:174743/1‑250 (MQ=255)
                                                                                                                                    aaaCCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGaa                                                                                                                                                                                                                     <  1:453960/146‑1 (MQ=255)
                                                                                                                                            ttcgatccccctacacccAAAGCGTCGGTATTAGCGTTGCAAACGGAGTTTGAAGCGCTGTGCGTTGAGCTGCCGATTTTGATCTCCGGCTCCATCGCCGACTCCATCTGGCTCACGCTCTCCGGGCCGACCACCGAAGCATTTTACAACTCATTGCGCCACGACGAAGCTCTTCCCTTTGGCCTGATCTGTGCGCTGGGGCCCGATTATCGGCGCCAGTACGTGCAGTAGTTGTCACGTATTGCTGAAt                                                                                                     <  2:426172/234‑1 (MQ=255)
                                                                                                                                            ttCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAAt                                                                                                     <  1:457402/250‑1 (MQ=255)
                                                                                                                                             tCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATg                                                                                                    <  2:239263/250‑1 (MQ=255)
                                                                                                                                                accgccTTAACGCCAAAGCTGCGGGATTTGCGGTGAAAGCGGAGGTTGATGCGCTGGGCGTTGAGCTGCCGATTATGGTCTCCGGCGCCATCACCGACGCCTTCGGGCGGACGCTCTCCGGGCAGACCACCGAAGCATTTTACGACTCATGGCGCCACGCCGGAGCTCTGACCTTTGGCCTGAGCTGTGCGCGGGGGCCCGATGAACTGCGCCAGTACGTGCAGGGGCTGGCACGGATTGCGGAATGCTAc                                                                                                <  2:69227/247‑1 (MQ=255)
                                                                                                                                                 cacCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCGTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCTCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACg                                                                                               <  2:303902/251‑1 (MQ=255)
                                                                                                                                                  acCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCCGGAGCTGTCACGGCTTGCGGAATGCTACGt                                                                                              >  1:445193/1‑251 (MQ=255)
                                                                                                                                                        aaCGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACcgc                                                                                        >  1:132200/1‑251 (MQ=255)
                                                                                                                                                          cGCCAAAGCGGCGGTATTTTCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACCCCGAAGCTCTGGCCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAAATGCGCCAGGACGTTCAGGAGCTGTACCGGATTGCGGAATGCTACCACACCgcgc                                                                                      >  2:99319/1‑251 (MQ=255)
                                                                                                                                                             ccaaGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGCGCTGCCGGTTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGAGCAGACCACCGAAGCATTTTAAAACTAATTGCGCCACGCCGAAGCTCTGCCCTTTGGCCTGAACTGAGCGCTGGGGCCCGAAGAACGGCGACAGTACGGGCAGGAGCTGTCCCGGAGTGCGGGCTGCTACGTCACCGCGAAcc                                                                                   >  2:126960/3‑251 (MQ=255)
                                                                                                                                                                  cggcggTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAAc                                                                              <  1:83570/251‑1 (MQ=255)
                                                                                                                                                                     cggTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACGGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGcc                                                                           <  2:457719/251‑1 (MQ=255)
                                                                                                                                                                       gTATTTGCGGTGAAAACGGAGTGTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCgg                                                                         <  1:442841/251‑1 (MQ=255)
                                                                                                                                                                                   acaaCGGAGATCGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCCCCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACAGAAGAATTTTACAACTCATTGCGCCACACCGCAGCTATGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGACAGTCCGAGCAGGAGCTGTCACGGCTTGCGGAATGCTAACTCACCGCGCACCCGAACGACGGGCTACCCAACGc                                                             >  2:276504/3‑251 (MQ=255)
                                                                                                                                                                                    aaaCGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGcc                                                            >  1:365305/1‑251 (MQ=255)
                                                                                                                                                                                      aCGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCtt                                                          <  1:105555/251‑1 (MQ=255)
                                                                                                                                                                                        ggAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTg                                                        >  1:320914/1‑251 (MQ=255)
                                                                                                                                                                                         gAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGACCGCCGGGCTACCCACCGCCTTTgg                                                       >  2:117431/1‑251 (MQ=255)
                                                                                                                                                                                          aGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGCGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGATACCCAACGCCTTTggg                                                      >  2:185256/1‑250 (MQ=255)
                                                                                                                                                                                               gAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTGCGGGCGCACGCTCTCCGGGCAGACAACCGAAGCATTTTACAACTCAGTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTTCGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCACCGCCTTTGGTGAGTa                                                 <  2:365707/251‑1 (MQ=255)
                                                                                                                                                                                               gAAGCGCTGGGCGTTGAGCTGCCGATTATGACCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTa                                                 <  1:411078/251‑1 (MQ=255)
                                                                                                                                                                                                   cgcTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCCCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAt                                             >  2:231168/1‑251 (MQ=255)
                                                                                                                                                                                                             tGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCACCGCCTTTGGTGCGTACGATCTCGACGCCg                                   >  1:411859/1‑251 (MQ=255)
                                                                                                                                                                                                             tGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGcc                                    <  1:252520/250‑1 (MQ=255)
                                                                                                                                                                                                               aGCTGCCGATTATGATCTCCGGATCCATCACCGACGCCTTCGGGCGCACTCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCAGTGCGCCACGCCGAAGCTCTGACCTTTGGCCGGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGGCACTGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCTGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGac                                 <  2:439204/251‑1 (MQ=255)
                                                                                                                                                                                                                 cTGCCGGTTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGGCGCCGccac                               >  2:322284/1‑251 (MQ=255)
                                                                                                                                                                                                                      gATTATGATCTCCGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGCCCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTAAGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAAAGCCTTTGGTGAGTACGATATCGACGCCGACACGATgg                          >  1:460703/1‑251 (MQ=255)
                                                                                                                                                                                                                       aTTATGATCTCCGGCACCATCCCCGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATg                           <  1:108043/249‑1 (MQ=255)
                                                                                                                                                                                                                                tcCGGCAACAACACCGACGCCATCGGGCGCACGCTCTCCGGGCAGCCCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCATGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGGACATGCAGGAGCTGTCACGGATTGCGGAATGCTATGTCACCGCGCACCCGAACGACGGGCTACCCAACGCCTTTGGTGAGAACGCTCTCGACGCCGaacagctggcacaac                    >  2:425691/1‑233 (MQ=255)
                                                                                                                                                                                                                                  cGGCACCATCACCGACGCCTTCGGGCGCACGCTCTCCGGGCATACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCGGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATAc              <  2:465986/251‑1 (MQ=255)
                                                                                                                                                                                                                                        cATCGCCGACGCCTTCGGGCGCCCGCCCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAAt        <  2:453601/251‑1 (MQ=255)
                                                                                                                                                                                                                                              cGACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGCCGCCGACACGAGGGCAAAACCGATACGTGAATGGgcgc  >  1:388736/1‑251 (MQ=255)
                                                                                                                                                                                                                                               gACGCCTTCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGgcgc  >  2:125809/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
ACCGCACGGCGTCTATTTCTCCGGACGTCAACGATCCGGCATTTCGTAATATCACTTTTGACGGGCTGGTGGCGGCTTATCGAGAGTCCACCAAAGCGCTGGTGGAAGGTGGCGCGGATCTGATCCTGATTGAAACCGTTTTCGACACCCTTAACGCCAAAGCGGCGGTATTTGCGGTGAAAACGGAGTTTGAAGCGCTGGGCGTTGAGCTGCCGATTATGATCTCCGGCACCATCACCGACGCCTCCGGGCGCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGC  >  NC_000913/4224237‑4224725

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