Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,271,122	A→G	V917A (GTC→GCC)	uvrA ←	excision nuclease subunit A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,271,122	0	A	G	100.0%	123.4 / NA	34	V917A (GTC→GCC)	uvrA	excision nuclease subunit A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/16);  total (18/16)

TAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTC  >  NC_000913/4270993‑4271255                                                                                                                                  |                                                                                                                                      tAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGAtctc                                                                                                                             <  1:250336/140‑1 (MQ=255)    aTTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTcgcc                                                                                                                          <  2:312045/140‑1 (MQ=255)         tGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGccac                                                                                                                     <  2:526917/140‑1 (MQ=255)                      aaaGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTc                                                                                                        <  1:326767/140‑1 (MQ=255)                       aaGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCt                                                                                                       <  1:244633/140‑1 (MQ=255)                           aaaaGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTc                                                                                                   >  1:384481/1‑140 (MQ=255)                                    cTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGa                                                                                          >  1:344473/1‑140 (MQ=255)                                              ggCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTc                                                                                >  2:502006/1‑140 (MQ=255)                                               gCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAAt                                                                                      <  2:442037/133‑1 (MQ=255)                                               gCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAAt                                                                                      >  1:442037/1‑133 (MQ=255)                                                       gATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGtttt                                                                       >  2:408072/1‑140 (MQ=255)                                                               cATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGt                                                               >  1:299668/1‑140 (MQ=255)                                                                     cTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGAt                                                         <  2:440351/140‑1 (MQ=255)                                                                              gCGTGCCGTGTGTTATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCaa                                                >  1:296619/1‑140 (MQ=255)                                                                              gCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCaa                                                >  2:533350/1‑140 (MQ=255)                                                                               cGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAAt                                               <  1:498978/140‑1 (MQ=255)                                                                               cGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAAt                                               <  1:408072/140‑1 (MQ=255)                                                                                  gCCGTGTGTTATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATcac                                            <  2:296619/140‑1 (MQ=255)                                                                                     gtgtgtGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATcaccac                                         >  1:143532/1‑140 (MQ=255)                                                                                            aTGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGgtg                                  <  2:398067/140‑1 (MQ=255)                                                                                                  cGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACaa                                       <  2:221655/129‑1 (MQ=255)                                                                                                  cGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACaa                                       >  1:221655/1‑129 (MQ=255)                                                                                                  cGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGccc                            >  1:92032/1‑140 (MQ=255)                                                                                                  cGCACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGccc                            >  2:268322/1‑140 (MQ=255)                                                                                                    cACTCCGCGACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTg                          >  1:23516/1‑140 (MQ=255)                                                                                                         cgcgACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATcgc                     <  2:92032/140‑1 (MQ=255)                                                                                                         cgcgACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATcgc                     <  1:509199/140‑1 (MQ=255)                                                                                                         cgcgACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATcgc                     <  1:502006/140‑1 (MQ=255)                                                                                                         cgcgACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATcgc                     <  1:449709/140‑1 (MQ=255)                                                                                                         cgcgACGGTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATcgc                     >  1:354972/1‑140 (MQ=255)                                                                                                               ggTTTCTGGCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGtt               >  2:288902/1‑140 (MQ=255)                                                                                                                       gCGTACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGt       >  1:492691/1‑140 (MQ=255)                                                                                                                          tACCGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACg    >  1:506991/1‑140 (MQ=255)                                                                                                                            ccGGAGGCGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTc  >  2:27229/1‑140 (MQ=255)                                                                                                                                  |                                                                                                                                      TAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGATCTCGCCGCCACCACTGCCGCCTTCTGGTCCCAGGTCGACAATCCAGTCAGCGGTTTTGATCACGTCGAGATTGTGCTCAATCACCACAATGGTGTTGCCCTGATCGCGCAGTTTATGCAGTACGTC  >  NC_000913/4270993‑4271255

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