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	96.1%	90.1 / NA	51	S170S (TCC→TCA)	LY180_09345	transposase IS609
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (27/22);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGTGAACCGGAGTGGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACAGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATGGCTGTGAGGGCGAATCT  >  CP006584/1964433‑1965021                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                       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                                                                                                                                                                                                                                                                            >  2:553907/1‑259 (MQ=17)                                                                 cTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGa                                                                                                                                                                                                                                                                            <  1:553907/259‑1 (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)                                                                                                                               tattgCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGaaaa                                                                                                                                                                       >  1:318858/3‑299 (MQ=14)                                                                                                                                        aTTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCCAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGTGTGAACCGGAGTTGATACTTCACTTACTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGCCCGTTACATTTTTACAAACACCCTGGACCTTACGGTAGTTCCGTTAGCACATCCCGCCCAGCTTTGCAAGAAAAATActcgtc                                                                                                                                                              >  2:296801/1‑294 (MQ=11)                                                                                                                                             aCCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACtt                                                                                                                                                                                                                                                                                            >  1:51857/1‑168 (MQ=14)                                                                                                                                             aCCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACtt                                                                                                                                                                                                                                                                                            <  2:51857/168‑1 (MQ=14)                                                                                                                                               cTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCCCTTACCGCAGGACTGGCTGACCGGGACATTTTTCACCACCCCCTGTACCTGCCGGCTATTCCGGTAACCACTCCCGCCCAGTTTCGCAACAAAAATACGGCTGTGTTCCt                                                                                                                                                       >  2:604678/1‑299 (MQ=11)                                                                                                                                                           cTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTa                                                                                                                                           <  1:453976/299‑1 (MQ=14)                                                                                                                                                                           gCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCg                                                                                                                           >  1:458942/1‑299 (MQ=17)                                                                                                                                                                              aGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGaa                                                                                                                        >  1:307263/1‑299 (MQ=17)                                                                                                                                                                                           ttttCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGt                                                                                                           <  1:444166/299‑1 (MQ=17)                                                                                                                                                                                             ttCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTcc                                                                                                         <  1:329916/299‑1 (MQ=17)                                                                                                                                                                                                tGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCTTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCg                                                                                                      >  1:93106/1‑299 (MQ=17)                                                                                                                                                                                                  aaaaCTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCt                                                                                                    <  1:328890/299‑1 (MQ=17)                                                                                                                                                                                                   aaaCTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCtt                                                                                                   <  1:447686/299‑1 (MQ=17)                                                                                                                                                                                                          ttACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAAt                                                                                            >  1:500233/1‑299 (MQ=17)                                                                                                                                                                                                                     aaGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAagc                                                                                 >  1:160925/1‑299 (MQ=14)                                                                                                                                                                                                                        aCTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCGTTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAagcagc                                                                              >  1:321973/1‑299 (MQ=14)                                                                                                                                                                                                                         cTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACTGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGcc                                                                             >  1:138113/1‑299 (MQ=14)                                                                                                                                                                                                                                ccgtgctaGCGTGGCGAGTTTGGCCACGCCAGCTGCCAGGCCGTTCATTGTTGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTcc                                                                      <  2:75852/292‑1 (MQ=11)                                                                                                                                                                                                                                       aGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCAgatcgg                                                               >  1:75852/1‑295 (MQ=9)                                                                                                                                                                                                                                                  gttAGCCACGCCAGCTTCCATCCCGATCATTGATGCTGAAGGGTGAACCGTAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAg                                                    <  2:189598/298‑1 (MQ=9)                                                                                                                                                                                                                                                          cGCCAGCATCCAGACGGAGTATTGATGCTGAAGGGTAAACCGGAGTTGATACTTCACTATCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGCCCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCa                                            <  1:152783/299‑1 (MQ=255)                                                                                                                                                                                                                                                                     aGCCCGATCATTGATGCTGAAGGGTGAACCTGATTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGTCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAg                                 <  1:50248/299‑1 (MQ=9)                                                                                                                                                                                                                                                                                tGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCa                      >  1:86324/1‑299 (MQ=14)                                                                                                                                                                                                                                                                                  aTGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGACAAATACGGCTGTTTTCCTGATCGAGCTGTACGCCCTGCGGGTAGCGCAATGCATCATTCTGTCCCCGCTTTTTGAATCCGGGAAAACCAGCCCGCTCCCGGAAGAAGTTTTTGTCAGCCCGCTCAAGGTCCTTCAGTTCCTGTTGCaac                    >  2:372218/1‑298 (MQ=11)                                                                                                                                                                                                                                                                                                   ggAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGG‑AAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATGGCTGTGAGGGCGAATCt  <  1:391746/299‑1 (MQ=17)                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                       ATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGT‑ATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGTGAACCGGAGTGGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACAGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATGGCTGTGAGGGCGAATCT  >  CP006584/1964433‑1965021

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