Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,501,974	C→T	T138I (ACC→ATC)	tehB →	tellurite methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,501,974	0	C	T	96.9%	92.7 / NA	32	T138I (ACC→ATC)	tehB	tellurite methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (12/19);  minor base A (1/0);  total (13/19)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTC  >  NC_000913/1501731‑1502217                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    cATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGCGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGCTTGTGGCGCCGATGGATATCGCTGAt                                                                                                                                                                                                                                              >  2:219269/1‑251 (MQ=255)            aaTGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCAATGAAAATCGGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTAc                                                                                                                                                                                                                                   <  2:197765/251‑1 (MQ=255)                                      gcgcATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGAAATCGCTGATTATCCATGTACCGGCGGCTTCCCGTTTGCCTTCaaa                                                                                                                                                                                                          >  2:462797/1‑250 (MQ=255)                                          aTTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAggg                                                                                                                                                                                                     <  2:432444/251‑1 (MQ=255)                                                    ttGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATAACGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAAAGGGGGAAATACGt                                                                                                                                                                                           >  2:22270/1‑251 (MQ=255)                                                            cTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTACCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATAttt                                                                                                                                                                                   >  1:254018/1‑250 (MQ=255)                                                              ggATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACg                                                                                                                                                                                 <  1:462797/251‑1 (MQ=255)                                                                         acacCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGgag                                                                                                                                                                      <  1:426773/251‑1 (MQ=255)                                                                         acacCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGgag                                                                                                                                                                      <  1:520769/251‑1 (MQ=255)                                                                          cacCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACATGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTTCCTTCAAAGAGGTAGAATTACGTCGATATTACGAAGGCTGGgaga                                                                                                                                                                     <  1:540375/251‑1 (MQ=255)                                                                                         tCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTTCGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATg                                                                                                                                                      >  2:220356/1‑251 (MQ=255)                                                                                                   cTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATTGAGATCGCTGATTATCCCTGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGGGAATTACGTCGATATTTCGAAGGCTGGGGGAGGGGGAAATACAATGAAggcgccgg                                                                                                                                            >  1:285775/1‑244 (MQ=255)                                                                                                          ttGATATACAGTACGATTTTATTCTTTCGTCTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATATCGCTGATTATCCAGGTACCGACGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGcg                                                                                                                                      <  2:440289/250‑2 (MQ=255)                                                                                                                   aGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTgcaccgcacc                                                                                                                            <  1:219269/251‑1 (MQ=255)                                                                                                                       cGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGTCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGGGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACg                                                                                                                        <  2:573742/251‑1 (MQ=255)                                                                                                                             tATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAAc                                                                                                                   <  1:22270/250‑1 (MQ=255)                                                                                                                                     cGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGt                                                                                                          <  1:321495/251‑1 (MQ=255)                                                                                                                                         tgtgGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGGGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAAAGGTAATCGAAtt                                                                                                       >  2:525656/1‑250 (MQ=255)                                                                                                                                                   aTGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGttt                                                                                            <  1:43027/251‑1 (MQ=255)                                                                                                                                                            gAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGAt                                                                                   <  2:118483/251‑1 (MQ=255)                                                                                                                                                                     aCCATTCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTGATCGTATTAAACTGCGTTTCGCCACGATGATGGCACg                                                                          >  1:552097/1‑251 (MQ=255)                                                                                                                                                                      ccATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGTTGAAATACAATGAAGACGTCTGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGt                                                                         <  1:585493/251‑1 (MQ=255)                                                                                                                                                                                           atatGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATGGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTGATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCggggatggt                                                    <  2:552097/251‑8 (MQ=255)                                                                                                                                                                                                 aaCGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAAcg                                              <  1:368584/251‑1 (MQ=255)                                                                                                                                                                                                                           aaCCTGATTGTGGCGGCGATGGATATCCCTGATTATCCATTTACCGTCGGCTTCCCGTTTGCCTTCCAAGAGGGGGAATTACGTCGCTATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCCCCGACGCCAACCGTAATCGTATTAAACTGCGTTTCGCCACGATGATGGCACGTAAAAAATGACCCGGTACACCCAAAACGCGTGAAAATTCCCCAcccttcgatgg                    >  2:556010/1‑241 (MQ=255)                                                                                                                                                                                                                            aCCTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAgatgat                   >  1:3352/1‑251 (MQ=255)                                                                                                                                                                                                                              cTGATTGTGGCGGCGATGGATATCTCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTTAAAGAGGGAGAATTACGGCGATATTACGAAGGCTGGGAGAGGGGGAAATACAATGAAGACGTCGGCGAGCTGCACCGCAACGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCCCGCTGCTGGCACGTAAAAAATGACCCGGGAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGAttt                 >  2:456640/1‑251 (MQ=255)                                                                                                                                                                                                                                gATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGAttg                 <  2:309318/249‑2 (MQ=255)                                                                                                                                                                                                                                   tgtgGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAAAGGTAATCGTATTAAACTGCGTTTCGCCACGATGCCGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTg             <  1:267807/250‑1 (MQ=255)                                                                                                                                                                                                                                    gtgGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGt            <  1:456640/250‑1 (MQ=255)                                                                                                                                                                                                                                           cgATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCttt    >  2:575655/1‑251 (MQ=255)                                                                                                                                                                                                                                             aTGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTc  >  2:130972/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    CATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTC  >  NC_000913/1501731‑1502217

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