Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,416,443	G→T	16.7%	V192V (GTC→GTA)	yfcI ←	transposase_31 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,416,443	0	G	T	16.7%	58.7 / 4.5	24	V192V (GTC→GTA)	yfcI	transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base G (9/11);  new base T (2/2);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

CTGTCTAAACCTCTATCCAGCATCTCCTGAGCAATACGCAGGGCTTCTTCGTGTTTGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTATCTCACCAATAAATGCACGAAAACGCTGGGCATCCCCTGTTTGTAATACGTAATTAAACAGGGCTTTTAGCTGTCTGTCATTAGTGTTCCCTGTAACTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCAACGGAAAAGCCGATGAATATATTTTGCGGGCTATAGCAGGCTCGGCAAATTCATCAAGCCAGCAGAGTGAATAAGGATAAGGACTTCTGCAACCATGATAAAACAGCATCGGGAGCACCAATGGAAGCTC  >  CP009273/2416189‑2416682                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                aggTCTAAACCTCTGTCCTGCATCGTCTGAGTAATTCGCAGGGCTTCTTCGTTTTTGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTATCTCACCAATAAATGCACGAAAACGCTGGGCATCCCCTGTTTGTAATACGTAATTAAACAGGGCTTTTAGCTGTCTGTCATTAGTGTTCCCTGTAACTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATcgcg                                                                                                                                                                                                                             <  1:120151/273‑1 (MQ=255)              tATCCAGCATCTCCTGAGCAATACGCAGGGCTTCTTCGTGTTTGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTATCTCACCAATAAATGCACGAAAACGCTGGGCATCCCCTGTTTGTAAAACGTAATTAAACAGGGCTTTTAGCTGTCTGTCATTAGTGTTCCCTGTAACTAGCAGCGAAACCAATTGGTCGACTAGACGCAACAGAttccttccccttg                                                                                                                                                                                                                     >  2:210393/1‑258 (MQ=255)              tATCCAGCATCTCCTGAGCAATACGCAGGGCTTCTTCGTGTTTGCCCTGCATTGCGCCTTCTCCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTATCTCACCAATAAATGCACGAAAACGCTGGGCCTCCCCTGTTTGTAATACGTAATTAAACAGAGCTTTTAGCTTTCTGTCATTAGTGTTCCCTGTAACTAGCATCGAAACAATTTGGTGGAATAATCCTAACatat                                                                                                                                                                                                                                 >  2:120147/1‑258 (MQ=255)                ttcatcATTTCCTTAGTAATTCGCAGGGCTTCTTCGTGTTTGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATTGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTATCTCACCAATAAATGCACGAAAACGCTGGGCATCCCCTGTTTGTAATACGTAATTAAACAGGGCTTTTAGCTGTCTGTCATTAGTGTTCCCTGTAACTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGttt                                                                                                                                                                                                              <  1:130409/270‑1 (MQ=255)                  cAGCATCTCCTGAGCAATACGCAGGGCTTCTTCGTGTTTGCCCTGCATTGCGCCTTCTTCACGTAAGCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTATCTCACCAATAAATGCACGAAAACGCTGGGCATCCCCTGTTTGTAATACGTAATTAAACAGGGCTTTTAGCTGTCTGTCATTAGTGTTCCCTGTAACTAGCAGCGAAACAATTTGGTCGACTAAt                                                                                                                                                                                                                                           <  1:210398/244‑1 (MQ=255)                                                                                                                                                       gCACGAAAACGCTGGGCATCCCCTGTTTGTACTACGTAATTAAACAGGGCTTTTAGCTGTCTGTCATTAGTGTTCCCTGTAACTAGCAGCTACACAATTTGGTCGACTAATCCCTACAGATCGCGCTGACGATTATGTTTCTGAATTTACTCCAATAACGCCATTTTCCGGTGTTGCATAAt                                                                                                                                                                    >  2:328371/1‑182 (MQ=255)                                                                                                                                                           gAAAACGCTGGGCATCCCCTGTTTGTAATACGTAATTAAACAGGGCTTTTAGCTGTCTGTCATTAGTGTTCCCTGTAACTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTCGCGGTGTTGCATAATCTCg                                                                                                                                                                <  1:328377/182‑1 (MQ=255)                                                                                                                                                                       cATCCCCTGTTTGTAATACGTAATTAAACAGGGCTTTTAGCTGTCTGTCATTAGTGTTCCCTGTAACTAGCAGCGAAACAATTTGGTCGa                                                                                                                                                                                                                                                <  2:28089/90‑1 (MQ=255)                                                                                                                                                                              tGTTTGTAATACGTAATTAAACAGGGCTTTTAGCTGTCTGTCATTAGTGTTCCCTGTAACTAGCAGCGAAACAATTTGGTCGa                                                                                                                                                                                                                                                <  2:33565/83‑1 (MQ=255)                                                                                                                                                                                                      ggCTTTTAGCTGTCTGTCATTAGTGTTCCCTGTAACTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTCATATCCACCAACGGAAAAGCTGATGAATATATTTTGCGGGCTATATCAGGCACGGCAAATTCATCAAGCCAGCAGAGTGAATAAGGAAAAGAACTACTGCAACCATGAAAAAACaacc                        >  1:232674/1‑272 (MQ=255)                                                                                                                                                                                                                  tctgtcATTAGTGTTCCCTGTAACTAGCAGCGAAACAATTTGGTCTACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTTCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCAAc                                                                                                                                  >  2:204519/1‑157 (MQ=255)                                                                                                                                                                                                                      tcATTAGTGTTCCCTGTAACTAGCAGCGAAACAATTTGGTCTACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAAGCTCGTCATCCGGCACCACGGTAATATCCACCAACGGaa                                                                                                                              <  1:204524/157‑1 (MQ=255)                                                                                                                                                                                                                            gtgtTCCCTGTAACTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCAACGGAAAAGCCGATGAATATATTTTGCGGGCTATAGCAGGCTCGGCAAATTCATCAAGCCAGCAGAGTGAATAAGGATAAGGACTTCTGCAACCATGATAAAACAGCATCGGGAGCACCAATGGAAGctc  >  2:179520/1‑275 (MQ=255)                                                                                                                                                                                                                                tCCCTGTAACTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATATTCTCGTCATcc                                                                                                                                                          >  2:222761/1‑119 (MQ=255)                                                                                                                                                                                                                                tCCCTGAAACTAGCAGCGAAACAATTTGGTCGACt                                                                                                                                                                                                                                              <  1:309125/35‑1 (MQ=255)                                                                                                                                                                                                                                 cccTGTAACTAGCAGCGAAACAATTTGGTCTACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCg                                                                                                                                                         >  2:12376/1‑119 (MQ=255)                                                                                                                                                                                                                                    tGTAACTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGca                                                                                                                                                      <  1:222766/119‑1 (MQ=255)                                                                                                                                                                                                                                     gTAACTAGCAGCGAAACAATTTGGTCTACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGcac                                                                                                                                                     <  1:12376/119‑1 (MQ=255)                                                                                                                                                                                                                                              agcGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCAACGGAAAAGCCGATGAatatat                                                                                                              >  1:33044/1‑149 (MQ=255)                                                                                                                                                                                                                                                cGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCAACGGAAAAGCCGATGAATATATTTTGCGGGCTATAGCAGGCTCGGCAAATTCATCAAGCcag                                                                      >  2:254566/1‑187 (MQ=255)                                                                                                                                                                                                                                                  aaaCAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCg                                                                                                                                                         >  2:350248/1‑102 (MQ=255)                                                                                                                                                                                                                                                  aaaCAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCAACGGAAAAGCCGATGAATATATTTTg                                                                                                          <  2:33043/149‑1 (MQ=255)                                                                                                                                                                                                                                                    aCAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCAACGGAAAAGCCGATGAATATATTTTGCGGGCTATAGCAGGCTCGGCAAATTCATCAAGCCAGCaga                                                                  <  1:254571/187‑1 (MQ=255)                                                                                                                                                                                                                                                      aaTTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGcac                                                                                                                                                     <  1:350254/102‑1 (MQ=255)                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                CTGTCTAAACCTCTATCCAGCATCTCCTGAGCAATACGCAGGGCTTCTTCGTGTTTGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTATCTCACCAATAAATGCACGAAAACGCTGGGCATCCCCTGTTTGTAATACGTAATTAAACAGGGCTTTTAGCTGTCTGTCATTAGTGTTCCCTGTAACTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCAACGGAAAAGCCGATGAATATATTTTGCGGGCTATAGCAGGCTCGGCAAATTCATCAAGCCAGCAGAGTGAATAAGGATAAGGACTTCTGCAACCATGATAAAACAGCATCGGGAGCACCAATGGAAGCTC  >  CP009273/2416189‑2416682

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