Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	882,742	A→G	intergenic (‑12423/‑14239)	unknown ← / → unknown	tRNA‑Gly‑GCC/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	882,742	0	A	G	100.0%	145.8 / NA	50	intergenic (‑12423/‑14239)	unknown/unknown	tRNA‑Gly‑GCC/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/23);  total (27/23)

TTGGAAATGGTAATCAAAAAAGATGACGGAAAAGAAATTGTTGTAGTGCAAAAAGCCGAACCTAACTTAGTGCCGGGTGCGCGGGTACGTATCGTCGGCGGTTCACGTTTGAACGTTTCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATACTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGCCTATTCCAAAGAATATGCACAACTTGAGGATGTTGTGAAAAGTTTTAACCGCTGGAATCAGCTCAATTCCAATCTGGAAGAA  >  1/882524‑882982                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 ttGGAAATGGTAATCAAAAAAGATGACGGAAAAGAAATTGTTGTAGTGCAAAAAGCCGAACCTAACTTAGTGCCGGGTGCGCGGGTACGTATCGTCGGCGGTTCACGTTTGAACGTTTCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTcc                                                                                                                                                                                                                  <  1:383421/251‑1 (MQ=255)                                   aaaTTGTTGTAGTGCAAAAAGCCGAACCTAACTTAGTGCCGGGTGCGCGGGTACGTATCGTCGGCGGTTCACGTTTGAATGTTTCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCtatta                                                                                                                                                                                  <  1:432277/249‑1 (MQ=255)                                   aaaTTGTTGTAGTGCAAAAAGCCGAACCTAACTTAGTGCCGGGTGCGCGGGTACGTATCGTCGGCGGTTCACGTTTGAACGTTTCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCtatta                                                                                                                                                                                  <  1:847828/249‑1 (MQ=255)                                    aaTTGTTGTAGTGCAAAAAGCCGAACCTAACTTAGTGCCGGGTGCGCGGGTACGTATCGTCGGCGGTTCACGTTTGAACGTTTCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTcc                                                                                                                                                                                                                  >  1:869322/1‑216 (MQ=255)                                               tGCAAAAAGCCGAACCTAACTTAGTGCCGGGTGCGCGGGTACGTATCGTCGGCGGTTCACGTTTGAACGTTTCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGAt                                                                                                                                                                    <  1:113142/251‑1 (MQ=255)                                               tGCAAAAAGCCGAACCTAACTTAGTGCCGGGTGCGCGGGTACGTATCGTCGGCGGTTCACGTTTGAACGTTTCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGAt                                                                                                                                                                    <  1:496284/251‑1 (MQ=255)                                               tGCAAAAAGCCGAACCTAACTTAGTGCCGGGTGCGCGGGTACGTATCGTCGGCGGTTCACGTTTGAACGTTTCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGAt                                                                                                                                                                    <  1:301990/251‑1 (MQ=255)                                                            aCCTAACTTAGTGCCGGGTGCGCGGGTACGTATCGTCGGCGGTTCACGTTTGAACGTTTCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGaa                                                                                                                                                        >  1:502983/1‑250 (MQ=255)                                                                 aCTTAGTGCCGGGTGCGCGGGTACGTATCGTCGGCGGTTCACGTTTGAACGTTTCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTAt                                                                                                                                                  >  1:425716/1‑251 (MQ=255)                                                                   ttAGTGCCGGGTGCGCGGGTACGTATCGTCGGCGGTTCACGTTTGAACGTTTCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAAtt                                                                                                                                                                        <  1:295844/227‑1 (MQ=255)                                                                                cgcgGGTACGTATCGTCGGCGGTTCACGTTTGAACGTTTCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGc                                                                                                                                    >  1:19057/1‑250 (MQ=255)                                                                                     gTACGTATCGTCGGCGGTTCACGTTTGAACGTTTCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGtt                                                                                                                              >  1:170925/1‑251 (MQ=255)                                                                                        cGTATCGTCGGCGGTTCACGTTTGAACGTTTCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAAcc                                                                                                                                                                                          <  1:863108/188‑1 (MQ=255)                                                                                                        cACGTTTGAACGTTTCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGAtt                                                                                                           <  1:561498/251‑1 (MQ=255)                                                                                                          cGTTTGAACGTTTCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCtatta                                                                                                                                                                                  >  1:587623/1‑178 (MQ=255)                                                                                                            tttGAACGTTTCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAg                                                                                                                                                          >  1:99818/1‑200 (MQ=255)                                                                                                            tttGAACGTGTCTGCACTTTAAGACTACCTATTCCTAAAAGGCTTTTGGCCAAATCCGGAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATGTTTTTATTTCCTTTTCTATTTTCCTTCCGAATGTCAAATAATGAAACCCTCTATTATCaa                                                                                                                                                                              <  1:700165/180‑1 (MQ=255)                                                                                                                    tttCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCt                                                                                                                                                                                                                 >  1:285279/1‑137 (MQ=255)                                                                                                                    tttCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAaataa                                                                                                                                                                                                 >  1:77479/1‑153 (MQ=255)                                                                                                                    tttCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGa                                                                                                >  1:232685/1‑250 (MQ=255)                                                                                                                            cTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAAt                                                                                                                                                                             <  1:590289/165‑1 (MQ=255)                                                                                                                               tAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGc                                                                                    <  1:29340/251‑1 (MQ=255)                                                                                                                                             ccTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGa                                                                                                                                                <  1:584475/177‑1 (MQ=255)                                                                                                                                              cTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGCCTATTCCAAAGAata                                                                     >  1:634927/1‑251 (MQ=255)                                                                                                                                                    ggCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGCCTATTCCAAAGAATATGcaca                                                               >  1:627523/1‑251 (MQ=255)                                                                                                                                                     gCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGCCTATTCCAAAGAATATGcac                                                                <  1:820230/249‑1 (MQ=255)                                                                                                                                                     gCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGCCTATTCCAAAGAATATGCACaa                                                              >  1:459562/1‑251 (MQ=255)                                                                                                                                                             cGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGCCTATTCCAAAGAATATGCACAACTTGAGGa                                                      <  1:111689/251‑1 (MQ=255)                                                                                                                                                                           cTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATg                                                                                                                                                                                               >  1:135582/1‑100 (MQ=255)                                                                                                                                                                             tttttACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGCCTATTCCAAAGAATATGCACAACTTGAGGATGTTGTGAAAAGtttt                                      >  1:136476/1‑251 (MQ=255)                                                                                                                                                                                     gCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAAtt                                                                                                                                                                        <  1:153810/113‑1 (MQ=255)                                                                                                                                                                                     gCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGCCTATTCCAAAGAATATGCACAACTTGAGGATGTTGTGAAAAGTTTTAACCg                                 >  1:365588/1‑248 (MQ=255)                                                                                                                                                                                                       ttttCATTCTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAag                                                                                                                                              >  1:445212/1‑121 (MQ=255)                                                                                                                                                                                                               cTGTGCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGa                                                                                                             >  1:222138/1‑146 (MQ=255)                                                                                                                                                                                                                 gtgCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAAccc                                                                                                                                                                                         >  1:3241/1‑68 (MQ=255)                                                                                                                                                                                                                 gtgCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCctct                                                                                                                                                                                      <  1:223036/71‑1 (MQ=255)                                                                                                                                                                                                                 gtgCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAAt                                                                                                                                                                             <  1:412642/80‑1 (MQ=255)                                                                                                                                                                                                                 gtgCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAg                                                                                                                                                          <  1:400481/99‑1 (MQ=255)                                                                                                                                                                                                                 gtgCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGa                                                                                                                        <  1:373002/133‑1 (MQ=255)                                                                                                                                                                                                                 gtgCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGAc                                                                                               >  1:196557/1‑158 (MQ=255)                                                                                                                                                                                                                 gtgCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGCCTATTCCAAAGAATATGCACAACTTGAGGa                                                      <  1:353481/199‑1 (MQ=255)                                                                                                                                                                                                                 gtgCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGCCTATTCCAAAGAATATGCACAACTTGAGGa                                                      >  1:351228/1‑199 (MQ=255)                                                                                                                                                                                                                 gtgCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGCCTATTCCAAAGAATATGCACAACTTGAGGATGTTGTGAAAAGTTTTAACCGCTGGaa                           <  1:727332/226‑1 (MQ=255)                                                                                                                                                                                                                 gtgCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGCCTATTCCAAAGAATATGCACAACTTGAGGATGTTGTGAAAAGTTTTAACCGCTGGAATCAGCTCAATTc               <  1:527973/238‑1 (MQ=255)                                                                                                                                                                                                                 gtgCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGCCTATTCCAAAGAATATGCACAACTTGAGGATGTTGTGAAAAGTTTTAACCGCTGGAATCAGCTCAATTCCAATCTGgaagaa  >  1:471632/1‑251 (MQ=255)                                                                                                                                                                                                                 gtgCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGCCTATTCCAAAGAATATGCACAACTTGAGGATGTTGTGAAAAGTTTTAACCGCTGGAATCAGCTCAATTCCAATCTGgaagaa  >  1:604174/1‑251 (MQ=255)                                                                                                                                                                                                                 gtgCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGCCTATTCCAAAGAATATGCACAACTTGAGGATGTTGTGAAAAGTTTTAACCGCTGGAATCAGCTCAATTCCAATCTGgaagaa  >  1:627488/1‑251 (MQ=255)                                                                                                                                                                                                                 gtgCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGCCTATTCCAAAGAATATGCACAACTTGAGGATGTTGTGAAAAGTTTTAACCGCTGGAATCAGCTCAATTCCAATCTGgaagaa  >  1:675337/1‑251 (MQ=255)                                                                                                                                                                                                                 gtgCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGCCTATTCCAAAGAATATGCACAACTTGAGGATGTTGTGAAAAGTTTTAACCGCTGGAATCAGCTCAATTCCAATCTGgaagaa  >  1:685836/1‑251 (MQ=255)                                                                                                                                                                                                                 gtgCTAGAATGCTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGCCTATTCCAAAGAATATGCACAACTTGAGGATGTTGTGAAAAGTTTTAACCGCTGGAATCAGCTCAATTCCAATCTGgaaga   >  1:360373/1‑250 (MQ=255)                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 TTGGAAATGGTAATCAAAAAAGATGACGGAAAAGAAATTGTTGTAGTGCAAAAAGCCGAACCTAACTTAGTGCCGGGTGCGCGGGTACGTATCGTCGGCGGTTCACGTTTGAACGTTTCTGCACTTTAAGTCTACCTATTCCTAAAAGGCTTTTGGCGAAATCCGAAAGCCTTTTTTACAGCCTAAACTTCAACAACATTTTCATTCTGTGCTAGAATACTCCGATTTTTTTATTTCCTTTTCTATTTTCCTTCTGAAGGTCAAATAATGAAACCCTCTATTATCAATAAATTAGATAGCTTAAAAGAACGTTATGAAGAGCTGGAAGCCTTGTTAGGCGATGCTTCCGTGATTAGCGATCAGGACAAATTCAGAGCCTATTCCAAAGAATATGCACAACTTGAGGATGTTGTGAAAAGTTTTAACCGCTGGAATCAGCTCAATTCCAATCTGGAAGAA  >  1/882524‑882982

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