Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	2,251,586	T→C	E104E (GAA→GAG)	LY180_10745 ←	major tail tube protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	2,251,586	0	T	C	94.4%	21.1 / ‑6.6	19	E104E (GAA→GAG)	LY180_10745	major tail tube protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (6/11);  total (7/12)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACTTCGACGGCCACCGTTTCGCCGGTATCGTCACGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATT  >  CP006584/2251334‑2251871                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                              ctcagAACCGTGGTTAACCGGTCAGCCGGGCTGTCGGGATGAGAACCACCGGCCGATTTTGCGGCGGTGTTTCTCCGGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCGCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCt                                                                                                                                                                                                                                                  <  2:424619/294‑1 (MQ=1)                                                          ggtGCGGGGGTGTTGCCGCGCCTGGTGGAGGCCGTTCACCTTCGCAATCATGTTGATGGTGTGGTTGTGGACAGCCTCCTTACCGGCCACCGTCAGCCGGAAATAGGTGCAGGCCCCGGAGTTTTTCGACTCGGTGTCTGCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAAcc                                                                                                                                                                                         <  2:123423/296‑1 (MQ=11)                                                             gcggcgGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGaa                                                                                                                                                                                      <  1:12017/298‑1 (MQ=14)                                                                                          cGTTCACCTTCTCATTCGTGTTGTGGGTGTCTATTTCGCCCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCAt                                                                                                                                                         <  2:130590/298‑1 (MQ=9)                                                                                           gTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTCTCGATTTCTTTCTGCCGTCCACGCATGACCACTTCGACGGCCACCGTTTAGCGGGCATCTTCGCGCTGGTCAGACCCGTGAAACCGAATCGGCACAGCATCCACCTCGGTCGCGGCTGAAAGCCCCCATAAAACCGCATCCGACAAGCCACCCGGCGGCCagcccact                                                                                                                                                        >  2:149185/1‑291 (MQ=11)                                                                                            ttCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTg                                                                                                                                                       >  1:107760/1‑298 (MQ=14)                                                                                                                ggtggggccaaTTGCGACCACCTCCTTCCCGGCCATCGTCGGCCGGAAATGGGTGCAGACCACGGAGAGTTTCGACTCGGTGTCTTCTCCCCGTTTGCCCTCGCCGGTGTCGATTTCGTGCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGGAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAgg                                                                                                                                   <  2:123435/288‑1 (MQ=11)                                                                                                                 ttggtgcggtgttCGACCCGATCCTTGCCGTCCATCGGCAGCCGGAAATAGGTGCATACCGCGGAGATTGTCGACTCGGGGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGc                                                                                                                                  <  2:423772/287‑1 (MQ=11)                                                                                                                                      tCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGc                                                                                                             >  1:135644/1‑298 (MQ=14)                                                                                                                                         ttACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCCTGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAACCGGCACAGCATCCACACCGGTTGCGGCGCAAAGCTCCCCGATAACCGACTCCGGGACGCCACCGCGCGACCACCCCATTGACCACGCATAGTCATCAAAGCCGCGGTCTACCGGTGCGCTGccct                                                                                                          >  2:304385/1‑296 (MQ=11)                                                                                                                                                    cGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCa                                                                                               >  1:415929/1‑298 (MQ=14)                                                                                                                                                     gTCCGCCGGAAAGAGGGGCGGGCCAGGGGGTTTTTCGGTTCGGGGTTGTTTCCCTGTTTGCCCTCGCGGGTGTCGATTTCTTTCTGGCGTCCACGCATGACCGCCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCAc                                                                                              <  2:107760/298‑1 (MQ=11)                                                                                                                                                          ccGGAAATGGGTGGAGACCAGGGGGCTTTTCGCATCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGGCGATTTCTTTCTGACGGCCACGCATGACCACCTCGCCGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGa                                                                                         <  2:138356/298‑1 (MQ=11)                                                                                                                                                                             aCGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCCTCCACACCCGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCCAGCGACCACTCCAATCACACCGCACCCTCACCAAGGCCGCCGCCTCCCGGTGCGCTGCCGTTCCTCCCCCCCCCGCGATAGttttgccgcctccggg                                                                      >  2:387613/1‑286 (MQ=11)                                                                                                                                                                                gagaTTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGTCGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGACACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCa                                                                   >  1:136398/1‑298 (MQ=14)                                                                                                                                                                                          gcttggtgctttctCCCTTTTTGCCCCCGCGGTGGTCGATTTCTTTCTGACGTCCCCGCATGACCCCCTCCACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAg                                                         <  2:414014/288‑1 (MQ=11)                                                                                                                                                                                                     tctcCCTGTTTGCCCTCGCCGTTGTCGATTTCTTTCTGCCGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACt                                              <  1:16771/298‑1 (MQ=14)                                                                                                                                                                                                                                    ctttctGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAg               <  1:387613/298‑1 (MQ=14)                                                                                                                                                                                                                                                 aCGCCTGCCCACGCCGACGGCCACCGTTTCGCCGGTATCGTCGGGCGGGGAAGAGACTGCAAAACGAATCGGCACAGCATCCCAACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCAtt  <  2:415929/298‑1 (MQ=11)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                              TTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACTTCGACGGCCACCGTTTCGCCGGTATCGTCACGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATT  >  CP006584/2251334‑2251871

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