Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,036,373	C→T	*578* (TAG→TAA)	recJ ←	ssDNA exonuclease, 5' ‑‑> 3'‑specific

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,036,373	0	C	T	100.0%	144.3 / NA	46	*578* (TAG→TAA)	recJ	ssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/21);  total (25/21)

TGAGGTCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGCTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTA  >  NC_000913/3036240‑3036502                                                                                                                                      |                                                                                                                                  tGAGGTCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAAt                                                                                                                              <  2:987182/139‑1 (MQ=255)      tCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGcc                                                                                                                         <  2:211361/139‑1 (MQ=255)        cTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGGGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAAt                                                                                                                              <  1:1192787/132‑1 (MQ=255)        cTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAAt                                                                                                                              >  2:1192787/1‑132 (MQ=255)                gATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTc                                                                                                               >  2:684136/1‑139 (MQ=255)                            gATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTg                                                                                                   >  1:57962/1‑139 (MQ=255)                                            ggTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCcg                                                                                   >  1:945742/1‑139 (MQ=255)                                                            ttGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGatat                                                                   >  1:770613/1‑139 (MQ=255)                                                              gAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCg                                                                 <  2:1490346/139‑1 (MQ=255)                                                                 cTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGa                                                                <  1:1318605/137‑1 (MQ=255)                                                                 cTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGa                                                                >  2:1318605/1‑137 (MQ=255)                                                                               gTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTaaa                                                                                                                               <  2:635594/60‑1 (MQ=255)                                                                               gTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTaaa                                                                                                                               >  1:635594/1‑60 (MQ=255)                                                                                          gggATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCa                                     <  2:57962/139‑1 (MQ=255)                                                                                                 tACCCAATCCACGCTCTTTTTTATAGAGCAGATGACGTTAAATTTGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCgcgc                                      >  2:977125/1‑131 (MQ=255)                                                                                                 tACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGt                                                  <  2:1265041/119‑1 (MQ=255)                                                                                                 tACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGt                                                  >  1:1265041/1‑119 (MQ=255)                                                                                                 tACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCa                                     <  1:962752/132‑1 (MQ=255)                                                                                                 tACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCa                                     >  2:962752/1‑132 (MQ=255)                                                                                                 tACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGAGAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCgcgc                                      <  1:977125/131‑1 (MQ=255)                                                                                                    ccAATCCACGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATc                           <  1:684136/139‑1 (MQ=255)                                                                                                           aCGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCa                                     >  2:980709/1‑122 (MQ=255)                                                                                                           aCGCTCTTTTTTATAGAGAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCa                                     <  1:980709/122‑1 (MQ=255)                                                                                                                         agagAAGATGACGTTAAATTGGCCAGATATTGTCGATGATAATTGGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGAc      >  1:1181847/1‑139 (MQ=255)                                                                                                                           agaagaTGACGTTAAATTGGCCAGATATTGTCGATGGTAATTTGCAGGCTGCGGTTGCCGCGAGACTCGTTGATATc                                                                  >  2:1061542/1‑77 (MQ=255)                                                                                                                           agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATc                                                                  <  1:1061542/77‑1 (MQ=255)                                                                                                                           agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCa                                     >  2:671185/1‑106 (MQ=255)                                                                                                                           agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCa                                     <  1:671185/106‑1 (MQ=255)                                                                                                                           agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGcc                                 >  1:1620568/1‑110 (MQ=255)                                                                                                                           agaagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGcc                                 <  2:1620568/110‑1 (MQ=255)                                                                                                                             aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGt                                                                                        <  1:1007346/53‑1 (MQ=255)                                                                                                                             aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGt                                                                                        >  2:1007346/1‑53 (MQ=255)                                                                                                                             aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCcgc                                                                                  <  2:1259431/59‑1 (MQ=255)                                                                                                                             aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCcgc                                                                                  >  1:1259431/1‑59 (MQ=255)                                                                                                                             aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTg                                                                       <  2:1023348/70‑1 (MQ=255)                                                                                                                             aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTg                                                                       >  1:1023348/1‑70 (MQ=255)                                                                                                                             aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGa                                                                >  2:1315915/1‑77 (MQ=255)                                                                                                                             aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGa                                                                <  1:1315915/77‑1 (MQ=255)                                                                                                                             aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGc                                                      <  2:1226959/87‑1 (MQ=255)                                                                                                                             aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGc                                                      >  1:1226959/1‑87 (MQ=255)                                                                                                                             aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATc                           >  2:458800/1‑114 (MQ=255)                                                                                                                             aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATc                           <  1:458800/114‑1 (MQ=255)                                                                                                                             aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTa  >  2:254957/1‑139 (MQ=255)                                                                                                                             aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTa  >  2:1235982/1‑139 (MQ=255)                                                                                                                             aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTa  >  2:755883/1‑139 (MQ=255)                                                                                                                             aagaTGACGTTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTa  >  1:1616784/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                  TGAGGTCCTGAATGCGATTATTTACCGGATTAATTTCAAACATGGTCTGATTTCTTTTATTGAGCTAGTCAAAATGCGGTGATAAGAGCGGGATTGTACCCAATCCACGCTCTTTTTTATAGAGAAGATGACGCTAAATTGGCCAGATATTGTCGATGATAATTTGCAGGCTGCGGTTGCCGCGAAACTCGTTGATATCGAGCTTATAAGCCAGTTGCACTTCGCGCACGCCGTTATCCGGCCAGAGGGCGGTATCGACATTA  >  NC_000913/3036240‑3036502

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