Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	1,563,237	T→C	V173A (GTT→GCT)	LY180_07445 →	transposase IS609

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	1,563,237	0	T	C	76.0%	43.1 / 26.8	52	V173A (GTT→GCT)	LY180_07445	transposase IS609
Reads supporting (aligned to +/- strand):  ref base T (7/5);  new base C (17/21);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.38e-01

AACAGCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGACCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAA  >  CP006584/1563089‑1563528                                                                                                                                                      |                                                                                                                                                                                                                                                                                                     aaCAGCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGCCAGATGGCACAGCCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAAACTTCAGCGACAGTTAAGCCGCAAGGTCa                                                                                                                                                 >  2:411826/1‑298 (MQ=14)  aCAGCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCaa                                                                                                                                                <  1:48775/298‑1 (MQ=14)       cGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCa                                                                                                                                          >  1:513277/1‑298 (MQ=14)       cGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGATTTCAGAAAAACAAGAACACGCTGGCGAGACTTCAGCACAAGTTAAGCCACAAGGTCAA‑ATTcc                                                                                                                                          >  2:71168/1‑297 (MQ=255)         tATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGc                                                                                                                                        >  1:289370/1‑298 (MQ=14)         tATTTTTCGGCCGAAGCTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGc                                                                                                                                        <  2:27776/298‑1 (MQ=14)           tttttCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGcaa                                                                                                                                      >  1:398222/1‑298 (MQ=31)               tCTGCCGAAACTGGGCTGGATGCGCTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCaacaac                                                                                                                                  >  1:410710/1‑298 (MQ=14)                cTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACt                                                                                                                                 <  1:413858/298‑1 (MQ=31)                cTGCCGAAACTGGGCTGGATGCGCTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACt                                                                                                                                 <  1:540265/298‑1 (MQ=14)                 ggCCGAAACGGGGCTGGATGCGCTACCGGAACAGCCGTCAGGGCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTg                                                                                                                                <  2:547905/297‑1 (MQ=14)                  gCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCGGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAACGTGAAGTATACACTCCGGTTCACCCTTCAGCATCAAAGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTGCAGCGACAGTTAAGCCGCAAGGCCAACATACA‑CAACAACTgg                                                                                                                               >  2:288390/1‑298 (MQ=11)                         cTGGGCTGGATGCGCTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGgcagaagc                                                                                                                        >  1:36702/1‑298 (MQ=14)                          tGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGCCAGTTAAGCCGCAAGGCCAA‑ATTCAGCAACAACTGgcagaagca                                                                                                                       >  1:395923/1‑298 (MQ=14)                          tGGGCTGGATGCGCTACAGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAACCTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAACCCAGAAGACCCTGGCGAAACTTCAGCGACACTTAAGCCGCAAGGCCAA‑ATTCCGCAACAACTGgcagacgca                                                                                                                       >  2:567614/1‑298 (MQ=255)                            ggCTGGATGCGCTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGgcagaagcaga                                                                                                                     >  1:288563/1‑298 (MQ=14)                                  aTGCGCGCCCGGAACAGCCGTCAGGTCACGGGGGTTGTGAAAAATGTCACTGTCAGCCATTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCa                                                                                                               <  2:155586/298‑1 (MQ=9)                                        tACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATAc                                                                                                         >  1:56226/1‑298 (MQ=14)                                           cGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGc                                                                                                      <  1:124571/298‑1 (MQ=31)                                             gAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGCGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCCAAAGCAGAACCGCAAAAtcaagagc                                                                                                    >  2:198740/1‑291 (MQ=9)                                               aTAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACt                                                                                                  <  1:483192/298‑1 (MQ=31)                                                caGCCGTCAGGGCTAGGGTGTTGTGAAAAATGACACTGTCAGGGAGTCCTGCGGTAAGTGGTACATCATTATTCAGA‑CAAAAAGTGAAGTATCCACTCCCGATAACCCTGCAGCATCAATGGCCGGGATGGATGCTGGCGTGGCTAAACTCGCCACGATGTCGGATGGCACAGGCTTTGAGCCTGTAAACAGTGTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTg                                                                                                 <  2:145086/297‑1 (MQ=11)                                                      tCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGCAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATCCAGCAACAACTGGCAGAAGCAGAAACCCAACAAACACCGACTGcctccc                                                                                           >  2:506084/1‑293 (MQ=9)                                                      tCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGAc                                                                                                   <  1:38931/290‑1 (MQ=31)                                                      tCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGAc                                                                                                   >  2:38931/1‑290 (MQ=31)                                                        aGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTg                                                                                         <  1:299727/298‑1 (MQ=14)                                                            cACGGGTGTTGTGAAAAATGTCACTTTCAGCCAGTCCTGCGGTAAGTTGTAAATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCTGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACTCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATc                                                                                     <  1:70945/298‑1 (MQ=12)                                                            cACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATc                                                                                     >  1:75173/1‑298 (MQ=17)                                                                    ttgtGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTACCGCAAataa                                                                             >  2:89953/1‑297 (MQ=17)                                                                                tcgcTGTCAGACAGTGCGGCGGTAAGTGGTACATCAGTATT‑AGACCGGAAAGGGAAGGATCCGCTCCGGATCACCCTTCAGCATCCATGGTCGGGCTGGATGCTGGCGTGGCTGAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTGCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTGAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTa                                                                 <  2:191299/295‑1 (MQ=11)                                                                                tcccTGTGAGACAGTCCTGGGGTAGGGGGTACAGCAGTATTCAGA‑AAGAAAGTGAAGTATCCACGCCGGCTCACCCTTACGCAGCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTa                                                                 <  2:368713/295‑1 (MQ=11)                                                                                 cACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTAc                                                                <  1:265860/298‑1 (MQ=17)                                                                                       cacccAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑AAGAAATGGCAGTATCCACTCCGGCTCACCCTTCAGCAGCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACGCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTcac                                                          <  2:234444/295‑1 (MQ=255)                                                                                          ccAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACaaa                                                       <  1:89953/298‑1 (MQ=17)                                                                                              tCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCTCTCCGGCTCACCCTTCGGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCGGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTca                                                   <  2:207862/298‑1 (MQ=17)                                                                                                 tGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACaa                                                <  1:376765/298‑1 (MQ=17)                                                                                                     gTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACggc                                            >  1:157636/1‑296 (MQ=17)                                                                                                        aGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCCCAAAGCCACAACGGCCgc                                         >  2:389203/1‑297 (MQ=17)                                                                                                        aGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGt                                         <  2:157636/298‑1 (MQ=17)                                                                                                           ggTACATCAGTATTCAGA‑CAGAAAGTTAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTGTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACggc                                            <  2:554458/292‑3 (MQ=17)                                                                                                           ggTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACggc                                            >  1:554458/1‑290 (MQ=17)                                                                                                                       ttCAGA‑CAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGc                          <  1:389203/298‑1 (MQ=17)                                                                                                                           ga‑cagaAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATg                      >  1:184502/1‑298 (MQ=17)                                                                                                                              cagaAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCCCAAAGTCACAACGGCCGTCAGCAAAACCCACCCAACGAt                    >  2:341865/1‑298 (MQ=12)                                                                                                                              cagaAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGAt                    <  1:341865/298‑1 (MQ=17)                                                                                                                               agaAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGATTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGCCTGCATTCCTGTATCGCAAATATCCGCAGGGACTACATCCACAAAGTCACAAAAGCCGCCAGCAAAAACCaaccaaaaata                   >  2:367268/1‑288 (MQ=11)                                                                                                                               agaAAGTGAAGTATCCACGACTGCTCACCCTGCAGCATAAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATa                   <  2:16835/298‑1 (MQ=12)                                                                                                                                     tGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCAtt             <  1:126023/298‑1 (MQ=17)                                                                                                                                        aGTATCCACTCCGGCTCACCCTTCAGCATCAATGGCCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACACTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAg          <  1:429551/298‑1 (MQ=17)                                                                                                                                           aTCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGAt       <  2:75173/298‑1 (MQ=17)                                                                                                                                              cACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTg    >  1:503200/1‑298 (MQ=17)                                                                                                                                                cTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGaa  <  1:426881/298‑1 (MQ=17)                                                                                                                                                      |                                                                                                                                                                                                                                                                                                     AACAGCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGA‑CAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAA‑ATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGACCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAA  >  CP006584/1563089‑1563528

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: ‑