Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	1,964,699	C→T	V180I (GTC→ATC)	LY180_09345 ←	transposase IS609

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	1,964,699	0	C	T	95.8%	39.6 / NA	24	V180I (GTC→ATC)	LY180_09345	transposase IS609
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (13/10);  total (14/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGTGAACCGGAGTGGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACAGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCA  >  CP006584/1964419‑1964989                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                   ttGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGt                                                                                                                                                                                                                                                                                   >  1:190370/1‑298 (MQ=17)           tCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGt                                                                                                                                                                                                                                                                         <  1:313330/298‑1 (MQ=17)                                               tGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGGGACCCGGAGTTTATACTTCACTTTCTGTCTGAATAAGTATGTAcaa                                                                                                                                                                                                                                     >  2:317268/1‑296 (MQ=255)                                                             aCTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACt                                                                                                                                                                                                                       <  1:72514/298‑1 (MQ=17)                                                               tttATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTgg                                                                                                                                                                                                                     >  1:228721/1‑298 (MQ=17)                                                                                      tatTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCacaac                                                                                                                                                                                              >  1:208197/1‑298 (MQ=17)                                                                                                     aaTGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGCCCTGCTCACCGGGCCTTTTTTCCCAACCCCCGTgccactacg                                                                                                                                                                               >  2:10135/1‑290 (MQ=11)                                                                                                     aaTGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACg                                                                                                                                                                               <  1:29688/298‑1 (MQ=14)                                                                                                                            ttctgcttctgcCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGc                                                                                                                                                        <  1:317268/298‑1 (MQ=14)                                                                                                                                gcttctgcCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAg                                                                                                                                                    <  1:329776/298‑1 (MQ=14)                                                                                                                                                                                  tccgatctcAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCg                                                                                                  <  1:149204/290‑1 (MQ=12)                                                                                                                                                                                                      gTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTc                                                                              >  1:60377/1‑298 (MQ=17)                                                                                                                                                                                                          ttCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTc                                                                          <  1:10135/298‑1 (MQ=17)                                                                                                                                                                                                             tGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTcccc                                                                       >  1:152067/1‑298 (MQ=17)                                                                                                                                                                                                                                   aGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAagc                                                 >  1:249329/1‑298 (MQ=14)                                                                                                                                                                                                                                       tgtgCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGcc                                             >  1:314623/1‑298 (MQ=14)                                                                                                                                                                                                                                         tgCCACCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCg                                           <  1:5149/298‑1 (MQ=14)                                                                                                                                                                                                                                            cATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGCTGTACCACTTACCGCAGGACTGGCTGACCGTGACATTGTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTGCGCAGAAAAATACGGCCGCTTTCCTGCTCGAGCTTTACGCACCTCGGGTACAGGAAAGACCCATTCTTTCCCCCCTTTTTGCATCGGGGTAAggcagcagccgc                                        >  2:299899/1‑292 (MQ=11)                                                                                                                                                                                                                                                  aCAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGgaa                                  >  1:360387/1‑298 (MQ=14)                                                                                                                                                                                                                                                  aCAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGgaa                                  <  1:180292/298‑1 (MQ=14)                                                                                                                                                                                                                                                             aGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGGAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCCTTCTGTCCACGCTTTTTGAATCCGGGAAAAGAAGCCCGCCTCCTGAAGCAGTTTTTGt                       >  2:62927/1‑298 (MQ=11)                                                                                                                                                                                                                                                                  aGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAATGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGcc                  <  1:294799/298‑1 (MQ=14)                                                                                                                                                                                                                                                                      aCGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCt              >  1:76212/1‑298 (MQ=14)                                                                                                                                                                                                                                                                                  aGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCCGGCCGCCTCCGGAAGAAGTTTTTTTAAGCCCGCTCAAGGCCTTTcc  >  2:115965/1‑297 (MQ=9)                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                   TTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGTGAACCGGAGTGGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACAGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCA  >  CP006584/1964419‑1964989

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