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	90.7%	79.2 / 0.8	54	E104E (GAA→GAG)	LY180_10745	major tail tube protein
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (28/21);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

TTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACTTCGACGGCCACCGTTTCGCCGGTATCGTCACGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATT  >  CP006584/2251291‑2251871                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                               gtggccAGGGTAATCACATTTTCTTTTTTCTTGGTTTTCTCACTTCTGAGTCGGGGTTAGCCGGGCAGCCAGGCTGACCGGAGGAACATCACAGGCCGTGATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACTTc                                                                                                                                                                                                                                                                                              <  2:702330/294‑1 (MQ=1)                                           ctgcTGAATCGGGGTAAACCGGTCAGCCAGGATGGCCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCGTCTCAAGCATGTGGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACTTCGACGGCAACCGTTTCGCCGGTATCGTCACGCTGGTAAGAGcc                                                                                                                                                                                                                                                    <  1:190109/295‑1 (MQ=0)                                                                                                                      ccAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCAcc                                                                                                                                                                        >  1:150473/1‑298 (MQ=14)                                                                                                                        aGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCGGATAAACGAATCCCGGAAGCCACCgc                                                                                                                                                                      >  1:801237/1‑297 (MQ=9)                                                                                                                            ggTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAAACGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGACACCGAGCGa                                                                                                                                                                  >  1:18713/1‑298 (MQ=14)                                                                                                                             gTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGAc                                                                                                                                                                 >  1:416567/1‑298 (MQ=14)                                                                                                                             gTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGAc                                                                                                                                                                 >  1:573857/1‑298 (MQ=14)                                                                                                                                 aCGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACt                                                                                                                                                             <  1:302400/298‑1 (MQ=14)                                                                                                                                 aCGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGAAAACCGAATCCGGGAAGCCCCCGAGGGACCACt                                                                                                                                                             >  1:589442/1‑298 (MQ=9)                                                                                                                                     cGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGCCGGCCACCGTTTCGCCGGTATCGTCCCGCTGGTCAGAGCCTGCAAAACGAATCGGCACAGCCTCCACACCGGTTGCGGCGTAA‑AGCTCCCAGACACCCCACTCCGGGAACCCACcgaccgacaccgccac                                                                                                                                                         >  2:3136/1‑290 (MQ=11)                                                                                                                                      gTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAACTAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGCCGCCCACGCATGACCACTTCGACGGCCACCGTTTCGCCGGTATCGTCACGCTGGTAAGAGCCCTCAAAACGAATCGGCCCAGCATCCACACCGGTTGCGCCGCAA‑AGCACCCATATAACCGCATCCGGCAACCCACCGCGCTCCCACTcaacg                                                                                                                                                        >  2:835339/1‑294 (MQ=11)                                                                                                                                       ttCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTg                                                                                                                                                       >  1:827603/1‑298 (MQ=14)                                                                                                                                          cccTTTTCAATTCGGTTGAGGGGGTCGATTTCGACACGCTCCGTACCGTCCGTGGTCGGCCGGAAATAGGTGCGGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTGCTGACGCCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACa                                                                                                                                                    <  2:342311/297‑1 (MQ=11)                                                                                                                                                gcaatatttgtgtgggTGTAGATTTCGACCCTCTTCTTGCCGTCCCTCGTCGGCCTGAAATAGTTGCAGACCCCAGTGATTTTCGATTCGGTTTCTTCTCCCGGTGTGCCCTC‑GCCGGTGTCGATTGCTTTCTGACGTCCACGCATGACCACCTCGACGGCCAGCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAATCGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCAt                                                                                                                                              <  2:79219/284‑1 (MQ=11)                                                                                                                                                       tGTTGATGGTGTCGATTTCGACCAGCGCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGTCTACGGAGCTTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGCCGATTCCTTTCTGACGTCCACGCATGACCCCCTCGCCGGCCACCGTGTCGCCGGTATCGTCTCGCCGGTAAGAGCCTGCAACCCGAATCGGCACAACATCCACCCCGGTCGCCGCGTAA‑AGCTCCCACAGCACCGCACCCGCCAAGCCACAGAGCGCCCCACCCATCGACAgccacgctccctt                                                                                                                                       >  2:266183/1‑287 (MQ=11)                                                                                                                                                           gATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAgg                                                                                                                                   <  1:99834/298‑1 (MQ=14)                                                                                                                                                                tgtCGATTTCGACCAGCTCCTTACCATCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAg                                                                                                                              >  1:887773/1‑298 (MQ=14)                                                                                                                                                                 gtggttttCGACCCGCTCCTTCCGGTCCCCCTGCCGCCCGAAATAGGTGCAGCCCACGGTGATTGTCGACTCGGTGTCTTCTCCCTGTTTGCCCTG‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAgg                                                                                                                             <  2:312969/293‑1 (MQ=11)                                                                                                                                                                      tttccgccctcTCGTTACCGTCCATCGTCGGCTGGAAAGTGGGGCAGGCCACGGGGGTTTGTGACTCGGGGTCTTGTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTGCTGCCGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGGTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTAc                                                                                                                        <  2:958503/288‑1 (MQ=11)                                                                                                                                                                      tttCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACCGCCACCGTTTCGCCGGTCTCGTCGCGCGGGTAAGAGCCCACCAAACTAAGCGGCACAGCATCCACACCGTTTGCCGCGTAA‑ACCTCCCACATAACCGAATCCGGGAAGACACCCACCCACCCCTCCATTAACAGCCCACCGCCCTCAAGGcccaagcaccc                                                                                                                        >  2:867418/1‑289 (MQ=11)                                                                                                                                                                      tttCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTATCCTCTCCCTGTTTGCCCGC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAACCGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTAc                                                                                                                        <  1:93324/298‑1 (MQ=9)                                                                                                                                                                       ttagACCAGCGCCCTTGCGTCCATCGGCGGGCGTGAAGAGGTCCAGACGGCGTGAATTTCCTACTCGGTGTCTTCTGCCTGTTTGCCCTC‑GCGGGTGTCGCTTTCTTTCTGACGTCCACGCAGGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAAGCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTAcc                                                                                                                       <  2:787825/295‑1 (MQ=11)                                                                                                                                                                          gACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCCGCAAACCGAATCGGCACAGCACCCACCCCGGTTGCGGCGTAA‑AGCTCCCCGATAACCGACTCCGGGAAGCCCCCGAGCCACAACTCCATTGACCGCCCATCACCCTCACGGCCGCGGTCTACCggc                                                                                                                    >  2:755507/1‑297 (MQ=11)                                                                                                                                                                            ccAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCTTCCCGCTGCTAACAGCCTGCAAACCGACTCCGCACAGCATCCACACCGGTTCCTGCGACC‑AGCCCCCAGCCCACCGCATCCGGGACGCCACCGCCCGACCACCACCTTGACACCCCCTCCTCATCAAGGCAGCGGTATCCCgccgc                                                                                                                  >  2:734320/1‑294 (MQ=11)                                                                                                                                                                               gCTACTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACTTCGACGGCCCACGTTTCCCCGGTATCGTCACGCTGGTAAGAACCTGAAAAACGAATCGCCACAGCACCCACACAGGTTGCGCAGTAACAGC‑GCAACATCACCGAATCCGCGAACCCACGGCGCGACCACCCCATCGACAGCGCATCGTCATCAAGTCAGAGGCCTCCCGGTgaccg                                                                                                               >  2:526325/1‑294 (MQ=11)                                                                                                                                                                                 tCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGc                                                                                                             >  1:630746/1‑298 (MQ=14)                                                                                                                                                                                    ttACCGTCCATCGTCAGCCGGAAATAGGTGCGGACCCCGGAGATTTGCGACTCGGTGGCTGCGCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGCCGTCCACGCATGACCACCTCGACGGCCACCGTTGCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGt                                                                                                          <  2:155095/298‑1 (MQ=11)                                                                                                                                                                                       ccGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTCA‑AGCCCCACGATCACCGACTCCGGGACGCCACCGAGCGACCCCTCCATAGACAGCGCATCGTCCGCACGGCCGAGGGCACCCGGTGCCCTGccctccg                                                                                                       >  2:628531/1‑293 (MQ=11)                                                                                                                                                                                            cATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATcccc                                                                                                  <  1:930734/298‑1 (MQ=14)                                                                                                                                                                                                gTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACCAATCGGCACAGCATCCACACCGGGTGACGCGTAA‑AGCACCCAGATAACCACATCCAGGAAGCCACCGAGCCACCACCCCCTTGACGACCCAGCGGCATCAAGGCCTAGGACTACCGGGCCGCCGCCGTCCATCCCCTCAc                                                                                              >  2:952744/1‑298 (MQ=11)                                                                                                                                                                                                      cGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCAAAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGCGCGACCAATCCATTGACAGCGCATCGTCATCCAGGCCGAGGTCTACCGGTGCGCTGGCGTTCATCCCCGCACCGAGAt                                                                                        >  1:855048/1‑298 (MQ=9)                                                                                                                                                                                                        gAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCCTTGACAGCGCATCGTCATCCAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACcgcgctcg                                                                                      >  1:537717/1‑294 (MQ=9)                                                                                                                                                                                                        gAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAg                                                                                      <  1:755507/298‑1 (MQ=14)                                                                                                                                                                                                           aTAGGTGCAGACCACGGTGATTTTCGCCTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGttt                                                                                   <  1:640727/298‑1 (MQ=14)                                                                                                                                                                                                                  gcgatcgcGGGGATTTCTGACGGGGTGGCTGCTCCCGGTTTGCCCTC‑GCCGGGGTCGCTTTCTTTCGGACGTCCACGCATGGCCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCt                                                                            <  2:627158/291‑1 (MQ=11)                                                                                                                                                                                                                   aggccctgtagaTTTTCGACTCGGTTGCTTCCCGCTGTTTGACCTT‑TCCTGTGTTGAGTTGTTTCCGGCTTCCACCAGTGCCCCGCGCGACGGCCACCGTTTCGCCGGTAGCGTCGCGCGGGTAAGAGCATGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCTGAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTAGCGGTGCGCTGCCGTGCATCCCCGCACCGCGATAGTTTTCGAGCtt                                                                           <  2:429595/289‑1 (MQ=11)                                                                                                                                                                                                                            ggaTTTTCGGTGCGGTGTCGTTTGCCTGTTTGCCCTC‑GCCGGTGTGGCTTTCTTTCTGCCTTCCACCCCGGACCGCCTCGACGGCCACCGGTTCGCCGGTATCGTCCCTGTGGGAGGTGCCGGCAAAACGAATAGGCACAGCTTCCACGCCGGTTGCGGCGGAA‑AGCTCCCAGATAACCGAAGCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGGGCTTACGGGTCAg                                                                  <  2:887773/297‑1 (MQ=11)                                                                                                                                                                                                                                       tCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCg                                                       >  1:642607/1‑298 (MQ=14)                                                                                                                                                                                                                                        cGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGt                                                      <  1:320871/298‑1 (MQ=14)                                                                                                                                                                                                                                        cGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGGTCACCGGTGCGCTGCCGTTCATCCCCGCACCGCGCGAGTTTTCGAGCTTCCGGGTCAGTTTTGGCCGCGt                                                      >  2:163318/1‑298 (MQ=255)                                                                                                                                                                                                                                                 ctcCCTGTTTGCCCGC‑GCCGGTGTCTTTTTTTTGCGGCCGTCCCCGCAGGACCCCCTCGACGGCCACCTTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCACAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACtt                                             <  2:303429/298‑1 (MQ=11)                                                                                                                                                                                                                                                  ttcgcggttGCCCCCGGCCGGTTGCGATTTGTTTTCGACGCCCACGCCTGCGCACCGCGACGGCCCCCGGTTGGCCGGTGTCGTCGCGCTGGGAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTGGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGCCCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACtt                                             <  2:630746/291‑1 (MQ=11)                                                                                                                                                                                                                                                      tGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCaa                                        >  1:293897/1‑298 (MQ=14)                                                                                                                                                                                                                                                        tttGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACg                                      >  1:198883/1‑298 (MQ=14)                                                                                                                                                                                                                                                           gCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGccc                                   <  1:64288/298‑1 (MQ=14)                                                                                                                                                                                                                                                                     ggTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGCCGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAATAGCCTGCAAAACGAATCGGCCAAGCATCCACACCGAGTGCGGCGTAA‑AACCCCCAGATAACCGAATCGGGGAAGCAACCTAGCGCCCACTCCATTGACAGCCGCTCCTCCTCACGGCCGAGGTCTACCGGTGCGCCGCCGTTCCTCCACGCACCACCACAGTTTGCGAGCTAACGCGGCAGGTTTTGCCGCCTGAGGAACTTTGCaaacgccagatcacga                          >  2:993268/1‑284 (MQ=2)                                                                                                                                                                                                                                                                            aTTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGt                   <  1:513612/298‑1 (MQ=14)                                                                                                                                                                                                                                                                              ttttttcgGCCGTCCCGGCATGCCGACTACGACGGCCACCGTGTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAAGCGGCACAGTAGCCACACCGGTTGCGGCGTAA‑AGCGCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATGGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTc                 <  2:416567/295‑1 (MQ=11)                                                                                                                                                                                                                                                                              ttctttctGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTc                 >  1:418272/1‑298 (MQ=14)                                                                                                                                                                                                                                                                                ctttctGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGCCTTTGCAACGCCCTGATAGCTGTAGCCGTTCAg               >  1:346679/1‑298 (MQ=14)                                                                                                                                                                                                                                                                                 tttctGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCCGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGACTCAGGGAAGCCACCCAGCGACCCCTCCATTGACAGCGCATCGACATCAAGGCCGAGGTCCACCGGGCCGGTGACGCTCATCCCCCCCCCGCCATCGTCTTCGAGTTTACGGGTCACCTTTTGCACCGAGACGGCCTTTTCTACGCCCCGATAACTGTAGTCGTTCAgc              >  2:574246/1‑297 (MQ=11)                                                                                                                                                                                                                                                                                          tCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTCGCCGTTCAGAAAGACGTTc     >  1:509148/1‑298 (MQ=14)                                                                                                                                                                                                                                                                                            cccGCATCACCACCTGGACGCCACCCGTTTCGCCGGGACCGTCGCGCGGGAAAGAGGCGGACAAAGGAATCGGCACAGCATGCACAACTGTTGCGGCGTAA‑CGCTCCCAGATAACCGGATCCGGGAAGCCACGGGGCGACCACTCCATTGCCAGCGCAGCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCAt   <  2:698324/296‑1 (MQ=11)                                                                                                                                                                                                                                                                                             aCGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCAtt  >  1:710758/1‑298 (MQ=14)                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                               TTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACTTCGACGGCCACCGTTTCGCCGGTATCGTCACGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATT  >  CP006584/2251291‑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: ‑