Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,955,564	A→T	L147Q (CTG→CAG)	recB ←	exodeoxyribonuclease V subunit RecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,955,564	0	A	T	100.0%	125.9 / NA	36	L147Q (CTG→CAG)	recB	exodeoxyribonuclease V subunit RecB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (25/11);  total (25/11)

TAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGTAAAGACTGCCGCTTC  >  NC_000913/2955486‑2955637                                                                               |                                                                          tAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACTGc                                                                         <  2:511355/81‑1 (MQ=255) tAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACTGc                                                                         <  2:511354/81‑1 (MQ=255)      gTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACTGCATGcc                                                                    <  2:511356/81‑1 (MQ=255)         gCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGa                                                                 >  2:511398/1‑81 (MQ=255)              gCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCa                                                             >  1:511379/1‑80 (MQ=255)                    aaTCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCAt                                                      >  1:511382/1‑81 (MQ=255)                          cgcAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGt                                                >  2:511417/1‑81 (MQ=255)                           gcAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGtt                                               <  1:511357/81‑1 (MQ=255)                           gcAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGtt                                               >  2:511397/1‑81 (MQ=255)                                  tGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATg                                        >  1:511407/1‑81 (MQ=255)                                      agcgtagcAGAGACTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCt                                    >  2:511416/1‑81 (MQ=255)                                      agcgtagcAGAGACTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCt                                    >  1:511383/1‑81 (MQ=255)                                      agcgtagcAGAGACTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCt                                    >  2:511389/1‑81 (MQ=255)                                        cgtagcAGAGACTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTgg                                  >  2:511412/1‑81 (MQ=255)                                           agcAGAGACTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCa                                >  1:511400/1‑80 (MQ=255)                                              agagACTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAg                            >  2:511424/1‑81 (MQ=255)                                                agaCTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGcc                          <  1:511358/81‑1 (MQ=255)                                                agaCTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGcc                          <  2:511359/81‑1 (MQ=255)                                                agaCTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGcc                          >  2:511387/1‑81 (MQ=255)                                                agaCTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGcc                          <  1:511360/81‑1 (MQ=255)                                                 gaCTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCg                         >  1:511408/1‑81 (MQ=255)                                                  aCTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGt                        >  1:511414/1‑81 (MQ=255)                                                   cTCATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTg                       <  1:511361/81‑1 (MQ=255)                                                     cATCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGaa                     >  2:511393/1‑81 (MQ=255)                                                       tCTTCAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGaa                     >  1:511390/1‑79 (MQ=255)                                                           cAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGt                 <  1:511362/79‑1 (MQ=255)                                                           cAATCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGTa                >  2:511419/1‑80 (MQ=255)                                                            aaTCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGTaaa              <  1:511363/81‑1 (MQ=255)                                                            aaTCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGTa                >  2:511411/1‑79 (MQ=255)                                                            aaTCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGTa                >  1:511388/1‑79 (MQ=255)                                                             aTCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGTa                >  2:511409/1‑78 (MQ=255)                                                             aTCAGCTGCTGCTCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGTa                >  1:511401/1‑78 (MQ=255)                                                                      tgctCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGGAAAAGCCGTGAATAGTAAAGACTGCCGCt    >  2:511413/1‑81 (MQ=255)                                                                      tgctCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGTAAAGACTGCCGCt    >  1:511420/1‑81 (MQ=255)                                                                       gctCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGTAAAGACTGCCGCtt   >  2:511410/1‑81 (MQ=255)                                                                        ctCAAACTGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGTAAAGACTGCCGCTTc  <  1:511364/81‑1 (MQ=255)                                                                               |                                                                          TAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGTAAAGACTGCCGCTTC  >  NC_000913/2955486‑2955637

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