Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,752,917	C→G	L375V (CTT→GTT)	recN →	DNA repair protein RecN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,752,917	0	C	G	100.0%	64.5 / NA	19	L375V (CTT→GTT)	recN	DNA repair protein RecN
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (3/16);  total (3/16)

ACAGCAGCAACTGGACGATCAGGCCGACTCACAAGAAACGCTTGCGCTGGCGGTAACGAAACATCATCAGCAGGCACTGGAAATCGCGCGCGCATTACACCAACAACGCCAGCAATATGCAGAAGAACTTGCACAGCTGATCACCGACAGTATGCATGCGCTCTCAATGCCGCATGGGCAGTTTACGATCGATGTTAAATTTGACGAGCATCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGT  >  NC_000913/2752790‑2753041                                                                                                                                |                                                                                                                             aCAGCAGCAACTGGACGATCAGGCCGACTCACAAGAAACGCTTGCGCTGGCGGTAACGAAACATCATCAGCAGGCACTGGAAATCGCGCGCGCATTACACCAACAACGCCAGCAATATGCAGAAGAAGTTGca                                                                                                                         <  2:617139/133‑1 (MQ=255)   agcagcAACTGGACGATCAGGCCGACTCACAAGAAACGCTTGCGCTGGCGGTAACGAAACATCATCAGCAGGCACTGGAAATCGCGCGCGCATTACACCAACAACGCCAGCAATATGCAGAAGAAGTTGcaca                                                                                                                       <  1:35793/133‑1 (MQ=255)     cagcaACTGGACGATCAGGCCGACTCACAAGAAACGCTTGCGCTGGCGGTAACGAAACATCATCAGCAGGCACTGGAAATCGCGCGCGCATTACACCAACAACGCCAGCAATATGCAGAAGAAGTTGCACAGc                                                                                                                     <  1:376171/133‑1 (MQ=255)         aaCTGGACGATCAGGCCGACTCACAAGAAACGCTTGCGCTGGCGGTAACGAAACATCATCAGCAGGCACTGGAAATCGCGCGCGCATTACACCAACAACGCCAGCAATATGCAGAAGAAGTTGCACAGCTGAt                                                                                                                 <  2:473016/133‑1 (MQ=255)             ggACGATCAGGCCGACTCACAAGAAACGCTTGCGCTGGCGGTAACGAAACATCATCAGCAGGCACTGGAAATCGCGCGCGCATTACACCAACAACGCCAGCAATATGCAGAAGAAGTTGCACAGCTGATCAcc                                                                                                             <  1:79416/133‑1 (MQ=255)               aCGATCAGGCCGACTCACAAGAAACGCTTGCGCTGGCGGTAACGAAACATCATCAGCAGGCACTGGAAATCGCGCGCGCATTACACCAACAACGCCAGCAATGTGCAGAAGAAGGTGGACAGCTGATCACCGa                                                                                                           >  1:160474/1‑133 (MQ=255)               aCGATCAGGCCGACTCACAAGAAACGCTTGCGCTGGCGGTAACGAAACATCATCAGCAGGCACTGGAAATCGCGCGCGCATTACACCAACAACGCCAGCAATATGCAGAAGAAGTTGCACAGCTGATCACCGa                                                                                                           <  1:501259/133‑1 (MQ=255)                  aTCAGGCCGACTCACAAGAAACGCTTGCGCTGGCGGTAACGAAACATCATCAGCAGGCACTGGAAATCGCGCGCGCATTACACCAACAACGCCAGCAATATGCAGAAGAAGTTGCACAGCTGATCACCGACAg                                                                                                        >  2:422257/1‑133 (MQ=255)                       gCCGACTCACAAGAAACGCTTGCGCTGGCGGTAACGAAACATCATCAGCAGGCACTGGAAATCGCGCGCGCATTACACCAACAACGCCAGCAATATGCAGAAGAAGTTGCACAGCTGATCACCGACAGTatgc                                                                                                   <  2:333780/133‑1 (MQ=255)                           aCTCACAATAAACGCCTGCGCTGGCGGTAACGAAACATCATCAGCAGGCACTGGAAATCGCGCGCGCATTACACCAACAACGCCAGCAATATGCAGAAGAAGTTGCACAGCTGATCACCGACAGTatgcatgc                                                                                               <  1:165071/133‑1 (MQ=255)                                           tGCGCTGGCGGTAACGAAACATCATCAGCAGGCACTGGAAATCGCGCGCGCATTACACCAACAACGCCAGCAATATGCAGAAGAAGTTGCACAGCTGATCACCGACAGTATGCATGCGCTCTCAATGCCGCAt                                                                               >  1:501666/1‑133 (MQ=255)                                              gcTGGCGGTATCGAATCATCATCAGCAGGCACTGGAAATCGAGCGCGCATTACACCAACAACGCCAGCAATATGCAGAAGAAGTTGCACAGCTGATCACCGACAGTATGCATGCGCTCTCAATGCCGCATggg                                                                            <  2:243289/133‑1 (MQ=255)                                              gcTGGCGGTAACGAAACATCATCAGCAGGCACTGGAAATCGCGCGCGCATTACACCAACAACGCCAGCAATATGCAGAAGAAGTTGCACAGCTGATCACCGACAGTATGCATGCGCTCTCAATGCCGCATggg                                                                            <  1:422257/133‑1 (MQ=255)                                              gcTGGCGGTAACGAAACATCATCAGCAGGCACTGGAAATCGCGCGCGCATTACACCAACAACGCCAGCAATATGCAGAAGAAGTTGCACAGCTGATCACCGACAGTATGCATGCGCTCTCAATGCCGCATggg                                                                            <  2:530930/133‑1 (MQ=255)                                              gcTGGCGGTAACGAAACATCATCAGCAGGCACTGGAAATCGCGCGCGCATTACACCAACAACGCCAGCAATATGCAGAAGAAGTTGCACAGCTGATCACCGACAGTATGCATGCGCTCTCAATGCCGCATggg                                                                            <  2:601188/133‑1 (MQ=255)                                                          gAAACATCATCAGCAGGCACTGGAAATCGCGCGCGCATTACACCAACAACGCCAGCAATATGCAGAAGAAGTTGCACAGCTGATCACCGACAGTATGCATGCGCTCTCAATGCCGCATGGGCAGTTTAcgatc                                                                <  1:474397/133‑1 (MQ=255)                                                                                                tACACCAACAACGCCAGCAATATGCAGAAGAAGTTGCACAGCTGATCACCGACAGTATGCATGCGCTCTCAATGCCGCATGGGCAGTTTACGATCGATGTTAAATTTGACGAGCATCACCTGGGCGCTGACgg                          <  1:448606/133‑1 (MQ=255)                                                                                                   acCAACAACGCCAGCCATATGCAGAAGAAGTTACACAGCTGTTCACCGACAGTATGCATGCGCTCTCAATGCCGCATGGGCAGTTTACGATCGATGTTAAATTTGACGAGCATCACCTGGGCGCTGACGGTGc                       <  2:233499/133‑1 (MQ=255)                                                                                                                        cAGAAGAAGTTGCACAGCTGATCACCGACAGTATGCATGCGCTCTCAATGCCGCATGGGCAGTTTACGATCGATGTTAAATTTGACGAGCATCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGt  <  2:501666/133‑1 (MQ=255)                                                                                                                                |                                                                                                                             ACAGCAGCAACTGGACGATCAGGCCGACTCACAAGAAACGCTTGCGCTGGCGGTAACGAAACATCATCAGCAGGCACTGGAAATCGCGCGCGCATTACACCAACAACGCCAGCAATATGCAGAAGAACTTGCACAGCTGATCACCGACAGTATGCATGCGCTCTCAATGCCGCATGGGCAGTTTACGATCGATGTTAAATTTGACGAGCATCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGT  >  NC_000913/2752790‑2753041

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