Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	566,304	G→T	intergenic (+202887/‑5715)	unknown → / → unknown	tRNA‑SeC(p)‑TCA/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	566,304	0	G	T	100.0%	126.0 / NA	43	intergenic (+202887/‑5715)	unknown/unknown	tRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (17/26);  total (17/26)

ACCTAAGAAATCTTTCGGCAGGGCGTTTGCAAGTAAAGTACCACCACCGGCCATGTAGGCGGCTGATAATGCATATAACAAAATGAAAAGTGAGAGTGTTGAAATAATACGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCGCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGAAGGTTTTATTTTTTAATAACATACATTAATCCTAAATAAAGTAGGCTGGTATTATCGCATAAAGTCGTCGGAGATTAAGTGCTAATGCGCCGATTTTGTTATAAAGTGCGGTCTTTTTTGAGAGCGTTTTTTAAATGGAAAGTGAATTTTA  >  1/566063‑566548                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     aCCTAAGAAATCTTTCGGCAGGGCGTTTGCAAGTAAAGTACCACCACCGGCCATGTAGGCGGCTGATAATGCATATAACAAAATGAAAAGTGAGAGTGTTGAAATAATACGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAgg                                                                                                                                                                                                                                             >  1:791744/1‑251 (MQ=255)             tttCGGCAGGGCGTTTGCAAGTAAAGTACCACCACCGGCCATGTAGGCGGCTGATAATGCATATAACAAAATGAAAAGTGAGAGTGTTGAAATAATACGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGc                                                                                                                                                                                                                                 <  1:667173/251‑1 (MQ=255)                      ggCGTTTGCAAGTAAAGTACCACCACCGGCCATGTAGGCGGCTGATAATGCATATAACAAAATGAAAAGTGAGAGTGTTGAAATAATACGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCcggcgg                                                                                                                                                                                                                           <  1:546618/248‑1 (MQ=255)                              cAAGTAAAGTACCACCACCGGCCATGTAGGCGGCTGATAATGCATATAACAAAATGAAAAGTGAGAGTGTTGAAATAATACGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCg                                                                                                                                                                                                                                           >  1:776656/1‑224 (MQ=255)                              cAAGTAAAGTACCACCACCGGCCATGTAGGCGGCTGATAATGCATATAACAAAATGAAAAGTGAGAGTGTTGAAATAATACGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACgg                                                                                                                                                                                                                 <  1:752732/250‑1 (MQ=255)                                     agtaCCACCACCGGCCATGTAGGCGGCTGATAATGCATATAACAAAATGAAAAGTGAGAGTGTTGAAATAATACGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGa                                                                                                                                                                                                         >  1:756825/1‑251 (MQ=255)                                         ccaccaccGGCCATGTAGGCGGCTGATAATGCATATAACAAAATGAAAAGTGAGAGTGTTGAAATAATACGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCAt                                                                                                                                                                                                     <  1:503386/251‑1 (MQ=255)                                           accaccGGCCATGTAGGCGGCTGATAATGCATATAACAAAATGAAAAGTGAGAGTGTTGAAATAATACGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGt                                                                                                                                                                                                                                                  <  1:896737/204‑1 (MQ=255)                                                  gCCATGTAGGCGGCTGATAATGCATATAACAAAATGAAAAGTGAGAGTGTTGAAATAATACGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGc                                                                                                                                                                                               <  1:5753/248‑1 (MQ=255)                                                           gcggcTGATAATGCATATAACAAAATGAAAAGTGAGAGTGTTGAAATAATACGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCgg                                                                                                                                                                                    <  1:322876/250‑1 (MQ=255)                                                                       gCATATAACAAAATGAAAAGTGAGAGTGTTGAAATAATACGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTTCGGCggtcaa                                                                                                                                                                       >  1:326609/1‑246 (MQ=255)                                                                       gCATATAACAAAATGAAAAGTGAGAGTGTTGAAATAATACGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGGTAAt                                                                                                                                                                       >  1:109306/1‑251 (MQ=255)                                                                                   aTGAAAAGTGAGAGTGTTGAAATAATACGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGgcgccgc                                                                                                                                                            >  1:790856/1‑247 (MQ=255)                                                                                           tGAGAGTGTTGAAATAATACGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGAAGGt                                                                                                                                                     <  1:699316/249‑1 (MQ=255)                                                                                                  gtTGAAATAATACGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGGTAATTAAGGCGCTGCCGAAGGGTTTAttttt                                                                                                                                            >  1:791438/1‑251 (MQ=255)                                                                                                      aaaTAATACGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGGAGGGTTTATTTTTtaat                                                                                                                                        >  1:328711/1‑251 (MQ=255)                                                                                                       aataatACGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCAt                                                                                                                                                                                                     >  1:815672/1‑189 (MQ=255)                                                                                                         taataCGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGgc                                                                                                                                                                 >  1:267545/1‑223 (MQ=255)                                                                                                         taataCGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGGAGGGTTTATTTTTtaataa                                                                                                                                      >  1:312087/1‑250 (MQ=255)                                                                                                              cGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGGTAATTAAGGCGCTGCCGAAGGTTTTATTTTTTAATAacataca                                                                                                                                >  1:105870/1‑251 (MQ=255)                                                                                                                tCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGAAGGTTTTATTTTTTAATAacatacat                                                                                                                               <  1:273558/250‑1 (MQ=255)                                                                                                                tCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGAAGGTTTTATTTTTTAATAacatacat                                                                                                                               <  1:393886/250‑1 (MQ=255)                                                                                                                tCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGAAGGTTTTATTTTTTAATAACATACAtt                                                                                                                              <  1:385762/251‑1 (MQ=255)                                                                                                                               atatTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGAAGGTTTTATTTTTTAATAACATACATTAATCCTAAATAAAGt                                                                                                               <  1:157230/251‑1 (MQ=255)                                                                                                                                 atTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGAAGGTTTTATTTTTTAATAACATACATTAATCCTAAATAAAGTAg                                                                                                             <  1:274285/251‑1 (MQ=255)                                                                                                                                            aaGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGAAGGTTTTATTTTTTAATAACATACATTAATCCTAAATAAAGTAGGCTGGtattat                                                                                                  <  1:895418/251‑1 (MQ=255)                                                                                                                                               gTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGAAGGTTTTATTTTTTAATAACATACATTAATCCTAAATAAAGTAGGCTGGTATTATCg                                                                                                <  1:577101/250‑1 (MQ=255)                                                                                                                                                  gCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGAAGGTTTTATTTTTTAATAACATACATTAATcc                                                                                                                         <  1:264769/222‑1 (MQ=255)                                                                                                                                                       gCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGAAGGTTTTATTTTTTAATAACATACATTAATCCTAAATAAAGTAGGCTGGTATTATCg                                                                                                <  1:348393/242‑1 (MQ=255)                                                                                                                                                              tcTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCAc                                                                                                                                                                                             >  1:517958/1‑142 (MQ=255)                                                                                                                                                                                               aaaaGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGAAGGTTTTATTTTTTAATAACATACATTAATCCTAAATAAAGTAGGCTGGTATTATCGCATAAAGTCGTCGGAGATTAAGTGCTAATGCGCCGATTTTGTTATAAAg                                               <  1:864763/251‑1 (MQ=255)                                                                                                                                                                                                          aCTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGc                                                                                                                                                                                                                <  1:895246/79‑1 (MQ=255)                                                                                                                                                                                                           cTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCAc                                                                                                                                                                                             >  1:483811/1‑97 (MQ=255)                                                                                                                                                                                                           cTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGAAGGTTTTATTTTTTAATAACATACATTAATCCTAAATAAAGTAGGCTGGTATTATc                                                                                                 <  1:6495/189‑1 (MQ=255)                                                                                                                                                                                                             gTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATcc                                                                                                                                                                                                                                      >  1:881389/1‑54 (MQ=255)                                                                                                                                                                                                             gTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAgg                                                                                                                                                                  <  1:829496/122‑1 (MQ=255)                                                                                                                                                                                                                      gAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGGTAATTAAGGCGCTGCCGAAGGGTTTATTTTTTAATAAAATACATTAATCCTAAATAAAGTAGGGTGGGATTATCGCATAAAGTCGTCGGAGATTaaa                                                                          >  1:822253/1‑200 (MQ=255)                                                                                                                                                                                                                      gAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGAAGGTTTTATTTTTTAATAACATACATTAATCCTAAATAAAGTAGGCTGGTATTATCGCATAAAGTCGTCGGAGATTAAGTGCTAATGCGCCGATTTTGTTATaa                                                 <  1:233317/226‑1 (MQ=255)                                                                                                                                                                                                                      gAGTAACCAAAGTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGAAGGTTTTATTTTTTAATAACATACATTAATCCTAAATAAAGTAGGCTGGTATTATCGCATAAAGTCGTCGGAGATTAAGTGCTAATGCGCCGATTTTGTTATAAAGTGCGGTCTTTTTTg                                 <  1:27902/242‑1 (MQ=255)                                                                                                                                                                                                                                 gTCCGATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGGCGGAGGGTTTATTTTTTAATAACATACATTAATcc                                                                                                                         >  1:190819/1‑143 (MQ=255)                                                                                                                                                                                                                                     gATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAgg                                                                                                                                                                  <  1:597453/98‑1 (MQ=255)                                                                                                                                                                                                                                     gATGAGAAGTGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGAAGGTTTTATTTTTTAATAACATACATTAATCCTAAATAAAGTAGGc                                                                                                           <  1:612284/153‑1 (MQ=255)                                                                                                                                                                                                                                              tGCCTCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGAAGGTTTTATTTTTTAATAACATACATTAATCCTAAATAAAGTAGGCTGGTATTATCGCATAAAGTCGTCGGAGATTAAGTGCTAATGCGCCGATTTTGTTATAAAGTGCGGTCTTTTTTGAGAGCGTTTTTTAAATGGAAAGTGAATTtta  <  1:801297/249‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     ACCTAAGAAATCTTTCGGCAGGGCGTTTGCAAGTAAAGTACCACCACCGGCCATGTAGGCGGCTGATAATGCATATAACAAAATGAAAAGTGAGAGTGTTGAAATAATACGTCCGGGAATACCGAAATATTTTTCCGCCAAGGTTGCTACGCCATCATCTTTGCGTTTGGCAGTTTGATAAACTTCAACGAAAAGTAGTCCACTGTAAGCAAGGAGTAACCAAAGTCCGATGAGAAGTGCCGCGGTAAAGGCGAATCCCATGCCGGCGGAGGTTAACGGCATGGCGAGCATTCCGGCACCAATGGCGGTTCCGGCGATAATTAAGGCGCTGCCGAAGGTTTTATTTTTTAATAACATACATTAATCCTAAATAAAGTAGGCTGGTATTATCGCATAAAGTCGTCGGAGATTAAGTGCTAATGCGCCGATTTTGTTATAAAGTGCGGTCTTTTTTGAGAGCGTTTTTTAAATGGAAAGTGAATTTTA  >  1/566063‑566548

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