Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,196,861	A→G	intergenic (+79/+6)	C0293 → / ← ymfD	small regulatory RNA C0293/putative SAM‑dependent methyltransferase YmfD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,196,861	0	A	G	100.0%	42.8 / NA	13	intergenic (+79/+6)	C0293/ymfD	small regulatory RNA C0293/putative SAM‑dependent methyltransferase YmfD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (4/9);  total (4/9)

ATTAGGACAATTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAG  >  NC_000913/1196640‑1197091                                                                                                                                                                                                                              |                                                                                                                                                                                                                                       attaGGACAATTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgctt                                                                                                                                                                                                            <  1:213315‑M1/250‑24 (MQ=255)              tGAATTACCACTCCTTCCACCCGCGCAAATATTGTTAAATTGCCGGTTTTGTATCCACTTCTCACCCGTGACGCGCCAGGGGACAGCCAACATGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGTCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccaacaacgattaaactgcgcgtcggctcagttttc                                                                                                                                                                                               <  2:366449‑M1/250‑37 (MQ=255)               gAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcg                                                                                                                                                                                              <  1:316821‑M1/250‑38 (MQ=255)                          ccTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatg                                                                                                                                                                                   <  1:352840‑M1/250‑49 (MQ=255)                                       cAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaa                                                                                                                                                                       <  1:309646‑M1/249‑61 (MQ=255)                                                      tGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaa                                                                                                                                                       >  1:290569‑M1/1‑174 (MQ=255)                                                       accgGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGAACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaa                                                                                                                                                       >  2:239952‑M1/2‑173 (MQ=255)                                                                   tCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttac                                                                                                                                         >  1:18511‑M1/1‑161 (MQ=255)                                                                                    aCTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcc                                                                                                                        <  2:332692‑M1/251‑108 (MQ=255)                                                                                                                 ggTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgc                                                                                           <  2:294240‑M1/251‑137 (MQ=255)                                                                                                                                            aTGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcggggaacaccatcccattggagaacaatgaataaataaattgcgggggggtcatctcccagtacatcca                                                                 >  1:382222‑M1/1‑88 (MQ=255)                                                                                                                                                                                      acatCTTTAGTACTGAGACGGTTTAACCTGGGGAATTATTGTCTAAccatcaacgattaaactgcgcttcggcgcagttttcgtttccaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggttatcatctcccagtatatccatactaacaataaggttatttactcaaccaggcataaacatttt                      <  2:69191‑M1/251‑206 (MQ=255)                                                                                                                                                                                                           tttAACCTAGGGAATTATTGTCGTAccataaacgattaaactgcgctgcggcgcagttttcgtttacaggatgttgaaaggaaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggttatcatctcccagtatatccatactaacaataaggttatttactcaaccaggcataaacattttgttttgtgcgtgggaacagc  <  2:289151‑M1/250‑226 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                                                       ATTAGGACAATTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAG  >  NC_000913/1196640‑1197091

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