Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,298,035 T→C E35G (GAA→GGA)  yraQ ← putative inner membrane permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,298,0350TC100.0% 117.3 / NA 34E35G (GAA→GGA) yraQputative inner membrane permease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/18);  total (16/18)

CGAGGAAATAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTTCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTG  >  NC_000913/3297913‑3298169
                                                                                                                          |                                                                                                                                      
cgAGGAAATAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCa                                                                                                                        <  1:144243/139‑1 (MQ=255)
cgAGGAAATAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCa                                                                                                                        >  2:144243/1‑139 (MQ=255)
 gAGGAAATAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAGGGCCCCCATTTCACAAAcaaa                                                                                                                       >  2:38047/1‑136 (MQ=255)
    gAAATAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCCATACTGGTGGTTTTGGCAGCAGTAAACGCCTTGCCGTAGTAAGGTCCCCATTTCACATACCAgaga                                                                                                                    <  1:38047/139‑1 (MQ=255)
    gAAATAAATCCCCGCTTAATCCAACGCCGCCTGCCATGGTTTAGCATCCGCCTGCCCAAGTATATATTTACCGATATTGTCGGTTTCGGCACCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAcaga                                                                                                                    <  2:586419/139‑1 (MQ=255)
               tcGCGTAATCAAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGAc                                                                                                         >  2:616686/1‑139 (MQ=255)
                cgcgTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCa                                                                                                                                    <  2:817632/111‑1 (MQ=255)
                cgcgTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCa                                                                                                                                    >  1:817632/1‑111 (MQ=255)
                                                cATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCa                                                                                                                                    >  1:386118/1‑79 (MQ=255)
                                                cATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCa                                                                                                                                    <  2:386118/79‑1 (MQ=255)
                                                    cGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATg                                                                    >  2:420141/1‑139 (MQ=255)
                                                        tGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCg                                                                <  2:634068/139‑1 (MQ=255)
                                                                 aTAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCt                                                       <  1:223823/139‑1 (MQ=255)
                                                                               aCTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGa                                         <  1:420141/139‑1 (MQ=255)
                                                                               aCTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGa                                         >  2:107743/1‑139 (MQ=255)
                                                                               aCTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGa                                         >  2:220882/1‑139 (MQ=255)
                                                                                 tgtgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGAcc                                       <  1:616686/139‑1 (MQ=255)
                                                                                 tgtgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGAcc                                       >  2:813837/1‑139 (MQ=255)
                                                                                 tgtgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGAcc                                       >  1:184662/1‑139 (MQ=255)
                                                                                    gggTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCAtt                                                                                                                                  <  1:320651/45‑1 (MQ=255)
                                                                                    gggTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCAtt                                                                                                                                  >  2:320651/1‑45 (MQ=255)
                                                                                       tttCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCt                                                       >  1:708003/1‑117 (MQ=255)
                                                                                       tttCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCt                                                       <  2:708003/117‑1 (MQ=255)
                                                                                        ttCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCgg                                                                                     <  1:12491/86‑1 (MQ=255)
                                                                                        ttCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCgg                                                                                     >  2:12491/1‑86 (MQ=255)
                                                                                         tCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAgcgc                                                                                       <  2:659467/83‑1 (MQ=255)
                                                                                         tCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAgcgc                                                                                       >  1:659467/1‑83 (MQ=255)
                                                                                                      aCGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGaa                  <  2:184662/139‑1 (MQ=255)
                                                                                                      aCGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGaa                  <  1:220882/139‑1 (MQ=255)
                                                                                                      aCGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGaa                  <  1:813837/139‑1 (MQ=255)
                                                                                                         cTTTGCCGTAGTAAGGTCCCCATATCACATACCAGAGACCGGCAATGACAACGAGAAAGACAAGCGCGGGTTTCCACCACTGCATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATaa               <  2:168213/139‑1 (MQ=255)
                                                                                                           ttGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAAAGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGt             >  1:789573/1‑139 (MQ=255)
                                                                                                             gCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTa           <  2:153826/139‑1 (MQ=255)
                                                                                                                      aGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTg  >  1:817348/1‑139 (MQ=255)
                                                                                                                          |                                                                                                                                      
CGAGGAAATAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTTCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTG  >  NC_000913/3297913‑3298169

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