Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NZ_LN832404	3,872,045	T→G	15.8%	H32Q (CAT→CAG)	EK12AG100_RS19225 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_LN832404	3,872,045	0	T	G	15.8%	111.8 / 3.0	54	H32Q (CAT→CAG)	EK12AG100_RS19225	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/15);  new base G (7/2);  total (36/18)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAA  >  NZ_LN832404/3871752‑3872331                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                               tgaggaCGTTCGCGTTTGCCGTGCGGGTAATGTAGTAGAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACAAGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACGCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCAt                                                                                                                                                                                                                                                                                                <  1:853558/294‑1 (MQ=255)   agtgatttCGGGTTTGCCGGGCGTGTAATTTAGATCAAACTTAAATTGTTGTACTAAAATTTAGATAACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATAc                                                                                                                                                                                                                                                                                              <  2:786699/288‑1 (MQ=255)      cctgtCGCGTTTGCCGTGCGTGGAATGTAATCAAAACGTATATTTATGTAATACAAGTTAAATTACAAAAAGAACAATGCATAAAAAAATGACATTCGTCGAGCAGAAATCTGAAAAGGGAAATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCGTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGc                                                                                                                                                                                                                                                                                           <  2:418398/290‑1 (MQ=255)              gTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTGGTTTCCCATACCACAAAGGGTa                                                                                                                                                                                                                                                                                   >  1:627173/1‑294 (MQ=255)                      tgcgtgTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTCTTCTGCCATACCGCAAAGGGTAGTTCAGGa                                                                                                                                                                                                                                                                           >  1:893893/1‑294 (MQ=255)                            aaaagtatggcAAACTTAAATTGTTGTACGAAAATTTAGATAACAAAAAGAACAAGGCATAAAAAAAGGAGATGCGTCGGGAAGAAAGCTGAAAAGGGATATCAGGGGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTAt                                                                                                                                                                                                                                                                     <  2:467540/284‑1 (MQ=255)                             aaTGTAAAACAAATGTATTATGTTGTACTACAATTTATAGTACAAAAAGAAAAAGGCAAAAAAAAATGACATGCGTCGGGAAGAAATCTGAAAAGGAAAATCAGGCGCTAAACAGGAAGGAAAGAAGAGTATGCGTTCAACGGCTGAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTGTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATc                                                                                                                                                                                                                                                                    <  2:656636/294‑1 (MQ=255)                                      cAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGTTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGt                                                                                                                                                                                                                                                           >  1:593788/1‑294 (MQ=255)                                         aCTTATATTGTTGTACTAAAATTTAGATCACAAAAAGGATAATTCAGAAGAAAATGACTGGCGGCGGGAAGAAATCTGAAAAGGGAGATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAg                                                                                                                                                                                                                                                        <  2:195564/294‑1 (MQ=255)                                                           cAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTGCTATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTAt                                                                                                                                                                                                                                      >  1:775056/1‑294 (MQ=255)                                                            aaTTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAAAAAACAGTAGTTCAGGCAAGTATCAATCGCAGGCAGACAACATCACGGGTCTATc                                                                                                                                                                                                                                     >  1:805273/1‑294 (MQ=255)                                                                tAGATGAAAAAAAGAACAAAGGATAAAAAAATGGCATGGGTCGGGAAGAAATCTGAAAATGGAGATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAg                                                                                                                                                                                                                                 <  2:457022/294‑1 (MQ=255)                                                                tAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAg                                                                                                                                                                                                                                 >  1:188233/1‑294 (MQ=255)                                                                 agagCACAAAAAGAACAATGCATAAAAAAAAGACATGCGGCGGGAAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGAAATGAAAGAAGAGTATGCGTTAAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTGTAAGGGATGTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAgg                                                                                                                                                                                                                                <  2:565756/294‑1 (MQ=255)                                                                    tCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGCCTATCGCAGGCaa                                                                                                                                                                                                                             >  1:12479/1‑294 (MQ=255)                                                                                              aTGACAGGTGTGGGGAAGAAATCTGAAAAGGGGTAGCAGGTGCTAAACAGGAGGGAAAGGAGAGTATTCTTTGAACGGATTAGATACTCGTTTAAAGGATTAATCATGAAGTTGAAGTTTAAGGGATTTTTTAAGGCTGCCGGTTAATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTAt                                                                                                                                                                                                   <  1:797824/294‑1 (MQ=255)                                                                                               tGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTCTGCTAGACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGAAGGTAAAAGATAGTAATGGAACAAAAGGCTATg                                                                                                                                                                                                  >  1:282802/1‑294 (MQ=255)                                                                                                    ggtgTCGGTCAAAAATCTTAAAAGTGTTATCAGGCGCTAAACAGGATGGAAAGATGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTGGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGttttt                                                                                                                                                                                             <  2:564645/291‑1 (MQ=255)                                                                                                      cGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATAAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACACACGGCCATGTCATTGt                                                                                                                                                                                           >  2:749799/1‑294 (MQ=255)                                                                                                            gCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAACGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAACGATAGTAATGGAACAAAAGCCTATGTTTTTGAAGGGGa                                                                                                                                                                                     >  2:341001/1‑294 (MQ=255)                                                                                                                 aaTCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGACCAAAAGGCTATGTTTTTGTAGGGGAATCGc                                                                                                                                                                                >  1:114145/1‑294 (MQ=255)                                                                                                                 aaTCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTCTGCTATACCACAAAAGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGc                                                                                                                                                                                >  1:770298/1‑294 (MQ=255)                                                                                                                   gcTGAAAAGGTGTGTCAGGGGGTAAACGGGAGGGGAAGAAGGGTAGGCGTTCAACGGTGTAGCTACTCGTTTAAAGTATTAATCAGGAAGTTGAATTTAAAGGGATTTTTTAAGGGTGCCGGTTTAGTCGCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATAGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTg                                                                                                                                                                              <  2:59983/293‑1 (MQ=255)                                                                                                                          aggggtgtgaGGAGCTAAACACGAGGGGAGGTAGTAGATGCTGTAAACGGTTGAGTGAGTCGATTAAAGTATTGTTCATGATGTGGAAGTTAAAGGGATTTTTTAAGGTGGCGGGTATATTCCCACTGGGGCTGATGATTTCAGGTGGTATATCGCATGATCTGGTTTCCCAGAGCGCAAATGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGGAGGCAAAAGATTGTAATGGAACAAAAGGCTATGTTCTTGTAGGGGAATCGCTGGATGAcc                                                                                                                                                                       <  2:316395/285‑1 (MQ=255)                                                                                                                            ggagaccgggtggtAAAAAGGGGGGAAAGAAGGAGATGCTTGCCAGGGGTTAGCTATCCGGTTAAAGTATTAATGATGAAGATGAGTTTTAAGGGATTGTTTAAGGCTGCCGGTATATTCCCACTGGCGCTGATGCTTTCAGGCGGTAGCTCGTATGCACTGGTTTCCCGAACCGCAAAGGGTAGTTGAGGAAAGTATCAATCGCAGTCATACGCCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCtt                                                                                                                                                                     <  2:85564/281‑1 (MQ=255)                                                                                                                                       gcTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTATGTAGGGGAATCGCTGGATTACATTATCAATTATgg                                                                                                                                                          >  1:649025/1‑294 (MQ=255)                                                                                                                                        atAAACTGGATGGAAATAAGCGTATGCTTTCAACGACCTATCTACTCGTTTTAATGATCAATCATTAAGTTGAAGTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCAAATGACGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCAGACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAATGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGt                                                                                                                                                         <  2:775056/293‑1 (MQ=255)                                                                                                                                                 tgggTAAGGAGAGATTGCTTTCAAAGGTTTAGTTACACGTTTAAAGGGTTTATCATTAATTTGAATTTTAAGGGAGTTTTTAATGCTTCCGGTTTACTCACACTGGCGCTGATGCTTTCAGGCTGTATCTCGTAGGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTAACAATCGCAGTCAGACACCATCACTGTGCGATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGTGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGAc                                                                                                                                                <  2:474253/293‑1 (MQ=255)                                                                                                                                                 aggGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTGCTAGACAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCAACACTGGGTCATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTATGTAGGGGAATAGCTGGATTACATTATCACTTATGGTGACGATGAc                                                                                                                                                >  1:160893/1‑294 (MQ=255)                                                                                                                                                         aagaGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGACAACAAAAGGAAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAACGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTaa                                                                                                                                        >  1:731510/1‑294 (MQ=255)                                                                                                                                                              tATGTTGTCAAAGGCTAGTCTACTCGTGTAAGGGGTGAATGATGAAGTTGAAGTTTAAGGGATGTTTGAAGGCTGCCGGTATATGCCCACTGGCGCAGATGGTATAAGGATGTAGCTAGTGTGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGCATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAGTGGAACAAAAGGCTATGTTTTTGTAGGGGTATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGc                                                                                                                                   <  2:430364/294‑1 (MQ=255)                                                                                                                                                                      cAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCCGGCAAAAGATAGTACTGGAACCAAAGGCTATGGTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGTTGGTTCCGATGACATCGTTACGCTGCTCAAGGac                                                                                                                           >  2:450041/1‑293 (MQ=255)                                                                                                                                                                             ttAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCCCTGGGCTATCGCAGGCAAAAGCTAGTAATGGAACACAAGGCTATGTTTTTGTAGGGGCATCGATGGATTCCCTTAACCCTGATGATGCCGCGGCCATAGTTACGATCATAACTGATCCAGcca                                                                                                                    >  2:524635/1‑292 (MQ=255)                                                                                                                                                                             ttAGCTACTCGTTCAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTCGTTCTGGAAAGTATCAATCGAAGTCCGACCCCATCACTGGACTAGCGCCGCCAAAAGAGAGCAATGGAACAAAAGGCTATGTTTTTGTAGGGGAACCGCTGGGTCACCTTATGACTGATGGTGCCGATGTATTGGTTAGGATGGTCAATGTACCAgaca                                                                                                                    >  2:152016/1‑291 (MQ=255)                                                                                                                                                                                 cTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTCGTCTCCCATACAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCCGGCAAAAGCTAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTTCCGATGACATCGTTAAGATGCCCAATGCTCACGCACTTaa                                                                                                                >  1:278508/1‑294 (MQ=255)                                                                                                                                                                                           aaaGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTTGGCTATCGCAGGCAAAAGATAGCAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGTTGGATTACCTTATCACTTGTGGTTCCGATGCATTCGTTAAGGTGCTCAATGATCCAGCACTTAACCGGCaaaac                                                                                                      >  2:247232/1‑293 (MQ=255)                                                                                                                                                                                               gATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGTTTTGCCTTACAGCAAAGGTTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCCCTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTc                                                                                                  >  1:37940/1‑294 (MQ=255)                                                                                                                                                                                                         gAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGACCTTAACCGGCACAATCTTAAGGTTGCAGa                                                                                        >  2:723368/1‑294 (MQ=255)                                                                                                                                                                                                             tggAATTTTAAGGGAGTTTTTAAGGCAGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGAc                                                                                    <  1:467247/292‑1 (MQ=255)                                                                                                                                                                                                                 aTTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCACACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCAGATTACATCGTTAAGATGCTCAATGATCACGCACTTAACAGGCCCAATATTCGTGTTACCGATGACTCaa                                                                                >  2:384892/1‑294 (MQ=255)                                                                                                                                                                                                                 aTTTGAAGGGATTTTTTACGGCTTCCTGTTTATTCCCACTGGCTCTGATGCTTTCAGGCTGTAACTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTTAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGTAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCaa                                                                                <  1:719217/294‑1 (MQ=255)                                                                                                                                                                                                                            ttttttAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTa                                                                     >  1:569289/1‑294 (MQ=255)                                                                                                                                                                                                                              ttttAAGGCTGCCGGTTTATGCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTaaa                                                                   <  1:723368/294‑1 (MQ=255)                                                                                                                                                                                                                               tttAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCTTTAAGATGCTCAATGATCCAGCACTTAAACGGCACAATAATACGGTTGCCGATGAAGCAAGATTTGTTTTACAt                                                                  >  2:839279/1‑294 (MQ=255)                                                                                                                                                                                                                                ttAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTTTGCCTTAACGCAAAAGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGACCAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATTGTGCCGATGACATCGTTAAGATGCTCATTGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGCCGCAAGATATGTTTTAAAta                                                                 >  1:742624/1‑293 (MQ=255)                                                                                                                                                                                                                                  aaGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCCGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGTTGCCGATGACAACGTTAAGATGCTCCATGATCCAGCACTTAATCTGCACCACAGTCAGGTTGCCGGAGAAGCACAAATTGTTTTgaaagat                                                               >  2:467704/1‑290 (MQ=255)                                                                                                                                                                                                                                        gCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTTCCGATGACATCGTTAAGATGCTCAATGATCAAACACTTAACCGGCCCCATATTCAGGTTGCCGATGACGCAAGATTTGCTTTCAATGCGGGGaaa                                                         >  2:448778/1‑294 (MQ=255)                                                                                                                                                                                                                                                aTTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACTATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGTATTACATTATAACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGCTTTGTTTTAAATGCGGGGAACAAGAAAtt                                                 >  1:867330/1‑294 (MQ=255)                                                                                                                                                                                                                                                    ccACTTGCGCTCATGGTGTCAGGCGGTATCTCGTATTCTAGTGTGTCCCATACCGCAAAGGGTAGTTTAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGATATCGAATGCAAAAGATAGTAAGGGAACAAAAGGCTATGTTTTTGTAGGGGAAGCACTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACTGGAACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAAGGCGGGGAAAAAGAAATTTAcc                                             <  2:296439/294‑1 (MQ=255)                                                                                                                                                                                                                                                        tGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTTTGCCAGACCGAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATAACGCAAGCATTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGca                                         >  1:682348/1‑294 (MQ=255)                                                                                                                                                                                                                                                                  gCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGAACCAGCACTTACCCGGCCCAATCTTCAGGTTGCCGATTCCGCAAGATTTGTTTTAAATACGGGGAAACAGAAATTTACCGG‑ACAatcactcat                              >  2:797397/1‑287 (MQ=255)                                                                                                                                                                                                                                                                       ccagcggtgtggcGTATGCTATGGTGTGCCATACCGCAAGTGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCAGCACGGGGCTAGCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTAGGGTTTGGTAGGGGAATCGCTGGATTACCTTATCACTGAGGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAGATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTtact                          <  2:114997/282‑1 (MQ=255)                                                                                                                                                                                                                                                                                 ctcGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATTACGCAAGCTTTGTTTTAAATGCGGGGAAAACGAACTTTACCGGCACAATATGGCATtactactcagggca                >  2:343195/1‑287 (MQ=255)                                                                                                                                                                                                                                                                                 ctcGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCAATGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCACGATTTGTTTTACATGCGGGGAAAAAGCAATTTACCGGCACCATCTCGCTTGCCTACTCCCGGac                >  2:439216/1‑293 (MQ=255)                                                                                                                                                                                                                                                                                      aTGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACg           >  2:77828/1‑294 (MQ=255)                                                                                                                                                                                                                                                                                      aTGCTCGGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACGGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTGTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGAAGAAATCGATAAGATGCCCAATGAGCATCGACTTAACCGGCACAATACGCAGGTTGCCGATGACGCAAGTGTTGTTATAAATGCGGGGAAAAAGAAATTTCCAGGCACACAACCGCTGTACTACTACGGCAAAAACg           >  2:100970/1‑294 (MQ=255)                                                                                                                                                                                                                                                                                             ggTTTCCCATACCGCGAAGGGTGGTTCAGGAAAGGATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTAGGTTTTTGTAGGGAAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGAACCAGCACTTGACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGAGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACgaagaaga    <  2:12479/294‑1 (MQ=255)                                                                                                                                                                                                                                                                                               tttCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGaaa  >  2:168774/1‑294 (MQ=255)                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                               TGATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAA  >  NZ_LN832404/3871752‑3872331

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