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	97.1%	212.2 / ‑2.0	70	*313* (TAG→TAA)	tdcA	tdc operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (39/29);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-01

AGCAACCGGCAGATCGTATGTAATATGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAACTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCATATTGTGCCA  >  NC_000913/3265999‑3266257                                                                                                                                 |                                                                                                                                   aGCAACCGGCAGATCGTATGTAATATGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTA‑TTATGTTTAAATAAAAACAAATAATTAACCAACTTc                                                                                                                         >  1:682773/1‑139 (MQ=255)   cAACCGGCAGATCGTATGTAATATGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCt                                                                                                                        >  2:55758/1‑139 (MQ=255)        ggCAGATCGTATGTAATATGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTaa                                                                                                                   >  2:950369/1‑139 (MQ=255)            gATCGTATGTAATATGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTA‑TTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGcc                                                                                                              <  2:682773/139‑1 (MQ=255)                        tatGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTacac                                                                                                   <  1:950369/139‑1 (MQ=255)                           gCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCtt                                                                                                            <  2:1607416/127‑1 (MQ=255)                           gCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCtt                                                                                                            >  1:1607416/1‑127 (MQ=255)                           gCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACTCa                                                                                                >  1:21486/1‑139 (MQ=255)                               tataAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACtt                                                                                                                          <  2:1634032/109‑1 (MQ=255)                               tataAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACtt                                                                                                                          >  1:1634032/1‑109 (MQ=255)                               tataAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATtat                                                                                            <  1:637103/139‑1 (MQ=255)                                  aaaaCCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAAtaattaa                                                                                                                                 <  1:630905/99‑1 (MQ=255)                                  aaaaCCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAAtaattaa                                                                                                                                 >  2:630905/1‑99 (MQ=255)                                  aaaaCCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAg                                                                                         >  2:1661989/1‑139 (MQ=255)                                    aaCCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAAcc                                                                                                                               <  1:576999/99‑1 (MQ=255)                                    aaCCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAAcc                                                                                                                               >  2:576999/1‑99 (MQ=255)                                    aaCCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGAt                                                                                       >  2:1468526/1‑139 (MQ=255)                                    aaCCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGAt                                                                                       >  1:1271922/1‑139 (MQ=255)                                    aaCCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGAt                                                                                       >  2:565006/1‑139 (MQ=255)                                         ctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATaa                                                                                          <  2:378914/131‑1 (MQ=255)                                         ctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATaa                                                                                          >  1:378914/1‑131 (MQ=255)                                         ctTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATa                                                                                  >  2:1433416/1‑139 (MQ=255)                                             gCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATActct                                                                              <  1:55758/139‑1 (MQ=255)                                                tGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAAcc                                                                                                                               >  2:1030365/1‑87 (MQ=255)                                                tGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAAcc                                                                                                                               <  1:1030365/87‑1 (MQ=255)                                                  aaTTAAAGATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTAAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTAGAAGATGAATACTCTTTGGc                                                                         >  1:1199762/1‑139 (MQ=255)                                                            gTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAAtaattaa                                                                                                                                 >  2:1662183/1‑73 (MQ=255)                                                            gTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAAtaattaa                                                                                                                                 <  1:1662183/73‑1 (MQ=255)                                                            gTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACtt                                                                                                                          >  1:462699/1‑80 (MQ=255)                                                            gTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACtt                                                                                                                          <  2:462699/80‑1 (MQ=255)                                                            gTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTc                                                                                                                         >  1:11909/1‑81 (MQ=255)                                                            gTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTc                                                                                                                         <  2:11909/81‑1 (MQ=255)                                                            gTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAACTAACCAACTTCTATTAATTGCCt                                                                                                             >  1:326177/1‑93 (MQ=255)                                                            gTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAACTAACCAACTTCTATTAATTGCCt                                                                                                             <  2:326177/93‑1 (MQ=255)                                                                 taaccgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTATTTGcc                                                                                                              <  2:920790/87‑1 (MQ=255)                                                                 taaccgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGcc                                                                                                              >  1:920790/1‑87 (MQ=255)                                                                   accgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATtat                                                                                            >  1:34875/1‑103 (MQ=255)                                                                   accgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATtat                                                                                            <  2:34875/103‑1 (MQ=255)                                                                     cgtaaccgACACCTGTAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGGCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGAt                                                      >  1:1223984/1‑139 (MQ=255)                                                                     cgtaaccgACACCTGCAAGACAGATAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTacac                                                                                                   <  2:1016921/94‑1 (MQ=255)                                                                     cgtaaccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTacac                                                                                                   >  1:1016921/1‑94 (MQ=255)                                                                         accgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCt                                                                                                       >  1:1311417/1‑86 (MQ=255)                                                                         accgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCt                                                                                                       <  2:1311417/86‑1 (MQ=255)                                                                          ccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCt                                                                                                             <  1:1761046/79‑1 (MQ=255)                                                                          ccgACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCt                                                                                                             >  2:1761046/1‑79 (MQ=255)                                                                              cacCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCtttt                                             >  2:216183/1‑139 (MQ=255)                                                                               acCTTCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCttttt                                            <  1:1433416/139‑1 (MQ=255)                                                                               acCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGcc                                                                                                              >  1:260559/1‑73 (MQ=255)                                                                               acCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGcc                                                                                                              <  2:260559/73‑1 (MQ=255)                                                                               acCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTacac                                                                                                   >  2:1560752/1‑84 (MQ=255)                                                                               acCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTacac                                                                                                   <  1:1560752/84‑1 (MQ=255)                                                                               acCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCttttt                                            >  2:256694/1‑139 (MQ=255)                                                                                    cAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTa                                                                                                      >  2:86854/1‑76 (MQ=255)                                                                                    cAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTa                                                                                                      <  1:86854/76‑1 (MQ=255)                                                                                     aaGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATAc                                      >  2:581086/1‑139 (MQ=255)                                                                                            tcagGTGGATTATTTATGTTTAAATAAAAACAAATAATTAAGCAACTTCTATTAATTGCCTTCGTCTACACCCGTTATAAGATGAAt                                                                                   >  2:1126475/2‑87 (MQ=255)                                                                                            acagGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGcc                                                                                                              >  1:1038475/1‑60 (MQ=255)                                                                                            acagGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGcc                                                                                                              <  2:1038475/60‑1 (MQ=255)                                                                                            acagGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAAt                                                                                   <  1:1126475/87‑1 (MQ=255)                                                                                            acagGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAAtt                               >  2:951111/1‑139 (MQ=255)                                                                                             cagGTGGATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAAttt                              <  1:1468526/139‑1 (MQ=255)                                                                                                 tGGATTATTTGTGTTTAAATAAAAACAAATAATTAACCAACTTCACTTAATTCCCTTCGTCTACACTCCTCATATGATGAATACTCTTTAGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGAAAATTATTCg                          >  1:1806382/1‑139 (MQ=255)                                                                                                   gATTATTTATGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGAc                        >  2:1437906/1‑139 (MQ=255)                                                                                                       atttatGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATa                    <  2:1087681/139‑1 (MQ=255)                                                                                                         ttatGTTTAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATAcc                  <  1:256694/139‑1 (MQ=255)                                                                                                              tttAAATAAAAACAAATAATTAACCAACTTCTATTAATTGcc                                                                                                              <  2:663447/42‑1 (MQ=255)                                                                                                              tttAAATAAAAACAAATAATTAACCAACTTCTATTAATTGcc                                                                                                              >  1:663447/1‑42 (MQ=255)                                                                                                              tttAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAAt                                >  1:1357391/1‑120 (MQ=255)                                                                                                              tttAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAAt                                <  2:1357391/120‑1 (MQ=255)                                                                                                              tttAAATAAAAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGc             >  2:1267483/1‑139 (MQ=255)                                                                                                                  aataaaAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCatat         >  2:723437/1‑139 (MQ=255)                                                                                                                  aataaaAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCatat         <  1:723437/139‑1 (MQ=255)                                                                                                                   ataaaAACAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCATAtt        >  2:234489/1‑139 (MQ=255)                                                                                                                         aCAAATAATTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCATATTGTGCca  <  1:234489/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                   AGCAACCGGCAGATCGTATGTAATATGCATTATAAAACCTCTTCGCCTGAATTAAATATGTAGGTAACCGTAACCGACACCTGCAAGACAGACAGGTGGATTATTTATGTTTAAATAAAAACAAATAACTAACCAACTTCTATTAATTGCCTTCGTCTACACCCATTATAAGATGAATACTCTTTGGCTAATTCCACCAAAACCGATGCTGCTTTTTTAATACGATAATTTTTCGACCATACCGCGGCATATTGTGCCA  >  NC_000913/3265999‑3266257

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