Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,652,331 T→C intergenic (+2346/+397) ydfD → / ← ynfP lysis protein/protein YnfP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,652,3310TC100.0% 23.1 / NA 9intergenic (+2346/+397)ydfD/ynfPlysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (5/4);  total (5/4)

GCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGC  >  NC_000913/1652126‑1652539
                                                                                                                                                                                                             |                                                                                                                                                                                                                
ncTGTGTGCCCTTGGATGGGAGGAAATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACCCAGGCGGGAACCACTCCGGTCATTCATCTTATTCAGCCAGCCATCTACGCTCTCCGGAGTCAGCTGAACTTAACGAGGCTGAGCAAAGAGCATATCATTGATGTTCACCTCCGCGCGTCTGGCAGAACAGACAAACAAAAATGCCCCTTTGttt                                                                                                                                                                     >  1:68952/2‑251 (MQ=255)
                               cTGAAACATGGAACAATGATGATAAGGAGAAATTTAGCACAAGACCGTTTCACGGTACCAAAAACACAGGAGGGAACCAATCGGGTCGTTCATCTTATTAAGTCAGCAATCGACGCTCTCCGGAGGAAGATGACATTACAGAGGCTGAGCAAAGAGCCTATTCTTGATGTTAAACTCAGAGAGTATGGCAGAACAGAAAAAAAAAAATGCACCTTTGTTTTTCAACATGAAATGTCCGCGAGAGGaaaaa                                                                                                                                       >  1:219847/1‑250 (MQ=255)
                                        ggAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTT‑AACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAc                                                                                                                            <  1:150171/251‑1 (MQ=255)
                                                        ggAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCCCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCCTATCATTGATGTTCACCTCAGAGAGTATGGCAGCACAGACACACACACATGCACCTTTGTTTTTAAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAACATTTTACCGTTGAc                                                                                                             >  1:196391/1‑251 (MQ=255)
                                                                                                      gggAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCGCTCCGGAGTCAGATGACAGTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTcc                                                               <  1:101487/251‑1 (MQ=255)
                                                                                                        gAACCAATCGGGTCATTCATCTTTTTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGATCAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTTACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCga                                                             <  1:60612/251‑2 (MQ=255)
                                                                                                                         catCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGAAAAGAGCATATCATTGATGTTAACCTCAGCGAGTCTGGCAGAACAGAAAAACAAAAATGCACCTTTGttttt                                                                                                                                                                   >  1:120430/1‑132 (MQ=255)
                                                                                                                                              gACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGCGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTACAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATg                           >  1:164857/1‑247 (MQ=255)
                                                                                                                                                      ccGGANGCAGANGNNATTAACGAGACTGAGCAAANAGCATATCATTGATGNTCACNTCAGANAGNa                                                                                                                                                                                                        >  2:121167/1‑66 (MQ=255)
                                                                                                                                                                   aCATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAAGTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGAGAAAACGTGCCGGACTCCGCCATCGAAAAGCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGc  <  1:121167/251‑1 (MQ=255)
                                                                                                                                                                                                             |                                                                                                                                                                                                                
GCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGC  >  NC_000913/1652126‑1652539

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