Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	306,497	T→C	intergenic (‑39137/+38108)	unknown ← / ← unknown	SSU rRNA ## 16S rRNA, small subunit ribosomal RNA/tRNA‑Arg‑ACG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	306,497	0	T	C	100.0%	227.8 / NA	79	intergenic (‑39137/+38108)	unknown/unknown	SSU rRNA ## 16S rRNA, small subunit ribosomal RNA/tRNA‑Arg‑ACG
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/46);  total (33/46)

ACGCTGTCGATTTGTTGCACCAGCTCCAACCATGCATGACGGGTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGTCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGCATGATACAACTATCATGCCTAAAATCAAAAGAGAAAGATTAAAGAATAAAGCGGCTTAAATCTTCGTTTTCGAT  >  1/306253‑306736                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                aCGCTGTCGATTTGTTGCACCAGCTCCAACCATGCATGACGGGTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTc                                                                                                                                                                                                                                           <  1:640087/251‑1 (MQ=255)    cTGTCGATTTGTTGCACCAGCTCCAACCATGCATGACGGGTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCAc                                                                                                                                                                                                                                         >  1:548962/1‑250 (MQ=255)      gTCGATTTGTTGCACCAGCTCCAACCATGCATGACGGGTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCg                                                                                                                                                                                                                                       >  1:832539/1‑250 (MQ=255)       tCGATTTGTTGCACCAGCTCCAACCATGCATGACGGGTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGc                                                                                                                                                                                                                                      <  1:541762/250‑1 (MQ=255)           tttGTTGCACCAGCTCCAACCATGCATGACGGGTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATAc                                                                                                                                                                                                                                 <  1:210228/251‑1 (MQ=255)             tgttgCACCAGCTCCAACCATGCATGACGGGTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACaa                                                                                                                                                                                                                               >  1:829875/1‑251 (MQ=255)                tgCACCAGCTCCAACCATGCATGACGGGTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATg                                                                                                                                                                                                                            <  1:116208/251‑1 (MQ=255)                  cACCAGCTCCAACCATGCATGACGGGTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCg                                                                                                                                                                                                                          >  1:710010/1‑251 (MQ=255)                      aGCTCCAACCATGCATGACGGGTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACaaa                                                                                                                                                                                                                              <  1:351543/243‑1 (MQ=255)                       gCTCCAACCATGCATGACGGGTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATTCGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACaaa                                                                                                                                                                                                                              <  1:433542/242‑1 (MQ=255)                       gCTCCAACCATGCATGACGGGTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTa                                                                                                                                                                                                                     >  1:844778/1‑251 (MQ=255)                        cTCCAACCATGCATGACGGGTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACaaa                                                                                                                                                                                                                              >  1:868431/1‑241 (MQ=255)                           cAACCATGCATGACGGGTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAatat                                                                                                                                                                                                                 <  1:344658/251‑1 (MQ=255)                                atgcatgACGGGTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGa                                                                                                                                                                                                            >  1:871793/1‑251 (MQ=255)                                        cGGGTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCttt                                                                                                                                                                                                    <  1:422389/251‑1 (MQ=255)                                          ggTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATg                                                                                                                                                                                                                            <  1:263288/225‑1 (MQ=255)                                          ggTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCg                                                                                                                                                                                                  >  1:47621/1‑251 (MQ=255)                                           gTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCg                                                                                                                                                                                                  <  1:264001/250‑1 (MQ=255)                                              taATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAatat                                                                                                                                                                                                                 <  1:651232/232‑1 (MQ=255)                                              taATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTg                                                                                                                                                                                                >  1:131498/1‑249 (MQ=255)                                              taATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTg                                                                                                                                                                                                >  1:383051/1‑249 (MQ=255)                                              taATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTg                                                                                                                                                                                                >  1:151384/1‑249 (MQ=255)                                                   aTCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATAc                                                                                                                                                                                                                                 <  1:661538/211‑1 (MQ=255)                                                   aTCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAatat                                                                                                                                                                                                                 >  1:684217/1‑227 (MQ=255)                                                   aTCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGa                                                                                                                                                                                                            <  1:697851/232‑1 (MQ=255)                                                   aTCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGt                                                                                                                                                                                         <  1:503483/251‑1 (MQ=255)                                                   aTCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGt                                                                                                                                                                                         >  1:309050/1‑251 (MQ=255)                                                         tACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATcc                                                                                                                                                                                    <  1:599244/250‑1 (MQ=255)                                                         tACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATc                                                                                                                                                                                     <  1:253766/249‑1 (MQ=255)                                                                           cTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCat                                                                                                                                                                   <  1:567333/249‑1 (MQ=255)                                                                            ttCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATAtt                                                                                                                                                                <  1:338260/251‑1 (MQ=255)                                                                                ccGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAatat                                                                                                                                                                                                                 <  1:4857/198‑1 (MQ=255)                                                                                 cGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATAc                                                                                                                                                                                                                                 >  1:690661/1‑181 (MQ=255)                                                                                 cGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGccc                                                                                                                                                           <  1:10057/251‑1 (MQ=255)                                                                                  gCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCa                                                                                                                                                          >  1:509409/1‑251 (MQ=255)                                                                                        gCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCt                                                                                                                                                                                                      >  1:689875/1‑201 (MQ=255)                                                                                        gCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTgg                                                                                                                                                    >  1:38970/1‑251 (MQ=255)                                                                                        gCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTgg                                                                                                                                                    >  1:472487/1‑251 (MQ=255)                                                                                         cATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGc                                                                                                                                                   <  1:808559/251‑1 (MQ=255)                                                                                          aTTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCa                                                                                                                                                  <  1:408988/251‑1 (MQ=255)                                                                                           ttGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCaa                                                                                                                                                 <  1:6085/251‑1 (MQ=255)                                                                                           ttGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCaa                                                                                                                                                 <  1:849226/251‑1 (MQ=255)                                                                                             gAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAgg                                                                                                                                               >  1:581201/1‑251 (MQ=255)                                                                                                     cGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGaaa                                                                                                                                       <  1:615521/251‑1 (MQ=255)                                                                                                     cGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGaaa                                                                                                                                       <  1:424514/251‑1 (MQ=255)                                                                                                               gAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCgg                                                                                                                                >  1:749044/1‑248 (MQ=255)                                                                                                                  tttGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGaa                                                                                                                          <  1:712946/251‑1 (MQ=255)                                                                                                                        aTCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAg                                                                                                                    <  1:702372/251‑1 (MQ=255)                                                                                                                        aTCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAg                                                                                                                    <  1:860249/251‑1 (MQ=255)                                                                                                                                    tttGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATg                                                                                                                                                                                                                            <  1:602893/135‑1 (MQ=255)                                                                                                                                           atatAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGc                                                                                                                                                   >  1:89597/1‑201 (MQ=255)                                                                                                                                                      gCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATg                                                                                                                                                                                                                            >  1:815781/1‑117 (MQ=255)                                                                                                                                                                 agaTACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCg                                                                                                                                                                                          >  1:508117/1‑140 (MQ=255)                                                                                                                                                                 agaTACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGc                                                                            <  1:723300/250‑1 (MQ=255)                                                                                                                                                                 agaTACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGc                                                                            <  1:830297/250‑1 (MQ=255)                                                                                                                                                                   aTACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGCATg                                                                         <  1:38751/251‑1 (MQ=255)                                                                                                                                                                    tACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAatat                                                                                                                                                                                                        <  1:786430/123‑1 (MQ=255)                                                                                                                                                                     aCTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGCATGAt                                                                       <  1:819749/251‑1 (MQ=255)                                                                                                                                                                          aaCTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGCATGATACAAc                                                                  <  1:14468/251‑1 (MQ=255)                                                                                                                                                                          aaCTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGCATGATACAAc                                                                  <  1:214966/251‑1 (MQ=255)                                                                                                                                                                            cTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAg                                                                              <  1:707554/237‑1 (MQ=255)                                                                                                                                                                            cTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGCATGATACAACTa                                                                <  1:527679/251‑1 (MQ=255)                                                                                                                                                                              ttCCGCCAAATACGGCTGTGGACCG‑ATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGCATGATACAACTATCa                                                             <  1:605382/251‑1 (MQ=255)                                                                                                                                                                                     aaaTACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAatat                                                                                                                                                                                                                 <  1:857543/97‑1 (MQ=255)                                                                                                                                                                                     aaaTACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGCATGATACAACTATCATGCCTa                                                       >  1:167163/1‑251 (MQ=255)                                                                                                                                                                                       aTACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTTAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGCATGATACAACTATCATGCCTaaa                                                     >  1:172918/1‑251 (MQ=255)                                                                                                                                                                                         aCGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACTCCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGCATGATACAACTATCATGCCTAAAAt                                                   <  1:796736/251‑1 (MQ=255)                                                                                                                                                                                                          tCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTTATGAATCTCTCCATCTATAAACAAAAGGCATGATACAACTATCATGCCTAAAATCAAAAGAGAAAGATTaa                                  >  1:774989/1‑251 (MQ=255)                                                                                                                                                                                                          tCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGCATGATACAACTa                                                                >  1:343679/1‑221 (MQ=255)                                                                                                                                                                                                           ccGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATg                                                                                                                                                                                                                            >  1:332465/1‑64 (MQ=255)                                                                                                                                                                                                                    gCACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGCATGATACAACTATCa                                                             <  1:565070/214‑1 (MQ=255)                                                                                                                                                                                                                     cACCAAGGGCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCatat                                                                                                                                                                 >  1:833314/1‑113 (MQ=255)                                                                                                                                                                                                                           gggCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGCATGATACAACTATCATGCCTAAAATCAAAAGAGAAAGATTAAAGAATAAAGCGGCTTaa                 >  1:844847/1‑251 (MQ=255)                                                                                                                                                                                                                            ggCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCta                                                                                        <  1:681627/179‑1 (MQ=255)                                                                                                                                                                                                                             gCTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGCATGATACAACTATCATGCCTAAAATCAAAAGAGAAAGATTAAAGAATAAAGCGGCTTAAAt               >  1:863221/1‑251 (MQ=255)                                                                                                                                                                                                                              cTCAAACACCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGCATGATACAACTATCATGCCTAAAATCAAAAGAGAAAGATTAAAGAATAAAGCGGCTTAAATc              <  1:648472/251‑1 (MQ=255)                                                                                                                                                                                                                                     acCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAg                                                                                                                    <  1:250309/142‑1 (MQ=255)                                                                                                                                                                                                                                     acCATATTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGCATGATACAACTATCATGCCTAAAATCAAAAGAGAAAGATTAAAGAATAAAGCGGCTTAAATCTTCGttt       >  1:658349/1‑251 (MQ=255)                                                                                                                                                                                                                                          atTCTGCGCTGCCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGCATGATACAACTATCATGCCTAAAATCAAAAGAGAAAGATTAAAGAATAAAGCGGCTTAAATCTTCGTTTTCGAt  >  1:741788/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                ACGCTGTCGATTTGTTGCACCAGCTCCAACCATGCATGACGGGTATAATTATCCAGTACGGCATCGATGTTCAGCTTCACCGCCTGAGCATTGAATGGTTCGCCGTCGCTGAATTTGACATCTTTGCGTAATTTGAAAATATAGGTTTTGCCATCCGGCGAGATACTCCAACTTTCCGCCAAATACGGCTGTGGACCGAATTCCGTATTCAGCACCAAGGGCTCAAACACCATATTCTGCGCTGTCATTTCACCGCTATACAAATGCGGATTAATATCACGAATATCTTTCGTGCTGGCGTAATCCAATTCGTTGTTTTGCATATTCGCCCAGCTTGGCAAGGCGGTGAAAAGTGCGGTCAGAATTAAAGTTATTTTTTTCATGAATCTCTCCATCTATAAACAAAAGGCATGATACAACTATCATGCCTAAAATCAAAAGAGAAAGATTAAAGAATAAAGCGGCTTAAATCTTCGTTTTCGAT  >  1/306253‑306736

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