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%	131.2 / NA	42	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 (19/23);  total (19/23)

AAACTACCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCT  >  NC_000913/1532910‑1533180                                                                                                                                        |                                                                                                                                        aaaCTACCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTa                                                                                                                                      >  1:925949/1‑139 (MQ=255)          ggCGTCAGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCg                                                                                                                             <  2:206753/139‑1 (MQ=255)                aGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATc                                                                                                                       >  1:294039/1‑139 (MQ=255)                 gTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCg                                                                                                                      <  1:285773/139‑1 (MQ=255)                                   cTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCcgc                                                                                                    >  1:683670/1‑139 (MQ=255)                                       cGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGa                                                                                                >  2:479480/1‑139 (MQ=255)                                                gCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGa                                                                                       >  2:315619/1‑139 (MQ=255)                                                      ggCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACg                                                                                        <  1:255894/132‑1 (MQ=255)                                                      ggCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACg                                                                                        >  2:255894/1‑132 (MQ=255)                                                      ggCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAacaaca                                                                                 <  1:479480/139‑1 (MQ=255)                                                             ttGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAagca                                                                          <  2:252060/139‑1 (MQ=255)                                                                      aaaCGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTtgct                                                                 <  2:334447/139‑1 (MQ=255)                                                                       aaCGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTtgctg                                                                >  1:150477/1‑139 (MQ=255)                                                                       aaCGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTtgctg                                                                >  2:184763/1‑139 (MQ=255)                                                                       aaCGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTtgctg                                                                <  2:150477/139‑1 (MQ=255)                                                                          gCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAg                                                             >  1:68491/1‑139 (MQ=255)                                                                              tACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGatat                                                         >  1:623417/1‑139 (MQ=255)                                                                                    acacATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGa                                                   <  2:641246/139‑1 (MQ=255)                                                                                        aTGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGACCTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGtt                                               >  2:92105/1‑139 (MQ=255)                                                                                        aTGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGtt                                               >  1:810353/1‑139 (MQ=255)                                                                                        aTGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGtt                                               >  2:449285/1‑139 (MQ=255)                                                                                          gACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGttat                                             <  1:184763/139‑1 (MQ=255)                                                                                          gACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGttat                                             >  1:616447/1‑139 (MQ=255)                                                                                               tcgcaATGTGTGTGAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTc                                        <  2:616447/136‑1 (MQ=255)                                                                                                cacaATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCg                                                                                                                      <  1:308613/60‑1 (MQ=255)                                                                                                cacaATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCg                                                                                                                      >  2:308613/1‑60 (MQ=255)                                                                                                    aTGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCAt                                   <  2:925949/139‑1 (MQ=255)                                                                                                    aTGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCAt                                   <  2:84575/139‑1 (MQ=255)                                                                                                       tgtgtTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTg                                <  1:92105/139‑1 (MQ=255)                                                                                                        gtgtTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGa                               <  2:623417/139‑1 (MQ=255)                                                                                                        gtgtTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGa                               <  1:449285/139‑1 (MQ=255)                                                                                                         tgtTACTCTTACCGACATCCTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAACGGTTTCCCCGTGAACTGGACGAACAACGAAAAGCAGCACTTGCTGCAGACATTACCGACGTTATTATTCGTGATCTGaa                              <  2:610842/139‑1 (MQ=255)                                                                                                         tgtTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACt                                                                                             >  2:663488/1‑76 (MQ=255)                                                                                                         tgtTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACt                                                                                             <  1:663488/76‑1 (MQ=255)                                                                                                         tgtTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGt                                      <  2:404239/131‑1 (MQ=255)                                                                                                         tgtTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGt                                      >  1:404239/1‑131 (MQ=255)                                                                                                         tgtTAATTTTACTGACATCTATCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGaa                              >  1:631033/1‑139 (MQ=255)                                                                                                                   aCTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGAc                    >  1:482768/1‑139 (MQ=255)                                                                                                                      gACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGt                 <  2:482768/139‑1 (MQ=255)                                                                                                                                aaaaGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCa       <  2:810353/139‑1 (MQ=255)                                                                                                                                   aGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTg    <  2:721915/139‑1 (MQ=255)                                                                                                                                     gAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGct  <  2:294039/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        AAACTACCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCATTGCT  >  NC_000913/1532910‑1533180

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: ‑