Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	1,563,748	A→G	E343E (GAA→GAG)	LY180_07445 →	transposase IS609

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	1,563,748	0	A	G	92.7%	85.0 / ‑0.5	41	E343E (GAA→GAG)	LY180_07445	transposase IS609
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (18/20);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTACACAAGCCAGCGTTGCGCGTGCTGTGGTCATACAGCGAAAGAAAATCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAGAACG  >  CP006584/1563590‑1564041                                                                                                                                                               |                                                                                                                                                                                                                                                                                                         gaaaaaTCAGGTTTAAACCGTTCGCAACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGGGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGt                                                                                                                                                            <  2:723137/296‑1 (MQ=9) gCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGt                                                                                                                                                            >  1:328821/1‑298 (MQ=9)  cAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGtt                                                                                                                                                           >  1:624421/1‑298 (MQ=9)   aaaaTCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTACACAAGCCAGCGTTGCGCGTGCTGTGGTCATACAGCGAAAGAAAATCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTg                                                                                                                                                          >  1:693931/1‑298 (MQ=31)   aaaaTCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTg                                                                                                                                                          >  1:752932/1‑298 (MQ=9)     aaTCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTACACAAGCCAGCGTTGCGCGTGCTGTGGTCATACAGCGCAAGAAAATCGCCTGTCACAAAGTAAATTACGATGCCAGGTATGTGGATATACAGCGAACGCCGCTGTAAATGGCACTCGTAACATTCTAGCGGCGGGGCCCGCCGTC‑CTT‑CCTGTG‑GGGTGGATGTTGCAGTCCGgcctatcgatgac                                                                                                                                                        >  2:246767/1‑288 (MQ=1)       caaGGTTTAAACGGTTTGATAATGGATCGGGGCTGGAATGAAAGGTACCGCCAGCTTGGGGACAAGCAGCTCTGGCGGGGCGGTCAGGTGCTTGCTGTTCCGCCAGGGTATACAAGCCAGCGTTGGGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCAAAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGc                                                                                                                                                      <  2:810642/296‑1 (MQ=11)          ggTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAgg                                                                                                                                                   >  1:9360/1‑298 (MQ=9)                     ttCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAACTGGCGCTCGTAACATTTTAGCGGCGGGGCACCCCCTTCCTCCCCTGTGCGAG‑AG‑TGGTCCGGCCAGCCCGCCCGTTG‑AGCAGGAAccccaccaa                                                                                                                                        >  2:797576/1‑292 (MQ=11)                       cGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAAt                                                                                                                                      <  2:332138/298‑1 (MQ=9)                                           tATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGaa                                                                                                                  >  1:633976/1‑298 (MQ=9)                                            aTGAAATGCGCCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTACACAAGCCAGCGTTGCGCGTGCTGTGGTCATACAGCGAAAGAAAATCGCCTGTCACAAAGTAAATTCAGATGCCAGGTAGGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTACCATTTTAGCGGCCGGGCACGCCGGT‑CTCTCCTGTG‑CAG‑AGATGGTGCAGTCAGGCCGCCCGTTGCAGCAGTAACCCACCCACATAGTTAAGGCGACAGCCTGaaa                                                                                                                 >  2:22220/1‑297 (MQ=11)                                              gAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGt                                                                                                               <  1:322359/298‑1 (MQ=9)                                               aaaTGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTa                                                                                                              >  1:215718/1‑298 (MQ=9)                                                aaTGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAg                                                                                                             <  1:454285/298‑1 (MQ=12)                                                  tGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCa                                                                                                           <  1:750889/298‑1 (MQ=12)                                                              ttGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGCACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTc                                                                                               >  1:386153/1‑298 (MQ=12)                                                                     cagggagcccgggcgtggggGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAgg                                                                                        <  2:806946/287‑1 (MQ=255)                                                                              cggggcgggcgGTCAGGTGCTTGCTGTCCCGCCAGCGTATACCAGCCAGCGTTGCGCGTAATGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTggagg                                                                               <  2:341019/291‑1 (MQ=255)                                                                                        ggTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAatat                                                                     <  1:603203/298‑1 (MQ=18)                                                                                             ggTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGGAGCCGGGATCCACGGCCTTCAGGGCGTGGAGGATGTCAATATCGGtt                                                                >  1:4476/1‑298 (MQ=255)                                                                                                    gctGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCa                                                         >  1:609156/1‑298 (MQ=255)                                                                                                         tCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGAt                                                    <  1:725776/298‑1 (MQ=255)                                                                                                         tCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGAt                                                    <  1:770630/298‑1 (MQ=255)                                                                                                                   tttacaAGCCGGGGTGGGGGGTACTGGGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAac                                          <  2:624421/295‑1 (MQ=14)                                                                                                                    atacaAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAaca                                         <  1:31520/296‑1 (MQ=255)                                                                                                                    atacaAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAAGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAaca                                         <  1:784663/296‑1 (MQ=255)                                                                                                                    atacaAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCCTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCCGGGATCCCCGTCCTTCAGGGCGGGGAGGAGGGCAATATCGGTTTCGGTCCGCGCTTAAATTCaaaa                                         >  2:564030/3‑296 (MQ=14)                                                                                                                       caAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAAcacca                                      >  1:55187/1‑298 (MQ=255)                                                                                                                       caAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAAcacca                                      <  1:310525/298‑1 (MQ=255)                                                                                                                                 ttGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTAAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCCGGGGGGATGTCAATCACGGTTTCTGTCAGCGATTCAATTCACAACAACGTAAAAAAc                            >  2:733299/1‑298 (MQ=14)                                                                                                                                    cgggaacTGTGGTTATACAGAGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGGCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGcc                         <  2:135207/292‑1 (MQ=255)                                                                                                                                     gcgTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCt                        <  1:67545/298‑1 (MQ=255)                                                                                                                                      cgTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTc                       <  1:129148/298‑1 (MQ=255)                                                                                                                                      cgTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTc                       <  1:309016/298‑1 (MQ=255)                                                                                                                                         acTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCgg                    >  1:534467/2‑298 (MQ=255)                                                                                                                                             tgGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACCCCGTT‑CTTGCCTGTG‑GAG‑AGATGGGTCAGTCAGGCCGCCCGTTTAAGCAGGAACCCACCGAAATGCTTCAGGCCACCACCTGAACGGCACAGGGACCCCCGTCCTTCAGGGCGGGGAGGCAGGAAATATCGGTTTCTGTCACAGATTAAATCAAACAACACGTAACAAAACCCTCCGGCAAt                >  2:224936/1‑298 (MQ=14)                                                                                                                                                tCATACAGCGAAAGGGAAACGCATGTAAAAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGTGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTaa             <  2:386153/298‑1 (MQ=255)                                                                                                                                                   tACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCCGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACCGAGCAGGGATCCACGTCCTTCAGGGCGTGTGGGGTGTCAATATCGGGTTCTGTCACCGATTCAATTTAAACCCACGAAACAAACACCCCCCGCAAATaacac          >  2:64122/1‑295 (MQ=14)                                                                                                                                                     cAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGGTGGTGCAGTCAGGCCCCCCGTTGAAGCAGGACCCCCCCGAAATGATTCAGGAGACCGCCGGAACGCAGCCGGGATCCACGGCCTTCAGGGCCGGGAGGATTTCAATATCCGTTTCCGTCCGGGATTAAATTCCACACCACCAAACCAACGCCCCCCGcccataaacccc        >  2:13214/1‑287 (MQ=14)                                                                                                                                                        cGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACAACGTAACAAACGCATCCGGCAATTAATCCCCAGa     >  1:785579/1‑298 (MQ=255)                                                                                                                                                           aagagaACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAGAACg  >  1:164922/1‑298 (MQ=255)                                                                                                                                                               |                                                                                                                                                                                                                                                                                                         GCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTACACAAGCCAGCGTTGCGCGTGCTGTGGTCATACAGCGAAAGAAAATCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAGAACG  >  CP006584/1563590‑1564041

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