Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	641,296	A→G	intergenic (+143743/‑58047)	unknown → / → unknown	tRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	641,296	0	A	G	98.6%	185.5 / NA	70	intergenic (+143743/‑58047)	unknown/unknown	tRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (49/20);  total (50/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTACAAGGTATGGCGCATTCTTTCCGTATGTTGAAGGGACATCAGGCAGAATACGAACATCGTATTGATATTGACCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCATCCGGTTACATGGCCGGGCTTGTAGGGGCGTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCAACGACAACCTGCAGGATTTGAAAACCGGCTTGTTGGTGGAAGC  >  1/641054‑641541                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      gTACAAGGTATGGCGCATTCTTTCCGTATGTTGAAGGGACATCAGGCAGAATACGAACATCGTATTGATATTGACCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTa                                                                                                                                                                                                                                               >  1:296503/1‑251 (MQ=255)      aGGTATGGCGCATTCTTTCCGTATGTTGAAGGGACATCAGGCAGAATACGAACATCGTATTGATATTGACCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATgg                                                                                                                                                                                                                                          >  1:267394/1‑251 (MQ=255)                  ttCTTTCCGTATGTTGAAGGGACATCAGGCAGAATACGAACATCGTATTGATATTGACCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAg                                                                                                                                                                                                                              >  1:12738/1‑251 (MQ=255)                              gTTGAAGGGACATCAGGCAGAATACGAACATCGTATTGATATTGACCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCa                                                                                                                                                                                                                  >  1:179233/1‑251 (MQ=255)                                                  aaTACGAACATCGTATTGATATTGACCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCggg                                                                                                                                                                                              >  1:808991/1‑251 (MQ=255)                                                          cATCGTATTGATATTGACCTATCGCCGAAAACCATGATTTACATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAg                                                                                                                                                                                                                              <  1:436610/211‑1 (MQ=255)                                                            tCGTATTGATATTGACCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCg                                                                                                                                                                                    >  1:709593/1‑251 (MQ=255)                                                              gTATTGATATTGACCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTc                                                                                                                                                                                      >  1:195809/1‑247 (MQ=255)                                                              gTATTGATATAGACCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGcc                                                                                                                                                                                                               >  1:66310/1‑222 (MQ=255)                                                                attgatattgaCCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTc                                                                                                                                                                                      >  1:326560/1‑245 (MQ=255)                                                                   gatattgaCCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATcg                                                                                                                                                                              >  1:301929/1‑250 (MQ=255)                                                                   gatattgaCCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATcg                                                                                                                                                                              >  1:35021/1‑250 (MQ=255)                                                                    atattgaCCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATcgcg                                                                                                                                                                            >  1:265304/1‑251 (MQ=255)                                                                          aCCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCATCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTAttt                                                                                                                                                                                       >  1:426085/1‑234 (MQ=255)                                                                                 gCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGcc                                                                                                                                                                                                               >  1:391051/1‑203 (MQ=255)                                                                                  ccGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGtt                                                                                                                                                              >  1:747084/1‑251 (MQ=255)                                                                                        aCCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCttctgtctctataca                                                                                                                                                        >  1:834400/1‑238 (MQ=255)                                                                                             gATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCg                                                                                                                                                   >  1:865477/1‑251 (MQ=255)                                                                                               tttATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACTGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGgtg                                                                                                                                                                <  1:435955/236‑1 (MQ=255)                                                                                                   tatTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATg                                                                                                                                             >  1:528561/1‑251 (MQ=255)                                                                                                                 cTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTAttt                                                                                                                               >  1:822086/1‑251 (MQ=255)                                                                                                                   gTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGa                                                                                                                             >  1:761527/1‑251 (MQ=255)                                                                                                                      aTCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAac                                                                                                                          <  1:579067/251‑1 (MQ=255)                                                                                                                       tCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAac                                                                                                                          >  1:393164/1‑250 (MQ=255)                                                                                                                           gATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCg                                                                                                                     >  1:473961/1‑251 (MQ=255)                                                                                                                               gtgATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATgg                                                                                                                 >  1:215345/1‑251 (MQ=255)                                                                                                                                                cTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCg                                                                                                >  1:638921/1‑251 (MQ=255)                                                                                                                                                        cTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGAcc                                                                                        <  1:682004/251‑1 (MQ=255)                                                                                                                                                        cTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGAcc                                                                                        <  1:824811/251‑1 (MQ=255)                                                                                                                                                              cACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTc                                                                                  <  1:117671/251‑1 (MQ=255)                                                                                                                                                              cACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTc                                                                                  <  1:559563/251‑1 (MQ=255)                                                                                                                                                              cACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTc                                                                                  <  1:706413/251‑1 (MQ=255)                                                                                                                                                                  gATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCg                                                                              >  1:509520/1‑251 (MQ=255)                                                                                                                                                                   aTTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCg                                                                              <  1:583124/250‑1 (MQ=255)                                                                                                                                                                   aTTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCg                                                                              <  1:286/250‑1 (MQ=255)                                                                                                                                                                       cACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTc                                                                          >  1:868166/1‑250 (MQ=255)                                                                                                                                                                           tGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCg                                                                     <  1:46182/251‑1 (MQ=255)                                                                                                                                                                            gAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGcc                                                                                                                                                                                                               >  1:585336/1‑112 (MQ=255)                                                                                                                                                                            gAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTAt                                                                                                                                                     >  1:37506/1‑170 (MQ=255)                                                                                                                                                                            gAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGt                                                                    >  1:781818/1‑251 (MQ=255)                                                                                                                                                                             aaTTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTc                                                                                                                                                                                  >  1:75168/1‑140 (MQ=255)                                                                                                                                                                                  acacTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTca                                                               >  1:868192/1‑250 (MQ=255)                                                                                                                                                                                   cacTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATGGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTcacc                                                             <  1:582714/251‑1 (MQ=255)                                                                                                                                                                                     ctctTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGCTTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGtt                                                                                                                                                          >  1:492782/1‑156 (MQ=255)                                                                                                                                                                                     ctctTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGaa                                                                                                                            <  1:484920/186‑1 (MQ=255)                                                                                                                                                                                     ctctTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCg                                                                                                                     >  1:660813/1‑193 (MQ=255)                                                                                                                                                                                     ctctTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCt                                                                                       <  1:118085/223‑1 (MQ=255)                                                                                                                                                                                     ctctTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGc                                                                             <  1:374893/233‑1 (MQ=255)                                                                                                                                                                                     ctctTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTcacca                                                            >  1:660792/1‑250 (MQ=255)                                                                                                                                                                                     ctctTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTcacca                                                            >  1:138565/1‑250 (MQ=255)                                                                                                                                                                                       ctTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTcacc                                                             >  1:695602/1‑247 (MQ=255)                                                                                                                                                                                         tATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCa                                                       >  1:806741/1‑251 (MQ=255)                                                                                                                                                                                         tATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCa                                                       >  1:330189/1‑251 (MQ=255)                                                                                                                                                                                              gtggtggtTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCAc                                                                                >  1:724755/1‑221 (MQ=255)                                                                                                                                                                                              gtggtggtTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCAt                                                  >  1:415446/1‑251 (MQ=255)                                                                                                                                                                                                   ggtTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCg                                                                     <  1:707118/227‑1 (MQ=255)                                                                                                                                                                                                   ggtTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCaa                                                      >  1:637382/1‑242 (MQ=255)                                                                                                                                                                                                   ggtTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCa                                             >  1:881216/1‑251 (MQ=255)                                                                                                                                                                                                   ggtTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCa                                             >  1:584704/1‑251 (MQ=255)                                                                                                                                                                                                   ggtTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGAGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGc                                                        >  1:682568/1‑240 (MQ=255)                                                                                                                                                                                                      ttGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCAACg                                          <  1:256072/251‑1 (MQ=255)                                                                                                                                                                                                        gTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATg                                                                                                                                             >  1:663552/1‑150 (MQ=255)                                                                                                                                                                                                        gTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCAACGa                                         >  1:483732/1‑250 (MQ=255)                                                                                                                                                                                                        gTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCAACGa                                         >  1:164784/1‑250 (MQ=255)                                                                                                                                                                                                            aTTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCAACGACAAcc                                    <  1:190657/251‑1 (MQ=255)                                                                                                                                                                                                                  tGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCAACGACAACCTGCAgg                              <  1:84932/251‑1 (MQ=255)                                                                                                                                                                                                                      gTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGc                                                        >  1:354235/1‑221 (MQ=255)                                                                                                                                                                                                                        ttttATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTcaccac                                                           >  1:645144/1‑216 (MQ=255)                                                                                                                                                                                                                                gTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCAACGACAACCTGCAGGATTTGAAAACCgg                 >  1:267325/1‑250 (MQ=255)                                                                                                                                                                                                                                         tAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGAc                                                                                                               <  1:882114/147‑1 (MQ=255)                                                                                                                                                                                                                                               cagcgTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCAACGACAACCTGCAGGATTTGAAAACCGGCTTGTTGGTGGAAGc  <  1:486802/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      GTACAAGGTATGGCGCATTCTTTCCGTATGTTGAAGGGACATCAGGCAGAATACGAACATCGTATTGATATTGACCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCATCCGGTTACATGGCCGGGCTTGTAGGGGCGTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCAACGACAACCTGCAGGATTTGAAAACCGGCTTGTTGGTGGAAGC  >  1/641054‑641541

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