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%	123.4 / NA	39	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 (20/19);  total (20/19)

AAACTACCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCAT  >  NC_000913/1532910‑1533176                                                                                                                                        |                                                                                                                                    aaaCTACCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTa                                                                                                                                  >  2:1161268/1‑139 (MQ=255)           gCGTCAGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGc                                                                                                                        >  2:617427/1‑139 (MQ=255)            cGTCAGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGc                                                                                                                           <  2:407612/135‑1 (MQ=255)            cGTCAGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGc                                                                                                                           >  1:407612/1‑135 (MQ=255)            cGTCAGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGca                                                                                                                       <  1:617427/139‑1 (MQ=255)                 gTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCg                                                                                                                  <  2:559037/139‑1 (MQ=255)                                                      ggCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAacaaca                                                                             >  2:789683/1‑139 (MQ=255)                                                                tCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGttt                                                                                                     <  2:943073/105‑1 (MQ=255)                                                                tCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGttt                                                                                                     >  1:943073/1‑105 (MQ=255)                                                                  gTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACt                                                                 <  2:190754/139‑1 (MQ=255)                                                                       aaCGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTtgctg                                                            >  2:1178394/1‑139 (MQ=255)                                                                        aCGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTtgctgc                                                           >  1:737675/1‑139 (MQ=255)                                                                             ttACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAacaaca                                                                             <  2:1198531/116‑1 (MQ=255)                                                                             ttACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAacaaca                                                                             >  1:1198531/1‑116 (MQ=255)                                                                             ttACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGata                                                      >  2:757094/1‑139 (MQ=255)                                                                               aCCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCa                                                                                                                              <  1:616905/65‑1 (MQ=255)                                                                               aCCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCa                                                                                                                              >  2:616905/1‑65 (MQ=255)                                                                                    acacATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGa                                               <  2:737675/139‑1 (MQ=255)                                                                                      acaTGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACg                                             <  1:1161268/139‑1 (MQ=255)                                                                                          gACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGttat                                         <  1:789683/139‑1 (MQ=255)                                                                                               tCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTc                                    <  1:757094/139‑1 (MQ=255)                                                                                                  caATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTc                                 >  1:434458/1‑139 (MQ=255)                                                                                                   aaTGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCa                                >  2:605865/1‑139 (MQ=255)                                                                                                    aTGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCAt                               >  2:315405/1‑139 (MQ=255)                                                                                                        gtgtTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGa                           <  1:605865/139‑1 (MQ=255)                                                                                                        gtgtTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGa                           <  1:315405/139‑1 (MQ=255)                                                                                                             aaTTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAagcagc                                                                    <  1:1260028/93‑1 (MQ=255)                                                                                                             aaTTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAagcagc                                                                    >  2:1260028/1‑93 (MQ=255)                                                                                                                 ttACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGat                          <  1:878859/131‑2 (MQ=255)                                                                                                                 ttACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGat                          >  2:878859/1‑130 (MQ=255)                                                                                                                 ttACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAg                  >  2:5902/1‑139 (MQ=255)                                                                                                                   aCTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGa                                                                                            <  2:502268/63‑1 (MQ=255)                                                                                                                   aCTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGa                                                                                            >  1:502268/1‑63 (MQ=255)                                                                                                                   aCTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACt                                                                 <  2:627420/90‑1 (MQ=255)                                                                                                                   aCTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACt                                                                 >  1:627420/1‑90 (MQ=255)                                                                                                                       aCATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAATAAGCATCACTTGCTGCAGATATTACCGACGTTATTATTCGTCAt                               <  1:1042972/120‑1 (MQ=255)                                                                                                                       aCATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCAt                               >  2:1042972/1‑120 (MQ=255)                                                                                                                            tttCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATa       >  2:886252/1‑139 (MQ=255)                                                                                                                                 aaaGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCAt  <  1:1178394/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                    AAACTACCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCAT  >  NC_000913/1532910‑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: ‑