Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	1,563,485	A→G	T256A (ACC→GCC)	LY180_07445 →	transposase IS609

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	1,563,485	0	A	G	98.4%	123.6 / ‑5.3	63	T256A (ACC→GCC)	LY180_07445	transposase IS609
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (30/32);  total (30/33)
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.95e-01

AGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGACCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTACACAAGCCAGCGTTGCGCGTGCTGTGG  >  CP006584/1563192‑1563732                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        aGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCCCAAAGCCACAACGGCCgc                                                                                                                                                                                                                                                     >  2:389203/1‑297 (MQ=17) aGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGt                                                                                                                                                                                                                                                     <  2:157636/298‑1 (MQ=17)                ttCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGc                                                                                                                                                                                                                                      <  1:389203/298‑1 (MQ=17)                    gacagaAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATg                                                                                                                                                                                                                                  >  1:184502/1‑298 (MQ=17)                      cagaAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCCCAAAGTCACAACGGCCGTCAGCAAAACCCACCCAACGAt                                                                                                                                                                                                                                >  2:341865/1‑298 (MQ=12)                      cagaAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGAt                                                                                                                                                                                                                                <  1:341865/298‑1 (MQ=17)                       agaAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGATTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGCCTGCATTCCTGTATCGCAAATATCCGCAGGGACTACATCCACAAAGTCACAAAAGCCGCCAGCAAAAACCaaccaaaaata                                                                                                                                                                                                                               >  2:367268/1‑288 (MQ=11)                       agaAAGTGAAGTATCCACGACTGCTCACCCTGCAGCATAAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATa                                                                                                                                                                                                                               <  2:16835/298‑1 (MQ=12)                             tGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCAtt                                                                                                                                                                                                                         <  1:126023/298‑1 (MQ=17)                                aGTATCCACTCCGGCTCACCCTTCAGCATCAATGGCCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACACTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAg                                                                                                                                                                                                                      <  1:429551/298‑1 (MQ=17)                                   aTCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGAt                                                                                                                                                                                                                   <  2:75173/298‑1 (MQ=17)                                      cACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTg                                                                                                                                                                                                                >  1:503200/1‑298 (MQ=17)                                        cTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGaa                                                                                                                                                                                                              <  1:426881/298‑1 (MQ=17)                                             gctCACCCTTCAGCATCAATGGTGGGGCTGGAGGGTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCa                                                                                                                                                                                                         <  2:412149/296‑1 (MQ=17)                                             gctCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGGGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGCCTTCAGCTACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGCCAGAAGCAGAAACGCAAAACACAGCGACTGCATCCCTGTATCGCAAAGAGCCGCCGGGACCACCATCACCAAGTCAAAACGGCCGCCAGCCACAACAACGCACTGATACTAATGGAGGATCTGAAGCTCa                                                                                                                                                                                                         >  2:475121/3‑298 (MQ=11)                                                         gCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACCGCCGCCAGCAAAAACCACCCAACGCGAGTCAGTAAGGATTTTAAGGTCAGCAACACGCCaa                                                                                                                                                                                             >  2:332504/1‑298 (MQ=11)                                                              aaTGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCTGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCCCGCAATGATAGTCATTGAGGATTTGAAGGCCAGAAACCTGTCAAAGTca                                                                                                                                                                                        >  2:115804/1‑298 (MQ=12)                                                              aaTGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTca                                                                                                                                                                                        <  1:584578/298‑1 (MQ=17)                                                                      ggCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTACCGCAAAGATCCGCAGGGCCCACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGCTAGTCATTGAGGATTTGAAGGTAAGAAAGCTGGCAAAGTCAgcggcgcg                                                                                                                                                                                >  2:121532/1‑296 (MQ=11)                                                                        cTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTa                                                                                                                                                                              <  2:410710/298‑1 (MQ=17)                                                                           gATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAACCCACGCAATGATAGTCATTGAGGATTGGAAGGTCAGCAACCTGTCAAAGTCAGCAGGGGGTACgg                                                                                                                                                                           >  2:232520/1‑298 (MQ=12)                                                                                  gCGGGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTcagcca                                                                                                                                                                    <  2:288563/298‑1 (MQ=17)                                                                                         tAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCaa                                                                                                                                                                                                     >  1:583333/1‑258 (MQ=17)                                                                                         tAAACGCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCaa                                                                                                                                                                                                     <  2:583333/258‑1 (MQ=17)                                                                                            aCTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGGAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGgcgc                                                                                                                                                          >  1:252975/1‑298 (MQ=17)                                                                                             cTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACCGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCATAAGCAGAAACGCAAATTACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACCACGGCCGTCAGCCAAAACACCGCAGTGAACGTCCTAGCGGCTTTGACGGTCCGCAACACGCCAAAAGCAGACACGGGTCCGGTCAGCACGCCGGGGCCCa                                                                                                                                                         >  2:174204/1‑298 (MQ=11)                                                                                                    cGCTGTCAGAGGGCACAGGCTTTGAGCCTGTAAAAAGTTTTCAGAAAAACCGGAAGACGCTGGCGAGACTGCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCg                                                                                                                                                  <  2:289370/298‑1 (MQ=12)                                                                                                                        tttGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAAcc                                                                                                                              <  2:503200/298‑1 (MQ=17)                                                                                                                          tGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGt                                                                                                                            <  1:198740/298‑1 (MQ=17)                                                                                                                                tGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATa                                                                                                                      >  1:12663/1‑298 (MQ=17)                                                                                                                                                    aCCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAACCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGa                                                                                                  >  2:208401/1‑298 (MQ=17)                                                                                                                                                        gaagaCGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATg                                                                                              <  1:71168/298‑1 (MQ=17)                                                                                                                                                          agaCGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATgca                                                                                            <  1:508126/298‑2 (MQ=17)                                                                                                                                                                  gCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCa                                                                                                                                                                                                                                         <  2:409933/149‑1 (MQ=21)                                                                                                                                                                  gCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCa                                                                                                                                                                                                                                         >  1:409933/1‑149 (MQ=21)                                                                                                                                                                   cGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCt                                                                                   <  1:454029/298‑1 (MQ=17)                                                                                                                                                                   cGAGACTGCAGCGACAGGAGCGCCGCAAGGTCAAATTCAACCACCACTGGCAGAAGCAGAAACGCAAAAGCCAGCGACTACATTCCTGGATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTGAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCt                                                                                   <  2:370879/298‑1 (MQ=11)                                                                                                                                                                               gACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAagca                                                                       >  1:237036/1‑298 (MQ=17)                                                                                                                                                                                 cAGTAAAGCCAAAAGGTCAAATTCAGCAAAAACTGGCAAAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAagcagc                                                                     <  2:184502/298‑1 (MQ=12)                                                                                                                                                                                   gTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGctc                                                                   <  1:183425/298‑1 (MQ=17)                                                                                                                                                                                   gTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGctc                                                                   >  1:373869/1‑298 (MQ=17)                                                                                                                                                                                         ccGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCtggcgt                                                             <  1:332504/298‑1 (MQ=17)                                                                                                                                                                                             aaGGTCAAATTACGCAACAACGAGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGCCTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCtggcgtggcg                                                         <  2:513277/298‑1 (MQ=9)                                                                                                                                                                                              aGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAg                                                                                                                                                                                           <  1:239084/167‑1 (MQ=17)                                                                                                                                                                                              aGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAg                                                                                                                                                                                           >  2:239084/1‑167 (MQ=17)                                                                                                                                                                                               ggTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGt                                                       >  1:326929/1‑298 (MQ=14)                                                                                                                                                                                               ggTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATAAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGctct                                                                  <  1:49506/287‑1 (MQ=14)                                                                                                                                                                                               ggTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATAAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGctct                                                                  >  2:49506/1‑287 (MQ=14)                                                                                                                                                                                                     aTTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGtg                                                 <  1:246895/298‑1 (MQ=14)                                                                                                                                                                                                         aGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGACATGCACCGCCCGCTTGAGGCCAAGCACCCCTGGCGGGGCGGGCAGGtgcttg                                             >  2:337737/1‑298 (MQ=11)                                                                                                                                                                                                          gCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGtgcttgc                                            >  1:17105/1‑298 (MQ=14)                                                                                                                                                                                                               aacTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTc                                       >  2:263865/1‑298 (MQ=14)                                                                                                                                                                                                               aacTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTc                                       <  1:263865/298‑1 (MQ=14)                                                                                                                                                                                                                   ggCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGtgcttagatcggaaga                                   >  1:44917/1‑287 (MQ=9)                                                                                                                                                                                                                    gaggaaccagaaACGCAAAAAACAGCGCTGGCATTCCTGTATCGCAAATACCCGCAGGGACTACCTTCAAAAAGTCCAAACGGCCGTCAGCAAAAACCACGCAATGATGGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCa                                  <  2:56226/295‑1 (MQ=11)                                                                                                                                                                                                                         agcagaaACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTa                             >  1:60682/1‑298 (MQ=14)                                                                                                                                                                                                                         agcagaaACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACGGGATCAGGGCTGGTATGAAATGCACCGCCAGCTCGAGTCCAAGCAGCTCTGCCGTGGCCGTACGGTTCTTGCTGTTTCGCCAGCGTa                             >  2:33300/1‑298 (MQ=11)                                                                                                                                                                                                                           aaaaaaCGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTata                           <  1:46658/295‑3 (MQ=14)                                                                                                                                                                                                                                      aaTACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGt                <  1:475121/298‑1 (MQ=14)                                                                                                                                                                                                                                               aCTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGtac       >  1:100890/1‑296 (MQ=14)                                                                                                                                                                                                                                                  gCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACCAGCCAACGTTTCGCGTACtgt    >  1:108189/1‑298 (MQ=9)                                                                                                                                                                                                                                                   cattgccgttacgaaaaTATCCGACGGGAAGACCTTCAGAAATGCACAACGACCGTCAGCCAAACCCACGCAAGGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCCGGCAAAATCAGGTTTAAACCGTTGGATACGGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTACACAAGCCAGCGTTGCGCGTGCtgtg   <  2:258551/284‑1 (MQ=1)                                                                                                                                                                                                                                                    aTTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTgg  >  2:157191/1‑298 (MQ=14)                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        AGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGACCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTACACAAGCCAGCGTTGCGCGTGCTGTGG  >  CP006584/1563192‑1563732

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