Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 335,016 C→T G3R (GGA→AGA)  yahC ← uncharacterized protein YahC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913335,0160CT100.0% 63.7 / NA 20G3R (GGA→AGA) yahCuncharacterized protein YahC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/12);  total (8/12)

AACCAGCAGTCAGCGCGGCGAAAAACGCCCACACCATCCCGGAATAATTAACCGCCAGACTTCTGATAAACCCCGATTTACCGCCGCCTGCGGCATAAAAACTGCAAAAGGCGACGAAGCCAATGGTTCCTAATAATTCCCATCCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGACGGTGACACCCGTTGCCGTGAGTCCGTTCATCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCAC  >  NC_000913/334791‑335260
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
aaCCAGCAGTCAGCGCGGCGAAAAACTCCCACACCATCCCGGAATAATTAACCGCCAGACTTCTGATAAACCCCGATTTACCGCCGCCTGCGGCATAAAAACTGCAAAAGGCGACGAAGCCAATGGTTCCTAATAATTCCCATCCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGACGGTGACACCCGTTGCCGTGAGTCTGTTCATCTTGTTCCCCTCTGGTTAg                                                                                                                                                                                                                             <  2:153559/251‑1 (MQ=255)
               cgacgAAAAACGCCACCCCCCTCCCGGAATAATTAACCGCCAGACTTCTGATAAAACCCGATTTACCTCCACCTGCGGCATAAAAACTGCAAAAGGCGACGAAGCCAATGGTTCCTAATAATTCCCATCCTTGAGAGAGCCCTACATGTTAAGACACTAACTGCCAAAGCTCTGCACAAATACCGACGGTGTCACCAGTTGCCGTGAGTCTGTTCATCGTGTTACCTCCTTGTTAGACATGGACGATAGCt                                                                                                                                                                                                              <  2:162937/248‑1 (MQ=255)
                             caccccATCCCGGAACAATTAACCGCCAGACTTCTGATAAACCCCGATTTCCCGCCGCCTGCGGCATAAAAACTGCAAAAGGCGACGAAGCCAATGGTTCCTAATAATCCCCATCCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGACGGTGACACCCTTTGCCGTGAGTCTGTTCATCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGc                                                                                                                                                                                                <  2:106507/247‑1 (MQ=255)
                                acCATCCCGGAATAATTAAACGCCAGACTTCTGATAAACCACGATTTACCGCCGCCTGCGGCATAAAAACTGCAAAAAGCGACGAAGCCAATGGTTCCTACTAATTCCCATCCTTGAGAGAGCCCGACATGTGAAGACACTATCTGCCAATGCCCTGCACAAATACCGACGGTGACACCCGTTGCCGTGAGTCTGTTCATCTTGTTCACCTCTGGTTAGACATCGACGAAAGCTACTCTTActtgatc                                                                                                                                                                                                >  2:75153/1‑242 (MQ=255)
                                           ataatTAACCGCCAGACTTCTGATAAACCCCGATTTACCGCCGCCTGCGGCATAAAAACTGCAAAAGGCGACGAAGCCAATGGTTCCTAATAATTCCCATCCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGACGGTGACACCCGTTGCCGTGAGTCTGTTCATCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTACTCTTACAGGTGCAGACTATAGAGaaa                                                                                                                                                                                  <  2:123204/251‑1 (MQ=255)
                                              atcagcggCCAGACTTCTGATAAACCCCGATTTACCGCCGCCTGCGGCATAAAAACTGCAAAAGGCGACGAAGCGAATGGTGCCTAATAAGACCCAACCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGTCGGTGACTCCCGTTGCCGTGAGTCTGTTCATCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGAGTATAGAGAAATaa                                                                                                                                                                               <  1:48657/244‑1 (MQ=255)
                                                    cGCCAGACTTCTGATAAACCCCGATTTACCGCCGCCTGCGGCATAAAAACTGCAAAAGGCGACGAAGCCAATGGTTCCTAATAATTCCCATCCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGACGGTGACACCCGTTGCCGTGAGTCTGTTCATCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTg                                                                                                                                                                         >  1:219895/1‑251 (MQ=255)
                                                    cGCCAGACTTCTGATAAACCCCGATTTACCGCCGCCGGCGGCATAAAAACTGCAAAAGGCGACGAAGCCAATGGTTCCTAATAATTCCCATCCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGACGGTGACACCCGTTGCCGTGAGTCTGTTCATCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTg                                                                                                                                                                         <  2:200398/251‑1 (MQ=255)
                                                     gCCAGACTTCTGATAAACCCCGATTTACCGCCGCCTGCGGCATAAAAACTGCAAAAGGCGACGAAGCCAATGGTTCCTAATAATTCCCATCCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGACGGTGACACCCGTTGCCGTGAGTCTGTTCATCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTg                                                                                                                                                                         >  1:157381/1‑250 (MQ=255)
                                                                             ttACCGCCGCCTGCGGCATAAAAACTGCAAAAGGCGACGAAGCCAATGGTTCCTAATAATTCCCATCCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGACGGTGACACCCGTTGCCGTGAGTCTGTTCATCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCa                                                                                                                                                 <  2:219895/250‑1 (MQ=255)
                                                                                         gCGGCATAAAAACTGCAACAGGCGACGAAGCCAATGGTTCATAATAATTCCCATCCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGACGGAGACCCCCGTTGCCGAGAGTCTGTTCCTCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTCCTCTTACATGTGCAGACTATAGAGAAATAAAACCTGCTTGAGTGAAAGGGTATCGGGTCAAAGAAACAAatat                                                                                                                                    >  2:107675/1‑251 (MQ=255)
                                                                                                       gCAAAAGGCGACGAAGCCAATGGTTCCTAATAATTCCCATCCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGACGGTGACACCCGTTGCCGTGAGTCTGTTCATCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGcac                                                                                                                       <  2:237908/250‑1 (MQ=255)
                                                                                                         aaaaGGCGACGAAGCCAATGGTTCCTAATAATTCCCCTCCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGACGGTGACACCCGTTGCCGTGAGTCTGTTCATCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACt                                                                                                                    >  2:100361/1‑251 (MQ=255)
                                                                                                                                          cccAACCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGACGGTGACACCCGTTGCCGTGAGTCTGTTCATCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACaaaaa                                                                                   <  2:157381/251‑1 (MQ=255)
                                                                                                                                                           ccGACATGGGAAGACACTACCTGCCAAAGCCCTGCACAAATACCGACGGTGACACCCGTTGCCGTGAGTATGTTCATCTTGTTCCCCTCTGGTTAGACATGGACGACAGCTACTCTTACATGTGCATACTATAGAGTAATAAAACCTGATTGAGTGACATGGTATCGGTTCAAAGAAACAAATATTGCACTACCGCACACTGCGACACGATTGTTGAATATTATTcatgacataggtgagataaggcgaac                                                                  >  2:67333/1‑227 (MQ=255)
                                                                                                                                                              aCATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGACGGTGACACCCGTTGCCGTGAGTCTGTTCATCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTa                                                                >  1:141203/1‑250 (MQ=255)
                                                                                                                                                                aggtgAAGACACTAACTGCCAAATCCCTGCACAAATACCGACGGTGACAGCCGTTGCCGTGAGTCTGTTCAGCTTGTTCCCCTCTTGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGTAAAGATTGTGGAATATGATACAACATAAATGTTAGTTACGGTGAATTAAAt                                                             <  2:234535/249‑1 (MQ=255)
                                                                                                                                                                                    aaaGCCCTGCACAAATACCGACGGTGACACCCGTTGCCGTGAGTCTGTTCACCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATc                                          >  1:236662/1‑250 (MQ=255)
                                                                                                                                                                                              acaAATACCGACGGTGACACCCGTTGCCGTGAGTCTGTTCATCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGtt                               <  2:127537/251‑1 (MQ=255)
                                                                                                                                                                                                                            gAGTCTGTTCATCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCcac  <  1:102047/250‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
AACCAGCAGTCAGCGCGGCGAAAAACGCCCACACCATCCCGGAATAATTAACCGCCAGACTTCTGATAAACCCCGATTTACCGCCGCCTGCGGCATAAAAACTGCAAAAGGCGACGAAGCCAATGGTTCCTAATAATTCCCATCCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGACGGTGACACCCGTTGCCGTGAGTCCGTTCATCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCAC  >  NC_000913/334791‑335260

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