Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,986,312	C→T	G334R (GGA→AGA)	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,986,312	0	C	T	100.0%	212.4 / NA	71	G334R (GGA→AGA)	unknown	16S rRNA (guanine(1207)‑N(2))‑methyltransferase (EC 2.1.1.172)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (39/32);  total (39/32)

ACACCCTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTTGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTTCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGGCGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCCGTGCACCGAATAAACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCATCATGAAACGGCGGGTTGGAAATGATCAGATCAAATTTGCCGTCAATATGGGAAAAC  >  1/1986073‑1986556                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     acacCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAAt                                                                                                                                                                                                                                           <  1:575442/251‑1 (MQ=255)     cttCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATatac                                                                                                                                                                                                                                       <  1:789987/249‑4 (MQ=255)               gtCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATat                                                                                                                                                                                                                                         <  1:327243/239‑2 (MQ=255)               gtCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTg                                                                                                                                                                                                                             >  1:824137/1‑251 (MQ=255)               gtCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTg                                                                                                                                                                                                                             >  1:644851/1‑251 (MQ=255)                                       tcgTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGa                                                                                                                                                                                                     >  1:849002/3‑251 (MQ=255)                                            aTTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATAcc                                                                                                                                                                                                                                      <  1:631294/213‑1 (MQ=255)                                            aTTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGAc                                                                                                                                                                                                    >  1:331687/1‑247 (MQ=255)                                                  aacaacGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGaaa                                                                                                                                                                                            <  1:55707/249‑1 (MQ=255)                                                         cgcCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGAtt                                                                                                                                                                                    >  1:556395/1‑250 (MQ=255)                                                           ggctgcacaggCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAg                                                                                                                                                                                 >  1:836921/10‑251 (MQ=255)                                                                          gTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACt                                                                                                                                                                                                                               <  1:325815/190‑1 (MQ=255)                                                                          gTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAgg                                                                                                                                                                  >  1:140673/1‑251 (MQ=255)                                                                           tGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGc                                                                                                                                                                 <  1:404751/251‑1 (MQ=255)                                                                           tGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGc                                                                                                                                                                 <  1:518554/251‑1 (MQ=255)                                                                                    gTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTAAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCa                                                                                                                                                        >  1:288100/1‑251 (MQ=255)                                                                                        caATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGaaa                                                                                                                                                                                            >  1:645939/1‑211 (MQ=255)                                                                                        caATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTg                                                                                                                                                    >  1:269025/1‑251 (MQ=255)                                                                                        caATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTg                                                                                                                                                    >  1:680419/1‑251 (MQ=255)                                                                                          atatTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGc                                                                                                                                                   >  1:61378/1‑250 (MQ=255)                                                                                                      tCGCCAAGTTTTACAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGc                                                                                                                                      <  1:278427/251‑1 (MQ=255)                                                                                                         ccAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTc                                                                                                                                    <  1:305890/250‑1 (MQ=255)                                                                                                                                    aTTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAAt                                                                                                                   <  1:155613/240‑1 (MQ=255)                                                                                                                                    aTTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGc                                                                                                         >  1:868576/1‑250 (MQ=255)                                                                                                                                    aTTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGc                                                                                                         >  1:205512/1‑250 (MQ=255)                                                                                                                                              gaTTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCa                                                                                              >  1:58342/1‑251 (MQ=255)                                                                                                                                               aTTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATc                                                                                               <  1:673408/249‑1 (MQ=255)                                                                                                                                               aTTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAg                                                                                             <  1:293125/251‑1 (MQ=255)                                                                                                                                                ttCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGcc                                                                                                                                                                                                                 <  1:173507/134‑1 (MQ=255)                                                                                                                                                   atatTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGc                                                                                                                                                                 >  1:14024/1‑179 (MQ=255)                                                                                                                                                   atatTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCg                                                                                         >  1:555162/1‑251 (MQ=255)                                                                                                                                                   atatTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCg                                                                                         >  1:15091/1‑251 (MQ=255)                                                                                                                                                        attCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAAc                                                                                    <  1:819720/251‑1 (MQ=255)                                                                                                                                                                    tgtCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGc                                                                        <  1:464344/251‑1 (MQ=255)                                                                                                                                                                     gtCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCg                                                                       >  1:818858/1‑251 (MQ=255)                                                                                                                                                                     gtCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGGAAGCg                                                                       >  1:770197/1‑251 (MQ=255)                                                                                                                                                                      tCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCg                                                                       <  1:768208/250‑1 (MQ=255)                                                                                                                                                                      tCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCg                                                                       <  1:25819/250‑1 (MQ=255)                                                                                                                                                                       cTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCAGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCggg                                                                     >  1:564581/1‑250 (MQ=255)                                                                                                                                                                         tAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTc                                                                                                                                    >  1:696893/1‑186 (MQ=255)                                                                                                                                                                              aataaGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATg                                                                                                                                            <  1:474490/173‑1 (MQ=255)                                                                                                                                                                                   gTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGcc                                                                                                        <  1:231117/204‑1 (MQ=255)                                                                                                                                                                                        aTTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAAc                                                                                                                                                           <  1:689851/148‑1 (MQ=255)                                                                                                                                                                                             aactaaAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAAcgg                                                                                  >  1:658888/1‑216 (MQ=255)                                                                                                                                                                                             aactaaAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGGCAAAGCCGTCATGaaaaagg                                               >  1:796712/1‑247 (MQ=255)                                                                                                                                                                                             aactaaAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGGTAGCGGTGTCAATGCCGTCATGaaaaagg                                               >  1:310541/1‑247 (MQ=255)                                                                                                                                                                                               ctaaAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGGCGTCATGaaaaggtgg                                             >  1:649459/1‑245 (MQ=255)                                                                                                                                                                                               ctaaAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGGAAGCGGGGGCAATGccggcatggaaaggggg                                             >  1:83038/1‑236 (MQ=255)                                                                                                                                                                                               ctaaAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGCGTACGCGGGGTCAATACCGTCCTGaaaaagcgg                                             >  1:315432/1‑245 (MQ=255)                                                                                                                                                                                                       cTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGaaa                                     <  1:18630/251‑1 (MQ=255)                                                                                                                                                                                                          cGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGcc                                                            <  1:311065/225‑1 (MQ=255)                                                                                                                                                                                                          cGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGGCAATGCCGGCATGaaaaagcgggatgggaaaga                                  >  1:731724/1‑234 (MQ=255)                                                                                                                                                                                                          cGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGGCAAAGCCCGCATGAAAAAGCGGGTTgggaaaga                                  >  1:578588/1‑245 (MQ=255)                                                                                                                                                                                                          cGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGGAAGCGGTGTCAATg                                                              >  1:625973/1‑223 (MQ=255)                                                                                                                                                                                                            cTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATgatc                                <  1:268844/251‑1 (MQ=255)                                                                                                                                                                                                            cTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATgatc                                <  1:309317/251‑1 (MQ=255)                                                                                                                                                                                                                    gCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGGAAGCGGGGGCCATGCCGTCATGAAAAAGCGGAGTGGGAATgatcagatagaa                        >  1:597404/1‑247 (MQ=255)                                                                                                                                                                                                                               tGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGCGTAAGCGGGGTCCATGCCGTCATGAAAAGGCGGGTTGGGAAAgaacagatc                           >  1:460007/1‑237 (MQ=255)                                                                                                                                                                                                                                   cacTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATc                                                                                               <  1:588644/165‑1 (MQ=255)                                                                                                                                                                                                                                   cacTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAAc                                                                                    >  1:559313/1‑176 (MQ=255)                                                                                                                                                                                                                                   cacTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGTCAATGCCGGCATGAAAAAGCGGGGTGGAAATgatc                                >  1:719455/1‑228 (MQ=255)                                                                                                                                                                                                                                   cacTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGGCAATGCCGGCCTGGAAAGGGGGGTTGGGAATGATTAGATCGAATTTggcggcaaaatg         >  1:675427/1‑240 (MQ=255)                                                                                                                                                                                                                                   cacTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGCACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGAGGGGGCCATGCCCTCCTGAAAAAGCGGATTTGAAATGATTAGATCGAATTTTCCGTCCATATg         >  1:623715/1‑251 (MQ=255)                                                                                                                                                                                                                                     cTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATATggg       <  1:829198/251‑1 (MQ=255)                                                                                                                                                                                                                                     cTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATATggg       <  1:66558/251‑1 (MQ=255)                                                                                                                                                                                                                                     cTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGTCAATGCCCTCACTAAAAGGCGGATTGGAAATGATCAGATCGAATTTTCCGTCAATATggg       >  1:702141/1‑251 (MQ=255)                                                                                                                                                                                                                                       tAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAAcggc                                                                                 >  1:674699/1‑175 (MQ=255)                                                                                                                                                                                                                                       tAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATgatc                                <  1:31835/224‑1 (MQ=255)                                                                                                                                                                                                                                       tAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATATGGGaa     <  1:729665/251‑1 (MQ=255)                                                                                                                                                                                                                                       tAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATATGGGaa     <  1:800969/251‑1 (MQ=255)                                                                                                                                                                                                                                          cTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATATGGGAAAac  <  1:739006/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     ACACCCTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTTGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTTCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGGCGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCCGTGCACCGAATAAACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCATCATGAAACGGCGGGTTGGAAATGATCAGATCAAATTTGCCGTCAATATGGGAAAAC  >  1/1986073‑1986556

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