Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	2,414,801	C→G	22.8%	intergenic (+343/‑2)	PP_2114 → / → PP_2115	transposase/conserved protein of unknown function

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	2,414,801	0	C	G	22.8%	49.0 / 13.7	22	intergenic (+343/‑2)	PP_2114/PP_2115	transposase/conserved protein of unknown function
Reads supporting (aligned to +/- strand):  ref base C (8/9);  new base G (2/3);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-03

TCAGCGAACAGTCCATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCG  >  AE015451/2414586‑2414984                                                                                                                                                                                                                        |                                                                                                                                                                                        tCAGCGAACAGTCCATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTc                                                                                                                                                                     <  1:116290/236‑1 (MQ=38)       aaCAGTCCATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATg                                                                                                                                                                                     >  1:429356/1‑214 (MQ=39)       aaCAGTCCATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAgg                                                                                                                                                               >  1:391473/1‑236 (MQ=40)            tCCATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGc                                                                                                                                                          <  1:470236/236‑1 (MQ=38)                             cacTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGtt                                                                                                                                                       <  1:550008/222‑1 (MQ=34)                                      aCCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTc                                                                                                                                                                     <  1:56070/199‑1 (MQ=32)                                                          aGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAg                                                                                                                                                                <  2:189818/184‑1 (MQ=31)                                                                                                  cGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGcctcct                                                                                                                                              <  1:19967/162‑1 (MQ=31)                                                                                                                    cGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGCTCGGCCTCCTGAAGCCGCGAAGGTCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgagatcaagatcaagat                                                  >  1:175238/1‑236 (MQ=255)                                                                                                                     ggCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTc                                                                                                                                        <  1:250480/149‑1 (MQ=32)                                                                                                                       ccccTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCC‑‑‑‑GGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCACCTCCCTTCGCAAGCCCTACGCTCGGCCTCCTGAAGTcgcaatct                                                                                                                                 <  1:358405/150‑6 (MQ=11)                                                                                                                       ccccTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCC‑‑‑‑GGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCACCTCCCTTCGCAAGCCCTACGCTCGGCCTCCTGAAGTcgcaatct                                                                                                                                 >  2:358404/1‑145 (MQ=11)                                                                                                                         ccTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCACCTGCCTTc                                                                                                                                                                     >  1:486087/1‑116 (MQ=9)                                                                                                                         ccTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCACCTGCCTTc                                                                                                                                                                     <  2:486086/116‑1 (MQ=9)                                                                                                                                          ggTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaa                                                                                                                             >  2:260436/1‑139 (MQ=33)                                                                                                                                               ccTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGATTGa                          >  1:71540/1‑233 (MQ=255)                                                                                                                                                                      cAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGaa                         <  2:688259/211‑1 (MQ=255)                                                                                                                                                                                 cGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGAGTcgag  >  1:61621/1‑221 (MQ=255)                                                                                                                                                                                          cATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGAGCGAGATCAAGAGCGCTGGGTTTGaa                         <  1:260437/191‑1 (MQ=255)                                                                                                                                                                                                           gcTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGGCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATGAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgagatga                                                            >  2:681538/1‑137 (MQ=255)                                                                                                                                                                                                           aggtGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGATTGAAAGCCggg                  >  2:100536/4‑179 (MQ=255)                                                                                                                                                                                                                  gggCATGGATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaag                                                                                              <  1:606557/98‑1 (MQ=255)                                                                                                                                                                                                                        |                                                                                                                                                                                        TCAGCGAACAGTCCATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCG  >  AE015451/2414586‑2414984

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