Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,533,045	C→T	intergenic (+93/‑7)	ydcC → / → pptA	H repeat‑associated putative transposase/4‑oxalocrotonate tautomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,533,045	0	C	T	100.0%	94.2 / NA	31	intergenic (+93/‑7)	ydcC/pptA	H repeat‑associated putative transposase/4‑oxalocrotonate tautomerase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/17);  total (14/17)

GCGTCAGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCAT  >  NC_000913/1532920‑1533176                                                                                                                              |                                                                                                                                    gCGTCAGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGc                                                                                                                        >  1:868007/1‑139 (MQ=255)      aGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATc                                                                                                                   <  2:245746/139‑1 (MQ=255)           ttACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACAt                                                                                                              <  1:872870/139‑1 (MQ=255)                         cTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCcgc                                                                                                <  2:741395/139‑1 (MQ=255)                                      gCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGa                                                                                   <  2:845553/139‑1 (MQ=255)                                          tGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGaacaa                                                                               >  1:777824/1‑139 (MQ=255)                                                        gTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACt                                                                 <  2:868007/139‑1 (MQ=255)                                                                    tACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGatat                                                     >  1:429639/1‑139 (MQ=255)                                                                    tACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGatat                                                     >  1:866245/1‑139 (MQ=255)                                                                        ggACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTAcc                                                 <  2:702010/139‑1 (MQ=255)                                                                           cacaTGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGAc                                              >  2:896926/1‑139 (MQ=255)                                                                                   aTTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTtattat                                      >  1:580928/1‑139 (MQ=255)                                                                                      cacaATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCg                                   >  2:787676/1‑139 (MQ=255)                                                                                        caATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTAc                                                  >  1:442607/1‑122 (MQ=255)                                                                                        caATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTAc                                                  <  2:442607/122‑1 (MQ=255)                                                                                          aTGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCAt                               <  1:896926/139‑1 (MQ=255)                                                                                          aTGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCAt                               <  2:429639/139‑1 (MQ=255)                                                                                              gtgtTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGa                           <  2:777824/139‑1 (MQ=255)                                                                                              gtgtTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGa                           <  2:866245/139‑1 (MQ=255)                                                                                               tgtTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATc                                                                                                                   >  1:277956/1‑50 (MQ=255)                                                                                               tgtTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATc                                                                                                                   <  2:277956/50‑1 (MQ=255)                                                                                               tgtTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGaa                          >  1:978064/1‑139 (MQ=255)                                                                                               tgtTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGaa                          >  1:312457/1‑139 (MQ=255)                                                                                                      tttACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTaaa                   <  1:896269/139‑1 (MQ=255)                                                                                                      tttACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGATTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGTACAGTaaa                   >  1:1051761/1‑139 (MQ=255)                                                                                                               aTCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGa                           >  2:308587/1‑122 (MQ=255)                                                                                                               aTCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGa                           <  1:308587/122‑1 (MQ=255)                                                                                                                  tttCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTTATCTGa                           >  2:42124/1‑119 (MQ=255)                                                                                                                  tttCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTTATCTGa                           <  1:42124/119‑1 (MQ=255)                                                                                                                      aaaaGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGTACAGTAAAGACAGTTCGATAAGCa   <  2:1051761/139‑1 (MQ=255)                                                                                                                       aaaGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCAt  <  2:924217/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                    GCGTCAGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCAT  >  NC_000913/1532920‑1533176

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