Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,093,886	C→G	intergenic (‑34333/‑38217)	unknown ← / → unknown	tRNA‑Ser‑GGA/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	1,093,886	0	C	G	97.9%	126.4 / ‑6.2	48	intergenic (‑34333/‑38217)	unknown/unknown	tRNA‑Ser‑GGA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (18/29);  total (19/29)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACATTCTGCAACGCTTGATAAATAAGCTGTGCCTGTTCGCCGGAGAGGGCGGTTTTATTTTCGTCGGAAATAACCACCGCACTTTGTTTACCTACAGCGGCAAGCTGCATTTGGTAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCCAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGTAAGCGGAATTCAATTCCCCCACCAATTGATTTAAGCGATCGGAGGCGTAGCGTTGCCTGTCCGTAAAATCCGGCGTGAAAATGGTGTCATCACGTTTCATTT  >  1/1093650‑1094125                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                aCATTCTGCAACGCTTGATAAATAAGCTGTGCCTGTTCGCCGGAGAGGGCGGTTTTATTTTCGTCGGAAATAACCACCGCACTTTGTTTACCTACAGCGGCAAGCTGCATTTGGTAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGc                                                                                                                                                                                                                                   <  1:401289/251‑1 (MQ=255)      cTGCAACGCTTGATAAATAAGCTGTGCCTGTTCGCCGGAGAGGGCGGTTTTATTTTCGTCGGAAATAACCACCGCACTTTGTTTACCTACAGCGGCAAGCTGCATTTGGTAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGaa                                                                                                                                                                                                                               >  1:170217/1‑250 (MQ=255)                  ataaataaGCTGTGCCTGTTCGCCGGAGAGGGCGGTTTTATTTTCGTCGGAAATAACCACCGCACTTTGTTTACCTACAGCGGCAAGCTGCATTTGGCAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGc                                                                                                                                                                                                                   >  1:888223/1‑250 (MQ=255)                   taaataaGCTGTGCCTGTTCGCCGGAGAGGGCGGTTTTATTTTCGTCGGAAATAACCACCGCACTTTGTTTACCTACAGCGGCAAGCTGCATTTGGTAAGTGCCTTTTTCCAAGTCCGGTAGCTGTATTCCTACACGGCGCCATCCTTCTTGCGTTAAAGTTTTGTACTCTAACCCGCGAATGCCTCGGCTGCGGGTTTTTTCTTGAATTTCAAAGCCCAGTTTCGGCAGGGCGGAACCAAAACCCTGccg                                                                                                                                                                                                                 >  1:339745/1‑250 (MQ=255)                     aataaGCTGTGCCTGTTCGCCGGAGAGGGCGGTTTTATTTTCGTCGGAAATAACCACCGCACTTTGTTTACCTACAGCGGCAAGCTGCATTTGGTAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAgg                                                                                                                                                                                                               >  1:765338/1‑251 (MQ=255)                              tgCCTGTTCGCCGGAGAGGGCGGTTTTATTTTCGTCGGAAATAACCACCGCACTTTGTTTACCTACAGCGGCAAGCTGCATTTGGTAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCa                                                                                                                                                                                                       <  1:711837/250‑1 (MQ=255)                                ccTGTTCGCCGGAGAGGGCGGTTTTATTTTCGTCGGAAATAACCACCGCACTTTGTTTACCTACAGCGGCAAGCTGCATTTGGTAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAat                                                                                                                                                                                                    <  1:208348/251‑1 (MQ=255)                                 cTGTTCGCCGGGGAGGGCGTTTTTATTTTCTTCGG‑AATAACCACCGCACTTTGTTTACCTACAGCGGCAAGCTGCATTTGGTAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAatga                                                                                                                                                                                                  <  1:89478/251‑1 (MQ=255)                                 cTGTTCGCCGGAGAGGGCGGTTTTATTTTCGTCGGAAATAACCACCGCACTTTGTTTACCTACAGCGGCAAGCTGCATTTGGTAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAatg                                                                                                                                                                                                   <  1:131386/251‑1 (MQ=255)                                  tGTTCGCCGGAGAGGGCGGTTTTATTTTCGTCGGAAATAACCACCGCACTTTGTTTACCTACAGCGGCAAGCTGCATTTGGTAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAatga                                                                                                                                                                                                  <  1:207638/251‑1 (MQ=255)                                    ttCGCCGGAGAGGGCGGTTTTATTTTCGTCGGAAATAACCACCGCACTTTGTTTACCTACAGCGGCAAGCTGCATTTGGTAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAatgatg                                                                                                                                                                                                <  1:744936/251‑1 (MQ=255)                                     tCGCCGGAGAGGGCGGTTTTATTTTCGTCGGAAATAACCACCGCACTTTGTTTACCTACAGCGGCAAGCTGCATTTGGTAAGTGCC‑TTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGt                                                                                                                                                                                               >  1:206510/1‑250 (MQ=255)                                        ccGGAGAGGGCGGTTTTATTTTCGTCGGAAATAACCACCGCACTTTGTTTACCTACAGCGGCAAGCTGCATTTGGTAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGt                                                                                                                                                                                               >  1:326312/1‑248 (MQ=255)                                                          ttttCGTCGGAAATAACCACCGCACTTTGTTTACCTACAGCGGCAAGCTGCATTTGGTAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGt                                                                                                                                                                          <  1:293221/251‑1 (MQ=255)                                                           tttCGTCGGAAATAACCACCGCACTTTGTTTACCTACAGCGGCAAGCTGCATTTGGTAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTc                                                                                                                                                                         <  1:757088/251‑1 (MQ=255)                                                                                    ttGTTTACCTACAGCGGCAAGCTGCATTTGGTAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGc                                                                                                                                                                                                                      >  1:747051/1‑181 (MQ=255)                                                                                        ttACCTACAGCGGCAAGCTGCATTTGGTAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTgagg                                                                                                                                            <  1:563029/251‑1 (MQ=255)                                                                                                           gCATTTGGTAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGc                                                                                                                                                                                                                      >  1:860925/1‑158 (MQ=255)                                                                                                           gCATTTGGTAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTgaggag                                                                                                                                          <  1:651662/234‑1 (MQ=255)                                                                                                           gCATTTGGTAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGAt                                                                                                                          >  1:62654/1‑250 (MQ=255)                                                                                                                         ccTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTa                                                                                                                              <  1:99777/232‑1 (MQ=255)                                                                                                                              ttCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGTaa                                                                                                      >  1:882258/1‑251 (MQ=255)                                                                                                                                             cTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGTAAGCGGaattcaat                                                                                          <  1:834461/248‑1 (MQ=255)                                                                                                                                                cATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGTAAGCGGAATTCAATTccccc                                                                                    <  1:63932/251‑1 (MQ=255)                                                                                                                                                  ttCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGc                                                                                                                                                                                                              <  1:704822/127‑1 (MQ=255)                                                                                                                                                  ttCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAatgatg                                                                                                                                                                                                >  1:867469/1‑141 (MQ=255)                                                                                                                                                  ttCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGc                                                                                                                                                                                   <  1:806969/154‑1 (MQ=255)                                                                                                                                                  ttCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCg                                                                                                                                                                                 <  1:226304/156‑1 (MQ=255)                                                                                                                                                  ttCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATc                                                                                                                         >  1:381407/1‑212 (MQ=255)                                                                                                                                                  ttCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGt                                                                                                        <  1:309955/229‑1 (MQ=255)                                                                                                                                                  ttCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGTAAGCGGaattc                                                                                             <  1:189550/240‑1 (MQ=255)                                                                                                                                                  ttCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGTAAGCGGAATTCAATTCCCcca                                                                                   >  1:602736/1‑250 (MQ=255)                                                                                                                                                  ttCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGTAAGCGGAATTCAATTCCCcca                                                                                   >  1:243987/1‑250 (MQ=255)                                                                                                                                                  ttCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGTAAGCGGAATTCAATTCCCcca                                                                                   >  1:236546/1‑250 (MQ=255)                                                                                                                                                       acacGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGTAAGCGGAATTCAATTCCCCCACCAAtt                                                                             <  1:584807/251‑1 (MQ=255)                                                                                                                                                                  ttcttcTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGa                                                                                                                           <  1:281850/194‑1 (MQ=255)                                                                                                                                                                  ttcttcTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGTAAGCGGAATTCAATTCCCCCACCAATTGATTTAAGCGa                                                                  >  1:170551/1‑251 (MQ=255)                                                                                                                                                                    cttcttGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGTAAGCGGAATTCAATTCCCCCACCAATTGATTTAAGCGAt                                                                 >  1:681569/1‑250 (MQ=255)                                                                                                                                                                           tgATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGt                                                                                                        <  1:781485/204‑1 (MQ=255)                                                                                                                                                                            gATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGc                                                                                                                                                                            <  1:528549/135‑1 (MQ=255)                                                                                                                                                                            gATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGTAAGCGGAATTCAATTCCCCCACCAATTGATTTAAGCGATCGGAGGCGt                                                        >  1:233925/1‑251 (MQ=255)                                                                                                                                                                              tAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGTAAGCGGAATTCAATTCCCCCACCAATTGATTTAAGCGATCGGAGGCGTAg                                                      <  1:605359/251‑1 (MQ=255)                                                                                                                                                                                  ggTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATc                                                                                                                         >  1:698633/1‑180 (MQ=255)                                                                                                                                                                                     ttGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGTAAGCGGAATTCAATTCCCCCACCAATTGATTTAAGCGATCGGAGGCGTAGCGTTGcc                                               <  1:744778/251‑1 (MQ=255)                                                                                                                                                                                       gTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGTAAGCGGAATTCAATTCCCCCACCAATTGATTTAAg                                                                     <  1:38056/227‑1 (MQ=255)                                                                                                                                                                                                    cgAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTg                                                                                                                                               <  1:450762/140‑1 (MQ=255)                                                                                                                                                                                                                cTGCCGGTTTCTTCCTGAATTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGTAAGCGGAATTCAATTCCCCCACCAATTGATTTAAGCGATCGGAGGCGTAGCGTTGCCTGTCCGTAAAATCCGGCGTGAAAATgg                    <  1:166611/251‑1 (MQ=255)                                                                                                                                                                                                                                  aTTTCAAAGCCGAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGTAAGCGGAATTCAATTCCCCCACCAATTGATTTAAGCGATCGGAGGCGTAGCGTTGCCTGTCCGTAAAATCCGGCGTGAAAATGGTGTCATCACGTTTCAttt  <  1:50587/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                ACATTCTGCAACGCTTGATAAATAAGCTGTGCCTGTTCGCCGGAGAGGGCGGTTTTATTTTCGTCGGAAATAACCACCGCACTTTGTTTACCTACAGCGGCAAGCTGCATTTGGTAAGTGCCTTTTTCCAAGTCCGGCAGCTGCATTCCTACACGGCGCCATTCTTCTTGTGATAAAGGTTTGTACTCTAACTCGCGAATGCCTCGGCTGCCGGTTTCTTCCTGAATTTCAAAGCCCAGTTTCGGCAGAGCGGAACCTAAACGCTGCCAGGCGTTATCAAATGATGTCTGCGGTTCTGCCGTTGGCGTCGGTGGCGATGGCGGACTGAATTGGTGAGGAGGCTGCGCCGTTAAGATCCTGCTGTTGTTTTTGGTAAGCGGAATTCAATTCCCCCACCAATTGATTTAAGCGATCGGAGGCGTAGCGTTGCCTGTCCGTAAAATCCGGCGTGAAAATGGTGTCATCACGTTTCATTT  >  1/1093650‑1094125

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