Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	1,563,751	T→C	N344N (AAT→AAC)	LY180_07445 →	transposase IS609

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	1,563,751	0	T	C	97.7%	103.5 / NA	43	N344N (AAT→AAC)	LY180_07445	transposase IS609
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (28/14);  total (29/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTACACAAGCCAGCGTTGCGCGTGCTGTGGTCATACAGCGAAAGAAAATCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAGAACGC  >  CP006584/1563596‑1564042                                                                                                                                                            |                                                                                                                                                                                                                                                                                                     tCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGc                                                                                                                                                       >  1:1010414/1‑298 (MQ=9)          aaccgttgagaaatGGAGCAGGGCTGGTAGGAAATGCACCGCCAGCTTGAGTACAAGACGCTCTGGGGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCa                                                                                                                                              <  2:962705/285‑1 (MQ=11)            ccGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCCGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATCTACAGCGAACGCCGATGTAACTGGCCCTCGTCACATTTTAGCGGCGGGGCACCCCGTTCTTGCCTTGGGTG‑GAGTTGTGCACTCAAGCCCCCCGCTGAGGCAGGGACCCAcc                                                                                                                                            >  2:42262/1‑298 (MQ=11)                   aTACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTACACAAGCCAGCGTTGCGCGTGCTGTGGTCATACAGCGAAAGAAAATCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGGGGAG‑AAGTGGGGGCATCAGGCCGTCCGTTGAAACGCGAACCCACAGAAATGa                                                                                                                                     >  2:1015609/1‑298 (MQ=11)                    tACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGAt                                                                                                                                    >  1:951115/1‑298 (MQ=9)                       tGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCa                                                                                                                                 >  1:872483/1‑298 (MQ=9)                         gATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATACACAGCGAACGCCGCTGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCCCGCCCTTCTTGAC‑GTGGGGTGGATGGTGCATCCCGGCCGCCCGGTGCAACAGGACACCACCAAAATggtcaaag                                                                                                                               >  2:18113/1‑291 (MQ=11)                                   ggTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTg                                                                                                                     >  1:857755/1‑298 (MQ=9)                                     tATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGaa                                                                                                                   <  1:644891/298‑1 (MQ=9)                                               accgccAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGGTCTTGCCTGTGGAG‑AGATGGTGCCGTCAGGCCGCCCGTTGAAGCGGGGACCCACCGAAATGATTCAGGCGACAGGCCGAACGGAGCAggg                                                                                                         >  2:255825/2‑298 (MQ=11)                                                     aGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCAc                                                                                                   <  1:255825/298‑1 (MQ=17)                                                      gCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACg                                                                                                  >  1:535691/1‑298 (MQ=17)                                                       cTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGt                                                                                                 >  1:91622/1‑298 (MQ=17)                                                        ttGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTc                                                                                                >  1:688686/1‑298 (MQ=17)                                                           aGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCCCGCCGTTCTTGCCTGTGGGG‑AGATGGTGCAGTCCGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCCGGCGACCGCCTGAACGTAGCAGGGATCCCCGTCCtc                                                                                             >  2:768896/1‑297 (MQ=11)                                                                 agcagcTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGc                                                                                       >  1:339799/1‑298 (MQ=17)                                                                  gcagcTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCg                                                                                      >  1:836082/1‑298 (MQ=17)                                                                      ctctGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTgga                                                                                  <  1:888451/298‑1 (MQ=17)                                                                             gtggcgGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGACGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCCACAGGCTGAACGTAGCAGGGAACCACGTCCTTCCGGGCGTGGGGGATGTc                                                                           >  1:1021172/1‑298 (MQ=11)                                                                                  ggTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAatat                                                                      <  1:412479/298‑1 (MQ=18)                                                                                    tCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGCCCTTCAGGGCGTGGAGGATGTCAATATCg                                                                    >  1:106954/1‑298 (MQ=16)                                                                                        gTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGGGGAG‑AGGTGGTGCAGGCATGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGGCAGCCTGAACGTAGCAGGGACCCAAGCCCTTCCGGGCCTGGGGGATGGCACTACCGGttt                                                                >  2:364430/1‑298 (MQ=11)                                                                                             tgctGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTc                                                           <  1:468331/298‑1 (MQ=255)                                                                                              gctGTGCCGCCAGCGTTTGCAAGCCGGCGTTGCGCGTACTGGGGTGATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGAGGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCa                                                          <  2:865995/298‑1 (MQ=11)                                                                                                  ttCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGa                                                      <  1:136436/298‑1 (MQ=255)                                                                                                    ccGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGAtt                                                    <  1:768896/298‑1 (MQ=255)                                                                                                          gCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGAGGTAAATATCGGTTTTTGTTAGCGATTAAATTc                                              >  2:315394/1‑298 (MQ=255)                                                                                                                   aGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAgg                                                                                                          >  1:85790/1‑229 (MQ=17)                                                                                                                   aGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAgg                                                                                                          <  2:85790/229‑1 (MQ=17)                                                                                                                    gCCCGCGTTGGGCGTTATGGTGTCATACAGCGAAAGAGAACCGCCTGTCAAAAAGTAAATTCAGAGGGCAGGTATGTGGATATACAGCGAACGCCGATGTAAAGGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCGGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGt                                    <  2:431093/298‑1 (MQ=14)                                                                                                                          gTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACaaa                              >  1:468171/1‑298 (MQ=255)                                                                                                                              cgcgTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGTGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGcc                          <  1:154591/298‑1 (MQ=255)                                                                                                                                       tgGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGCG‑AGATGGTGCAGTGAGGCCGCCCGTTGAAGCAGGAACCCACCGAACTGATTCAGGCGACAGCCCGCACTCAGCAGGGAGCCACGTCCTTCAGGGCGTGGAGGATGTCAATAACGGGGTCTGGCAGCGAGTAAATTACACACCCAGTAACAAAACCCACCGGccaa                 >  2:528270/1‑295 (MQ=11)                                                                                                                                         gTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTa               >  1:40016/1‑298 (MQ=255)                                                                                                                                          tCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTaa              >  1:557865/1‑298 (MQ=255)                                                                                                                                           cATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAAt             <  1:315394/298‑1 (MQ=255)                                                                                                                                             tACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATcc           >  1:157340/1‑298 (MQ=255)                                                                                                                                                aGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTCCCGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACCTAACAACCCCCTCCCGCAATTAATccccc        >  2:60227/1‑297 (MQ=255)                                                                                                                                                 gCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGGAACAAACGCCTCCGGCGATTAATCCCCAg       >  1:122965/1‑298 (MQ=255)                                                                                                                                                   gaaagagaACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAGaa     >  1:943049/1‑298 (MQ=255)                                                                                                                                                    aaagagaACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAGAAc    <  1:382014/298‑1 (MQ=255)                                                                                                                                                     aagagaACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAGAACg   >  1:469980/1‑298 (MQ=255)                                                                                                                                                      agagaACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTGGACTGTGGCG‑AGATCGTGCCGTCCGGCCGCACGGTGAAGCAGGAACCCACCAAAATGATTAAGCAGCCAGCCTGAACGGAGCAGGGCTACACGTCCTACAAGGCATGGAGGATGTCAATAGCGCGTTCTGTCACAGATTAACTTCAACACCACGTAACAAACGCCTCCGGCAATAAATCCCAACAACGc  >  2:580539/1‑298 (MQ=14)                                                                                                                                                            |                                                                                                                                                                                                                                                                                                     TCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTACACAAGCCAGCGTTGCGCGTGCTGTGGTCATACAGCGAAAGAAAATCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAGAACGC  >  CP006584/1563596‑1564042

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