Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	2,251,616	A→G	R94R (CGT→CGC)	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,616	0	A	G	94.6%	89.6 / 0.0	60	R94R (CGT→CGC)	LY180_10745	major tail tube protein
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (27/26);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

CTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACTTCGACGGCCACCGTTTCGCCGGTATCGTCACGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATTAACTTGAGTTTGCGCGGCATTGCCATCGGTCAG  >  CP006584/2251333‑2251904                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                  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)                                                                                                                                                                                                                                                                        cACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATTAACTTGAGTTTGCGCGGCAtt              <  1:860996/298‑1 (MQ=14)                                                                                                                                                                                                                                                                         aCCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATTAACTTGAGTTTGCGCGGCATTg             <  1:880846/298‑1 (MQ=14)                                                                                                                                                                                                                                                                           cGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATTAACTTGAGTTTGCGCGGCATTGcc           <  1:989683/298‑1 (MQ=14)                                                                                                                                                                                                                                                                            gTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATTAACTTGAGTTTGCGCGGCATTGCCa          >  1:144188/1‑298 (MQ=14)                                                                                                                                                                                                                                                                             tttCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATTAACTTGAGTTTGCGCGGCATTGCCAt         >  1:716124/1‑298 (MQ=14)                                                                                                                                                                                                                                                                                  ccGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATTAACTTGAGTTTGCGCGGCATTGCCATCGGTc    <  1:879965/298‑1 (MQ=14)                                                                                                                                                                                                                                                                                    ggTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATTAACTTGAGTTTGCGCGGCATTGCCATCGGTCAg  <  1:526647/298‑1 (MQ=14)                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                  CTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTC‑GCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACTTCGACGGCCACCGTTTCGCCGGTATCGTCACGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAA‑AGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATTAACTTGAGTTTGCGCGGCATTGCCATCGGTCAG  >  CP006584/2251333‑2251904

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