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%	100.6 / ‑6.1	54	S170S (TCC→TCA)	LY180_09345	transposase IS609
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (26/27);  total (27/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.75e-01

TCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGTGAACCGGAGTGGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACAGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATG  >  CP006584/1964437‑1965004                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                       tCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGCGTTTAGCCACGCCAGCATCCAGCCCCATCATTGATGCTGAAGGGTGAACCGGGGTTGGTACtt                                                                                                                                                                                                                                                                                 >  2:454149/1‑298 (MQ=12)  tccgatcTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTc                                                                                                                                                                                                                                                                                <  1:416893/292‑1 (MQ=12)     tGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACCCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTAATACTTCACt                                                                                                                                                                                                                                                                             >  1:3408/1‑298 (MQ=12)           tCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCCTCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGAGGCTGAAGGGGGAAACGGGGTTGATACTTCagagacgga                                                                                                                                                                                                                                                                       >  2:416893/1‑290 (MQ=11)           tCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTtctgt                                                                                                                                                                                                                                                                       >  1:66010/1‑298 (MQ=17)                gCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGaa                                                                                                                                                                                                                                                                  <  1:454149/298‑1 (MQ=17)                 cGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAAt                                                                                                                                                                                                                                                                 <  1:357806/298‑1 (MQ=17)                   tGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGACTCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATCATTCCCTTTCTGTCTGAAtcc                                                                                                                                                                                                                                                               >  2:467452/1‑296 (MQ=255)                       ttgtgcttgcGGTGGTTGTGTCTTTTTGAAGGTGTTCCCTGGGGATAGTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCGTGCTGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTCAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑At                                                                                                                                                                                                                                                          <  2:556085/289‑1 (MQ=11)                             tGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCa                                                                                                                                                                                                                                                    >  1:584648/1‑298 (MQ=17)                             tGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCa                                                                                                                                                                                                                                                    >  1:561820/1‑298 (MQ=17)                                   tCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTAcc                                                                                                                                                                                                                                              >  2:68017/1‑298 (MQ=17)                                      ttGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGACTGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTGATTGCTGAATTTGACCTTGCGGCTTACCTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTCTTCTGAACACTGTTTACAGGCCCAAAGACTGTGCACTCTGACATCGTGGCGAGTTTAGCCACGCCAGCAGCCAGCCCGAACATTGATGCCGAAGGTGGAACCGGAGTTGATACTTCACTTTCTGTCTGCAT‑CTGTATGTACCACTTTCCGCa                                                                                                                                                                                                                                           >  2:228481/1‑298 (MQ=11)                                               tATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCt                                                                                                                                                                                                                                  >  1:103675/1‑298 (MQ=17)                                                 tGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGATTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGa                                                                                                                                                                                                                                <  1:508752/298‑1 (MQ=17)                                                  gAAGGTAGTCCCTGCGGTTATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGAc                                                                                                                                                                                                                               <  2:325408/298‑1 (MQ=17)                                                            ccTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGGTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTCGCGAGTTTAGCAACGCCCGCATCCAGCCCGATCATTTATGCTGCAGGGTGACCCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑GTGTACCACTTACCGCCGGACTGGCTGACCGTGCCAttt                                                                                                                                                                                                                     >  2:206405/1‑298 (MQ=11)                                                                           gATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTg                                                                                                                                                                                                      >  2:174977/1‑298 (MQ=17)                                                                             tACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACCCCCGCGAc                                                                                                                                                                                                    >  1:100799/1‑298 (MQ=12)                                                                             tACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGAc                                                                                                                                                                                                    >  1:125349/1‑298 (MQ=17)                                                                                 ggAATGCTGTGGCTGTTTTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGa                                                                                                                                                                                                <  2:388613/298‑1 (MQ=9)                                                                                 ggAAGGCAGTGGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGa                                                                                                                                                                                                <  2:287905/298‑1 (MQ=9)                                                                                   aaTGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACg                                                                                                                                                                                              <  2:471293/298‑1 (MQ=14)                                                                                       cagtgcttGTATTTGTCGTTTCTGCTTCGGCCAGTTCTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTa                                                                                                                                                                                          <  2:261925/291‑1 (MQ=255)                                                                                             cTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTCCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCgg                                                                                                                                                                                    <  1:171956/298‑1 (MQ=9)                                                                                                               cttctgcCAGTTATTGCTGAATTCGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGt                                                                                                                                                                  <  1:275741/298‑1 (MQ=9)                                                                                                                              gCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATa                                                                                                                                                   <  1:345926/298‑1 (MQ=14)                                                                                                                               cTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGTGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATAc                                                                                                                                                  >  1:516812/1‑298 (MQ=14)                                                                                                                                tGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACg                                                                                                                                                 >  1:301746/1‑298 (MQ=14)                                                                                                                                tGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACg                                                                                                                                                 <  1:297909/298‑1 (MQ=14)                                                                                                                                      tgtcgttGCGGTTTAACTGTCGTTGAAGCCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGtt                                                                                                                                           <  2:420514/293‑1 (MQ=9)                                                                                                                                      tGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGtt                                                                                                                                           <  1:467452/298‑1 (MQ=14)                                                                                                                                       gACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCAGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGGCATTTTTCACAACACCCGTGCCCTGACGGCTCTTCCGGTAGCGCATCCAGCCCAGTTTCGGCAAAAAAATACAGGTGttt                                                                                                                                          >  2:29559/1‑298 (MQ=11)                                                                                                                                              cGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGAt                                                                                                                                   <  1:218083/298‑1 (MQ=14)                                                                                                                                                          cGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGc                                                                                                                       <  1:564331/298‑1 (MQ=17)                                                                                                                                                            cTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGccc                                                                                                                     <  2:131953/298‑1 (MQ=17)                                                                                                                                                                             tcttctGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGc                                                                                                    >  1:398269/1‑298 (MQ=17)                                                                                                                                                                              cttctGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGca                                                                                                   <  1:68017/298‑1 (MQ=17)                                                                                                                                                                                           tGAAAACTGTTTACAGGCCCAAAGATTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTcccc                                                                                      >  1:192781/1‑298 (MQ=17)                                                                                                                                                                                              aaaCTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCt                                                                                   >  1:57302/1‑298 (MQ=17)                                                                                                                                                                                                      tACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAAt                                                                           <  2:204063/298‑1 (MQ=17)                                                                                                                                                                                                         aGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCgg                                                                        >  1:189973/1‑298 (MQ=17)                                                                                                                                                                                                                aaGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAag                                                                 >  1:378027/1‑298 (MQ=14)                                                                                                                                                                                                                 aGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAagc                                                                >  1:562504/1‑298 (MQ=14)                                                                                                                                                                                                                  gACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAagca                                                               >  2:223752/1‑298 (MQ=14)                                                                                                                                                                                                                     tgtgCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGTGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGcc                                                            >  2:143167/1‑298 (MQ=14)                                                                                                                                                                                                                        gCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGc                                                         <  1:398558/298‑1 (MQ=14)                                                                                                                                                                                                                                       ggCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGTGAGCCGGAGTGGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACAGTGACATTTTTCACAACACCCGTGACCTGACGGCTGTTCCGGTAGCGCATCCAGCCCAGTTTCGGAAGAAAAATACTGCTGTTTTCATGCTCGACCTTTACCCCACACGGGTTATGGAATGCCTCATTCCGTCCCCCCCTTTTGTCATGGGGAACAACATCCTGCCTCCGGacgaaggta                                          >  2:336820/1‑295 (MQ=11)                                                                                                                                                                                                                                              ttAGCACCGCCAGCGTCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAg                                   <  2:18604/298‑1 (MQ=9)                                                                                                                                                                                                                                               tAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGc                                  >  2:507628/1‑298 (MQ=14)                                                                                                                                                                                                                                                 tccccGCCAGCCTCCAGCCCTATCATTGATGCTGAACGGGGATCCGGAGTTGATACTACCCTTTCTGTGTGAATACTG‑TTGTAACACTTACAGCATGACTGGCTGACCGTGACATGTTTCACAGCACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCTTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGccc                                <  2:3408/294‑1 (MQ=11)                                                                                                                                                                                                                                                        cAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGGACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAg                         <  2:584648/298‑1 (MQ=14)                                                                                                                                                                                                                                                                ccccctgtccttGATGCTGAAGGGTAAACCGTAGTTGTGACTTTACTTTCTGTCTGAATACTG‑ATGTCCCACTTACCCCAGGACGGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGCCGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCa                 <  2:301746/288‑1 (MQ=11)                                                                                                                                                                                                                                                                               gCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGTGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATg  <  1:143167/298‑1 (MQ=14)                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                       TCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGTGAACCGGAGTGGATACTTCACTTTCTGTCTGAATACTG‑ATGTACCACTTACCGCAGGACTGGCTGACAGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATG  >  CP006584/1964437‑1965004

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