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%	98.8 / NA	27	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 (12/15);  total (12/15)

TATTCACTGTTATTAGGACAATTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTT  >  NC_000913/1196629‑1197055                                                                                                                                                                                                                                         |                                                                                                                                                                                                   tATTCACTGTTATTAGGACAATTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgatt                                                                                                                                                                                   >  1:160902‑M1/1‑238 (MQ=255)                  aCAATGATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTAGCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTGAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAgcatcaacgattaaactgcgcttcggcgca                                                                                                                                                                 <  2:25849‑M1/251‑31 (MQ=255)                             ttACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGCCTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgttt                                                                                                                                                       >  1:98639‑M1/1‑210 (MQ=255)                                    tCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcag                                                                                                                                                                >  2:559915‑M1/1‑203 (MQ=255)                                    tCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCGTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcag                                                                                                                                                                <  1:559915‑M1/234‑32 (MQ=255)                                     ccTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgt                                                                                                                                              >  1:139242‑M1/1‑202 (MQ=255)                                            cccGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagg                                                                                                                                       <  1:28763‑M1/251‑57 (MQ=255)                                                cTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaa                                                                                                                                    >  1:209932‑M1/1‑191 (MQ=255)                                                            tAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggca                                                                                                                        >  1:272370‑M1/1‑179 (MQ=255)                                                                     ggTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagccc                                                                                                              >  1:625017‑M1/1‑170 (MQ=255)                                                                         ttGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgcc                                                                                                           >  2:638388‑M1/1‑166 (MQ=255)                                                                             aTCAACTCCTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagtta                                                                                                      <  1:13928‑M1/251‑90 (MQ=255)                                                                                       cACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcggga                                                                                              <  1:251126‑M1/249‑98 (MQ=255)                                                                                                 ccgCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatccta                                                                                  <  1:21628‑M1/250‑110 (MQ=255)                                                                                                     cAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcctcctcttgg                                                                              >  1:28290‑M1/1‑138 (MQ=255)                                                                                                            cAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGGCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaa                                                                       <  1:603370‑M1/251‑121 (MQ=255)                                                                                                                   cAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaaactcataccttgtgggaacaactaaaa                                                                 >  2:223665‑M1/1‑124 (MQ=255)                                                                                                                       cATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaa                                                            <  1:628371‑M1/251‑132 (MQ=255)                                                                                                                          tGGGTGCAATCACCTTAGCGTTCAGGGACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggaaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaat                                                          <  2:303767‑M1/250‑134 (MQ=255)                                                                                                                                                cAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggttatcatctcccagt                                    <  2:39769‑M1/250‑156 (MQ=255)                                                                                                                                                   gTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggttatcatctcccagtatat                                <  2:28290‑M1/251‑160 (MQ=255)                                                                                                                                                          cGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggttatcatcgcccagtatatccatact                         <  2:98639‑M1/251‑167 (MQ=255)                                                                                                                                                                gTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTGTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggttatcatctcccagtatatccatactaacaat                   <  2:231316‑M1/251‑173 (MQ=255)                                                                                                                                                                               aGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaaactcctcccatgggagaacaaagaataaataaaatttggggttttcctctccccgtttatccatactaacaataaagttattttctc     >  1:468714‑M1/1‑64 (MQ=255)                                                                                                                                                                                gcgcGGCCCCTTCACATACATCTTTAGTACTGAGAATGTTTAACCTAGGGAATTATTGTCTTAccatcaacgattaaactgcgtttcggcgcagtttgcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggttatcatctcccagtatatccatactaacaataaggttatttactcaa   <  2:80287‑M1/251‑189 (MQ=255)                                                                                                                                                                                 cgcgGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggggaacaaagaataatggaaatgcgggggtatcatctccccgtatatccctactaacaaaaagggtatttactcaac  >  1:178887‑M1/1‑62 (MQ=255)                                                                                                                                                                                                                   cTGTTTAACCTAGGGAATTATTGTCGTAccatcaacgattaaactgcgcttcggcgcagttttcgtttacaggatgttgaaagggaaaattctggggcaaaaaaagcccgccagttacggcgggaaacctcatcctatgggagaacaatgaataatgaaattgcggggttatcatcgcccagtatatccatactaaca                     <  2:339707‑M1/198‑171 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                   TATTCACTGTTATTAGGACAATTATGAATTACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTT  >  NC_000913/1196629‑1197055

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