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	95.8%	33.6 / ‑3.9	24	A155A (GCC→GCA)	LY180_10860	terminase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (14/9);  total (14/10)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

GGCGGTGACTTCAAGGAAATTGACCTGCTCGGACGCCAGATTGAACGACTGGCACGGGTCAACCGTTACAGTCAGACCGGCAACGAGGCAGACCTTAATCCGAACGTCGCTAACCGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCCGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAAT‑ACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCAACACCGCGCAGAGCCATA  >  CP006584/2265652‑2266159                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                    gggggTAACATCAAGGAAAATCACACGCGCGGGCGCCAGGTTGAACGCCTGGAACGGCGGAACCTTTAGCGTCAGACCGGCAACGAGGGAGCCCTGAATCCGAACGTCGCTAACCGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGGGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTGGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTc                                                                                                                                                                                                                    <  2:607513/296‑1 (MQ=1)                                             cagccgggcgctggTCAACCGGTACGCTCAGAGCGGCAACGGGGGTGGCCGTAGTCCGAACGGCGCTATCGGGAACAAGGGCGGGCGTCGCAAACCGAAAAAGAATGGTTTCAGTGACTAGGCCATCGAAAAGCGGGGGCAGATTTTCGTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCCGGGCTTGAGCACCGCATCCGCGATATCCTGAGATCCCGCCAGATTGGCGCGACGTTTTATTTTGCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGa                                                                                                                                                                        <  2:129151/287‑1 (MQ=1)                                                                                                 aaTCCGAACGTCGCTAACCGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGCTTTTTCTGTCTGCCAGTAAAACGCAGGCGTATGCGGTCCGCGAAT‑Acatcat                                                                                                                   >  1:480390/1‑299 (MQ=9)                                                                                                    ccGAACGTCGCTAACCGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAAT‑ACATCATCGc                                                                                                                >  1:133456/1‑299 (MQ=14)                                                                                                        aCGTCGCTAACCGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAAT‑ACATCATCGCCttt                                                                                                            >  1:244428/1‑299 (MQ=14)                                                                                                             cgtgaCCGCAAAAAGGGCGGGCGTGGGAAACCGTAAAAGGATTTTTTGTGGGACGGGGCCATCGAAAAGCGGGAGGAGATTTTTTTTGGGCAGTCTTTCGAATATCTGTTGCACTGGTATCGCGCAGGGCTTGAGCACGGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAAT‑ACATCATCGCCTTTGCCCg                                                                                                       <  2:8602/295‑1 (MQ=11)                                                                                                               tAACCGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTAGCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCCCGACGTCTTATGTGTCCCGCTAGGCGCTGCTGCGCGCCCTTAACACCGGTCATAACCAGATTTTTCTGTCGCCCAGTACAACGGAGACGTACGTGTTTCGCGAATCACA‑CATTCCCTTTGCCCGGc                                                                                                     >  2:186580/1‑299 (MQ=11)                                                                                                                    gCAAAAAAGGCGGGCGTCGCAAACCGATAACGAATTTTTTGAGTGACGAGGCCATTGAAAAGCGGGAGCAGATTTTTGTTGAGGCGGCTTTCGAATATCGGTTTCACTGGTATCTGGCAGTGCTTGAGCAACGCATCCGCGTTATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAAT‑ACATCATCGCCTTTGCCCGGCTGgtt                                                                                                <  2:408837/299‑1 (MQ=11)                                                                                                                           aGGCGGGCTTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGTGGCGATGCTGCGCGCCCTGAAAACCGGTAATAACCCGATTTTTCTGTCGGCCAGTAAAACGCAAGCGTATGTGTTCCGCGAAT‑ACATGATCGCCTTTGCCCGGCTGgttgacattg                                                                                         >  1:113615/1‑299 (MQ=255)                                                                                                                            ggcgggcgTCGCAAACCGAAAAAGAATTTTTACAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTAGCACTGGTATCGCGCAGGGCTTGAGCACCGCCTCCGCGATGTCCTGAAATACCGCACGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCTCGACCTGAAAACCGGTCATTACCAGCTTTTTCTGTCGGCCCGTCAAACTACGACGTATGTTTTCTGCGAAG‑ACATCATCACCTTTTCTCGGCTAgttgatggtga                                                                                        >  2:126527/1‑299 (MQ=11)                                                                                                                            ggcggactTCGCAAAACCAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGGGCAGATTTTCTTTGAGCAGTCTTTCGAATATCTGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAAT‑ACATCATCGCCTTTGCCCGGCTGgttgacgttga                                                                                        <  1:34214/299‑1 (MQ=9)                                                                                                                               gggcgTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTCTTATTTTTCCCGCCAGGACCTGCTGCCCGCCCTGAAAACCGGTACTAACCCGATTTTTCTGTTCGCCAGTCAAACCCAGGCGTACGTGTTCCGCCAAT‑ACCACCTCGCCTTTGCCCGCCTTgttcccgttggctt                                                                                     >  2:412947/1‑295 (MQ=11)                                                                                                                                         aaCCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGCAT‑ACATCCTCGCCTTTGCCCGGCTGGTTGACGTTTAACTGACCGGTGAc                                                                           >  1:186955/1‑299 (MQ=9)                                                                                                                                          aCCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAAT‑ACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGAcc                                                                          >  1:22938/1‑299 (MQ=14)                                                                                                                                              aaaaaGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAAT‑ACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAAt                                                                      >  1:209613/1‑299 (MQ=17)                                                                                                                                               aaaaGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTAGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAAT‑ACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAAtt                                                                     <  1:231275/299‑1 (MQ=17)                                                                                                                                                              gggACGAGGCTAGCGTAGAGCTGGGGGACATTGCCGTTGAGGAGCTTTTCGAATATCAGTTGCACCGGTATGGCCAAGGGCTTGAGCACGCCGTCCGCGATATCCTGAAAGCCCGCCAGATTGGCGCGACGTGTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAGCCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAAT‑ACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAAc                                                      <  2:313467/297‑1 (MQ=11)                                                                                                                                                                 aCGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAAT‑ACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGgc                                                   <  1:38515/299‑1 (MQ=17)                                                                                                                                                                    aGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGCAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGACGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCACTAACACGCAGCCCTAGTTGCTCCGCGAAT‑CCATCATCCCCTTTGCCCGCCTGGTTGCCGTTGACCTGACCGGCGCCCCAATCTTCCTGGCCACTCACGGCCCa                                                >  2:647045/1‑299 (MQ=11)                                                                                                                                                                                                 tgtGAGCAGTATTTCGAATATCAGGTGCACTGGTATCGCGCATGGCTTGAGCACCTCATCCGCGATACCCTGCAATCCCGCCAGATTGTCGCGACGTTTTATTTTTCCCGCGGGGCGCTGCTGCGCGCCCTGAAAACCGTTCATAACCAGATTTTTCTGTCGGCCAGTCAAACGCAGGCGTATGTGTTCCGCGAAT‑ACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCaa                   <  1:268386/297‑1 (MQ=11)                                                                                                                                                                                                   tGAGCAGTCTTTCGAATATCAGTTGCACTGGTCTCGCGCAGGGCTTGAGCACCGCATACGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAACCGCAGGCGTATGTGTTCCGCGAAT‑ACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCAaca                 <  1:131670/299‑1 (MQ=17)                                                                                                                                                                                                        aGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCCAAGCGCTGCTGCGCGCCCGGAAAACCGGTCATAACCAGACTCTTCTGTCGGCCAGTCAAACGCACGCGTATGTGTTCCGCGAAT‑ACACCATCGCCCTTGCCCGCCTGGTCTACGATCACCTGACCATTGACCCACTTGTCCCCGGCAGTCCCAGCGCGACACTGACTTTACCCGTCTCCACCACCCACACctcg            >  2:441581/1‑299 (MQ=11)                                                                                                                                                                                                                  atatCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAAT‑ACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCTCCAAATCCAACACCGCGCAGAGCCATa  >  1:362064/1‑299 (MQ=18)                                                                                                                                                                                                                  atatCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAAT‑ACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCAACACCGCGCAGAGCCATa  >  1:479624/1‑299 (MQ=18)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                    GGCGGTGACTTCAAGGAAATTGACCTGCTCGGACGCCAGATTGAACGACTGGCACGGGTCAACCGTTACAGTCAGACCGGCAACGAGGCAGACCTTAATCCGAACGTCGCTAACCGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCCGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAAT‑ACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCAACACCGCGCAGAGCCATA  >  CP006584/2265652‑2266159

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