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

GCAGATCGTATGTAATATGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAACTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCATATTGTGCCACAGG  >  NC_000913/3266007‑3266261                                                                                                                         |                                                                                                                                       gCAGATCGTATGTAATATGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAAc                                                                                                                                <  2:1446084/129‑1 (MQ=255) gCAGATCGTATGTAATATGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAAc                                                                                                                                >  1:1446084/1‑129 (MQ=255)               atatGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTaca                                                                                                        >  2:937509/1‑139 (MQ=255)                tatGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGcc                                                                                                                  >  1:113808/1‑128 (MQ=255)                tatGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGcc                                                                                                                  <  2:113808/128‑1 (MQ=255)                tatGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCt                                                                                                                 <  2:1378340/129‑1 (MQ=255)                tatGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCt                                                                                                                 >  1:1378340/1‑129 (MQ=255)                tatGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTacac                                                                                                       <  2:658581/139‑1 (MQ=255)                 atGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACAcc                                                                                                      <  2:784705/139‑1 (MQ=255)                    cATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCAt                                                                                                   >  1:529078/1‑139 (MQ=255)                               ctctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGaa                                                                                        <  2:529078/139‑1 (MQ=255)                                 ctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAAcc                                                                                                                                   <  1:1151433/94‑1 (MQ=255)                                 ctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAAcc                                                                                                                                   >  2:1151433/1‑94 (MQ=255)                                 ctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCt                                                                                                                 <  2:947877/112‑1 (MQ=255)                                 ctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATa                                                                                      >  2:271343/1‑139 (MQ=255)                                 ctTCGCCTGAATTAAATATGCAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCt                                                                                                                 >  1:947877/1‑112 (MQ=255)                                   tCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACAcc                                                                                                      <  2:710737/121‑1 (MQ=255)                                   tCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACAcc                                                                                                      >  1:710737/1‑121 (MQ=255)                                          aaTTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGcc                                                                                                                  <  1:490805/102‑1 (MQ=255)                                          aaTTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGcc                                                                                                                  >  2:490805/1‑102 (MQ=255)                                          aaTTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACAcc                                                                                                      >  1:839948/1‑114 (MQ=255)                                          aaTTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACAcc                                                                                                      <  2:839948/114‑1 (MQ=255)                                          aaTTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGc                                                                             >  2:830534/1‑139 (MQ=255)                                                    gTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTccacc                                                                   >  2:1367725/1‑139 (MQ=255)                                                        gtaaccgtaaccGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAAcc                                                                                                                                   >  2:163061/1‑71 (MQ=255)                                                        gtaaccgtaaccGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAAcc                                                                                                                                   <  1:163061/71‑1 (MQ=255)                                                        gtaaccgtaaccGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTac                                                                                                         <  2:740888/97‑1 (MQ=255)                                                        gtaaccgtaaccGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTac                                                                                                         >  1:740888/1‑97 (MQ=255)                                                           accgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATtat                                                                                                >  2:1178031/1‑103 (MQ=255)                                                           accgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATtat                                                                                                <  1:1178031/103‑1 (MQ=255)                                                             cgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCtt                                                                                                                <  1:767370/85‑1 (MQ=255)                                                             cgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCtt                                                                                                                >  2:767370/1‑85 (MQ=255)                                                               taaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTaa                                                                                                                       >  2:1233280/1‑76 (MQ=255)                                                               taaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTaa                                                                                                                       <  1:1233280/76‑1 (MQ=255)                                                                                acagacagGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGAt                                       <  1:1367725/139‑1 (MQ=255)                                                                                   gacagGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAAt                                    >  1:64957/1‑139 (MQ=255)                                                                                    acagGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTacac                                                                                                       >  2:362988/1‑71 (MQ=255)                                                                                    acagGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTacac                                                                                                       <  1:362988/71‑1 (MQ=255)                                                                                    acagGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAAtt                                   >  2:167833/1‑139 (MQ=255)                                                                                    acagGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAAtt                                   >  1:854367/1‑139 (MQ=255)                                                                                    acagGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAAtt                                   >  1:594866/1‑139 (MQ=255)                                                                                         tGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATActct                                                                                  <  2:54585/87‑1 (MQ=255)                                                                                         tGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATActct                                                                                  >  1:54585/1‑87 (MQ=255)                                                                                           gATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGAt                                       >  2:996082/1‑128 (MQ=255)                                                                                           gATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGAt                                       <  1:996082/128‑1 (MQ=255)                                                                                           gATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGAc                            >  2:741410/1‑139 (MQ=255)                                                                                                  tatGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACcg                     >  1:313022/1‑139 (MQ=255)                                                                                                          aataaaAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATtat                                                                                                <  1:271880/56‑1 (MQ=255)                                                                                                          aataaaAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATtat                                                                                                >  2:271880/1‑56 (MQ=255)                                                                                                                     aTAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCATATTGTGCCACAgg  <  2:594866/139‑1 (MQ=255)                                                                                                                     aTAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCATATTGTGCCACAgg  >  2:1192295/1‑139 (MQ=255)                                                                                                                         |                                                                                                                                       GCAGATCGTATGTAATATGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAACTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCATATTGTGCCACAGG  >  NC_000913/3266007‑3266261

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