Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,037,599	C→T	A170T (GCA→ACA)	recJ ←	ssDNA‑specific exonuclease RecJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,037,599	0	C	T	100.0%	44.9 / NA	20	A170T (GCA→ACA)	recJ	ssDNA‑specific exonuclease RecJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/10);  total (10/10)

CAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTAGGGTTAATGATCGCTTCCGCTGCGGGTAATGTGTCGCCTGGCAAATGGTGATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTAC  >  NC_000913/3037366‑3037788                                                                                                                                                                                                                                          |                                                                                                                                                                                              cAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATTTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTAGGGTTAATGATCGCTTCCGCTGTGGGTAATGTGTCGCCTg                                                                                                                                                                              >  1:37324/1‑251 (MQ=255)      ngCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGATCTGCCAGGTTAGGCATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCTCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTAGGGTTAATGATCGCTTCCGCTGTGGGTAATGTGTCGCCTGGCaa                                                                                                                                                                          >  1:43679/2‑250 (MQ=255)                      cAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTAGGGTTAATGATCGCTTCCGCTGTGGGTAATGTGTCGCCAGGCAAATGGTGATCGGTAACa                                                                                                                                                          >  1:465782/1‑250 (MQ=255)                              cgccgTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCATTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTTCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCGGCCAGTGATTTCGAGGGGAAATTTCAGTCGCGCAAGTTAGGGTTAATGATCGCTTCCGCTGTGGGTAATGTGTCGCCTGGCAAATGGTGATCGGTAACAATAACCggg                                                                                                                                                 <  1:293379/248‑1 (MQ=255)                                                 cAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATACAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGACGCGCAAGTTAGGGTTAATGATCGCTTCCGCTGTGGGTAATGTGTCGCCTgg                                                                                                                                                                             >  1:36065/1‑204 (MQ=255)                                                       gACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGAGCCAAAAGGTGCGCAGCGCAAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTAGGGATAATGATCGCTTCCGCTGTGGGTAATGTGTCGCCTGGCAAATGGTGATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCa                                                                                                                         >  1:193388/1‑250 (MQ=255)                                                         ccAGATCCAGCAGTGCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGAGCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGTAAATTACAGTCGCGCAAGTTAGGGTTAATGATCGCTTCCGCTGTGGGTAATGTGTCGCCTGGCAAATGGTGATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAAcc                                                                                                                      <  1:284051/251‑1 (MQ=255)                                                                    agTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTAGGGTTAATGATCGCTTCCGCTGTGGGTAATGTGTCGCCTGGCAAATGGTGATCGGTAACAATAACCGGGATGCCCAAcgag                                                                                                                                     <  1:425921/225‑1 (MQ=255)                                                                                            aTGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACCCCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTAGGGTTAATGATCGCTTCCGCTGTGGGTAATGTGTCGCCTGGCAAATGGTGATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCTTGCTCAACCCCCGCATGggagga                                                                                                        <  1:412332/230‑1 (MQ=255)                                                                                                  cGCTCATCAAACCAGCCCTGAGCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGGGATTTCGACGGGAAATTACAGTCGCGCAAGTGATGGTTAATGATCGCTTCCGCTGTGGGGAATGTGTCGCCTGGCAAATGGTGATCGGTAACAATAACCGGGATGCCCAAcgagag                                                                                                                                   <  1:472367/197‑3 (MQ=255)                                                                                                             ccAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTGAGGGTTAATGATCGCTTCCGCTGTGGGTAATGTGTCGCCTGGCAAATGGTGATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCGCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTg                                                                                     <  1:218631/232‑1 (MQ=255)                                                                                                                            cAAAAAGGTGCGCAGCGCCAGCATCAGATCAACCGCCACAACCACGCCTGCCAGTGATTTCGACGGGACATTACAGTAGCGCAAGTTGGGGTTCATGATCGCTTCCGCTGTGGGTAATGTGTCGCCTGGCAAATGGTGATCGGTACCACTAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTgct                                                   >  1:408146/1‑249 (MQ=255)                                                                                                                                                        tAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTAGGGTTAATGATCGCTTCCGCTGTGGGTAATGTGTCGc                                                                                                                                                                                 >  1:6484/1‑97 (MQ=255)                                                                                                                                                        naaaaCGNNACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTAGGGTTAATGATCGCTTCCGCTGTGGGTAATGTGTCGc                                                                                                                                                                                 <  2:6484/96‑1 (MQ=255)                                                                                                                                                                           cTGCCAGTGATTTGGACGGGAAATGACAGTCGCGCAAGTTAGGGTTAATGATCGCTTCCGCTGTGGGTAATGTTTCGCCTGGCAAATGGTGATCGGTAACAATAACCGGGATGCCCAAcgan                                                                                                                                     <  1:55562/122‑2 (MQ=255)                                                                                                                                                                             gCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTAGGGTTAATGAGCGCTTCCGCTGTGGGTAATGTGTCGCCTGGCAAATGGTGATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATggg                                                                                                            <  1:356225/145‑1 (MQ=255)                                                                                                                                                                             gCCAGTGANTNNNANNGGAAATTACAGTCGCGCAAGTTAGGGTTAATGATCGCTTCCGCTNTGGGTAATGTGTCGCNTGNCAAANGGTNATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATggg                                                                                                            >  2:356225/1‑145 (MQ=255)                                                                                                                                                                             gCCAGGGATTTCGACGGGAAATTACAGTCGCGCAAGTTAGGGTTAATGATCTTTTCCGCTGTGGGGAATGTGTCGCCTGGCAAATGGTGATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATTCCGTTATCCGCCGTGACAATGAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTGCCGGGCTTAAGCCGTTACCGTCTTCGAAACGGTTTGGTAc  <  1:472392/251‑1 (MQ=255)                                                                                                                                                                              ccAGTGATTTCGACGGGAAATTATAGTCGCGCAAGTTCGGGTTAATGATCGCTTCCGCTGTGGGTAATGTGTCGCCTGGCAAATGGTGATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGATATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATg                                                           >  1:486644/1‑193 (MQ=255)                                                                                                                                                                              ccAGTGATTTCGACGGGAAATTATAGTCGCGCAAGTTAGGGTTAATGATCGCTTCCGCTGTGGGTAATGTGTCGCCTGGAACATGGTGATCGGGAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCAGCATGGGAGGATATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATg                                                           >  1:486638/1‑193 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                              CAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTAGGGTTAATGATCGCTTCCGCTGCGGGTAATGTGTCGCCTGGCAAATGGTGATCGGTAACAATAACCGGGATGCCCAACGAGCGAGCGTGCTCAACCCCCGCATGGGAGGAAATACCGTTATCCACCGTGACAATTAACTGCGCGCCACGGGCATGGGCCTGATCGACCACTTCCGGGCTTAAGCCGTAACCGTCTTCGAAACGGTTTGGTAC  >  NC_000913/3037366‑3037788

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