Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	235,589	A→G	R138R (CGT→CGC)	rnhA ←	ribonuclease HI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	235,589	0	A	G	100.0%	49.6 / NA	16	R138R (CGT→CGC)	rnhA	ribonuclease HI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (7/9);  total (7/9)

AATCCGATGAAACAGAGTAAAAACAAGCCACGAATTCGCCAGGCGGTTGGAGCCACCCGGCAATGTCGTAAACCACAGGCTTAAACTTCAACTTGGTAGCCTGTATCTTCCAGTGTGGGATTCATCGCCGCGGCACGAGCCAGTTCATCACAGCGTTCGTTTTCCGGGTGTCCGGCATGGCCTTTAACCCATTCCCATTTGATTTGATGCTGCCCCAATGCAGCATCAAGACGTTGCCAGAGATCGACATTTT  >  NC_000913/235455‑235707                                                                                                                                       |                                                                                                                       aaTCCGATGAAACAGAGTAAAAACAAGCCACGAATTCGCCAGGCGGTTGGAGCCACCCGGCAATGTCGTAAACCACAGGCTTAAACTTCAACTTGGTAGCCTGTATCTTCCAGTGTGGGATTCATCGCCGCGGCGCGAGc                                                                                                                   <  1:259636/140‑1 (MQ=255)     cGATGAAACAGAGTAAAAACAAGCCACGAATTCGCCAGGCGGTTGGAGCCACCCGGCAATGTCGTAAACCACAGGCTTAAACTTCAACTTGGTAGCCTGTATCTTCCAGTGTGGGATTCATCGCCGCGGCGCGAGCCAGt                                                                                                               <  1:108826/140‑1 (MQ=255)       aTGAAACAGAGTAAAAACAAGCCACGAATTCGCCAGGCGGTTGGAGCCACCCGGCAATGTCGTAAACCACAGGCTTAAACTTCAACTTGGTAGCCTGTATCTTCCAGTGTGGGATTCATCGCCGCGGCGCGAGCCAGTtc                                                                                                             >  1:350164/1‑140 (MQ=255)              agagTAAAAACAAGCCACGAATTCGCCAGGCGGTTGGAGCCACCCGGCAATGTCGTAAACCACAGGCTTAAACTTCAACTTGGTAGCCTGTATCTTCCAGTGTGGGATTCATCGCCGCGGCGCGAGCCAGTTCATCACAg                                                                                                      <  2:350164/140‑1 (MQ=255)                  aaaaaaaaaGCCACGAATCCGCCAGGCGGTTGGAGCCACCCGGCAATGTCGTAAACCCCAGGCTAAAACTCCAACTTGGTAGCCTGTATTTTCCAGTGTGGGATTCATCGCCGCGGCGCGAGCCAGTTCATCACAGcgtt                                                                                                  <  2:350466/139‑1 (MQ=255)                            ccACGAATTCGCCAGGCGGTTGGAGCCACCCGGCAATGTCGTAAACCACAGGCTTAAACTTCAACTTGGTAGCCTGTATGTTCCAGTGTTGGATTCATCGCCGCGGCGCGAGCCAGTTCATCACAGCGTTCGTTTTCCgg                                                                                        >  2:318059/1‑140 (MQ=255)                                            cGGTTGGAGCCACCCGGCAATGTCGTAAACCACAGGCTTAAACTTCAACTTGGTAGCCTGTATCTTCCAGTGTGGGATTCATCGCCGCGGCGCGAGCCAGTTCATCACAGCGTTCGTTTTCCGGGTGTCCGGCa                                                                              <  1:271388/134‑1 (MQ=255)                                            cGGTTGGAGCCACCCGGCAATGTCGTAAACCACAGGCTTAAACTTCAACTTGGTAGCCTGTATCTTCCAGTGTGGGATTCATCGCCGCGGCGCGAGCCAGTTCATCACAGCGTTCGTTTTCCGGGTGTCCGGCa                                                                              >  2:271388/1‑134 (MQ=255)                                                 ggAGCCACCCGGCAATGTCGTAAACCACAGGCTTAAACTTCAACTTGGTAGCCTGTATCTTCCAGTGTGGGATTCATCGCCGCGGCGCGAGCCAGTTCATCACAGCGTTCGTTTTCCGGGTGTCCGGCATGGCCTTTAAc                                                                   <  2:127943/140‑1 (MQ=255)                                                         ccGGCAATGTCGTAAACCACAGGCTTAAACTTCAACTTGGTAGCCTGTATCTTCCAGTGTGGGATTCATCGCCGCGGCGCGAGCCAGTTCATCACAGCGTTCGTTTTCCGGGTGTCCGGCATGGCCTTTAAcccattccc                                                           >  1:69301/1‑140 (MQ=255)                                                         ccGGCAATGTCGTAAACCACAGGCTTAAACTTCAACTTGGTAGCCTGTATCTTCCAGTGTGGGATTCATCGCCGCGGCGCGAGCCAGTTCATCACAGCGTTCGTTTTCCGGGTGTCCGGCATGGCCTTTAAcccattccc                                                           <  2:27155/140‑1 (MQ=255)                                                                       aaCCACAGGCTTAAACTTCAACTTGGTAGCCTGTATCTTCCAGTGTGGGATTCATCGCCGCGGCGCGAGCCAGTTCATCACAGCGTTCGTTTTCCGGGTGTCCGGCATGGCCTTTAACCCATTCCCAGGTGATTTGAtgc                                             <  2:69301/140‑1 (MQ=255)                                                                                  tAAACTTCAACTTGGTAGCCTGTATCTTCCAGTGTGGGATTCATCGCCGCGGCGCGAGCCAGTTCATCACAGCGTTCGTTTTCCGGGTGTCCGGCATGGCCTTTAACCCATTCCCATTTGATTTGATGCTGCCCCAATgc                                  <  1:66122/140‑1 (MQ=255)                                                                                         cAACTTGGTAGCCTGTATCTTCCAGTGTGGGATTCATCGCCGCGGCGCGAGCCAGTTCATCACAGCGTTCGTTTTCCGGGTGTCCGGCATGGCCTTTAACCCATTCCCATTTGATTTGATGCTGccc                                        >  2:138777/1‑127 (MQ=255)                                                                                         cAACTTGGTAGCCTGTATCTTCCAGTGTGGGATTCATCGCCGCGGCGCGAGCCAGTTCATCACAGCGTTCGTTTTCCGGGTGTCCGGCATGGCCTTTAACCCATTCCCATTTGATTTGATGCTGccc                                        <  1:138777/127‑1 (MQ=255)                                                                                                      tGTATCTTCCAGTGTGGGATTCATCGCCGCGGCGCGAGCCAGTTCATCACAGCGTTCGTTTTCCGGGTGTCCGGCATGGCCTTTAACCCATTCCCATTTGATTTGATGCTGCCCCAATGCAGCATCAAGACGTTGCCaga              >  1:258383/1‑140 (MQ=255)                                                                                                                  tgtgGGATTCATCGCCGCGGCGCGAGCCAGTTCATCACAGCGTTCGTTTTCCGGGTGTCCGGCATGGCCTTTAACCCATTCCCATTTGATTTGATGCTGCCCCAATGCAGCATCAAGACGTTGCCAGAGATCGACAtttt  >  2:262863/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                       AATCCGATGAAACAGAGTAAAAACAAGCCACGAATTCGCCAGGCGGTTGGAGCCACCCGGCAATGTCGTAAACCACAGGCTTAAACTTCAACTTGGTAGCCTGTATCTTCCAGTGTGGGATTCATCGCCGCGGCACGAGCCAGTTCATCACAGCGTTCGTTTTCCGGGTGTCCGGCATGGCCTTTAACCCATTCCCATTTGATTTGATGCTGCCCCAATGCAGCATCAAGACGTTGCCAGAGATCGACATTTT  >  NC_000913/235455‑235707

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