Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	2,265,885	C→A	A155A (GCC→GCA)	LY180_10860 →	terminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	2,265,885	0	C	A	94.7%	35.5 / NA	19	A155A (GCC→GCA)	LY180_10860	terminase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (11/7);  minor base T (0/1);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCTTAATCCGAACGTCGCTAACCGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCCGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCAACACCGCGCAGAGCCATAACGGCGACCTGTA  >  CP006584/2265742‑2266172                                                                                                                                                |                                                                                                                                                                                                                                                                                                gACCTTAATCCGAACGTCGCTAACCGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAAt                                                                                                                                       >  1:369594/1‑298 (MQ=14)         accGGACGTCGCTAACCGAAACAAAGGCGAGAGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGGAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACatcatc                                                                                                                               <  2:415790/297‑1 (MQ=9)              aCGTCGCTAACCGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCtt                                                                                                                          <  1:316378/298‑1 (MQ=14)                                  ggcgggcgTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGgttgacgttg                                                                                                      >  1:166431/1‑298 (MQ=14)                                     gggcgTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTtgacc                                                                                                   <  1:68995/298‑1 (MQ=14)                                         gTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGGCATAACCAGATTTTTCTGTCGGCCAATAAAACCCCGGCGTATGTGTTCCGCCAATCCATCATGGCCTTTGCCCCGCTGGTTGACGTtgcccttcc                                                                                               >  2:119247/1‑295 (MQ=11)                                                aCCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGAc                                                                                        >  1:258841/1‑298 (MQ=14)                                                        aGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGt                                                                                >  1:67845/1‑298 (MQ=17)                                                        aGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACCCAGGCGTATGTGTTCCGCTAATACATCATCGCCTTTGCCCGGCTGGTTTACGTTCCCCTTACCCGTGCCCCAATTGt                                                                                >  2:333163/1‑298 (MQ=255)                                                                    ggggagggggacaggAAAAGGGGGAGCAGAATTTGGTGGGGCAGGCTGTCGAATAGGAGGTGGAGAGGAATCGGGATGGGGTTGGGAACCGAGTCCGAGATATCCTGAAAGCCCGCCCGATTGGCGCGACGGTTTATGTTTCCGGCGAGGCGCTGCTGCGGGCCCTGAAAACCGGTCAGAACCAGATGTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTGGACCTGACCGGTGACCCAATGGTCCTGGGCAATaa                                                                    <  2:397643/284‑1 (MQ=2)                                                                     tGACGCAGCCATCGAAAAGTTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCATTTGCACTGGTATCGCGCAGGGCTTAAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCAAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAAc                                                                   <  1:203461/298‑1 (MQ=12)                                                                                               aGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCCGATTTTTCTGTCGGCCCGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTAACGTTGACCTGACCGGTTACCCAATTGCCCTGGGCAATAACGGCGCACAACTGATTTTTCCCgcacc                                         >  2:332284/1‑294 (MQ=255)                                                                                               aGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGCCCCCATTGTCCTGGGCACTAACGGGCCAAAACTGAATTTTTCCGGccc                                         >  2:65391/1‑296 (MQ=255)                                                                                                                      gAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCAACACCGCGCAGAGcc                  >  1:410053/1‑298 (MQ=18)                                                                                                                         tatCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCAACACCGCGCAGAGCCATa               <  1:332284/298‑1 (MQ=18)                                                                                                                         tatCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTCAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAACCGCACGCGTCTCTGTTCCGCGACTACACCATCGCCTTTGCCCCGCTGGTTGACGCTGCACTGACCCGTGCCCCAATTGTCCCGGGCAAACACCGCCCCACCCCGCATTTTCTCGCACCCaacaccaacacccccccagaccaca               >  2:189469/1‑285 (MQ=11)                                                                                                                              gTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCAACACCGCGCAGAGCCATAACGGc          >  1:253601/1‑298 (MQ=18)                                                                                                                                  cACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCAACACCGCGCAGAGCCATAACGGCGAcc      <  1:119247/298‑1 (MQ=18)                                                                                                                                      ggTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGTGCAAAACTGATTTTTCTCGGCACCAACTCCAACACCGCGCAGAGCCATAACGGCGACCTGTa  <  1:379011/298‑1 (MQ=18)                                                                                                                                                |                                                                                                                                                                                                                                                                                                GACCTTAATCCGAACGTCGCTAACCGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCCGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCAACACCGCGCAGAGCCATAACGGCGACCTGTA  >  CP006584/2265742‑2266172

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