Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,954,158	G→A	L616L (CTG→TTG)	recB ←	exodeoxyribonuclease V subunit RecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,954,158	0	G	A	100.0%	50.2 / NA	20	L616L (CTG→TTG)	recB	exodeoxyribonuclease V subunit RecB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/6);  total (14/6)

CAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCAGAGTTTCAAAAACACTGTCGCGGTTCGAAAGGTAAACGGAAGGGATTTCCAGCAACGTTAAGGCATCGCGCACCTGGGCGGCCTCCTGGCGGCTGCGCACCAGCACACTGATGTCCGAAGCACGCACCGGACGCGCGTCGTCGCCGTT  >  NC_000913/2953940‑2954380                                                                                                                                                                                                                           |                                                                                                                                                                                                                               cAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAACGCACTACGCAGGGTGTTCTCACGTTCGGGCg                                                                                                                                                                                                >  1:374205/1‑251 (MQ=255)    gATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCACCAAGTTTTCAGCAATGTTAACCGCCGAAATACGCGCCCGAAGCATCGGCATAACGACACCATTGCGCCAGATATGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTATTCAGCGTTTCGCTATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAACGCACTACGCAg                                                                                                                                                                                                                    >  1:384186/1‑228 (MQ=255)                            gCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGCTTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAACGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGc                                                                                                                                                                      >  1:288248/1‑250 (MQ=255)                                                 ttCAGCATTGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATACCGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTGCCACATCCCCCGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTGGCCAACGCACTACGCAGGGTTTTCTCACGTGCGGGCGTCATCGCCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCAGAGtt                                                                                                                                                 <  1:250631/250‑1 (MQ=255)                                                     gCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGCTTTGCCGATAACCATCGAACTCTTCGACTACCACATACCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAACGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCAGAGTTTCaa                                                                                                                                             >  1:234119/1‑250 (MQ=255)                                                      cAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATCACGCCACGTTTGCGCCCGATTTGCCGATAACCATCCAACTCTTGCACTACCTAATCCGACTCATGTTCGTCAATGTCCACCGGTGCGATATCCTGCACGCTCATCCATATCATTTACGATGCCCACGCGCAACGCCCGGTGTTCTCACGGTAGTGCGTCAACACCGCCTGCACCAGCCAAAGCATTTCCTGCGCTTccgggatttcgtc                                                                                                                                            >  1:439344/1‑239 (MQ=255)                                                        aTGTTACGCGCCGACATCAGCGCCCGCAGCATCGTCATAACGCCACGTTTTCGCCAGATTGGCCGATAACCAGCGAACTCGTCGGCGTCCACATCCCACGCATTTTCGTCATGGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATCGACGTTGCCAACGCACTACGCAGGGTGTTCTCACGTTCg                                                                                                                                                                                                    <  1:294719/192‑1 (MQ=255)                                                                          aGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCACCCCCATCATTGACGTTGCCAACGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCAGAGTTTCAAAAACACTGTcgc                                                                                                                                 <  1:63901/241‑1 (MQ=255)                                                                              cccGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTATTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAACGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCaacaac                                                                                                                                                                            >  1:218168/1‑194 (MQ=255)                                                                                    gcaTCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAACGCACTACGCAGGGTGTTCTCACGt                                                                                                                                                                                                       >  1:30430/1‑161 (MQ=255)                                                                                           naTAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAACGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCAGAGTTTCAAAAACACTGTCGCGGTTCGAAAg                                                                                                                       >  1:109398/2‑234 (MQ=255)                                                                                                  ccACGTTTGCGCCAGATTTGCCGATANCCATCGNNCTNTTCGACTACCCCATCCCACNCANNNTNNNCNATGTTNAGNNTTTCGATATNCAGCGCGTTCAGCCCCATCATTGACGTTGCCAACGCACTACGCAGGGTGTTCTNACGTTCgggn                                                                                                                                                                                                 >  2:100842/1‑152 (MQ=255)                                                                                                  ccACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAACGCACTACGCAGGGTGTGCTCACGTTCGGGc                                                                                                                                                                                                 <  1:100842/153‑1 (MQ=255)                                                                                                                                         tCGACTACCACACCCCACGCATGTGCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAACGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCAGAGTTTCAAAAACACTGTCGCGGTTCGAAAGGTAAACGGAAGGGATTTCCAGCAACGGTAAGGCAGCGCGCACCTGGGCGGCCTCCTggcggc                                                         <  1:140924/250‑1 (MQ=255)                                                                                                                                                ccACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAACGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCAGAGTTTCAACAACACTGTCGCGGTTCGAAAGGTAAACGGAAGGGATTTCCAGc                                                                                                 >  1:400532/1‑203 (MQ=255)                                                                                                                                                                     cATTTTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGCCGTTGCCAACGCCCTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCAGAGTTTCAAAAACACTGTCGCGGTTCGAAAGGTAAACGGAAGGGATTTCCAGCAACGTTAAGGCATcg                                                                                  <  1:21280/197‑1 (MQ=255)                                                                                                                                                                                gTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAACGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCAGAGTTTCAAAAACACTGTCGCGGTTCGAAAGGTAAACGGAAGGGATTTCCAGCAACGTTAAGGCATCGCGCACCTGGGCGGCCTCCTGGCGGCTGCGCACCAGCACACTGATGTCCGAAGCACGCACCGGAcg                 >  1:11373/1‑251 (MQ=255)                                                                                                                                                                                               gcgTTCAGCCCCATCATTGCCGTTGCCAACGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCAGAGTTTCAAAAACACTGTCGCGGTTCGAAAGGTAAACGGAAGGGATTTCCAGCAACGTTAAGGCATCGCGCACCTGGGCGGCCTCCTGGCGGCTGCGCACCAGCACACTGATGTCCGAAGCACGCACCGGACGCGCGTCGTCGCCGtt  >  1:373381/1‑251 (MQ=255)                                                                                                                                                                                               gcgTTCAGCCCCATCATTGACGTTGCCAACGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCAGAGTTTCAAAAACACTGTCGCGGTTCGAAAGGTAAACGGAAGGGATTTCCAGCAACGTTAAGGCATCGCGCACCTGGGCGGCCTCCTGGCGGCTGCGCACCAGCACACTGATGTCCGAAGCACGCACCGGACGCGCGTCGTCGCCGtt  >  1:373378/1‑251 (MQ=255)                                                                                                                                                                                                              atTGACGTTGCCAACGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCANNGTTTCanna                                                                                                                                           >  2:15765/1‑99 (MQ=255)                                                                                                                                                                                                                           |                                                                                                                                                                                                                               CAAGATATCGGTAAGACGCCGCTCACCGCCTGCCGTTGCCAGCAAGTTTTCAGCAATGTTACGCGCCGACATCAGCGCCCGCAGCATCGGCATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCAGAGTTTCAAAAACACTGTCGCGGTTCGAAAGGTAAACGGAAGGGATTTCCAGCAACGTTAAGGCATCGCGCACCTGGGCGGCCTCCTGGCGGCTGCGCACCAGCACACTGATGTCCGAAGCACGCACCGGACGCGCGTCGTCGCCGTT  >  NC_000913/2953940‑2954380

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