Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,953,004	T→C	A1000A (GCA→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,953,004	0	T	C	100.0%	71.1 / NA	21	A1000A (GCA→GCG)	recB	exodeoxyribonuclease V subunit RecB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/11);  total (10/11)

AAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTTCCCACTGCGATTCAAAGCCGCCGAGTTCCAGTTTTTCCCGCACCCAGTTCGGGTCAACCGGCTGGGTAAAATCCAGGTCTTCAAA  >  NC_000913/2952875‑2953129                                                                                                                                  |                                                                                                                              aaaCTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGCGCCTggagga                                                                                                                     <  2:110509/140‑1 (MQ=255)  aaCTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGCGCCTGGAGGAc                                                                                                                    <  1:66878/140‑1 (MQ=255)        cGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGCGcc                                                                                                                            <  2:39620/126‑1 (MQ=255)        cGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGCGcc                                                                                                                            >  1:39620/1‑126 (MQ=255)             tAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGCGCCTGGAGGACAGCCGTGATcc                                                                                                         >  2:445901/1‑140 (MQ=255)                cGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGCGCCTGGAGGACAGCCGTGATCCACt                                                                                                      <  1:527633/140‑1 (MQ=255)                             cTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGCGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCgg                                                                                         <  2:334944/140‑1 (MQ=255)                             cTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGCGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCgg                                                                                         <  1:203785/140‑1 (MQ=255)                                       gCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGCGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTTCCCACTg                                                                                >  1:26927/1‑139 (MQ=255)                                       gCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGCGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTTCCCACTg                                                                                <  2:26927/139‑1 (MQ=255)                                          gTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGCGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCg                                                                                          <  1:323901/126‑1 (MQ=255)                                          gTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGCGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCg                                                                                          >  2:323901/1‑126 (MQ=255)                                                       cAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGGTTCATTGAGAGGCGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTTCCCACTGCGATTCAAAGCCGCCGa                                                               >  1:142340/1‑140 (MQ=255)                                                                         cACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGCGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTTCCCACTGCGATTCAAAGCCGCCGAGTTCCAGTTTTTCCCGCa                                             >  2:340089/1‑140 (MQ=255)                                                                                   ttGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGCGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTTCCCACTGCGATTCAAAGCCGCCGAGTTCCAGTTTTTCCCGCACCCAGTTCgg                                   >  1:495862/1‑140 (MQ=255)                                                                                         gcggAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGCGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTTCCCACTGCGATTCAAAGCCGCCGAGTTCCAGTTTTTCCCGCACCCAGTTCGGGTCAAc                             >  2:242185/1‑140 (MQ=255)                                                                                              aaGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGCGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTTCCCACTGCGATTCAAAGCCGCCGAGTTCCAGTTTTTCCCGCACCCAGTTCGGGTCAACCGGCt                        <  1:242185/140‑1 (MQ=255)                                                                                                           cTTACGCCGGTTTCATTGAGAGGCGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTTCCCACTGCGATTCAAAGCCGCCGAGTTCCAGTTTTTCCCGCACCCAGTTCGGGTCAACCGGCTGGGTAAAATCCAg           >  2:305260/1‑140 (MQ=255)                                                                                                                  cGGTTTCATTGAGAGGCGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTTCCCACTGCGATTCAAAGCCGCCGAGTTCCAGTTTTTCCCGCACCCAGTTCGGGTCAACCGGCTGGGTAAAATCCAGGTCTTca    <  1:445901/140‑1 (MQ=255)                                                                                                                    gTTTCATTGAGAGGCGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTTCCCACTGCGATTCAAAGCCGCCGAGTTCCAGTTTTTCCCGCACCCAGTTCGGGTCAACCGGCTGGGTAAAATCCAGGTCTTcaaa  >  1:432867/1‑140 (MQ=255)                                                                                                                               agGCGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTTCCCACTGCGATTCAAAGCCGCCGAGTTCCAGTTTTTCCCGCACCCAGTTCGGGTCAACCGGCTGGGTAAAATCCAGGTc        <  1:387648/123‑1 (MQ=255)                                                                                                                                  |                                                                                                                              AAACTGGCGGATTAACGTATCAAGCTGACTGGCGATAAGCGGTTCACTAATCGGCAGATAAAACTCCATCTCCACCTGTTTATTGCGGGCGGAAAGTTGACTCAGGCTTACGCCGGTTTCATTGAGAGGTGCCTGGAGGACAGCCGTGATCCACTCGGTCAATACCGGTTCCCACTGCGATTCAAAGCCGCCGAGTTCCAGTTTTTCCCGCACCCAGTTCGGGTCAACCGGCTGGGTAAAATCCAGGTCTTCAAA  >  NC_000913/2952875‑2953129

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