Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	1,645,397	C→T	intergenic (+1191/+69)	LY180_07765 → / ← LY180_07770	transposase/sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	1,645,397	0	C	T	91.9%	176.0 / ‑2.1	64	intergenic (+1191/+69)	LY180_07765/LY180_07770	transposase/sulfatase
Reads supporting (aligned to +/- strand):  ref base C (5/0);  new base T (31/26);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 6.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GGTATACAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTTGTCTGATTATTGATTTTTCGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGAG‑TTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTT  >  CP006584/1645162‑1645688                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                       ggTATACAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGGGAGTCGTGTAAAATATCGAGTTCGCATATCTTGTTGTCTGATTATCGATTTTTCGCTAAACCATTTTATCATATGAAAAGATGTG‑TATCTAc                                                                                                                                                                                                                                         >  2:62620/1‑298 (MQ=9)   tataCAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCCAATCTTGTTGTCTGATTATTGATTTTTCGCGAAACCATTTGATCATATTACAAGATTTG‑TATCTCCct                                                                                                                                                                                                                                       >  2:976811/1‑298 (MQ=2)     taCAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTTGTCTGATTATTCATTTTTAGCGAAACCATTTGATCATATTACAAGAtgattttccatctta                                                                                                                                                                                                                                      >  2:265314/1‑285 (MQ=0)         gCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGTTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGTCAGGTTACTGCCTAGTTCGGTAATAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAACCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGACAAAGAGACTCGTGTAAAATATTGAGTTCGACCCTATTGTTGTCTGAATATTTATTTTTCTCCAAACCCTTTACTCATATGCCAAGATGTT‑Tagcctacttacctt                                                                                                                                                                                                                                 >  2:773443/1‑285 (MQ=1)         gCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTTGTCTGCTTATTGATTTTTCGCGAAACGCTTTGATCATAGGACAAGATGTT‑TATCTACcttaactt                                                                                                                                                                                                                                 >  2:3519/1‑298 (MQ=9)                gCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTTGTTTGATTATTTATTTTTCGCGAAACCATTTGATTATATGACAAGATGTGTTACCTA‑CTTAAATTAattatt                                                                                                                                                                                                                          >  2:212881/1‑298 (MQ=0)                                                          gTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTCACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATaaa                                                                                                                                                                                >  1:505868/1‑298 (MQ=25)                                                                 gttatgttTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTgg                                                                                                                                                                         >  1:57620/1‑298 (MQ=255)                                                                  ttattgtttgtttgCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTGTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGt                                                                                                                                                                        <  2:730896/292‑1 (MQ=255)                                                                    atgttTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCa                                                                                                                                                                      <  1:916389/298‑1 (MQ=255)                                                                                cgttttGGTTGGAAGGTTACGGCCAAGTTCGGTAAAGGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCa                                                                                                                                                          <  2:70849/295‑1 (MQ=255)                                                                                        ttggCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATcaac                                                                                                                                                  <  1:413444/298‑1 (MQ=255)                                                                                         tggCGGGGTAGGGCCAAGTTAGGTAAGGGTGAGAGTTTTACAGTGAAGTAAGGCGTGGCAAGCCAACGTTAAGTTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATcaaca                                                                                                                                                 <  2:457987/298‑1 (MQ=17)                                                                                         tggCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTACGTGTAATTCGGGGCAGCATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTAGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATATACTTTACTTTAATGATTTTGATAAAAATCATTAGTGGCTTCATCAGGGTAAAGCATAACAGGCTGGACACCTGtttggcaaagatcacaa                                                                                                                                                 >  2:90235/1‑284 (MQ=11)                                                                                          ggCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATcaacaa                                                                                                                                                >  1:723043/1‑298 (MQ=255)                                                                                                   cGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTg                                                                                                                                       >  2:110454/1‑298 (MQ=255)                                                                                                    ggCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGc                                                                                                                                      <  1:17312/298‑1 (MQ=255)                                                                                                       cAAGTTGGGTAAGAGTGAGAGTTTTACGGTCAAGTAAGGCGTGGCAAGACAAAGTTAAGCTGTTGAGGTGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTGTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTc                                                                                                                                   <  2:671217/298‑1 (MQ=17)                                                                                                        aaGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCa                                                                                                                                  >  1:638754/1‑298 (MQ=255)                                                                                                        aaGATCGGTAAGAGTGAGAGGTTTACAGTCAAGTAAGGCGTGGCAAGCAAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGTGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCa                                                                                                                                  <  2:505181/298‑1 (MQ=255)                                                                                                              ggTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGTATTCACAAGGGTAAAGCCTAAAAGGCTGGTAAAGTGTTTTGGCACGATCAACAACGCCAATACCCTCATCGGtt                                                                                                                            >  2:464699/1‑298 (MQ=17)                                                                                                                         gagTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCt                                                                                                                 >  1:997062/1‑298 (MQ=255)                                                                                                                            ttttACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGat                                                                                                              <  1:31955/298‑1 (MQ=255)                                                                                                                                  aGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTAGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTAGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATCACAATTCTTAGGGGATTCAACAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATACACGAAGCAAATTCCCTTCTCGGTTAACTAATCGCAtatatacac                                                                                                        >  2:256235/1‑295 (MQ=14)                                                                                                                                      aaGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATca                                                                                                                                                    <  2:507785/250‑1 (MQ=255)                                                                                                                                      aaGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATca                                                                                                                                                    >  1:507785/1‑250 (MQ=255)                                                                                                                                      aaGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATca                                                                                                                                                    <  2:303529/250‑1 (MQ=255)                                                                                                                                      aaGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATca                                                                                                                                                    >  1:303529/1‑250 (MQ=255)                                                                                                                                       aGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCTAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTg                                                                                                   <  1:228911/298‑1 (MQ=255)                                                                                                                                        gtaagggtgggCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGc                                                                                                  <  2:141060/289‑1 (MQ=255)                                                                                                                                         tAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTATCTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCg                                                                                                 <  2:360936/298‑1 (MQ=255)                                                                                                                                          aaggCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTATGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGATCGATCAACAACGCCAATTGCCCCATCGGTTAACTGATCGCAGATATAACCttctttgt                                                                                                >  2:87326/4‑292 (MQ=18)                                                                                                                                                  gCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGt                                                                                        <  2:502117/298‑1 (MQ=255)                                                                                                                                                   cAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGt‑t                                                                                      <  1:787805/298‑1 (MQ=255)                                                                                                                                                     aGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATCGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCa                                                                                    >  1:860553/1‑298 (MQ=255)                                                                                                                                                       ccAACGTTAAGCTGTTGAGACGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACg                                                                                  >  1:434058/1‑298 (MQ=255)                                                                                                                                                               aaGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAACTATCGAGTTCGCACATGTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGATCAAGTGTTTTGGCCCGATCAACAACGCCAATTGACACATCGGTTAACTGAACGCGTATACTACCTTTTTCGGGTCCACGTCTCA‑GATTTTTGa                                                                          >  2:533205/1‑298 (MQ=18)                                                                                                                                                               aaGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAACGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCCCTGATATAACCTTTTGCGGGGACCCGT‑TCACGTTTTTTGa                                                                          >  2:894950/1‑298 (MQ=255)                                                                                                                                                                  cTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTACAGTATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACa                                                                       >  1:387415/1‑298 (MQ=255)                                                                                                                                                                      tGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGa                                                                   <  1:44043/298‑1 (MQ=255)                                                                                                                                                                        atggttgttAAGTTGAAGTTGGGGCAGAATGGGTAAAGGGAGTTGTGTAAAATGTCGAGTTCGATCATTTTGTGTTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TAGCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAg                                                                 <  2:562044/291‑1 (MQ=25)                                                                                                                                                                            gTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGa                                                             >  1:223818/1‑298 (MQ=255)                                                                                                                                                                               ttAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGAG‑tt                                                         >  2:88087/1‑298 (MQ=255)                                                                                                                                                                                 aaGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGAG‑TTGt                                                       >  2:364807/1‑298 (MQ=255)                                                                                                                                                                                   gtgtAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGAG‑TTGTaa                                                     >  1:960149/1‑298 (MQ=255)                                                                                                                                                                                         ttCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGAG‑TTGTAATATgtc                                               >  2:124860/1‑296 (MQ=255)                                                                                                                                                                                               gtaGAATTGGTAAAGAGAGTGTTTTAAAATATCTACTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTGAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGAG‑TTGTAATATGCCGTTCCt                                         <  1:306837/296‑1 (MQ=255)                                                                                                                                                                                                  gAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGTTATAACTCTTTGCGGGGTACCTT‑TCACGCTTTTTGAACTGTG‑AGGGGAGTTTTTAATATgtatgtcgttaa                                      >  2:669420/1‑287 (MQ=21)                                                                                                                                                                                                   aaTTGGTAAAGAGATTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGAG‑TTGTAATATGCCGTTCCTGCAg                                     <  2:456805/298‑1 (MQ=255)                                                                                                                                                                                                     ttGGTAAAGTGGGTAGTGTAAATATTTGCGTGCGTTCATTTTGTTTGCTAATTGTTGCATTTTGGCGAAACCAGTTGTTCATATGACAAGATGTG‑TATCTAACTTAACTTAATGATGTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGCTATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGAG‑TTGTAATATGCCGTTCCTGCAGCg                                   <  2:761654/298‑1 (MQ=255)                                                                                                                                                                                                      tgttAAAGAGGGTCTTGTAAAATATCGAGTTCGAACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATTTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGAG‑TTGTAATATGCCGTTCCTGCAGCGt                                  <  2:960149/295‑1 (MQ=255)                                                                                                                                                                                                        gTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGAG‑TTGTAATATGCCGTTCCTGCAGCGTgg                                <  1:384406/298‑1 (MQ=255)                                                                                                                                                                                                               gggtgttgtAAAGTATCGAGTTCGATAATTGTGTTGTCTGATTATTGATTTTTGGCGAAACCTTTTGATCATATGACAAGATGTG‑TATCTAACTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGTTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGAG‑TTGTAATATGCCGTTCCTGCAGCGTGGAATCCCa                         <  2:111898/292‑1 (MQ=255)                                                                                                                                                                                                               gagTAGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATGATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTAGCTTAACTTAATGATTTTGATAAAAATAATTAGGGGATTCATCAGGGTAAAGCATAAACGGCTGGTCAAGTGTTTTGGAACGATCAACAACGCCAATTGCCTCATCGGTTAACTGGTCGCTGATATAACCTATTGCGGGGACAAGT‑TGACGATTTGTGATCAGTGAAGGGGAG‑TTGTAATGTGCCGTTCCTGCAGCGGGGAATCCCa                         >  1:95161/1‑298 (MQ=25)                                                                                                                                                                                                                 gTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGAG‑TTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATa                       >  1:498441/1‑298 (MQ=255)                                                                                                                                                                                                                  tCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGAG‑TTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATaa                      >  1:709015/1‑298 (MQ=255)                                                                                                                                                                                                                     tgtAAAATTGCGAGTTCGCACATTTTGTTGTCTGATTTTTAATTTTTGGCGAGACCATTTGATGTTATGACCAGATGTG‑GTGCTACCTTATCTTAAGGATTTTGATAAAAATAATTATGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGTCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGAGGAT‑TTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGt                   <  2:173389/298‑1 (MQ=255)                                                                                                                                                                                                                      gtAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGAG‑TTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTa                  >  1:472150/1‑298 (MQ=255)                                                                                                                                                                                                                      gtAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGCCAAGTGTTTTGGCACGATCAACAACGCAAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAATAGTGAAGGGGAG‑TTGGACGATGCCGTTCCTGTGGCGTGGAATGCCGTAAATTa                  >  2:791001/1‑298 (MQ=37)                                                                                                                                                                                                                            tatCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCCGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGAG‑TTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCaaa            >  1:123854/1‑298 (MQ=255)                                                                                                                                                                                                                                  gTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATCACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGAG‑TTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACg      >  1:255004/1‑298 (MQ=255)                                                                                                                                                                                                                                   ttCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGCTATAACCTTTTGCGGGGACACGC‑TCACGATTTTTGAACAGTGAAGGGGAG‑TTGTTATATGCCGTTCCTGCAGCGGGAGCTCCACAAAAGTAATCAAAACCACgg     >  2:258994/1‑298 (MQ=37)                                                                                                                                                                                                                                     ctcacaTTTTTTTGTTTTATTTTTTATTTTTGCCGAAACCATTTTATCATATGACAAGATGTG‑TATTTACTTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGAG‑TTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGtt   <  2:787268/296‑1 (MQ=255)                                                                                                                                                                                                                                      gCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGAG‑TTGTAATATTCCGTTCCTGCAGCGTGGAATCCAATAAAATAATCAAAGCCACGGttt  >  2:448242/1‑298 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                       GGTATACAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTTGTCTGATTATTGATTTTTCGCGAAACCATTTGATCATATGACAAGATGTG‑TATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGT‑TCACGATTTTTGAACAGTGAAGGGGAG‑TTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTT  >  CP006584/1645162‑1645688

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