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

AGGCAGACCTTAATCCGAACGTCGCTAAC‑CGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCCGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGT‑GACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCAACACCGCGCAGAGCCAT  >  CP006584/2265737‑2266158                                                                                                                                                      |                                                                                                                                                                                                                                                                                   aggcgtaCCTTAATCCGAACGTCGCTAAC‑CGCAAAAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCCGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCcg                                                                                                                               <  2:314503/292‑1 (MQ=1)         cTTAAACCGAACGACGCTAAC‑CGGAACAAAGGCGGGGCGGGCATACCGAAAAAGAATTTGTGCAGTGACGAGGCCATCGAAAAGCTGGGGCATATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCTGAGGGCTTGAGCACCGCAGCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCGGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATAca                                                                                                                       <  2:632845/298‑1 (MQ=11)             gaccGAACGT‑GCCAACACGCAAGAACGGGTGGCGACGCAAACCGAAAACGAATTCGTGACGTTACGAGGCCATCAAAAATCTGGAGCACATTTTCTTTAAGCAGTCTTGCAAATATCAGCGGCACCGGGAGCGCGGAGGGCTTGACCGCCGCATCCGCGCTATCCTGGAGTCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTTCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATAcatcat                                                                                                                   <  2:230320/296‑1 (MQ=11)              tCCGAACGTCGCTAAC‑CGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACatcatc                                                                                                                  <  1:359496/298‑1 (MQ=14)                   aCGTCGCTAAC‑CGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGCTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCCTAACCAGATTTTTCTGTCGGCCAGGAAAACGCAGGCGTATGTGTTCCGCGAAAACATCCTCCCCtt                                                                                                             >  2:760574/1‑298 (MQ=9)                                gCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGgt                                                                                                 >  1:646257/1‑298 (MQ=14)                                    caaAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGACTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGgtttgc                                                                                             >  1:615670/1‑295 (MQ=14)                                      aaGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCGGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGgttgacgt                                                                                           <  1:64851/298‑1 (MQ=14)                                       gagcgggggTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAAGCCATCGAAAAGCTGGTGCAGATTTTCTTTGAGCAGTCTTTCGAAGATCAGTTGGACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGGGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGgttgacgtt                                                                                          <  2:223017/296‑1 (MQ=9)                                        ggcgggcgTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGgttgacgttg                                                                                         >  1:802324/1‑298 (MQ=14)                                          cgggcgTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCCGATTTTTTTGTCGGCCAGTAAAACGCAAGCGTAGGTGTTCCGCCAATCAATCATCCACTTTTCCCGGCCGggtgaccttgac                                                                                       >  2:360213/1‑298 (MQ=11)                                              cgTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTtgacctga                                                                                   <  1:341356/298‑1 (MQ=14)                                                          aaaaaGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGT‑GACCCaa                                                                      >  1:458898/1‑298 (MQ=14)                                                           aaaaGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGT‑GACCCAAt                                                                     <  1:66674/298‑1 (MQ=17)                                                                 aTTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGT‑GACCCAATTGTCCt                                                               <  1:106685/298‑1 (MQ=17)                                                                 aTTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGT‑GACCCAATTGTCCt                                                               <  1:803908/298‑1 (MQ=17)                                                                  atgtttCAGTGCTGGGACCAGTGAAAAGGTTGAGCGGTTTTTCTTTGAGCAGTTTTTCGAATTTGAGTTGCACTTATCTCGCGGAGGGCTTGAGCACGGCAGCCGCGATATCCTGAAATCCGGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGT‑GACCCAATTGTCCTg                                                              <  2:595071/295‑1 (MQ=11)                                                                    ttttCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGGTGACGTTGCCCTGCCCCGG‑TACCCAATTGCCCTggg                                                            >  2:331356/1‑298 (MQ=255)                                                                       tCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGT‑GACCCAATTGTCCTGGGCaa                                                         >  1:702929/1‑298 (MQ=17)                                                                                 ggCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGT‑GACCCAATTGTCCTGGGCAATAACGGCGCa                                               >  1:231696/1‑298 (MQ=17)                                                                                  gCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGCAAAACGCAGGCGTATGTCTTCCGCGAATACATCATCGCCTTTTCCCGGCTGGGTGACGTTCACCTGACCGGT‑CACCCAACTGTCCTGGGCaactacacggcaa                                              >  2:577114/1‑287 (MQ=11)                                                                                              cTGGAGCAGTTTTTCTTGTAGCCTGCGTGTGGATATAGGTGGCACTGGTCTCGCGCAGGGCTTGAGGACCGCATCCGCGATATCCTGAAATCCGGCCAGATTGGCGCGACGGTTTATTTTTCCCGCGGGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGT‑GACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTtc                                  <  2:44893/298‑1 (MQ=11)                                                                                               tGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGT‑GACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTtct                                 >  1:634754/1‑298 (MQ=17)                                                                                               tGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGT‑GACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTtct                                 >  1:194631/1‑298 (MQ=17)                                                                                                gaatcAGAGTTTAATTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTAAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGT‑GACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTtctc                                <  1:154838/294‑1 (MQ=12)                                                                                                     aGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGCAAACCGGTCATAACCAGATTTTTCTGTCTGCCAGTAAAACGCAGGCGTATGCGTTCCGCGAATCCATCATCGCCTTTGCCCGGCCGGTTCACGTTG‑CCATACCAGTGGACCCAACTGTCCCGGGCAATAACGGCGCCAACCCGGttttctcctgccc                           >  2:339031/1‑289 (MQ=11)                                                                                                            cTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCATCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGT‑GACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTcc                    <  1:344589/298‑1 (MQ=17)                                                                                                                     gTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTTACCGGG‑TACCCCATTGTCCCGGGCAATAACGGCGCAAAAATGATTTTTCCCGGCCCCAACTCCAACAccccc           >  2:77400/1‑295 (MQ=11)                                                                                                                             aaTATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCCTAACCCGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCCTCGCCTTTGCCCGGCTTGTTGACGTTGCCCTGACCGGT‑GACCCAATTGTCCTGGGCAATCCAGGCGCCAAACTGCTTTCTCTCGGCACCAACCCCAACACCCCGCAGAcccc   >  2:576046/1‑297 (MQ=11)                                                                                                                              atatCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGCCGTTGCCCTGACCGGT‑GACCCCATTGGCCTGGGCAATAAAGGCGCAAAACTGTTTTTTTCCCGGACCAACCCCAAAAcccccccccgccaa  >  2:757888/1‑286 (MQ=11)                                                                                                                                                      |                                                                                                                                                                                                                                                                                   AGGCAGACCTTAATCCGAACGTCGCTAAC‑CGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCCGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGT‑GACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCAACACCGCGCAGAGCCAT  >  CP006584/2265737‑2266158

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