Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,987,232	G→A	A27V (GCG→GTG)	unknown ←	16S rRNA (guanine(1207)‑N(2))‑methyltransferase (EC 2.1.1.172)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,987,232	0	G	A	100.0%	156.0 / NA	57	A27V (GCG→GTG)	unknown	16S rRNA (guanine(1207)‑N(2))‑methyltransferase (EC 2.1.1.172)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/33);  total (24/33)

CCAGCAAAACCTCCTGCCCGGCGGTGCATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGATGTTGTCTGTTGGCATAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGACGTTGTAAATCGGCGAGAAAATCATCGTTCGCCGCCCCGCTGAGCAGCAGGGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAAAAACGCGCCGACACACTTCCCGCCGAGATTTATTATTGCAGCAAATGGGCATCAGCCAATGGCAATTGG  >  1/1986993‑1987454                                                                                                                                                                                                                                                |                                                                                                                                                                                                                               ccaGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCAcccc                                                                                                                                                                                                                         >  1:748261/1‑247 (MQ=255)  cagtaaaaCTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGa                                                                                                                                                                                                                    >  1:70827/5‑251 (MQ=255)    gtaaaaCTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGa                                                                                                                                                                                                                    >  1:69385/3‑249 (MQ=255)                     gcggTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCg                                                                                                                                                                                                   >  1:374102/1‑249 (MQ=255)                     gcggTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCg                                                                                                                                                                                                   >  1:621028/1‑249 (MQ=255)                             ttGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGa                                                                                                                                                                                                                    <  1:53771/224‑1 (MQ=255)                                 gagagaAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCaa                                                                                                                                                                                     >  1:161844/1‑251 (MQ=255)                                     gaAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCgg                                                                                                                                                                                                                       <  1:500441/213‑2 (MQ=255)                                          tGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTc                                                                                                                                                                            <  1:111770/251‑1 (MQ=255)                                                  cTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAacac                                                                                                                                                                      <  1:696085/249‑1 (MQ=255)                                                  cTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACtt                                                                                                                                                                    <  1:833122/251‑1 (MQ=255)                                                    ggAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCg                                                                                                                                                                  >  1:227864/1‑251 (MQ=255)                                                                 cTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACaaa                                                                                                                                                     <  1:798096/251‑1 (MQ=255)                                                                  tGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAAt                                                                                                                                                    >  1:22942/1‑251 (MQ=255)                                                                  tGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAAt                                                                                                                                                    >  1:161682/1‑251 (MQ=255)                                                                   gTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGcagcagag                                                                                                                                                                                                           <  1:17806/195‑3 (MQ=255)                                                                               tCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATc                                                                                                                                       <  1:657975/251‑1 (MQ=255)                                                                               tCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATc                                                                                                                                       <  1:288961/251‑1 (MQ=255)                                                                               tCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATc                                                                                                                                       <  1:380476/251‑1 (MQ=255)                                                                                 cAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCt                                                                                                                                     <  1:284892/251‑1 (MQ=255)                                                                                 cAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCt                                                                                                                                     <  1:117471/251‑1 (MQ=255)                                                                                     aaaaCAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGc                                                                                                                                                                                 <  1:790000/203‑1 (MQ=255)                                                                                     aaaaCAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATgg                                                                                                                                 <  1:842067/251‑1 (MQ=255)                                                                                      aaaCAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGc                                                                                                                                <  1:268589/251‑1 (MQ=255)                                                                                       aaCAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCt                                                                                                                               <  1:282305/251‑1 (MQ=255)                                                                                            gATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGAc                                                                                                                                                                                                         <  1:753793/172‑1 (MQ=255)                                                                                             aTCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCgg                                                                                                                                                                                                                       >  1:102828/1‑156 (MQ=255)                                                                                                           gCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAACAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGc                                                                                                                                >  1:442945/1‑230 (MQ=255)                                                                                                              aTGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCAt                                                                                                        <  1:205206/251‑1 (MQ=255)                                                                                                                 aaaTTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATtat                                                                                                     <  1:325650/251‑1 (MQ=255)                                                                                                                   aTTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATAc                                                                                                   <  1:560422/251‑1 (MQ=255)                                                                                                                   aTTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATAc                                                                                                   <  1:594207/251‑1 (MQ=255)                                                                                                                                  ccACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACaaa                                                                                                                                                     >  1:481550/1‑186 (MQ=255)                                                                                                                                        cACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAac                                                                                                                                                                        >  1:462963/1‑161 (MQ=255)                                                                                                                                        cACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCt                                                                                                                                     <  1:687501/196‑1 (MQ=255)                                                                                                                                          cTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGc                                                                                                                                                                                 <  1:742007/150‑1 (MQ=255)                                                                                                                                          cTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGc                                                                                                                                >  1:382188/1‑199 (MQ=255)                                                                                                                                            ttGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATc                                                                                                                                       <  1:393692/190‑1 (MQ=255)                                                                                                                                              gtgGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTg                                                                                                 >  1:549478/1‑226 (MQ=255)                                                                                                                                                   gTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATAc                                                                                                   >  1:656655/1‑219 (MQ=255)                                                                                                                                                    tttGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAAc                                                                   <  1:673120/250‑1 (MQ=255)                                                                                                                                                    tttGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAAc                                                                   <  1:503310/250‑1 (MQ=255)                                                                                                                                                               cggTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTgg                                                                                                                         >  1:739912/3‑185 (MQ=255)                                                                                                                                                               cggTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAACGCGCCGacacac                                                       <  1:564288/249‑1 (MQ=255)                                                                                                                                                                    ggcgtAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAACGCGCCGACACACTTcc                                                   >  1:10092/5‑250 (MQ=255)                                                                                                                                                                          gTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGAtt                                                                                    >  1:598067/1‑211 (MQ=255)                                                                                                                                                                            cAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATcc                                                                                                                                      <  1:753925/159‑1 (MQ=255)                                                                                                                                                                                      gCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATc                                                                                                                                       >  1:666166/1‑148 (MQ=255)                                                                                                                                                                                       cTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATc                                                                                                                                       >  1:717476/1‑147 (MQ=255)                                                                                                                                                                                             gcgcGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATAc                                                                                                   <  1:583273/177‑1 (MQ=255)                                                                                                                                                                                                  ggCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAACGCGCCGACACACTTCCCGCCGAGATTTATTATTGCAGCAAATGGGCa                    >  1:450767/1‑251 (MQ=255)                                                                                                                                                                                                        gAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAACGCGCCGACACACTTCCCGCCGAGATTTATTATTGCAGCAAATGGGCATCAGc               >  1:536005/1‑250 (MQ=255)                                                                                                                                                                                                          aTCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAACGCGCCGACACACTTCCCGCCGAGATTTATTATTGCAGCAAATGGGCATCAGCCaa            >  1:185177/1‑251 (MQ=255)                                                                                                                                                                                                            cccTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGata                                                                             <  1:133275/184‑1 (MQ=255)                                                                                                                                                                                                              ctggcGTTGTAAATCGGCGAGAAAATCATCATTCACCCCCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAACGCGCCGACACACTTCCCGCCGAGATTTATTATTGCAGCAAATGGGCATCAGCCAATGGc        <  1:262236/247‑1 (MQ=255)                                                                                                                                                                                                              ctggcGTTGTAAATCGGCGAGAAAATCATCATTCACCCCCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAACGCGCCGACACACTTCCCGCCGAGATTTATTATTGCAGCAAATGGGCATCAGCCAATGGc        <  1:81522/247‑1 (MQ=255)                                                                                                                                                                                                                    ctGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAACGCGCCGACACACTTCCCGCCGAGATTTATTATTGCAGCAAATGGGCATCAGCCAATGGCAATTgg  <  1:894461/250‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                               CCAGCAAAACCTCCTGCCCGGCGGTGCATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGATGTTGTCTGTTGGCATAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGACGTTGTAAATCGGCGAGAAAATCATCGTTCGCCGCCCCGCTGAGCAGCAGGGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAAAAACGCGCCGACACACTTCCCGCCGAGATTTATTATTGCAGCAAATGGGCATCAGCCAATGGCAATTGG  >  1/1986993‑1987454

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