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%	122.3 / NA	40	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/20);  total (20/20)

AAACTACCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCA  >  NC_000913/1532910‑1533175                                                                                                                                        |                                                                                                                                   aaaCTACCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTa                                                                                                                                 >  2:544835/1‑139 (MQ=255)      aCCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCAt                                                                                                                            >  2:380385/1‑139 (MQ=255)         tGGCGTCAGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACAGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGcc                                                                                                                         >  1:16419/1‑139 (MQ=255)                agaccTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTCTAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACAt                                                                                                                   >  1:465579/4‑138 (MQ=255)                aGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACAt                                                                                                                   <  2:465579/138‑1 (MQ=255)                              gCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGttt                                                                                                    >  1:437460/1‑139 (MQ=255)                                                gCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTgg                                                                                      <  2:803292/135‑1 (MQ=255)                                                gCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTgg                                                                                      >  1:803292/1‑135 (MQ=255)                                                        cATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACaaa                                                                          >  2:477071/1‑139 (MQ=255)                                                                tCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAgcagca                                                                  >  2:500214/1‑139 (MQ=255)                                                                tCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAgcagca                                                                  >  2:296404/1‑139 (MQ=255)                                                                     taAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTtgc                                                             >  1:226867/1‑139 (MQ=255)                                                                      aaaCGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTtgct                                                            <  1:639175/139‑1 (MQ=255)                                                                      aaaCGCTTTACCGGACACATGACATTCAAAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTtgct                                                            <  1:380385/139‑1 (MQ=255)                                                                             ttACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAGTCATGCCGCACATCGACATTAAATGTTTTCCGCGTGaa                                                                                          >  1:522938/1‑102 (MQ=255)                                                                             ttACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGaa                                                                                          <  2:522938/102‑1 (MQ=255)                                                                              tACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGatat                                                    <  1:296404/139‑1 (MQ=255)                                                                                 cGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTAc                                                 <  2:194482/139‑1 (MQ=255)                                                                                    acacATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGa                                              >  1:79288/1‑139 (MQ=255)                                                                                        aTGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAagcagc                                                                   >  2:324890/1‑114 (MQ=255)                                                                                        aTGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAagcagc                                                                   <  1:324890/114‑1 (MQ=255)                                                                                          gACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGttat                                        <  2:226867/139‑1 (MQ=255)                                                                                               tCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTc                                   <  2:79288/139‑1 (MQ=255)                                                                                                  caATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAAGAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTc                                >  1:243118/1‑139 (MQ=255)                                                                                                  caATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTc                                <  2:437460/139‑1 (MQ=255)                                                                                                    aTGTGTGTTAATTTTACTGCCATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTCAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCAt                              <  1:526173/139‑1 (MQ=255)                                                                                                     tgtgtgTTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGTCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATc                             >  1:286013/1‑139 (MQ=255)                                                                                                        gtgtTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGa                          <  1:434611/139‑1 (MQ=255)                                                                                                         tgtTAATTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGaa                         >  2:358200/1‑139 (MQ=255)                                                                                                             aaTTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGa                          <  2:450884/134‑1 (MQ=255)                                                                                                             aaTTTTACTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGa                          >  1:450884/1‑134 (MQ=255)                                                                                                                   aCTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCAt                              >  2:713196/1‑124 (MQ=255)                                                                                                                   aCTGACATCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCAt                              <  1:713196/124‑1 (MQ=255)                                                                                                                          tCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGa                          <  2:397829/121‑1 (MQ=255)                                                                                                                          tCTTTCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGa                          >  1:397829/1‑121 (MQ=255)                                                                                                                              tCAAAAGGAGTGTAATCATGTCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAg    <  2:286013/139‑1 (MQ=255)                                                                                                                              tCAAAAGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAg    <  2:614635/139‑1 (MQ=255)                                                                                                                                aaaaGGCGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCGCTTGCTGCAGACATTACCGACGTTATTATGCGTCATCTGAACAGTAAACACAGTTCGATAAGCa  <  1:41778/139‑1 (MQ=255)                                                                                                                                   aGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGa                                              >  2:429025/1‑92 (MQ=255)                                                                                                                                   aGGAGTGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGa                                              <  1:429025/92‑1 (MQ=255)                                                                                                                                        |                                                                                                                                   AAACTACCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTCGTAATCTTGCCCTGGGCATTATTGTCGTATAAACGCTTTACCGGACACATGACATTCACAATGTGTGTTAATTTTACTGACATCTTTCAAAAGGAGCGTAATCATGCCGCACATCGACATTAAATGTTTTCCGCGTGAACTGGACGAACAACAAAAAGCAGCACTTGCTGCAGATATTACCGACGTTATTATTCGTCATCTGAACAGTAAAGACAGTTCGATAAGCA  >  NC_000913/1532910‑1533175

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