Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,986,910	G→A	R134R (CGC→CGT)	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,910	0	G	A	98.6%	205.9 / ‑4.8	70	R134R (CGC→CGT)	unknown	16S rRNA (guanine(1207)‑N(2))‑methyltransferase (EC 2.1.1.172)
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (46/23);  total (46/24)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GTTGAATATAACTGCCGATCACTCCCGCGCCGCAGCCCACGTCCAACACCCCGCCTGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGTTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAATTTGTGGCGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAACTGAAAATGATACAGACTGCAACGGCGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACGCCGCTACGATTTTCGCCCACCAGCAAAACCTCCTGCCCGGCGGTGCATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGATGTTGTCTGT  >  1/1986672‑1987154                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     gTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCg                                                                                                                                                                                                                                          >  1:151337/1‑251 (MQ=255) gTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGAGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCg                                                                                                                                                                                                                                          >  1:510831/1‑251 (MQ=255)    gAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGAtt                                                                                                                                                                                                                                       >  1:658031/1‑251 (MQ=255)    gAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGAtt                                                                                                                                                                                                                                       >  1:180815/1‑251 (MQ=255)      atatAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCGGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGAtttt                                                                                                                                                                                                                                     >  1:615956/1‑251 (MQ=255)              gccaaTCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAatat                                                                                                                                                                                                                             >  1:791005/5‑251 (MQ=255)               ccGATCACTCCCGCGCCGCAGCCCACGTCCAGCACCTCGCCTGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGTTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAATTTGTGGCGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAACTGAAAATGATACAGACTGCAACGGCGGGCGGAGTCGATTTTGTGAatat                                                                                                                                                                                                                             <  1:11417/250‑1 (MQ=255)                 aaTCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCg                                                                                                                                                                                                                           >  1:393609/2‑250 (MQ=255)                    cacccctgcgcCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGt                                                                                                                                                                                                                       >  1:839009/8‑251 (MQ=255)                                gcAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCa                                                                                                                                                                                                           >  1:374079/1‑251 (MQ=255)                                   gCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAatat                                                                                                                                                                                                                             >  1:381684/1‑230 (MQ=255)                                         gTCCAACGCCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCAtttttt                                                                                                                                                                                                     <  1:751717/248‑1 (MQ=255)                                         gTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGa                                                                                                                                                                                                                                 <  1:726661/220‑1 (MQ=255)                                         gTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAata                                                                                                                                                                                                                              <  1:716407/223‑1 (MQ=255)                                         gTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTcgg                                                                                                                                                                                                  >  1:364997/1‑251 (MQ=255)                                         gTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTcgg                                                                                                                                                                                                  >  1:406534/1‑251 (MQ=255)                                         gTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTcgg                                                                                                                                                                                                  >  1:874023/1‑251 (MQ=255)                                               cacCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGTCGGaa                                                                                                                                                                                             >  1:276473/1‑250 (MQ=255)                                                   ccgcccgcGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAacgcac                                                                                                                                                                                        <  1:782685/251‑1 (MQ=255)                                                      cccgCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAacgcacacc                                                                                                                                                                                     <  1:235148/248‑4 (MQ=255)                                                      cccgCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAacgcacacc                                                                                                                                                                                     <  1:685270/248‑4 (MQ=255)                                                      cccgCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAacgcacac                                                                                                                                                                                      <  1:831075/247‑3 (MQ=255)                                                              gTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTaa                                                                                                                                                                                                                     >  1:109320/1‑211 (MQ=255)                                                              gTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGa                                                                                                                                                                              >  1:645805/1‑250 (MQ=255)                                                               ttttttAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGAtt                                                                                                                                                                            <  1:7695/251‑1 (MQ=255)                                                                tttttAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCAtttttt                                                                                                                                                                                                     <  1:391161/225‑1 (MQ=255)                                                                                cAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGTTAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTaaaa                                                                                                                                                           >  1:327340/1‑251 (MQ=255)                                                                                cAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTaaaa                                                                                                                                                           >  1:267767/1‑251 (MQ=255)                                                                                             ccgccATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCcgg                                                                                                                                               >  1:726689/1‑250 (MQ=255)                                                                                                  aTCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGa                                                                                                                                                                              >  1:366719/1‑214 (MQ=255)                                                                                                  aTCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCCCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGta                                                                                                                                          >  1:768600/1‑249 (MQ=255)                                                                                                    ccAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACAt                                                                                                                                       <  1:452063/251‑1 (MQ=255)                                                                                                    ccAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACAt                                                                                                                                       <  1:899841/251‑1 (MQ=255)                                                                                                    ccAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACAt                                                                                                                                       <  1:385221/251‑1 (MQ=255)                                                                                                                  cTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAg                                                                                                                            >  1:150425/1‑248 (MQ=255)                                                                                                                    gAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACAt                                                                                                                                       <  1:599727/235‑1 (MQ=255)                                                                                                                    gAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCa                                                                                                                       >  1:699636/1‑251 (MQ=255)                                                                                                                    gAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCa                                                                                                                       >  1:194344/1‑251 (MQ=255)                                                                                                                    gAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCa                                                                                                                       >  1:110471/1‑251 (MQ=255)                                                                                                                               gCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTgag                                                                                                                                 <  1:275620/230‑1 (MQ=255)                                                                                                                                  aCGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATgg                                                                                                         >  1:99561/1‑251 (MQ=255)                                                                                                                                   cGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGa                                                                                                                                                                              >  1:352346/1‑181 (MQ=255)                                                                                                                                   cGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATgg                                                                                                         >  1:845888/1‑250 (MQ=255)                                                                                                                                   cGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATgg                                                                                                         >  1:749738/1‑250 (MQ=255)                                                                                                                                   cGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATgg                                                                                                         >  1:212149/1‑250 (MQ=255)                                                                                                                                    gCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCa                                                                                                       >  1:169108/1‑251 (MQ=255)                                                                                                                                          atagTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAacgcgcaccactacg                                                                                                                                                                               >  1:217985/1‑162 (MQ=255)                                                                                                                                           tagTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTc                                                                                                    >  1:529807/1‑247 (MQ=255)                                                                                                                                           tagTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTc                                                                                                    >  1:783633/1‑247 (MQ=255)                                                                                                                                                 aaTCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTAttc                                                                                          <  1:244789/251‑1 (MQ=255)                                                                                                                                                                  gATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACAt                                                                                                                                       >  1:34984/1‑189 (MQ=255)                                                                                                                                                                  gATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAg                                                                         >  1:381789/1‑251 (MQ=255)                                                                                                                                                                        tAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATcac                                                                   <  1:193855/251‑1 (MQ=255)                                                                                                                                                                            aTTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATc                                                                <  1:6723/250‑1 (MQ=255)                                                                                                                                                                            aTTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCg                                                               <  1:688114/251‑1 (MQ=255)                                                                                                                                                                               tCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGgcg                                                            <  1:794608/251‑1 (MQ=255)                                                                                                                                                                               tCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGgcg                                                            >  1:55890/1‑251 (MQ=255)                                                                                                                                                                                          ttCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGaaa                                                 <  1:821885/251‑1 (MQ=255)                                                                                                                                                                                                                         gAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTgg                   >  1:53459/1‑250 (MQ=255)                                                                                                                                                                                                                           aaaTGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGaaa                                                                                                            >  1:772392/1‑159 (MQ=255)                                                                                                                                                                                                                             aTGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACAt                                                                                                                                       >  1:392336/1‑130 (MQ=255)                                                                                                                                                                                                                             aTGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTgagag                                                                                                                               >  1:105046/1‑138 (MQ=255)                                                                                                                                                                                                                             aTGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAggg                                         >  1:105341/1‑224 (MQ=255)                                                                                                                                                                                                                             aTGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAAt                                   <  1:520717/230‑1 (MQ=255)                                                                                                                                                                                                                               gATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTAttc                                                                                          >  1:311964/1‑173 (MQ=255)                                                                                                                                                                                                                                aTACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTggt           >  1:414259/1‑249 (MQ=255)                                                                                                                                                                                                                                       cTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATcaca                                                                  <  1:807440/189‑1 (MQ=255)                                                                                                                                                                                                                                       cTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACttt                       >  1:233494/1‑232 (MQ=255)                                                                                                                                                                                                                                       cTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTggtg          <  1:741319/245‑4 (MQ=255)                                                                                                                                                                                                                                         gCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTcggt  >  1:369840/1‑248 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     GTTGAATATAACTGCCGATCACTCCCGCGCCGCAGCCCACGTCCAACACCCCGCCTGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGTTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAATTTGTGGCGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAACTGAAAATGATACAGACTGCAACGGCGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACGCCGCTACGATTTTCGCCCACCAGCAAAACCTCCTGCCCGGCGGTGCATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGATGTTGTCTGT  >  1/1986672‑1987154

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