Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,266,127	C→T	*313* (TAG→TAA)	tdcA ←	tdc operon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,266,127	0	C	T	100.0%	157.3 / NA	49	*313* (TAG→TAA)	tdcA	tdc operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/22);  total (27/22)

GCAACCGGCAGATCGTATGTAATATGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAACTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCATATTGTGCCA  >  NC_000913/3266000‑3266257                                                                                                                                |                                                                                                                                   gCAACCGGCAGATCGTATGTAATATGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACtt                                                                                                                          >  2:1497649/1‑138 (MQ=255) gCAACCGGCAGATCGTATGTAATATGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACtt                                                                                                                          <  1:1497649/138‑1 (MQ=255)    aCCGGCAGATCGTATGTAATATGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAAcc                                                                                                                               <  2:138381/130‑1 (MQ=255)    aCCGGCAGATCGTATGTAATATGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAAcc                                                                                                                               >  1:138381/1‑130 (MQ=255)                           cATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAAcc                                                                                                                               <  1:1067568/107‑1 (MQ=255)                           cATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAAcc                                                                                                                               >  2:1067568/1‑107 (MQ=255)                              tataAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAAc                                                                                                                            >  1:286373/1‑107 (MQ=255)                              tataAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAAc                                                                                                                            <  2:286373/107‑1 (MQ=255)                              tataAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATtat                                                                                            >  1:778030/1‑139 (MQ=255)                              tataAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATtat                                                                                            >  2:1544471/1‑139 (MQ=255)                                taAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATaa                                                                                          <  2:778030/139‑1 (MQ=255)                                   aaCCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTacac                                                                                                   <  2:474340/127‑1 (MQ=255)                                   aaCCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTacac                                                                                                   >  1:474340/1‑127 (MQ=255)                                        ctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAAtaattaa                                                                                                                                 <  2:1148430/92‑1 (MQ=255)                                        ctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAAtaattaa                                                                                                                                 <  1:1377927/92‑1 (MQ=255)                                        ctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAAtaattaa                                                                                                                                 >  1:1148430/1‑92 (MQ=255)                                        ctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAAtaattaa                                                                                                                                 >  2:1377927/1‑92 (MQ=255)                                        ctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTATACCCATTATaa                                                                                          <  1:670886/131‑1 (MQ=255)                                        ctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTATACCCATTATaa                                                                                          >  2:670886/1‑131 (MQ=255)                                          tCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAAcc                                                                                                                               >  2:1051724/1‑92 (MQ=255)                                          tCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAAcc                                                                                                                               <  1:1051724/92‑1 (MQ=255)                                                           gTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCtt                                                                                                            >  1:2037241/1‑94 (MQ=255)                                                           gTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCtt                                                                                                            <  2:2037241/94‑1 (MQ=255)                                                           gTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTccacc                                                               >  2:738518/1‑139 (MQ=255)                                                              ggTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTaa                                                                                                                   <  1:307699/84‑1 (MQ=255)                                                              ggTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTaa                                                                                                                   >  2:307699/1‑84 (MQ=255)                                                                taaccgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTaa                                                                                                                   >  2:815493/1‑82 (MQ=255)                                                                taaccgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTaa                                                                                                                   <  1:815493/82‑1 (MQ=255)                                                                      taaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGAtgc                                                    >  1:296030/1‑139 (MQ=255)                                                                        accgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGAtgctg                                                  >  2:1365213/1‑139 (MQ=255)                                                                        accgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGAtgctg                                                  >  2:178305/1‑139 (MQ=255)                                                                             cacCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACAcc                                                                                                  <  1:1511465/86‑1 (MQ=255)                                                                             cacCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACAcc                                                                                                  >  2:1511465/1‑86 (MQ=255)                                                                              acCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGcc                                                                                                              <  2:1078923/73‑1 (MQ=255)                                                                              acCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGcc                                                                                                              >  1:1078923/1‑73 (MQ=255)                                                                              acCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCt                                                                                                             <  1:228423/74‑1 (MQ=255)                                                                              acCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCt                                                                                                             >  2:228423/1‑74 (MQ=255)                                                                              acCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCttttt                                            >  1:1231162/1‑139 (MQ=255)                                                                                           acagGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCt                                                                                                       >  2:170037/1‑67 (MQ=255)                                                                                           acagGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCt                                                                                                       <  1:170037/67‑1 (MQ=255)                                                                                                  gATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGAc                        >  2:295943/1‑139 (MQ=255)                                                                                                   attattTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGAcc                       <  1:178305/139‑1 (MQ=255)                                                                                                     tatttatGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCAt                     >  1:711049/1‑139 (MQ=255)                                                                                                         tatGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACcg                 >  1:1055827/1‑139 (MQ=255)                                                                                                         tatGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACcg                 <  1:1544471/139‑1 (MQ=255)                                                                                                            gTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGAt                                   <  1:1772843/118‑1 (MQ=255)                                                                                                            gTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGAt                                   >  2:1772843/1‑118 (MQ=255)                                                                                                               taaataaaAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCat           >  2:1476820/1‑139 (MQ=255)                                                                                                                 aataaaAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCatat         >  1:696156/1‑139 (MQ=255)                                                                                                                   taaaAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCATATtg       <  1:1476820/139‑1 (MQ=255)                                                                                                                      aaaCAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCATATTGTGc    <  2:1231162/139‑1 (MQ=255)                                                                                                                       aaCAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCATATTGTGcc   <  1:295943/139‑1 (MQ=255)                                                                                                                        aCAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCATATTGTGCca  <  2:696156/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                   GCAACCGGCAGATCGTATGTAATATGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAACTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCATATTGTGCCA  >  NC_000913/3266000‑3266257

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