Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,326,337	C→T	I286I (ATC→ATT)	unknown →	Chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,326,337	0	C	T	98.6%	193.5 / ‑3.3	73	coding (858/1902 nt)	unknown	Chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (40/32);  minor base . (1/0);  total (41/32)
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.91e-01

ATGAAGTAGGCGGCGAAAAAACCTTTGAAGTATTAGCAACCAACGGTGATACACACTTAGGTGGTGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATCACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGATTTTAGTCGGCGGTCAAACCCGTATGCCGCTTGTTCAACAAGAAGTTGAAAAATTCTTCGGCAAAGCACCACGTAAAGACGTGAATCC  >  1/1326102‑1326579                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   aTGAAGTAGGCGGCGAAAAAACCTTTGAAGTATTAGCAACCAACGGTGATACACACTTAGGTGGTGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAAc                                                                                                                                                                                                                                      >  1:581854/1‑250 (MQ=255)   gAAGTAGGCGGCGAAAAAACCTTTGAAGTATTAGCAACCAACGGTGATACACACTTAGGTGGTGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCgg                                                                                                                                                                                                                                   >  1:496244/1‑251 (MQ=255)           cggcgAAAAAACCTTTGAAGTATTAGCAACCAACGGTGATACACACTTAGGTGGTGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAac                                                                                                                                                                                                                           <  1:511774/251‑1 (MQ=255)                   aaaCCTTTGAAGTATTAGCAACCAACGGTGATACACACTTAGGTGGTGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAAc                                                                                                                                                                                                                   <  1:276774/251‑1 (MQ=255)                     aCCTTTGAAGTATTAGCAACCAACGGTGATACACACTTAGGTGGTGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACATTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACAt                                                                                                                                                                                                                 <  1:284793/251‑1 (MQ=255)                      ccTTTGAAGTATTAGCAACCAACGGTGATACACTCTTAGGTGGAGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAAc                                                                                                                                                                                                                   >  1:283366/1‑248 (MQ=255)                         ttGAAGTATTAGCAACCAACGGTGATACACACTTAGGTGGTGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCaaa                                                                                                                                                                                                             >  1:219959/1‑251 (MQ=255)                          tGAAGTATTAGCAACCAACGGTGATACACACTTAGGTGGTGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAAt                                                                                                                                                                                                            <  1:201842/251‑1 (MQ=255)                                aTTAGCAACCAACGGTGATACACACTTAGGTGGTGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACt                                                                                                                                                                                                       >  1:326456/1‑250 (MQ=255)                                               tGATACACACTTAGGTGGTGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTagag                                                                                                                                                                                        <  1:463393/250‑1 (MQ=255)                                                gATACACACTTAGGTGGTGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTa                                                                                                                                                                                                                      <  1:862455/219‑1 (MQ=255)                                                      cacTTAGGTGGTGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGt                                                                                                                                                                                >  1:726450/1‑251 (MQ=255)                                                      cacTTAGGTGGTGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGt                                                                                                                                                                                >  1:484523/1‑251 (MQ=255)                                                      cacTTAGGTGGTGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACAT‑ACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGtt                                                                                                                                                                               >  1:738342/1‑251 (MQ=255)                                                        cTTAGGTGGTGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAg                                                                                                                                                                                                                                                >  1:882060/1‑185 (MQ=255)                                                        cTTAGGTGGTGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGt                                                                                                                                                                                >  1:297628/1‑249 (MQ=255)                                                        cTTAGGTGGTGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTg                                                                                                                                                                              <  1:672463/251‑1 (MQ=255)                                                            ggtggtGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAAGCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAac                                                                                                                                                                                                                           <  1:374374/202‑1 (MQ=255)                                                              tggtgAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGAt                                                                                                                                                                         <  1:11631/250‑1 (MQ=255)                                                                gtgAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTa                                                                                                                                                                      >  1:582206/1‑251 (MQ=255)                                                                gtgAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTa                                                                                                                                                                      >  1:384713/1‑251 (MQ=255)                                                                       cTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGAt                                                                                                                                                                         <  1:113394/241‑1 (MQ=255)                                                                         ttGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGAt                                                                                                                                                                         <  1:900576/239‑1 (MQ=255)                                                                                ccGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTg                                                                                                                                                                              <  1:88024/227‑1 (MQ=255)                                                                                ccGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTg                                                                                                                                                      >  1:528540/1‑251 (MQ=255)                                                                                   tATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGc                                                                                                                                                   <  1:309048/251‑1 (MQ=255)                                                                                    atcatcAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGcc                                                                                                                                                  <  1:80635/251‑1 (MQ=255)                                                                                    atcatcAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGc                                                                                                                                                   >  1:750624/1‑250 (MQ=255)                                                                                    atcatcAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGc                                                                                                                                                   >  1:742959/1‑250 (MQ=255)                                                                                       atcaACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGc                                                                                                                                                   >  1:834477/1‑247 (MQ=255)                                                                                         caACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTa                                                                                                                                              >  1:242090/1‑250 (MQ=255)                                                                                         caACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTAc                                                                                                                                             <  1:745482/251‑1 (MQ=255)                                                                                           actactTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGt                                                                                                                                           >  1:380063/1‑251 (MQ=255)                                                                                            ctactTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAAc                                                                                                                                                                                                                                      >  1:694296/1‑159 (MQ=255)                                                                                             tactTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTAt                                                                                                                                         <  1:270849/251‑1 (MQ=255)                                                                                              actTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCaaaa                                                                                                                                                                                              >  1:334740/1‑197 (MQ=255)                                                                                              actTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTAtt                                                                                                                                        >  1:653263/1‑251 (MQ=255)                                                                                              actTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTAtt                                                                                                                                        >  1:529976/1‑251 (MQ=255)                                                                                                ttAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGc                                                                                                                                      <  1:541108/251‑1 (MQ=255)                                                                                                 tAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTAc                                                                                                                                             >  1:301013/1‑243 (MQ=255)                                                                                                   gTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGt                                                                                                                                                                    <  1:13491/218‑1 (MQ=255)                                                                                                          aGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGa                                                                                                                             <  1:64054/250‑1 (MQ=255)                                                                                                           gTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGt                                                                                                                                           <  1:226131/235‑1 (MQ=255)                                                                                                           gTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATg                                                                                                                           >  1:857293/1‑251 (MQ=255)                                                                                                           gTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATg                                                                                                                           >  1:813268/1‑251 (MQ=255)                                                                                                                     gAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAAt                                                                                                                  >  1:603209/1‑250 (MQ=255)                                                                                                                          gggCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGc                                                                                                                                                                                                                                               <  1:659936/120‑1 (MQ=255)                                                                                                                            gCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGt                                                                                                            >  1:77017/1‑249 (MQ=255)                                                                                                                                 gACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTaaaa                                                                                                                               <  1:350099/225‑1 (MQ=255)                                                                                                                                     tACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATcgac                                                                                                 <  1:290244/251‑1 (MQ=255)                                                                                                                                      aCGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATcgacg                                                                                                >  1:566563/1‑251 (MQ=255)                                                                                                                                        gCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGt                                                                                                                                           <  1:676422/206‑1 (MQ=255)                                                                                                                                              aCCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGAtt                                                                                        >  1:436278/1‑251 (MQ=255)                                                                                                                                               cccGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAAGCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGAttt                                                                                       <  1:430410/251‑1 (MQ=255)                                                                                                                                                       caatgcaaCGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGATTTTAGTcgg                                                                                >  1:857323/1‑250 (MQ=255)                                                                                                                                                                gTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAg                                                                                                                                                                                                                                                >  1:101861/1‑81 (MQ=255)                                                                                                                                                                gTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGt                                                                                                                                                                                <  1:270795/145‑1 (MQ=255)                                                                                                                                                                gTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGATTTTAGTCGGCGGTCAAAc                                                                       >  1:562263/1‑250 (MQ=255)                                                                                                                                                                gTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGATTTTAGTCGGCGGTCAAAc                                                                       >  1:816198/1‑250 (MQ=255)                                                                                                                                                                         aagCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGATTTTAGTCGGCGGTCAAACCCGTATGc                                                               <  1:633972/249‑1 (MQ=255)                                                                                                                                                                                          aaaaTCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGATTTTAGTCGGCGGTCAAACCCGTATGCCGCTTGTTCAACAAGAAGt                                            <  1:298629/251‑1 (MQ=255)                                                                                                                                                                                             aTCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGATTTTAGTCGGCGGTCAAACCCGTATGCCGCTTGTTCAACAAGAAGTTGa                                         >  1:816299/1‑251 (MQ=255)                                                                                                                                                                                                  gCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAac                                                                                                                                                                                                                           <  1:696523/68‑1 (MQ=255)                                                                                                                                                                                                  gCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGATTTTAGTCGGCGGTCAAACCCGTATGCCGCTTGTTCAACAAGAAGTTGAAAAAt                                    >  1:835022/1‑251 (MQ=255)                                                                                                                                                                                                     ttcttcCGCCCAACAAACCGATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGATTTTAGTCGGCGGTCAAACCCGTATGCCGCTTGTTCAACAAGAAGTTGAAAAAtt                                   >  1:304170/1‑249 (MQ=255)                                                                                                                                                                                                                        gATGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGATTTTAGTCGGCGGTCAAACCCGTATGCCGCTTGTTCAACAAGAAGTTGAAAAATTCTTCGGCAAAGCACCACGTaa              >  1:359143/1‑251 (MQ=255)                                                                                                                                                                                                                         aTGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGATTTTAGTCGGCGGTCaaa                                                                        >  1:617883/1‑192 (MQ=255)                                                                                                                                                                                                                         aTGTAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGATTTTAGTCGGCGGTCAAACCCGTATGCCGCTTGTTCAACAAGAAGTTGAAAAATTCTTCGGCAAAGCACCACGTaaa             >  1:277256/1‑251 (MQ=255)                                                                                                                                                                                                                            tAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGGGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGATTTTAGTCGGCGGTCAAACCCGTATGCCGCTTGTTCAACAAGAAGTTGAAAAATTCTTCGGCAAAGCACCACGTAAAGa           <  1:347055/250‑1 (MQ=255)                                                                                                                                                                                                                            tAAATCTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGATTTTAGTCGGCGGTCAAACCCGTATGCCGCTTGTTCAACAAGAAGTTGAAAAATTCTTCGGCAAAGCACCACGTAAAGa           <  1:651470/250‑1 (MQ=255)                                                                                                                                                                                                                                 cTTCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGATTTTAGTCGGCGGTCAAACCCGTATGCCGCTTGTTCAACAAGAAGTTGAAAAATTCTTCGGCAAAGCACCACGTAAAGACGTGaa     >  1:181041/1‑251 (MQ=255)                                                                                                                                                                                                                                   tCCTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTAtt                                                                                                                                        >  1:565157/1‑118 (MQ=255)                                                                                                                                                                                                                                    ccTTACATTACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGATTTTAGTCGGCGGTCAAACCCGTATGCCGCTTGTTCAACAAGAAGTTGAAAAATTCTTCGGCAAAGCACCACGTAAAGACGTGAATcc  <  1:424162/251‑1 (MQ=255)                                                                                                                                                                                                                                      ttacattacaGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGATTTTAGTCGGCGGTCAAACCCGTATGCCGCTTGTTCAACAAGAAGTTGAAAAATTCTTCGGCAAAGcac                     >  1:654479/1‑230 (MQ=255)                                                                                                                                                                                                                                      ttacattacaGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTACAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGATTTTAGTCGGCGGTCAAACCCGTATGCCGCTTGTTCAACAAGAAGt                                            >  1:196177/1‑207 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   ATGAAGTAGGCGGCGAAAAAACCTTTGAAGTATTAGCAACCAACGGTGATACACACTTAGGTGGTGAAGACTTTGATAACCGTATCATCAACTACTTAGTGGATGAGTTCAAAAAAGAACAGGGCATCGACCTACGCAACGACCCGCTTGCAATGCAACGTTTAAAAGAAGCGGCAGAAAAAGCCAAAATCGAGCTTTCTTCCGCCCAACAAACCGATGTAAATCTTCCTTACATCACAGCGGATGCAACCGGTCCTAAACACTTAAACATCAAATTGACTCGTGCAAAATTAGAGTCTTTGGTTGAGGATTTAGTGGCGAAATCCCTTGAGCCTGTACGTATTGCATTAAAAGATGCAGGCAAATCCCCGTCTGAAATCGACGATGTGATTTTAGTCGGCGGTCAAACCCGTATGCCGCTTGTTCAACAAGAAGTTGAAAAATTCTTCGGCAAAGCACCACGTAAAGACGTGAATCC  >  1/1326102‑1326579

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