Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	1,964,559	G→A	N226N (AAC→AAT)	LY180_09345 ←	transposase IS609

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	1,964,559	0	G	A	97.9%	84.9 / NA	49	N226N (AAC→AAT)	LY180_09345	transposase IS609
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (25/22);  minor base C (1/0);  total (27/22)
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.86e-01

GCCACGCCAGAGCTGCTTGTACTCAAGCTGGCGGCGCATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGTGAACCGGAGTGGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACAGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAG  >  CP006584/1964280‑1964852                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                      gccacgccaGAGCTGCTTGTACTCAAGCTGGCGGCGCATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGc                                                                                                                                                                                                                                                                                    <  1:98499/299‑1 (MQ=14)           agCTGCTTGTACTCAAGCTGGCGGCGCATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTg                                                                                                                                                                                                                                                                          <  1:147864/299‑1 (MQ=14)            gctgctTGTACTCAAGCTGGCGGCGCATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGTGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGt                                                                                                                                                                                                                                                                         >  1:491223/1‑299 (MQ=14)            gctgctTGTACTCAAGCTGGCGGCGCATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGt                                                                                                                                                                                                                                                                         >  1:173684/1‑299 (MQ=14)            gctgctTGTACTCAAGCTGGCGGCGCATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGAGGATATTTGCGATACGGGAACGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCCGTTATTGTTTAATTTGACCTTGCGCCTTTACTGt                                                                                                                                                                                                                                                                         >  2:137649/1‑299 (MQ=9)                ctTGTACTCAAGCTGGCGGCGCATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGTGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCt                                                                                                                                                                                                                                                                     <  1:572909/299‑1 (MQ=14)                    tccTCAAGCTGGCGGCGCATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAg                                                                                                                                                                                                                                                                 <  1:136473/297‑1 (MQ=14)                        cAAGCTGGCGGCGCATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCTGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGtctc                                                                                                                                                                                                                                                             <  1:168267/299‑1 (MQ=14)                            cTGGCGGCGCATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCa                                                                                                                                                                                                                                                         >  1:397742/1‑299 (MQ=14)                             cgggggaccgtttCCTACAGGCCCTGGTCCAGTGTCGACGTGTTTAAACCGGCGTTTGCCCGGACTGTGGCGCCCGCCTGGCTGACCGGACCCGCTGCTGACTTGGACATGTTGCTGCCCTTCAAATCCGCAATGACTATCATTGCGTGGTGTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAg                                                                                                                                                                                                                                                        <  2:542772/289‑1 (MQ=11)                                     cATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGtcttct                                                                                                                                                                                                                                                <  1:156833/299‑1 (MQ=14)                                          cATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGttt                                                                                                                                                                                                                                           >  1:132118/1‑299 (MQ=14)                                                     tGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAACGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACt                                                                                                                                                                                                                                <  1:137649/299‑1 (MQ=14)                                                              aTCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGACTGCAGTCGCTGT‑ATTTTGCCTTTCTGCTTCTGCCCGTTATTGCTGAATTTGCCCTTGCTGCTTACCTGTCGCTGTAGTCTCGCCAGCCTCCTTTGGTTTTTCTTATAACCGTTTACAgt                                                                                                                                                                                                                       >  2:328890/1‑298 (MQ=11)                                                                     ggTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAg                                                                                                                                                                                                                >  1:278824/1‑299 (MQ=17)                                                                                        ccGGTCATTGCGCCCCGTCTGGCTGCCCGTCCCCGCTGCTGACGTTGACATGTTGCTGACCTTCAAATCCTCTATGACTATCATTGCGTGGTGTTTGCTGACGGTCGGTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCg                                                                                                                                                                                             <  2:632477/299‑1 (MQ=255)                                                                                                    ccccGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAg                                                                                                                                                                                 >  1:598427/1‑299 (MQ=17)                                                                                                                    atACCCGCTGTTACCTTTGACCTGTTTCTGCACTTAAAATTCTCAATGACTTTCCTTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGGCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAg                                                                                                                                                                 <  2:569173/298‑1 (MQ=11)                                                                                                                      aCCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATAAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGcc                                                                                                                                                               >  1:174626/1‑299 (MQ=17)                                                                                                                           ctggtgACTTTGACCTGTTGCTGACCTTCAAGGCCGACGTGACGCTCATTGCGTGGGTTTTGCTGACGGTCGTTGTTTCTTTATGAAGGTAGTCCCTGCGGCTATTTGCGAGACGGGAATGCAGTCGGTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAGCTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGatc                                                                                                                                                          <  2:561211/299‑3 (MQ=11)                                                                                                                              ctgACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGCATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTCACTGCCGCTGAAGTCTCGCCACCGTCTTCTGGTTGTTCTGAAACCTCTTTCCAGGCCCAACAACAGTCCCCTCTAACAGAGCGGCGAGTTTAGCCACGCCAGCTGCCCGCCCGAccgct                                                                                                                                                       >  2:329916/1‑296 (MQ=11)                                                                                                                                 aCTCTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTTGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGAt                                                                                                                                                    >  1:2493/1‑299 (MQ=17)                                                                                                                                                       cAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGCTACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTCTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCCGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTCACAGCGTGGCGATTGTAGCCACGCCAGAACCCCGCCCGACCATATATCTTGACGTGTGAAcacgggagg                                                                                                                              >  2:447686/1‑291 (MQ=11)                                                                                                                                                          gtcttccatggtgaTCATTGCGTGTGTTTGGTCGAGGGTCGTTGTGACTTTATGAAGGTAGTTCCTGCGGAGATTTGCGTTACGGGATTGCAGTCGCGGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGTCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTGTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATa                                                                                                                           <  2:173684/286‑1 (MQ=11)                                                                                                                                                                       aTCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTTAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTTGGCCATCTGACAGCGTGGCGAGTTTCGCCACGCCCACCTCCAGCCCGATCATCGATCCTCAAGGGTGAACCGGACTTAGAACTTCACATtcttt                                                                                                              >  2:284372/1‑297 (MQ=11)                                                                                                                                                                         cATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTGTGTGAAGGTAGTCCCTGCTGATAGTTGCGATACGGGAATGCAGTCGCTGT‑TTTTTGCGTATCTGCTTCTGCCAGTCGGCTCTGCATTTGACCTTGCGGCTGAACTGTCGCTGAAGTCTCGCCAGCGTCTTTTGGTTTTTCTGACAACAGTTTCCAGGCCCTAAGCCTGTGCCATCTGACACCGTGGCGCGGTATGTCACCCCACCATCCCGTCAGACCCATAATACTGAACGGGGCACAGCGGTGGACAATTCACtgactgact                                                                                                            >  2:161337/1‑291 (MQ=2)                                                                                                                                                                               gTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATAc                                                                                                      <  1:33989/299‑1 (MQ=17)                                                                                                                                                                               gTGGTTTTTGCTGACGGTCGTTGTGACTGTATGAAGGTAGTCCCTGCGGATATTTGCGATAAGGGAATGAAGTCGCTGG‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTAGCCGGCGTCTTCTGGTTGTTCTGCAAACTGTTTACAGGCCAGAAGACTGTGCCCTCTGAACGCGTGGCGGGTTTAGCCACCCCAGCCGCCACCCACAGCACTGATGCGGACTGGTGCACAGGAGTGGATGTTCCACTTtcttacgtcatac                                                                                                      >  2:33989/1‑289 (MQ=11)                                                                                                                                                                                     tttGCTGACGGTCGTTGTGACTTCATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTTCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCCAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGt                                                                                                <  1:361102/299‑1 (MQ=12)                                                                                                                                                                                           gACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACt                                                                                          <  1:240342/299‑1 (MQ=17)                                                                                                                                                                                               gTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTAcc                                                                                      <  1:89618/299‑1 (MQ=17)                                                                                                                                                                                                           ttATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGAAGGACTGGCt                                                                          >  1:313521/1‑299 (MQ=17)                                                                                                                                                                                                           gttggAAGGTGGTCTCTGGTGGTTTTTTCGGTACGGGAATGTAGT‑GTTGTGATTTTGCGTTTTTGCTTTGGCTAGTTATTGCTGAATTTGGCCTTGCGGCGTAACTGTCGCTGAAGTATCGCCAGCGTCTTCTGGTTTTGCTGAAAACTTTTTACAGGCCCAACGACTGTGCCATCTGACGGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCt                                                                          <  2:364056/295‑1 (MQ=11)                                                                                                                                                                                                                   gTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGa                                                                  >  1:225186/1‑299 (MQ=17)                                                                                                                                                                                                                     aGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGATTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTTACa                                                                >  1:272278/1‑299 (MQ=12)                                                                                                                                                                                                                          cTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGa                                                                                                   <  1:553907/259‑1 (MQ=17)                                                                                                                                                                                                                          cTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGa                                                                                                   >  2:553907/1‑259 (MQ=17)                                                                                                                                                                                                                            gCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTca                                                         >  1:163322/1‑299 (MQ=17)                                                                                                                                                                                                                                             gggAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTCACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGAGTGTTTCGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTTATCCGGAGTTGCTACTTCACTTTTTGTCTGAACACTGATGTACCACTTACGCCGGCACTGGCTGACCTTGACATTTTTAACCACAATCGTCAccatc                                        >  2:494565/1‑296 (MQ=11)                                                                                                                                                                                                                                              tgAATGCAGTTGCTGT‑TTTTTTCGGTTCGGCTTCTGCACTTTATTGATGGATTTGACCTTGTGGCTTAAGTGGCGCTGAAGTCGCGCCAGCGGCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGTCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGAc                                       <  2:598427/298‑1 (MQ=11)                                                                                                                                                                                                                                              ggAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGGTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGAATCCAGCCCGCTCATTGATGCTGAAGGGTGCACCGGAGTTGATACTTCACTTTCTGTCTGATTACTGATGTACCACTTACCGCCGGATTGGCTGCCCGTGACATATTTCACAACCACCGTGTCCTGAc                                       >  2:439868/1‑299 (MQ=11)                                                                                                                                                                                                                                                 aTGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGc                                    >  1:606812/1‑299 (MQ=14)                                                                                                                                                                                                                                                    cAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATCCTTCCCTTTCTGTCTGAATATTGATCTACCACTTCCCGCAGGACCGGCTGACCGTCACATTTTTCACAAAACCCGTTACCTGACgccgat                                 >  2:444166/1‑294 (MQ=11)                                                                                                                                                                                                                                                       tCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTcc                              <  1:284372/299‑1 (MQ=14)                                                                                                                                                                                                                                                             t‑ATTTTGCTTTTCTTCTTCGGCGAGTTATTGCTGAATTTGGCCTTGCGGCTTAACTGTCGCTGAAGTCTCGGCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGAGCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACGTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAgg                        <  2:174626/299‑2 (MQ=11)                                                                                                                                                                                                                                                                          ctgcttctgcCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCa            >  1:189598/1‑299 (MQ=14)                                                                                                                                                                                                                                                                             cttctgcCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGtt         >  1:81237/1‑299 (MQ=14)                                                                                                                                                                                                                                                                                ctgcCTGTTATTGCTTAATTTGACCTTGCGGCTTAACTGTCTCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCg      <  2:81237/299‑1 (MQ=9)                                                                                                                                                                                                                                                                                    cAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCCGCCCGATCATTTATGCGGAAGGGTGAACCGGATTTGATACTTCACTTTCTGTCTGCATACTGATGTACCACTTACCGCAGGACTCGCGGACCGTGACATTTTTCACAACCCCCGTGGCCGTCCGGCTAGTCCGGGCCTGCACCCAGCCCACTCTCgggct  >  2:358555/1‑296 (MQ=11)                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                      GCCACGCCAGAGCTGCTTGTACTCAAGCTGGCGGCGCATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGTGAACCGGAGTGGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACAGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAG  >  CP006584/1964280‑1964852

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