Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,088,010 2 bp→TT intergenic (+150/‑273) metK → / → galP methionine adenosyltransferase/galactose:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,088,0100CT100.0% 101.9 / NA 29intergenic (+150/‑274)metK/galPmethionine adenosyltransferase/galactose:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/13);  total (16/13)
*NC_0009133,088,0110CT100.0% 100.8 / NA 29intergenic (+151/‑273)metK/galPmethionine adenosyltransferase/galactose:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/13);  total (16/13)

AACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTA‑ATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATC  >  NC_000913/3087767‑3088245
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                           
aaCCGAAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTac                                                                                                                                                                                                                                       >  2:30446/1‑251 (MQ=255)
     cagcaTACGGTCAATTTGGTCGTGAACATTTCCCGTGGGAATAAACCGACAAAGCGCAGCTGCTGAGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGTCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCCGATTATCATCACAACTGAAATTGATTACacca                                                                                                                                                                                                                                   >  2:241148/1‑250 (MQ=255)
              gTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAg                                                                                                                                                                                                                         <  1:102850/251‑1 (MQ=255)
                  ctttgcacGTGCACATTTCACGCGGGAAAAACCCGGCACCTCGCAGCTGCTGCGCAATGCTGCAGGTCTGAAGTACTCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGACGCCTCCGGCAAAAAAACCGCCCGCAGAATAACATCATACATCCTGCTCACGGTTCACCGCAGATTATCATAACAACTGAAATTGATTCCACCAACCACAACAGACa                                                                                                                                                                                                                      >  2:63278/8‑250 (MQ=255)
                     tGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCCCCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAg                                                                                                                                                                                                                   <  1:289140/250‑1 (MQ=255)
                      ggTCGTGAACATTTCCCGTGGGAAAAAACCGACAAACCGCAGCTGCTGCGCGATGCTGCCGGTCTGAATTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGAGAGGCGTTCACGCCGCATCCTGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGa                                                                                                                                                                                                                  <  1:210459/250‑1 (MQ=255)
                         cGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTg                                                                                                                                                                                                              >  1:5678/1‑251 (MQ=255)
                                    cccGGGGGAAAAAACCGCCAAAGCGCAGCTTCTTCGCGATGCGGCCGGTCTGAAGTAATCTTTCTTCCCCTGCGTGCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTAGCATCACAACTGAAATTGATTACACCAACCACAACGGACAAAGAGTTTTAATATTTTCa                                                                                                                                                                                                   <  2:266845/251‑1 (MQ=255)
                                         gggAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCAtatta                                                                                                                                                                                              <  1:85688/251‑1 (MQ=255)
                                                                    tgcGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACAtt                                                                                                                                                                   >  1:205226/1‑251 (MQ=255)
                                                                           tgctgcCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTc                                                                                                                                                            >  1:224242/1‑251 (MQ=255)
                                                                                   gTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAACAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGt                                                                                                                                                    >  1:255447/1‑251 (MQ=255)
                                                                                              tctttcttCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGACTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGAAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAAttac                                                                                                                                           >  2:97544/1‑247 (MQ=255)
                                                                                                 ttcttcACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGa                                                                                                                                        <  1:119873/249‑1 (MQ=255)
                                                                                                 ttcttcACCTGCGTTCAAAGGCCAGCCGCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCa                                                                                                                                          <  2:99889/247‑1 (MQ=255)
                                                                                                      cACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGc                                                                                                                                 >  1:102419/1‑251 (MQ=255)
                                                                                                                  aaGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTA‑ATCTatg                                                                                                                    >  1:102455/1‑251 (MQ=255)
                                                                                                                                  tGGCCTTGTTCTCTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATACCATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACACCAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGATGTTACATTTCTTTTCAGTAAAGACTTAATTGCAGATAACAGCGTTTA‑ATCTATGATGATATAACTCAatt                                                                                                    >  2:257795/1‑251 (MQ=255)
                                                                                                                                      cTTTTTCTTTTGGATAGGCTTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCACAGCTTATCATCACAACTGAAATTGATTACACCAACCACAACAGAACAAGATTTGTACTATTTTCATATTATTATTCGGTTTTCACAGTTATTACCTTTCTTTTCAGTACAGTCTTAATGGCAGATAACAGC‑TTTACATCTATGATGATATAACTCAATTCttt                                                                                                >  2:231105/1‑251 (MQ=255)
                                                                                                                                                 ggATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTA‑ATCTATGATGATATAACTCAATTATTTTCATGCACTTa                                                                                     <  1:425828/251‑1 (MQ=255)
                                                                                                                                                                              aaaaCCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTA‑ATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTa                                                         >  1:327693/1‑250 (MQ=255)
                                                                                                                                                                                aaCCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTA‑ATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGc                                                      >  1:245389/1‑251 (MQ=255)
                                                                                                                                                                                         cacaATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTA‑ATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACt                                             >  2:40445/1‑251 (MQ=255)
                                                                                                                                                                                                    catTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTA‑ATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTg                                  <  1:30446/251‑1 (MQ=255)
                                                                                                                                                                                                             tGATCACGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTA‑ATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGCTTACACTGATGTGATTTGCTTCACATc                         >  2:101042/1‑251 (MQ=255)
                                                                                                                                                                                                                  ccGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTA‑ATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCttttt                    <  1:321603/250‑1 (MQ=255)
                                                                                                                                                                                                                   cGTTTCACCGCAGATTATCATCACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTA‑ATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTa                   <  1:47500/251‑1 (MQ=255)
                                                                                                                                                                                                                                   atcatcACAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTA‑ATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTTATTTTCTTCACAGCTTTTTACGTCGTACTCACCTAt   <  1:159544/251‑1 (MQ=255)
                                                                                                                                                                                                                                     catcaCAACTGAAATTGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTA‑ATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATc  <  1:14409/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                           
AACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTA‑ATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATC  >  NC_000913/3087767‑3088245

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