Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,427,059	G→T	20.0%	R322M (AGG→ATG)	PP_3042 →	terminase large subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,427,059	0	G	T	20.0%	33.9 / 2.9	20	R322M (AGG→ATG)	PP_3042	terminase large subunit
Reads supporting (aligned to +/- strand):  ref base G (9/7);  new base T (2/2);  total (11/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

CTCATCGAGGAGCACCACAAGGCCGAAATGCTGGCTCAAGGCGAATGGCGCGCTCATGCGCCGGGCGATGGCGAGACGGTGAGTTTTCACCTCAATGCGCTCTATGCCCCGCTGGGATGGCAGTCATGGGTCACGCTGGCCAGGGAATACGAAGAAGCTAAGCGCGCCCAGAGCCG‑TGGCAAATTGTACCCCATGCAGGTGTTCTACAACACGCGCCTGGCCGAGGTGTGGGATTGTGCGATTGAGCAGACCAAGGCCGAAGTGCTGCAGGC  >  NC_002947/3426917‑3427188                                                                                                                                               |                                                                                                                                   cTCATCGAGGAGCACCACAAGGCCGAAATGCTGGCTCAAGGCGAATGGCGCGCTCATGCGCCGGGCGATGGCGAGACGGTGAGTTTTCACCTCAATGCGCTCTATGCCCCGCTGGGATGGCAGTCATGGGTCACGCTGGCCAGGGAATAc                                                                                                                             >  2:149726/1‑150 (MQ=255)                         gAAATGCTGGCTCAAGGCGAATGGCGCGCTCATGCGCCGGGCGATGGCGAGACGGTGAGTTTTCACCTCAATGCGCTCTATGCCCCGCTGGGATGGCAGTCATGGGTCACGCTGGCCAGGGAATACGAAGAAGCTAAGCGCGCCCAGAGc                                                                                                     >  1:260261/1‑150 (MQ=255)                               cTGGCTCAAGGCGAATGGCGCGCTCATGCGCCGGGCGATGGCGAGACGGTGAGTTTTCACCTCAATGCGCTCTATGCCCCGCTGGGATGGCAGTCATGGGTCACGCTGGCCAGGGAATACGAAGAAGCTAAGCGCGCCCAGAGCCG‑TGGc                                                                                              >  1:307463/1‑150 (MQ=255)                                  gCTCAAGGCGAATGGCGCGCTCATGCGCCGGGCGATGGCGAGACGGTGAGTTTTCACCTCAATGCGCTCTATGCCCCGCTGGGATGGCAGTCATGGGTCACGCTGGCCATGGAATACGAAGAAGCTAAGCGCGCCCAGAGCCG‑TGGCaaa                                                                                           >  2:186369/1‑150 (MQ=255)                                  gCTCAAGGCGAATGGCGCGCTCATGCGCCGGGCGATGGCGAGACGGTGAGTTTTCACCTCAATGCGCTCTATGCCCCGCTGGGATGGCAGTCATGGGTCACGCTGGCCATGGAATACGAAGAAGCTAAGCGCGCCCAGAGCCG‑TGGCaaa                                                                                           >  2:186266/1‑150 (MQ=255)                                       aGGCGAATGGAGCGCTCATGCGCCGGGCGATGGCGAGACGGTGAGTTTTCACCTCAATGAGCCCGGAGACCCGCTGGGACGGCAGTCATGGGTCACGCTGGCCAGGGAATACGAAGGAGCTAAGCGCGCCCAGAGCCG‑TGGCAAATTGTa                                                                                      <  1:144793/150‑1 (MQ=255)                                                       cATGCGCCGGGCGATGGCGAGACGGTGAGTTTTCACCTCAATGCGCTCTATGCCCCGCTGGGATGGCAGTCATGGGTCACGCTGGCCATGGAATACGAAGAAGCTAAGCGCGCCCAGAGCCG‑TGGCAAATTGTACCCCATGCAGGTGTTc                                                                      <  1:186266/150‑1 (MQ=255)                                                       cATGCGCCGGGCGATGGCGAGACGGTGAGTTTTCACCTCAATGCGCTCTATGCCCCGCTGGGATGGCAGTCATGGGTCACGCTGGCCATGGAATACGAAGAAGCTAAGCGCGCCCAGAGCCG‑TGGCAAATTGTACCCCATGCAGGTGTTc                                                                      <  1:186369/150‑1 (MQ=255)                                                         tGCGCCGGGCGATGGCGAGACGGTGAGTTTTCACCTCAATGCGCTCTATGCCCCGCTGGGATGGCAGTCATGGGTCACGCTGGCCAGGGAATACGAAGAAGCTAAGCGCGCCCAGAGCCG‑TGGCAAATTGTACCCCATGCAGGTGTTCTa                                                                    <  2:307463/150‑1 (MQ=255)                                                                 gCGATGGCGAGACGGTGAGTTTTCACCTCAATGCGCTCTATGCCCCGCTGGGATGGCAGTCATGGGTCACGCTGGCCAGGGAATACGAAGAAGCTAAGCGCGCCCAGAGCCG‑TGGCAAATTGTACCCCATGCAGGTGTTCTACAACAcgc                                                            <  1:149726/150‑1 (MQ=255)                                                                   gATGGCGAGACGGTGAGTTTTCACCTCAATGCGCTCTATGCCCCGCTGGGATGGCAGTCATGGGTCACGCTGGCCAGGGAATACGAAGAAGCTAAGCGCGCCCAGAGCCG‑TGGCAAATTGTACCCCATGCAGGTGTTCTACAACACgcgc                                                          >  1:109954/1‑150 (MQ=255)                                                                                    ttttCACCTCAATGCGCTCTATGCCCCGCTGGGATGGCAGTCATGGGTCACGCTGGCCAGGGAATACGAAAAAGCTCAG‑GGGAGAAGGGGCGTTGGGAAATGGGAGCCCAGGCAGGGGTTCTACAACACGCGCCTGGCCGAGGTGTGGGa                                         >  1:164942/1‑150 (MQ=255)                                                                                                 gcgcTCTATGCCCCGCTGGGATGGCAGTCATGGGTCACGCTGGCCAGGGAATACGAAGAAGCTAAGCGCGCCCAGAGCCG‑TGGCAAATTGTACCCCATGCAGGTGTTCTACAACACGCGCCTGGCCGAGGTGTGGTATTGTGCGATTGAg                            <  2:109954/150‑1 (MQ=255)                                                                                                     tctATGCCCCGCTGGGATGGCAGTCATGGGTCACGCTGGCCAGGGAATACGAAGAAGCTAAGCGCGCCCAGAGCCG‑TGGCAAATTGTACCCCATGCAGGTGTTCTACAACACGCGCCTGGCCGAGGTGTGGGATTGTGCGATTGAGCAGa                        >  2:78840/1‑150 (MQ=255)                                                                                                           ccccGCTGGGATGGCAGTCATGGGTCACGGTCGCCAGGGAATGCGAAGAAGCTAAGCGCGCCCAGAGCCG‑TGGCAAATTGTACCCCATGCAGGTGTTCTACAACACGCGCCTGGCCGAGGTGTGGGATTGTGCGATTGAGCAGACCaatg                  <  2:164942/150‑3 (MQ=255)                                                                                                            cccGCTGGGATGGCAGTCATGGGTCACGCTGGCCAGGGAATACGAAGAAGCTAAGCGCGCCCAGAGCCG‑TGGCAAATTGTACCCCATGCAGGTGTTCTACAACACGCGCCTGGCCGAGGTGTGGGATTGTGCGATTGAGCAGACCAAGGc                 >  2:385681/1‑150 (MQ=255)                                                                                                            cccGCTGGGATGGCAGTCATGGGTCACGCTGGCCAGGGAATACGAAGAAGCTAAGCGCGCCCAGAGCCG‑TGGCAAATTGTACCCCATGCAGGTGTTCTACAACACGCGCCTGGCCGAGGTGTGGGATTGTGCGATTGAGCAGACCAAGGc                 >  2:46923/1‑150 (MQ=255)                                                                                                                      tGGCAGTCATGGGTCACGCTGGCCAGGGAATACGAAGAAGCTAAGCGCGCCCAGAGCCG‑TGGCAAATTGTACCCCATGCAGGTGTTCTACAACACGCGCCTGGCCGAGGTGTGGGATTGTGCGATTGAGCAGACCAAGGCCGAAGtgctg       <  1:78840/150‑1 (MQ=255)                                                                                                                        gCAGTCATGGGTCACGCTGGCCAGGGAATACAAAGAAGCTAAGCGCGCCCAGAGCCG‑TGGCAAATTGTACCCCATGCAGGTGTTCTACAACACGCGCCTGGCCGAGGTGTGGGATTGTGCGATTGAGCAGACCAAGGCCGAAGTGCTGCa     >  2:330892/1‑150 (MQ=255)                                                                                                                           gTCATGGGTCACGCTGGCCAGGGAATACAAAGAAGCTAAGCGCGCCCAGAGCCG‑TGGCAAATTGTACCCCATGCAGGTGTTCTACAACACGCGCCTGGCCGAGGTGTGGGATTGTGCGATTGAGCAGACCAAGGCCGAAGTGCTGCAGgc  <  1:330892/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                   CTCATCGAGGAGCACCACAAGGCCGAAATGCTGGCTCAAGGCGAATGGCGCGCTCATGCGCCGGGCGATGGCGAGACGGTGAGTTTTCACCTCAATGCGCTCTATGCCCCGCTGGGATGGCAGTCATGGGTCACGCTGGCCAGGGAATACGAAGAAGCTAAGCGCGCCCAGAGCCG‑TGGCAAATTGTACCCCATGCAGGTGTTCTACAACACGCGCCTGGCCGAGGTGTGGGATTGTGCGATTGAGCAGACCAAGGCCGAAGTGCTGCAGGC  >  NC_002947/3426917‑3427188

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