Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	641,326	G→C	intergenic (+143773/‑58017)	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,326	0	G	C	100.0%	216.6 / NA	77	intergenic (+143773/‑58017)	unknown/unknown	tRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (51/26);  total (51/26)

GTTGAAGGGACATCAGGCAGAATACGAACATCGTATTGATATTGACCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCATCCGGTTACATGGCCGGGCTTGTAGGGGCGTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCAACGACAACCTGCAGGATTTGAAAACCGGCTTGTTGGTGGAAGCCACCCCTTGGCGTCAACAGGTGGCATT  >  1/641083‑641568                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   gTTGAAGGGACATCAGGCAGAATACGAACATCGTATTGATATTGACCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCa                                                                                                                                                                                                                                             >  1:179233/1‑251 (MQ=255)                     aaTACGAACATCGTATTGATATTGACCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCggg                                                                                                                                                                                                                         >  1:808991/1‑251 (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:35021/1‑250 (MQ=255)                                      gatattgaCCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATcg                                                                                                                                                                                                         >  1:301929/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:706413/251‑1 (MQ=255)                                                                                                                                 cACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTc                                                                                                             <  1:559563/251‑1 (MQ=255)                                                                                                                                     gATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCg                                                                                                         >  1:509520/1‑251 (MQ=255)                                                                                                                                      aTTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCg                                                                                                         <  1:286/250‑1 (MQ=255)                                                                                                                                      aTTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCg                                                                                                         <  1:583124/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:584704/1‑251 (MQ=255)                                                                                                                                                                      ggtTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCa                                                                        >  1:881216/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:164784/1‑250 (MQ=255)                                                                                                                                                                           gTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCGTCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCAACGa                                                                    >  1:483732/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)                                                                                                                                                                                                                       tCCGGTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCCCAGCAACCATTTCCAACGACAACCTGCAGGATTTGAAAACCGGCTTGTTGGTGGAAGCCAccc                        <  1:61681/250‑1 (MQ=255)                                                                                                                                                                                                                          ggTTACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGAcc                                                                                                                   <  1:277278/156‑1 (MQ=255)                                                                                                                                                                                                                            ttACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCAACGACAACCTGCAGGATTTGAAAACCGGCTTGTTGGTGGAAGCCACCCCTTGGc                  <  1:538599/251‑1 (MQ=255)                                                                                                                                                                                                                             tACATGGCCGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCAACGACAACCTGCAGGATTTGAAAACCGGCTTGTTGGTGGAAGCCACCCCTTGGCg                 >  1:158384/1‑251 (MQ=255)                                                                                                                                                                                                                                     cGGGCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCAACGACAACCTGCAGGATTTGAAAACCGGCTTGTTGGTGGAAGCCACCCCTTGGCGTCAACAgg         >  1:803944/1‑251 (MQ=255)                                                                                                                                                                                                                                       ggCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGtt                                                                                                              <  1:671233/148‑1 (MQ=255)                                                                                                                                                                                                                                       ggCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCAACGACAAcc                                                               <  1:371416/195‑1 (MQ=255)                                                                                                                                                                                                                                       ggCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCAACGACAACCTGCAGGATTTGAAAACCGGCTTGTTGGTGGAAGCCACCCCTTGGCGTCAACAGGTg       >  1:376895/1‑251 (MQ=255)                                                                                                                                                                                                                                        gCTTGTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCAACGACAACCTGCAGGATTTGAAAACCGGCTTGTTGGTGGAAGCCACCCCTTGGCGTCAACAGGTg       <  1:777443/250‑1 (MQ=255)                                                                                                                                                                                                                                            gTAGGGGCCTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCAACGACAACCTGCAGGATTTGAAAACCGGCTTGTTGGTGGAAGCCACCCCTTGGCGTCAACAGGTGGCAtt  >  1:442378/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   GTTGAAGGGACATCAGGCAGAATACGAACATCGTATTGATATTGACCTATCGCCGAAAACCATGATTTATATTTTGATCGCGACTGTTATCTTGATTGTGATTTCTTTACATCACTTTATCGCTGCCGCACCGATTTCACCTGAATTAACACTCTTATTGGTGGTGGTTTGTACATTCCTTGCAGTTTTTATCGGTTTTTTCATAGCGGCAGCATCCGGTTACATGGCCGGGCTTGTAGGGGCGTCTTCCAGCCCGATTTCCGGTATCGGGATTATTTCCGTCATCGCGATTTCTTTGGTGTTGGTTTCTATCGGTAATGCATCCGGTTTATTTGAAACCACCGATGGACAGAAATTTTTGACCGCACTGACCTTGTTCACCGCTTCCATCGTGCTCACCACAGCAACCATTTCCAACGACAACCTGCAGGATTTGAAAACCGGCTTGTTGGTGGAAGCCACCCCTTGGCGTCAACAGGTGGCATT  >  1/641083‑641568

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: ‑