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%	203.7 / NA	63	*313* (TAG→TAA)	tdcA	tdc operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/28);  total (35/28)

ATATGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAACTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCATATTGTGCCACAGG  >  NC_000913/3266021‑3266261                                                                                                           |                                                                                                                                       atatGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTaca                                                                                                        >  2:70871/1‑139 (MQ=255)  tatGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGcc                                                                                                                  <  2:185249/128‑1 (MQ=255)  tatGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGcc                                                                                                                  >  1:185249/1‑128 (MQ=255)  tatGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTacac                                                                                                       >  2:1004355/1‑139 (MQ=255)  tatGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTacac                                                                                                       >  1:190095/1‑139 (MQ=255)   atGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACAcc                                                                                                      <  2:190095/139‑1 (MQ=255)    tGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACAccc                                                                                                     <  2:384257/139‑1 (MQ=255)       aTTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTaca                                                                                                        >  1:247130/1‑133 (MQ=255)       aTTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTaca                                                                                                        <  2:247130/133‑1 (MQ=255)         tataAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACAcc                                                                                                      >  1:194775/1‑133 (MQ=255)         tataAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACAcc                                                                                                      <  2:194775/133‑1 (MQ=255)                 ctctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGaa                                                                                        <  1:1004355/139‑1 (MQ=255)                   ctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGCCACCAGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCaa                                                                                                                                 <  2:887086/96‑1 (MQ=255)                   ctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCaa                                                                                                                                 >  1:887086/1‑96 (MQ=255)                   ctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCt                                                                                                                 <  1:1011928/112‑1 (MQ=255)                   ctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCt                                                                                                                 >  2:1011928/1‑112 (MQ=255)                   ctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATa                                                                                      >  1:1301164/1‑139 (MQ=255)                   ctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATa                                                                                      >  1:680332/1‑139 (MQ=255)                     tCGCCTGAATTAAATATGTCGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTcca                                                                                                          <  2:806713/117‑3 (MQ=255)                     tCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCt                                                                                                                 >  2:486815/1‑110 (MQ=255)                     tCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCt                                                                                                                 <  1:486815/110‑1 (MQ=255)                     tCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTa                                                                                                          >  1:806713/1‑117 (MQ=255)                           gAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTgg                                                                              >  2:860983/1‑139 (MQ=255)                              ttAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCa                                                                                                                                  <  2:885013/84‑1 (MQ=255)                              ttAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCa                                                                                                                                  >  1:885013/1‑84 (MQ=255)                                      gTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTccacc                                                                   >  2:1112468/1‑139 (MQ=255)                                      gTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGGTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAg                                                                                             <  2:275474/113‑1 (MQ=255)                                      gTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGGTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAg                                                                                             >  1:275474/1‑113 (MQ=255)                                        aGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCaa                                                                 >  1:1236214/1‑139 (MQ=255)                                        aGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCaa                                                                 >  1:87421/1‑139 (MQ=255)                                           taaccgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGcc                                                                                                                  <  2:539137/87‑1 (MQ=255)                                           taaccgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGcc                                                                                                                  >  1:539137/1‑87 (MQ=255)                                              ccgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAAcc                                                                                                                                   >  1:766217/1‑67 (MQ=255)                                              ccgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAAcc                                                                                                                                   <  2:766217/67‑1 (MQ=255)                                              ccgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGa                                                           >  1:1107264/1‑139 (MQ=255)                                               cgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCt                                                                                                           <  1:725095/90‑1 (MQ=255)                                               cgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCt                                                                                                           >  2:725095/1‑90 (MQ=255)                                               cgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTacac                                                                                                       >  1:781024/1‑94 (MQ=255)                                               cgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTacac                                                                                                       <  2:781024/94‑1 (MQ=255)                                               cgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGAt                                                          >  1:561012/1‑139 (MQ=255)                                                 taaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGAtgc                                                        >  2:604582/1‑139 (MQ=255)                                                   accgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCt                                                                                                                 >  1:382163/1‑80 (MQ=255)                                                   accgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCt                                                                                                                 <  2:382163/80‑1 (MQ=255)                                                    ccgACACCTACAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGAtgctgc                                                     <  1:860983/139‑1 (MQ=255)                                                         acCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAAcc                                                                                                                                   <  2:235806/56‑1 (MQ=255)                                                         acCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAAcc                                                                                                                                   >  1:235806/1‑56 (MQ=255)                                                         acCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATaa                                                                                              >  2:642818/1‑93 (MQ=255)                                                         acCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATaa                                                                                              <  1:642818/93‑1 (MQ=255)                                                                   cagacagGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTacac                                                                                                       >  1:1130659/1‑74 (MQ=255)                                                                   cagacagGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTacac                                                                                                       <  2:1130659/74‑1 (MQ=255)                                                                      acagGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAAtt                                   >  2:464201/1‑139 (MQ=255)                                                                        agGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCt                                                                                                           <  1:1016317/65‑1 (MQ=255)                                                                        agGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCt                                                                                                           >  2:1016317/1‑65 (MQ=255)                                                                            ggATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGa                             <  1:604582/139‑1 (MQ=255)                                                                              attattTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGAcc                           >  1:557931/1‑139 (MQ=255)                                                                               ttatttatGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCa                          <  1:464201/139‑1 (MQ=255)                                                                                 atttatGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATa                        <  1:1112468/139‑1 (MQ=255)                                                                                        tttAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGc                 >  1:1296801/1‑139 (MQ=255)                                                                                            aataaaAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTAGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATa                                           <  2:276219/109‑1 (MQ=255)                                                                                            aataaaAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTAGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATa                                           >  1:276219/1‑109 (MQ=255)                                                                                                 aaaCAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCATATTGTGc        <  2:943009/139‑1 (MQ=255)                                                                                                  aaCAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCATATTGTGcc       <  2:680332/139‑1 (MQ=255)                                                                                                       aTAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCATATTGTGCCACAgg  >  2:1196015/1‑139 (MQ=255)                                                                                                           |                                                                                                                                       ATATGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAACTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCATATTGTGCCACAGG  >  NC_000913/3266021‑3266261

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