Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,954,152	T→C	T618A (ACG→GCG)	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,152	0	T	C	100.0%	47.0 / NA	15	T618A (ACG→GCG)	recB	exodeoxyribonuclease V subunit RecB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (6/9);  total (6/9)

CATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCAGAGTTTCAAAAACACTGTCGCGGTTCGAAAGGTAAACGG  >  NC_000913/2954030‑2954271                                                                                                                           |                                                                                                                        cATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGCTGCCAGCGCa                                                                                                               >  1:404175/1‑133 (MQ=255)  aTAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGCTGCCAGCGCAc                                                                                                              <  1:109075/133‑1 (MQ=255)     aCGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGCTGCCAGCGCACTAc                                                                                                           <  1:356685/133‑1 (MQ=255)      cGCCACGTTTGCGCCTGATTTGACGATAACTATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGCTGCCAGGGCACTACg                                                                                                          <  1:196179/133‑1 (MQ=255)       gCCACGTTTGAGCCATATTTGCCGCTTACCATCGAACTCTTAGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGCTGCCAGCGCACTACGc                                                                                                         <  2:230591/133‑1 (MQ=255)          aCGTTTGCCCCAGATATGCCGACAACCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGCTGCCAGCGCACTACGCAgg                                                                                                      <  1:141438/133‑1 (MQ=255)                       aTTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGCTGCCAGCGCACTACGCAGGGTGTTCTCACGtt                                                                                         <  1:385012/133‑1 (MQ=255)                         ttGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGCTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCg                                                                                       <  2:25932/133‑1 (MQ=255)                                 aaCCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGCTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGtcat                                                                               >  1:122187/1‑133 (MQ=255)                                                                cGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGCTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTgc                                                >  2:115744/1‑133 (MQ=255)                                                                   aTGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGCTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCt                                             <  2:404175/133‑1 (MQ=255)                                                                            aTTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGCTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCAGAGtt                                    <  2:97744/133‑1 (MQ=255)                                                                                        ttCGATATCCAGCGCGTTCAGCCCCATCATTGACGCTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCAGAGTTTCAAAAACACTg                        >  2:513792/1‑133 (MQ=255)                                                                                          cGATATCCAGCGCGTTCAGCCCCATCATTGACGCTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAAGCAAAGCATATCCAGCGTTTCCAGAGGTTCAAAAACACTGTc                      >  2:156905/1‑133 (MQ=255)                                                                                                   gcgcGTTCACCCCCATCATTGCCGCGGCCAGCGCACTCCGCGGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCAGAGTTTCAAAAACACTGTCGCGGTTCGa             <  2:122187/133‑1 (MQ=255)                                                                                                    cgcgTTCAGCCCCATCATTGACGCTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCAGAGTTTCAAAAACACTGTCGCGGTTCGaa            >  2:405472/1‑133 (MQ=255)                                                                                                              cccATCATAGACGCTCCCCGCGCACTCCGCTGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCAGAGTTTCAAAAACACTGTCGCGGTTCGAAAGGTAAACgg  <  1:143075/133‑1 (MQ=255)                                                                                                                           |                                                                                                                        CATAACGCCACGTTTGCGCCAGATTTGCCGATAACCATCGAACTCTTCGACTACCACATCCCACGCATGTTCGTCATTGTTCAGCGTTTCGATATCCAGCGCGTTCAGCCCCATCATTGACGTTGCCAGCGCACTACGCAGGGTGTTCTCACGTTCGGGCGTCATCACCGCCTGCAACAACCAAAGCATTTCCTGCGCTTCCAGAGTTTCAAAAACACTGTCGCGGTTCGAAAGGTAAACGG  >  NC_000913/2954030‑2954271

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