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.3%	70.3 / NA	28	T138I (ACC→ATC)	tehB	tellurite methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (8/18);  total (8/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCT  >  NC_000913/1501727‑1502210                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             gatgCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACTATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAAGATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCTATGGATATCGc                                                                                                                                                                                                                                           <  2:42111/248‑1 (MQ=255)  aCGCATGGGATAAAAATGCCAGTAGTATCGACAACGTCGAGCGCAGTAAATCCATTGAAAATCTGGATAAGTTACACACCCGAGTCTTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGTCGATGGATATCGCt                                                                                                                                                                                                                                          <  1:256317/251‑1 (MQ=255)                   gacattaGTATCGCCACCGTCGAGAGCATTAAATCCATTGAAAATCTGTATTATTTACACTCCCGAGTCGTTGATCTGAATAACCTCACAATTGATAGACAGTACGATTTTATTCTTTTGACTGTTGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACCAC‑CTGATTGTGGCGGCGATGGATATCGCTTATTATCCATGTACCGt                                                                                                                                                                                                                         <  2:550124/245‑1 (MQ=255)                    ccATGAGTTTCGCCAACGTCGATCGCATTATATCCATTGAAAATCTGGATAATTTACACTCCCGAGTCGTTGATCTTAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTTCTGATGTTCCTTGAGGCTAAAACCAGCCCCGGGTTGATTGCCAATATGCAACTTTGCACTAAACCTGGTGGTTACATC‑CTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTc                                                                                                                                                                                                                        <  1:434777/251‑1 (MQ=255)                                   acggcGAGCGCATTAAAATCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTTCGATTGTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAGCCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGcc                                                                                                                                                                                                         <  1:335778/247‑1 (MQ=255)                                       cGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTc                                                                                                                                                                                                      <  1:374722/250‑1 (MQ=255)                                         aGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCaa                                                                                                                                                                                                    >  1:520831/1‑250 (MQ=255)                                                 aaaTCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAggg                                                                                                                                                                                              >  1:428910/1‑248 (MQ=255)                                                                       aTTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGc                                                                                                                                                                     <  1:309706/251‑1 (MQ=255)                                                                           acacacCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCGGTGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGGTATTACGAAGGCTggg                                                                                                                                                                 >  1:499633/1‑251 (MQ=255)                                                                                      tCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGGGAGGGGGaa                                                                                                                                                       >  1:127608/1‑250 (MQ=255)                                                                                                gAATAACCTCACATTTGATAGACAGTACGATTTTATTCGTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGAGTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAg                                                                                                                                            <  1:208788/251‑1 (MQ=255)                                                                                                   tAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGt                                                                                                                                                                                                                              >  1:379200/1‑166 (MQ=255)                                                                                                   tAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACCAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGt                                                                                                                                                                                                                              <  2:379200/166‑1 (MQ=255)                                                                                                     aCCTCCCATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGCCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAAAGAAGACGTc                                                                                                                                       <  2:236713/251‑1 (MQ=255)                                                                                                     aCCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTc                                                                                                                                       <  2:347107/251‑1 (MQ=255)                                                                                                                 aTAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAagag                                                                                                                                                                                                <  1:230120/182‑1 (MQ=255)                                                                                                                 aTAGACAGTACGATGTTATTCTTTCGACTGTGGTTATGATGTTCCGTGAGGCTAAACCCATCCGCGGGGTGATTGCCCATATGTAACGTTGCACTAAACCTTGTGGTTACAAC‑CTGATTGTGGCGGTGATGGATATCGCTGATTATCCATGTACCGTCGGCGTCCCGTTTGCCTTCAAAGAGGGAGAAGTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATGCAATGAAGACGTCGGCGCGCTgcac                                                                                                                           <  2:157170/251‑1 (MQ=255)                                                                                                                   aGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAACCCAGCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCGATGGATAACGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATGACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTgcatc                                                                                                                          <  2:423268/250‑3 (MQ=255)                                                                                                                                       ttCGACTGTGGTGCGGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGGACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTaa                                                                                                       <  2:499633/249‑1 (MQ=255)                                                                                                                                                     tGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGGGAAATACAATGAAGACGGCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGt                                                                                       >  1:124289/1‑251 (MQ=255)                                                                                                                                                                    acgcttgtccgatctggTTGATTGTCAATATGCAACGTTGCACTTAACCTGGTGTTTACAAC‑CTGATTGTGGCGGCGATGGATATCTCTGATTATCCATGTACCGTCGGTTGCACTTTTGCCTTCAAAGAGGGAGAATGACGGTGATAAGCCGAAGGCTGGGATAGGTTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTAGTAAACTGCGTTTCGCCACGATGc                                                                          <  2:83603/234‑1 (MQ=255)                                                                                                                                                                      aaaaCCATCCCCGGGTTGATTGCCAATATGCAACGTTGCTCTAAACCTGGTGGTTACTACTCT‑ATTGGGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCTTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCTAGCTGCACCGCACCGACGCCAGCGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGc                                                                      <  2:451332/251‑1 (MQ=255)                                                                                                                                                                                   ggTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCGATGGATATCGCTGATTATCCCTGTAACGTCGGCTTCCCGCTTGCCTTCACAGAGGGAGAATTACGTCGATATTACGGAGGCTGGGAGAGGGGGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGCATTAAACTGCGTTTCGCCACGCTGcagatcggaagagccaac                                                         >  1:83603/1‑234 (MQ=255)                                                                                                                                                                                           gCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCTGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTaa                                                 <  2:127608/251‑1 (MQ=255)                                                                                                                                                                                                 atggaaCGTTGCTCTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACaa                                           <  1:298109/247‑1 (MQ=255)                                                                                                                                                                                                                                         gtgGCGGCGATGGATATCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGtt    >  1:163278/1‑251 (MQ=255)                                                                                                                                                                                                                                           ggcggcGATGGATATCTCTGATTCTACATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGTCAACGGTAATCGTATTAAACAGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTTCCCACGCTGAGATGATTTACTGttct  <  1:348709/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             GACGCATGGGATAAAAATGCCATGAGTATCGCCAACGTCGAGCGCATTAAATCCATTGAAAATCTGGATAATTTACACACCCGAGTCGTTGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAAC‑CTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCT  >  NC_000913/1501727‑1502210

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