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	85.0%	27.1 / ‑1.6	20	A155A (GCC→GCA)	LY180_10860	terminase
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base A (10/7);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GGTGACTTCAAGGAAATTGACCTGCTCGGACGCCAGATT‑GAACGACTGGCACGGGTCAACCGTTACAGTCAGACCGGCAACGAGGCAGACCTTAATCCGAACGTCGCTAACCGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCCGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCAACACCGCGCAGAGCCATAACGGCGAC  >  CP006584/2265655‑2266167                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                           cgtgagggtgggggaaTTGAGCTGTTGGGACG‑CAGATTGGATGGCCGGGCATGGGTCATCCTGTACGTTCAGCCCGGCACAGGGGCAGACATTAGTCTGAACGTGGCTAACCGCAACAAAGGCGGGCGTCGGAAACCGAAAAAGAATTTTTTCAGGGACGAGGCCATCGAAAAGCTGGAGCAGATTTTGTTTGAGCGGTCTTGCGAATATCAGTTGCACTGGTATCGCGCCGGGCTTGAGCTCCGCATTCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTccc                                                                                                                                                                                                                         <  2:188629/284‑1 (MQ=1)                                                          cacccGTTCAAGACAAACCGGCAACAGGGCAGACCTTAAACCGACCGTCGCTAGCCGCAACCAAGGCCGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAATCTGGAGCAGATTTTCTTTGAGCAGTCGTTCGAATATCAGTTGCACTGGTATCGCGCCGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGcc                                                                                                                                                                 <  2:463523/295‑1 (MQ=1)                                                                                               aaTCCGAACGTCGCTAACCGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATAcatca                                                                                                                            <  1:447537/298‑1 (MQ=14)                                                                                                 tccagaCGTCGCTAACCCCAACAAAGGCGGGCGTAGCAAACCGAAAAAGAATTTTTTCAGGGACGAGGCCATCGAAAAGCTGGGGCAGGTTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACatcatc                                                                                                                          <  2:81702/293‑1 (MQ=9)                                                                                                      aCGTCGCTAACCGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCtt                                                                                                                     >  1:365883/1‑298 (MQ=14)                                                                                                           gCTAACCGAAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGccc                                                                                                                <  1:17444/298‑1 (MQ=14)                                                                                                               aCCGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCt                                                                                                            >  1:316944/1‑298 (MQ=14)                                                                                                                      caaAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGCCGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGACAACCGGTCATACCCAGATTTTCCTGTCGGCCAGTACAACGCAGGCGTATGTGTTCCGCGATTCCCTCCTCGCCCTTGCCCGGCTgggtacc                                                                                                     >  2:546980/1‑293 (MQ=11)                                                                                                                              ggcgTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTtgacct                                                                                             >  1:655464/1‑298 (MQ=14)                                                                                                                                             aaaaGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAAt                                                                              <  1:443978/298‑1 (MQ=17)                                                                                                                                              aaaGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAAtt                                                                             >  1:216005/1‑298 (MQ=17)                                                                                                                                                                   ggCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCCGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTATATTTTTCCCGCGAGGCGCTGCTGCGCACCCTGAACACCGGTCATAACCAGATTTTTCTGTCGGTCAGTAAAACGCAGCAGTAGTTGCTCCGCCAATACCTCATCTCCTTTGCACGGCTGGTTGACGTTTACCTAACCGGTGCCCCAATTCTCCTGGGcctgaaatgtcca                                                        >  2:176721/1‑286 (MQ=11)                                                                                                                                                                    gCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGGCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCaa                                                       >  1:15654/1‑298 (MQ=17)                                                                                                                                                                      cATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCaaaa                                                     <  1:175959/298‑1 (MQ=17)                                                                                                                                                                       aaagAAAAGGTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACGGGTATCGCGCCGGGATTGAGCACCGCATCCGGGCTATCCTGAAATCCCGCCAGATTGGCGCGACGTTTGATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAAc                                                    <  2:416630/295‑1 (MQ=12)                                                                                                                                                                              aGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGCGGCGCTGCTGCGCGCCCTGACAACCGGTCATAACCAGATTTTTCCGTCGGCCAGTAAAACCCCGGCGTATGTGTTCCGCGACTACCTCACCCCCTTTGCCCGCCCGGTTGACGTTGACCTAACGGGGCACCCAATCGTCCGGGGCAATACCGgccaacacctcatttt                                             >  2:304852/1‑284 (MQ=11)                                                                                                                                                                                  ggAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTtctc                                         <  1:318599/298‑1 (MQ=17)                                                                                                                                                                                                    gCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCACCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCAACAcc                       >  1:410115/1‑298 (MQ=17)                                                                                                                                                                                                       gTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCAACAACTCCAACACcgcg                    >  1:586392/1‑298 (MQ=17)                                                                                                                                                                                                                         atcttGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCAACACCGCGCAGAGCCATAACGGCGAc  <  1:527215/294‑1 (MQ=18)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                           GGTGACTTCAAGGAAATTGACCTGCTCGGACGCCAGATT‑GAACGACTGGCACGGGTCAACCGTTACAGTCAGACCGGCAACGAGGCAGACCTTAATCCGAACGTCGCTAACCGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCCGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCAACACCGCGCAGAGCCATAACGGCGAC  >  CP006584/2265655‑2266167

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