Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	1,964,727	G→T	S170S (TCC→TCA)	LY180_09345 ←	transposase IS609

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	1,964,727	0	G	T	98.2%	98.0 / NA	55	S170S (TCC→TCA)	LY180_09345	transposase IS609
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (27/27);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGTGAACCGGAGTGGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACAGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAG‑CAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATGGCTGTGAGGGCGAAT  >  CP006584/1964437‑1965019                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                      tCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACtt                                                                                                                                                                                                                                                                                                >  1:82024/1‑298 (MQ=17)    aTGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGCTACTTCAc                                                                                                                                                                                                                                                                                             >  1:153183/1‑298 (MQ=17)          ttcttttggttgTTTGTGCTGCCGGGCGTGTTTACTTTATGAATGGAGTCCCTGGGTGTATTTGCGTTACGGGAATGCGGTGTCTTTATTTTTCGTTTCTGCTTCTTCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTtctg                                                                                                                                                                                                                                                                                       <  2:380992/287‑1 (MQ=11)          aTCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTtctg                                                                                                                                                                                                                                                                                       >  1:406986/1‑298 (MQ=17)            cATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTtctgtc                                                                                                                                                                                                                                                                                     <  1:632377/298‑1 (MQ=17)                               aCGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACt                                                                                                                                                                                                                                                                  >  1:493556/1‑298 (MQ=17)                               aCGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACt                                                                                                                                                                                                                                                                  >  1:218217/1‑298 (MQ=17)                                      ttGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACACGCCCAAAGACTGTGCCATCTGAAAGCGTGGCGAGGTTAGCAACGCCAGAATCCAGCCCGACCACTGATGCTGAAGGGTGACCCGGAGTGGATAATTCACTCTTTGTCTCCATCTGGATCTACAACTTTCCGcc                                                                                                                                                                                                                                                           >  2:238668/1‑297 (MQ=11)                                            gttgtttgtggTAGTCCCTCCGGAGTTTTGGGATCCGGGTGTGCAGGCGCTGTATTTTGGGTTTTTGCTTCTGCCAGTTAGTGCTGAATGTGTCCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCATCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGCTGTACCACTTACCGCAGGACTg                                                                                                                                                                                                                                                     <  2:381656/289‑1 (MQ=11)                                               tgtGAAGGTAGGTCCTTTGGGTATTTGCGGTACGGGAATGGAGGTGCTGTGTTGTGCGTTTCTGCTTCTGCCACTTATTGCGTAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCt                                                                                                                                                                                                                                                  <  2:469195/296‑1 (MQ=11)                                                       tAGTCCCTTGTGATATTTGCGAGACGGGAATGCAGTCGTCTTTGTTTGAGTTTCTGCGTCTGCCAGTTATTGCTGACTTTGACCTTGCGGCTTAACTGTCGCTGCAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAGAACTGTTTACAGGCCCCAAGACTGTGCCATCTGACAGCGTTGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCCTTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGa                                                                                                                                                                                                                                          <  2:529028/298‑1 (MQ=11)                                                                   atatTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCacaa                                                                                                                                                                                                                              >  1:475751/1‑298 (MQ=17)                                                                       ttGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACAcc                                                                                                                                                                                                                          >  1:633572/1‑298 (MQ=17)                                                                            aTACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGa                                                                                                                                                                                                                     >  1:477332/1‑298 (MQ=17)                                                                                      gcgtacgttGTATTTCGGATTTCTGTTTCTGCCAGTTATTGCTGCATTTGGCCTTACGGCTTAGCGTTCGCTGCAGTCTCGCCTGCTGCTTCTTGTTTTTCTGCAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCt                                                                                                                                                                                                           <  2:258011/290‑1 (MQ=11)                                                                                      gCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCAAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCt                                                                                                                                                                                                           <  1:141362/298‑1 (MQ=9)                                                                                         gTCTCTGTATTTTGCGTTTCTGATTCGGCCAGTTATTGCGGGATTTGACCTTGGGGCTTAATGGTCGGTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGGTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTAAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTAtt                                                                                                                                                                                                        <  2:43871/298‑1 (MQ=11)                                                                                         gTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACCGCGTGGCGAGTTTAGCCACCCCAGCATCCAGCCCGACAATTCATCCTGAAGGGTGAACCGGAGTTGATACCTCACTTTCTGGCTGAATACCTAAGTACCACTTACCGCAGGACTGGCCTACCGTGCCATTTTTTGACACACCCGTCCCCTCACggttctt                                                                                                                                                                                                        >  2:328908/1‑293 (MQ=11)                                                                                                tATTTTGCGTTTCTGCTTCTGGCAGTTTTGGCTGATTTTGACCTGGCGGCTTAACTGTCGCTGAAGTCTCGCCATCGTCTTCGGGTTTTTCTGAAAGCTTTTTACAGGCCCCAAGACTGTGCCATCTGACAGCGTGGCGAGTTGGGCCACGCCAGCATACAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAGTACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAg                                                                                                                                                                                                 <  2:243625/298‑1 (MQ=11)                                                                                                                     ccAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGc                                                                                                                                                                            <  1:238668/298‑1 (MQ=14)                                                                                                                      cAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCGAGCCCGATCATTGATGCTGAAGGGGGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTAACACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTAACCTGACGGCTATTCCGGTATCGCCTCCAGCCCAGTTTCGGCa                                                                                                                                                                           >  1:68034/1‑298 (MQ=11)                                                                                                                        gttattGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGa                                                                                                                                                                         >  1:487036/1‑298 (MQ=14)                                                                                                                                    tttGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAACCTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGAATCCAGCCCGATCATTGATGCTGAAGGGGGAACCGGAGTTGATAATTCACTTTCTGTCTGAATACTGATTTACCACTTACCGCAGGACTGGCTGACCGGCACATTTTTGACAACACCCGTGACCTGACGGCTCTTCCGGTAGCCGAGCCAGCCCATTTTGTGCAGCACAATACGGCtt                                                                                                                                                             >  2:310163/1‑297 (MQ=11)                                                                                                                                              cGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGAt                                                                                                                                                   <  1:555863/298‑1 (MQ=14)                                                                                                                                                  ttAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAACATACGGCTGTTTTctggtccgga                                                                                                                                               >  2:214979/1‑289 (MQ=9)                                                                                                                                                  ttAACTGTCGCTAAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAg                                                                                                                                               >  1:173123/1‑298 (MQ=14)                                                                                                                                                   tAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGc                                                                                                                                              >  1:93185/1‑298 (MQ=14)                                                                                                                                                            cTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGccc                                                                                                                                     <  1:214979/298‑1 (MQ=17)                                                                                                                                                              gAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGTCCGCTATTCCGGGAGCGCATCCACCCCAGTTTCGGCAGAAAAAACCGGATGTTTTCCTGACCGAGCTTTACccccct                                                                                                                                   >  2:88745/1‑296 (MQ=11)                                                                                                                                                                  tctcGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCggg                                                                                                                               >  1:642355/1‑298 (MQ=17)                                                                                                                                                                           cGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCGGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAAt                                                                                                                      <  1:484959/298‑1 (MQ=17)                                                                                                                                                                                          cTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTccc                                                                                                       <  1:8345/298‑1 (MQ=17)                                                                                                                                                                                                acggTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGAt                                                                                                                                                   <  2:445174/245‑1 (MQ=14)                                                                                                                                                                                                aCTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGAt                                                                                                                                                   >  1:445174/1‑248 (MQ=14)                                                                                                                                                                                                                  gACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAag‑ca                                                                              <  1:326598/298‑1 (MQ=14)                                                                                                                                                                                                                   aCTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCAAAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGACAGAAAAATACGGCTGCTTTCCTGATCGAGCTTTACGCCCTGCCGGTCGCGGGTAGCATCCTCCTTTCCCCGCCTTTTGAATCGGGGAAaag‑aag                                                                             >  2:281108/1‑298 (MQ=11)                                                                                                                                                                                                                    cTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAag‑cagc                                                                            <  1:524990/298‑1 (MQ=14)                                                                                                                                                                                                                                   gCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAG‑CAGCCCGCTTCCGgaagaa                                                             <  1:364856/298‑1 (MQ=14)                                                                                                                                                                                                                                        cggattttAGCCATGGCAACGTCTGGCCCGTGCATTGATGCTGAAGGATGAACCGCAGTTGATACTTCACTTTCTTTCTGAATTATTATGTACCCCTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCACAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAG‑CAGCCCGCTTCCGGAAGAAGtttt                                                        <  2:320636/293‑1 (MQ=11)                                                                                                                                                                                                                                             tgttgcccccgcAGCATCCAGCCAGGTCATTGATGTTGAAGGGTGAACCGGAGTTTATAATTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTAACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAG‑CAGCCCGCTTCCGGAAGAAGTTTTTGTaa                                                   <  2:410190/287‑1 (MQ=11)                                                                                                                                                                                                                                                      gCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGCTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCAGCTTTTTGAATTGGGGAAAGG‑CAGCCCGCTTCCGCGAGAATGTTTTGGAAGCCCGCCCa                                          >  2:320063/1‑298 (MQ=11)                                                                                                                                                                                                                                                                aGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCAGGTAGCGCATCCAGCCCAGCTTCGGCAGAAAAACACGGCTGTTTTCCTGATCGAGCTGTACGCCCTGCGGGTAGAGGAATGCATTATTCTGTCCCCGCTTTTTGAATCGGGGAAAAG‑CAGCCCGCTTCCAGAAGAAGGTTGTGTAAGCCCGCTCAAGGTttttga                                >  2:442767/1‑296 (MQ=11)                                                                                                                                                                                                                                                                   ccGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAG‑CAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTg                             >  1:335266/1‑298 (MQ=14)                                                                                                                                                                                                                                                                            gATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGCCTGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGACTGTTTTCCTGATCGAGCTTTACGCCCTACTCGATGCGGCGTGCATCATTTTTTCCTCGTTGTTTTATTCGGGGACAAGACCACCCG‑TCTCGAACTAAGTTTTTGTCAGCCCCCTCAAGGCGTTTTAGTGTGTGttacg                    >  2:556977/1‑295 (MQ=11)                                                                                                                                                                                                                                                                              tGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAG‑CAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAg                                        >  1:200665/1‑276 (MQ=14)                                                                                                                                                                                                                                                                              tGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAG‑CAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAg                                        <  2:200665/276‑1 (MQ=14)                                                                                                                                                                                                                                                                                   aaGGGTTAACCGGAGTTGATACTTCCCTTTCTGTTTGAATACTGATGTACCCCTTACCGCAGGACTGGCTGACCGTGACCTTTTTCACAACACCCGTGACCTGACTGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAG‑CAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATGGCtg             <  2:218217/298‑1 (MQ=9)                                                                                                                                                                                                                                                                                     gggTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCAtt                                                                                                              <  2:507384/200‑1 (MQ=14)                                                                                                                                                                                                                                                                                     gggTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCAtt                                                                                                              >  1:507384/1‑200 (MQ=14)                                                                                                                                                                                                                                                                                     gggTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAG‑CAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCCCAGGGTTTTTCCGTTCCTTTTTCAATGGGtggg           >  2:3729/1‑296 (MQ=11)                                                                                                                                                                                                                                                                                        tGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCCGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGCAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAG‑CAGCCCGCTTCCGGAAGCAGTTTTTTTACGCCCGGTCAAGGTCTCTACGTGACTGTTGCAAAGGCTGTGAgg        >  2:505566/1‑298 (MQ=11)                                                                                                                                                                                                                                                                                        tGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCACTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAG‑CAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATGGCTGTGAgg        <  1:367853/298‑1 (MQ=14)                                                                                                                                                                                                                                                                                         gAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAG‑CAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATGGCTGTGAggg       <  1:470739/298‑1 (MQ=14)                                                                                                                                                                                                                                                                                              ggAGTTGATACTTCACTTTCTGTCTGAATCCTGCTGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAG‑CAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATGGCTGTGAGGGCGAAt  <  1:392480/298‑1 (MQ=17)                                                                                                                                                                                                                                                                                              ggAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAG‑CAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATGGCTGTGAGGGCGAAt  <  1:442767/298‑1 (MQ=17)                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                      TCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGTGAACCGGAGTGGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACAGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAG‑CAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATGGCTGTGAGGGCGAAT  >  CP006584/1964437‑1965019

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