Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,205,324	C→T	*311* (TAG→TAA)	ttdR ←	transcriptional activator of ttdABT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,205,324	0	C	T	98.4%	194.8 / NA	62	*311* (TAG→TAA)	ttdR	transcriptional activator of ttdABT
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (27/34);  minor base G (1/0);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCGATGCTCTCTTGCCTGATCCTGCTGCGTCATCATGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATCTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCAGACACGTAATTTCATGCTGCGATAGAGCGGCTCCCGGTAAACAGCCCAGATATTGGCGCTT  >  NC_000913/3205190‑3205456                                                                                                                                       |                                                                                                                                     tttCGATGCTCTCTTGCCTGATCCTGCTGCGTCATCATGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACa                                                                                                                                  <  2:1044179/139‑1 (MQ=255) tttCGATGCTCTCTTGCCTGATCCTGCTGCGTCATCATGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACa                                                                                                                                  <  2:219666/139‑1 (MQ=255) tttCGATGCTCCCTTGCCTGATCCTGCTGCGTCATCATGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACa                                                                                                                                  <  2:1080739/139‑1 (MQ=255)                   tGATCCTGCTGCGTCATCATGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCAt                                                                                                                <  2:1204289/139‑1 (MQ=255)                        ctgctgCGTCATCATGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGc                                                                                                           >  2:1070276/1‑139 (MQ=255)                          gctgcGTCATCATGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTa                                                                                                                                    >  2:261118/1‑112 (MQ=255)                             gcGTCATCATGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGcc                                                                                                      <  2:402344/139‑1 (MQ=255)                              cGTCATCATGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGccc                                                                                                     <  1:631839/139‑1 (MQ=255)                              cGTCATCATGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGccc                                                                                                     <  1:1070276/139‑1 (MQ=255)                                tcatcaTGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAg                                                                                                   <  2:1472980/139‑1 (MQ=255)                                tcatcaTGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAg                                                                                                   >  1:1472980/1‑139 (MQ=255)                                tcatcaTGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAg                                                                                                   >  2:1382393/1‑139 (MQ=255)                                 catcatGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACa                                                                                                                                  <  1:1504527/107‑1 (MQ=255)                                 catcatGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACa                                                                                                                                  >  2:1504527/1‑107 (MQ=255)                                 catcatGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGc                                                                                                  <  1:1382393/139‑1 (MQ=255)                                  atcatGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGcc                                                                                                 <  1:896170/139‑1 (MQ=255)                                  atcatGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGcc                                                                                                 <  1:763312/139‑1 (MQ=255)                                   tcatGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCg                                                                                                >  1:717833/1‑139 (MQ=255)                                           acaTCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGc                                                                                        <  2:162085/139‑1 (MQ=255)                                            caTCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCa                                                                                       >  2:328902/1‑139 (MQ=255)                                                 aaCGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGt                                                                                               >  2:393247/1‑126 (MQ=255)                                                 aaCGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGt                                                                                               <  1:393247/126‑1 (MQ=255)                                                 aaCGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCACAGCCGTTGCTGGCACCATg                                                                                  >  2:1192472/1‑139 (MQ=255)                                                   cGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATgcc                                                                                >  1:362364/1‑139 (MQ=255)                                                       tgtgGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCa                                                                            >  1:1166907/1‑139 (MQ=255)                                                       tgtgGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCa                                                                            >  1:257962/1‑139 (MQ=255)                                                       tgtgGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCa                                                                            <  2:717833/139‑1 (MQ=255)                                                       tgtgGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCa                                                                            >  2:110774/1‑139 (MQ=255)                                                         tgGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGa                                                                          >  1:1168573/1‑139 (MQ=255)                                                              tcgtcAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGc                                                                                                       <  2:770581/105‑1 (MQ=255)                                                              tcgtcAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGc                                                                                                       >  1:770581/1‑105 (MQ=255)                                                              tcgtcAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATGTACATGACCTGATAGCCTTGATCGGGCTTG‑CCAGGCGTTGCTGGCACCATGCCGCCAGAAATTcc                                                                    >  2:262680/1‑139 (MQ=255)                                                                gtcAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCa                                                                   >  2:631176/1‑139 (MQ=255)                                                                gtcAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCa                                                                   >  1:1441901/1‑139 (MQ=255)                                                                gtcAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCa                                                                   <  2:1508321/139‑1 (MQ=255)                                                                gtcAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCa                                                                   <  1:110774/139‑1 (MQ=255)                                                                      agaCAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCAGa                                                             <  2:257962/139‑1 (MQ=255)                                                                          aaaaaTCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCAGACACg                                                         >  1:967440/1‑139 (MQ=255)                                                                            aaaTCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCAGACACGTa                                                       <  2:1166907/139‑1 (MQ=255)                                                                                    aCGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATATCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCAGACACGTAATTTCAtg                                               <  2:967440/139‑1 (MQ=255)                                                                                               aagaaagCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCAGACACGTAATTTCATGCTGCGATagag                                    <  2:362364/139‑1 (MQ=255)                                                                                                    agCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCAGACACGTAATTTCATGCTGCGATAGAGCGGCt                               >  2:904443/1‑139 (MQ=255)                                                                                                      cgcgAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCAGACACGTAATTTCATGCTGCGATAGAGCGGc                                >  2:476009/1‑136 (MQ=255)                                                                                                      cgcgAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCAGACACGTAATTTCATGCTGCGATAGAGCGGc                                <  1:476009/136‑1 (MQ=255)                                                                                                         gAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGc                                                                                                                      <  2:865345/47‑1 (MQ=255)                                                                                                         gAAAAGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGc                                                                                                                      >  1:865345/1‑47 (MQ=255)                                                                                                             aGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCAGACACGTAATTTCATGCTGCGATAGAGCGGCTCCCGGTaaa                      <  1:631176/139‑1 (MQ=255)                                                                                                             aGGATCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCAGACACGTAATGTCATGCTGCGATAGAGCGGCTCCCGGTaaa                      <  2:1441901/139‑1 (MQ=255)                                                                                                                aTCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGt                                                                                               >  1:270923/1‑63 (MQ=255)                                                                                                                aTCCCGAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGt                                                                                               <  2:270923/63‑1 (MQ=255)                                                                                                                     gAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCAGACACGTAATTTCATGCTGCGATAGAGCGGCTCCCGGTAAACAGCCCAg              >  1:557384/1‑139 (MQ=255)                                                                                                                     gAGTGATTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCAGACACGTAATTTCATGCTGCGATAGAGCGGCTCCCGGTAAACAGCCCAg              >  2:1354748/1‑139 (MQ=255)                                                                                                                          aTTTCTGGTGGATTTACATGACCTGATAGCCTTCATCggg                                                                                                            <  2:281945/40‑1 (MQ=255)                                                                                                                          aTTTCTGGTGGATTTACATGACCTGATAGCCTTCATCggg                                                                                                            >  1:281945/1‑40 (MQ=255)                                                                                                                          aTTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGc                                                                                                           >  2:1417653/1‑41 (MQ=255)                                                                                                                          aTTTCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGc                                                                                                           <  1:1417653/41‑1 (MQ=255)                                                                                                                            ttCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCAGACACGTAATTTCATGCTGCGATAGAGCGGCTCCCGGTAAACAGCCCAGATATTgg       <  2:557384/139‑1 (MQ=255)                                                                                                                             tCTGGTGGATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCAGACACGTAATTTCATGCTGCGATAGAGCGGCTCCCGGTAAACAGCCCAGATATTGgc      <  2:1168573/139‑1 (MQ=255)                                                                                                                                ggtggATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCAGACACGTAATTTCATGCTGCGATAGAGCGGCTCCCGGTAAACAGCCCAGATATTGGCGCt   <  1:328902/139‑1 (MQ=255)                                                                                                                                 gtggATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCAGACACGTAATTTCATGCTGCGATAGAGCGGCTCCCGGTAAACAGCCCAGATATTGGCGCtt  <  2:1442342/139‑1 (MQ=255)                                                                                                                                  tggATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCa                                                               >  1:77879/1‑77 (MQ=255)                                                                                                                                  tggATTTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCa                                                               <  2:77879/77‑1 (MQ=255)                                                                                                                                       |                                                                                                                                     TTTCGATGCTCTCTTGCCTGATCCTGCTGCGTCATCATGACAACATCCAACGTCTGTGGCGTCGTCAGGAGACAAAAATCTGGACGAAATTCAAAAGAAAGCGCGAAAAGGATCCCGAGTGATTTCTGGTGGATCTACATGACCTGATAGCCTTCATCGGGCTTGCCCAGCCGTTGCTGGCACCATGCCGCCAGAAATTCCACGCAGACACGTAATTTCATGCTGCGATAGAGCGGCTCCCGGTAAACAGCCCAGATATTGGCGCTT  >  NC_000913/3205190‑3205456

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