Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	494,419	G→A	P5P (CCG→CCA)	recR →	gap repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	494,419	0	G	A	100.0%	88.2 / NA	28	P5P (CCG→CCA)	recR	gap repair protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/12);  total (16/12)

AAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCC  >  NC_000913/494284‑494549                                                                                                                                        |                                                                                                                                   aaGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTg                                                                                                                                 >  1:1090000/1‑139 (MQ=255)     ccTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTGTTaa                                                                                                                             >  2:1008234/1‑139 (MQ=255)       tggtggCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTGTTAaca                                                                                                                           <  1:838223/139‑1 (MQ=255)              tGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTGTTAACACAGCTTa                                                                                                                    >  1:1003671/1‑139 (MQ=255)                     ttCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTGTTAACACAGCTTATGGAAGc                                                                                                             <  2:1090000/139‑1 (MQ=255)                                       cgtATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCggg                                                                                           <  2:428833/139‑1 (MQ=255)                                                gaaACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTGt                                                                                                                                <  2:951784/93‑1 (MQ=255)                                                gaaACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTGt                                                                                                                                >  1:951784/1‑93 (MQ=255)                                                   aCGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGaa                                                                               >  2:895026/1‑139 (MQ=255)                                                                          tGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGt                                                                                        <  1:1011985/107‑1 (MQ=255)                                                                          tGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGt                                                                                        >  2:1011985/1‑107 (MQ=255)                                                                                   cGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTGTTAAcaca                                                                                                                         >  1:398080/1‑65 (MQ=255)                                                                                   cGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTGTTAAcaca                                                                                                                         <  2:398080/65‑1 (MQ=255)                                                                                   cGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCg                                                                                         >  1:730718/1‑97 (MQ=255)                                                                                   cGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCg                                                                                         <  2:730718/97‑1 (MQ=255)                                                                                    ggAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGt                                                                             >  1:1039532/1‑108 (MQ=255)                                                                                    ggAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGt                                                                             <  2:1039532/108‑1 (MQ=255)                                                                                      aaTGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGtt                                                          <  1:105657/125‑1 (MQ=255)                                                                                      aaTGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGtt                                                          >  2:105657/1‑125 (MQ=255)                                                                                             cTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCACTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCg                                     <  1:1008234/139‑1 (MQ=255)                                                                                                                 gCCGTTCTGATGCAAACCAGCCCACTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTgg                 >  1:918331/1‑139 (MQ=255)                                                                                                                     ttCTGATGCAAACCAGCCCACTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCa             >  2:10286/1‑139 (MQ=255)                                                                                                                               aaCCGGCCCACTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGcgcg                                         <  1:963547/101‑1 (MQ=255)                                                                                                                               aaCCGGCCCACTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGcgcg                                         >  2:963547/1‑101 (MQ=255)                                                                                                                               aaCCAGCCCACTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCAcc   >  1:406381/1‑139 (MQ=255)                                                                                                                                aCCAGCCCACTGTTAACACAGCTTATGGAAGCACTGCGCAGTCTGCCGGGCGTTCGCCCGAAGTCGGCGCAGCGTATGGCGTTCTCGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCAccc  >  2:747903/1‑139 (MQ=255)                                                                                                                                    gCCCACTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTc                    >  2:461435/1‑117 (MQ=255)                                                                                                                                    gCCCACTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTc                    <  1:461435/117‑1 (MQ=255)                                                                                                                                        |                                                                                                                                   AAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTTCTGATGCAAACCAGCCCGCTGTTAACACAGCTTATGGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCC  >  NC_000913/494284‑494549

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