Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 610,269 C→T intergenic (+112716/‑89074) 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
*1610,2690CT100.0% 219.7 / NA 81intergenic (+112716/‑89074)unknown/unknowntRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (37/44);  total (37/44)

ATCGCATTCGTGCCGCCTTGTTAGCCGGTATTCGCGCGGCGATTTTATGGCAACAGGTGGGCGGCACCAAATGGCAGTTGTTGTTCTTTCGAGGCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGCGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAGACATCGAAAAAACCACCAATCACGACGTAAAAGCGGTGGAATATTTCCTCAAAGAAAAAAGCGCAGCGTTGCTG  >  1/610026‑610509
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
aTCGCATTCGTGCCGCCTTGTTAGCCGGTATTCGCGCGGCGATTTTATGGCAACAGGTGGGCGGCACCAAATGGCAGTTGTTGTTCTTTCGAGGCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCg                                                                                                                                                                                                                                           >  1:733024/1‑251 (MQ=255)
 tCGCATTCGTGCCGCCTTGTTAGCCGGTATTCGCGCGGCGATTTTATGGCAACAGGTGGGCGGCACCAAATGGCAGTTGTTGTTCTTTCGAGGCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCgg                                                                                                                                                                                                                                          <  1:848892/251‑1 (MQ=255)
          tGCCGCCTTGTTAGCCGGTATTCGCGCGGCGATTTTATGGCAACAGGTGGGCGGCACCAAATGGCAGTTGTTGTTCTTTCGAGGCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTaa                                                                                                                                                                                                                                  <  1:695456/250‑1 (MQ=255)
          tGCCGCCTTGTTAGCCGGTATTCGCGCGGCGATTTTATGGCAACAGGTGGGCGGCACCAAATGGCAGTTGTTGTTCTTTCGAGGCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTaa                                                                                                                                                                                                                                  <  1:7737/250‑1 (MQ=255)
          tGCCGCCTTGTTAGCCGGTATTCGCGCGGCGATTTTATGGCAACAGGTGGGCGGCACCAAATGGCAGTTGTTGTTCTTTCGAGGCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTaa                                                                                                                                                                                                                                  <  1:781353/250‑1 (MQ=255)
              gccTTGTTAGCCGGTATTCGCGCGGCGATTTTATGGCAACAGGTGGGCGGCACCAAATGGCAGTTGTTGTTCTTTCGAGGCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAAtt                                                                                                                                                                                                                              >  1:502881/1‑250 (MQ=255)
               ccTTGTTAGCCGGTATTCGCGCGGCGATTTTATGGCAACAGGTGGGCGGCACCAAATGGCAGTTGTTGTTCTTTCGAGGCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGTCGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTc                                                                                                                                                                                                                             <  1:733640/250‑1 (MQ=255)
               ccTTGTTAGCCGGTATTCGCGCGGCGATTTTATGGCAACAGGTGGGCGGCACCAAATGGCAGTTGTTGTTCTTTCGAGGCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTCTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTc                                                                                                                                                                                                                             <  1:893128/250‑1 (MQ=255)
                   gTTAGCCGGTATTCTCGCGGCGATTTTATGGCAACAGGTGGGCGGCACCAAATGGCAGTTGTTGTTCTTTCGAGGCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCgtgt                                                                                                                                                                                                                        >  1:317885/1‑251 (MQ=255)
                           gTATTCGCGCGGCGATTTTATGGCAACAGGTGGGCGGCACCAAATGGCAGTTGTTGTTCTTTCGAGGCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAg                                                                                                                                                                                                                >  1:18848/1‑251 (MQ=255)
                           gTATTCGCGCGGCGATTTTATGGCAACAGGTGGGCGGCACCAAATGGCAGTTGTTGTTCTTTCGAGGCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTACACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAg                                                                                                                                                                                                                >  1:721912/1‑251 (MQ=255)
                                cgcgcgGCGATTTTATGGCAACAGGTGGGCGGCACCAAATGGCAGTTGTTGTTCTTTCGAGGCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGtt                                                                                                                                                                                                           <  1:752746/251‑1 (MQ=255)
                                    cggcgATTTTATGGCAACAGGTGGGCGGCACCAAATGGCAGTTGTTGTTCTTTCGAGGCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGt                                                                                                                                                                                                        <  1:252735/250‑1 (MQ=255)
                                                     cAGGTGGGCGGCCCCAAATGGCAGTTGTTGTTCTTTCGAGGCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTg                                                                                                                                                                                                      <  1:348397/235‑1 (MQ=255)
                                                        gTGGGCGGCACCAAATGGCAGTTGTTGTTCTTTCGAGGCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGgcggcg                                                                                                                                                                                   >  1:860509/1‑251 (MQ=255)
                                                                ccacAAATGGCATTTGTTGTTCTTTCGAGGCAAAATCGCCGCCCCCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTTATGGAGCTTCACCCTATGCACCTCAGCGCATTAACCGCACTCTCCCCCATTGACGGCCGTTATAAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGcac                                                                                                                                                                            <  1:451401/247‑1 (MQ=255)
                                                                          cAGTTGTTGTTCTTTCGAGGCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGa                                                                                                                                                                  <  1:634726/250‑1 (MQ=255)
                                                                                     cTTTCGAGGCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCttt                                                                                                                                                      <  1:18614/251‑1 (MQ=255)
                                                                                           gggCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTc                                                                                                                                                                                                                             <  1:611277/173‑1 (MQ=255)
                                                                                                    cgccgcCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCaaa                                                                                                                                       >  1:763867/1‑251 (MQ=255)
                                                                                                    cgccgcCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCaaa                                                                                                                                       >  1:537331/1‑251 (MQ=255)
                                                                                                                 caacaaATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGa                                                                                                                                                                  >  1:444283/1‑211 (MQ=255)
                                                                                                                 caacaaATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAg                                                                                                                                                                <  1:382111/213‑1 (MQ=255)
                                                                                                                  aacaaATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACg                                                                                                                         <  1:597446/251‑1 (MQ=255)
                                                                                                                   aaaaaTTCTCTCTTCTCTTTAATCTAACTCATGGAGCTTCACGCGATGAAACTCAGCGCATTAACCGAACTCTCCCCCATTGACGGACGTTATCTGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAAGTCCGTGTAGCGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCCTCTTTGTCGCAGGAAGCAAACGATTACCTTGACgg                                                                                                                        <  1:872420/249‑1 (MQ=255)
                                                                                                                     aaaTTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACgg                                                                                                                                                                                                                   <  1:669011/158‑1 (MQ=255)
                                                                                                                       aTTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGACTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTAtt                                                                                                                    >  1:443231/1‑251 (MQ=255)
                                                                                                                           tctcttctctTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAg                                                                                                                <  1:573359/251‑1 (MQ=255)
                                                                                                                           tctcttctctTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAg                                                                                                                <  1:492297/251‑1 (MQ=255)
                                                                                                                              cttctcttTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGa                                                                                                                                                                  >  1:542312/1‑198 (MQ=255)
                                                                                                                                tctcttTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAAt                                                                                                           >  1:413080/1‑251 (MQ=255)
                                                                                                                                   cttTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGa                                                                                                                                                                  >  1:801197/1‑193 (MQ=255)
                                                                                                                                      taatttaaCTCATGGAGCTTTACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGa                                                                                                                                                                  <  1:425039/190‑1 (MQ=255)
                                                                                                                                      taatttaaCTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCACTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTg                                                                                                                            <  1:47177/228‑1 (MQ=255)
                                                                                                                                       aatttaaCTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCa                                                                                                    >  1:557871/1‑251 (MQ=255)
                                                                                                                                              cTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTACTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTg                                                                                              <  1:795119/250‑1 (MQ=255)
                                                                                                                                                cATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAg                                                                                                                <  1:203307/230‑1 (MQ=255)
                                                                                                                                                cATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATg                                                                                           >  1:565027/1‑251 (MQ=255)
                                                                                                                                                    gAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGcac                                                                                        <  1:174853/250‑1 (MQ=255)
                                                                                                                                                      gCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGa                                                                                                                                                                  >  1:645005/1‑174 (MQ=255)
                                                                                                                                                      gCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACaa                                                                                      >  1:830456/1‑250 (MQ=255)
                                                                                                                                                      gCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACaa                                                                                      >  1:869394/1‑250 (MQ=255)
                                                                                                                                                        ttCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGc                                                                                   <  1:270335/251‑1 (MQ=255)
                                                                                                                                                           aCGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCGGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGcac                                                                                        <  1:590456/243‑1 (MQ=255)
                                                                                                                                                           aCGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAAc                                                                                     >  1:110641/1‑246 (MQ=255)
                                                                                                                                                            cGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGGTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAg                                                                                                                <  1:437225/218‑1 (MQ=255)
                                                                                                                                                              cTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTaaa                                                                             >  1:236944/1‑251 (MQ=255)
                                                                                                                                                               tATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAg                                                                            <  1:53078/251‑1 (MQ=255)
                                                                                                                                                               tATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAg                                                                            <  1:535711/251‑1 (MQ=255)
                                                                                                                                                                aTGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCt                                                                                                                              <  1:888386/200‑1 (MQ=255)
                                                                                                                                                                aTGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAg                                                                                                                <  1:552812/214‑1 (MQ=255)
                                                                                                                                                                aTGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGcac                                                                                        <  1:832506/238‑1 (MQ=255)
                                                                                                                                                                aTGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAAc                                                                                     >  1:792929/1‑241 (MQ=255)
                                                                                                                                                                aTGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTaa                                                                              >  1:455031/1‑248 (MQ=255)
                                                                                                                                                                          gcgcATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGAtt                                                                                                                                  >  1:374228/1‑186 (MQ=255)
                                                                                                                                                                             cATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAGACATCGAAAAAac                                                               >  1:501342/1‑250 (MQ=255)
                                                                                                                                                                               ttAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATg                                                                                           <  1:60973/220‑1 (MQ=255)
                                                                                                                                                                               ttAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAGACATCGAAAAAacca                                                             >  1:396218/1‑250 (MQ=255)
                                                                                                                                                                                  aCCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAGACATCGAAAAAACCACCaa                                                         >  1:589862/1‑251 (MQ=255)
                                                                                                                                                                                        ctctCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTc                                                                                                                                                                                                                             >  1:737173/1‑81 (MQ=255)
                                                                                                                                                                                          ctcCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGATTTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACg                                                                                                                                                                    >  1:745421/1‑136 (MQ=255)
                                                                                                                                                                                          ctcCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTc                                                                                                                                                                                                                             >  1:468495/1‑79 (MQ=255)
                                                                                                                                                                                          ctcCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAGACATCGAAAAAACCACCAATCacgacg                                                 >  1:596817/1‑251 (MQ=255)
                                                                                                                                                                                            cccccATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGa                                                                                                                                                                  <  1:666021/136‑1 (MQ=255)
                                                                                                                                                                                            cccccATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCcttctt                                                                                                                                                       <  1:885059/147‑1 (MQ=255)
                                                                                                                                                                                              cccATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAg                                                                                                                                                                <  1:663995/136‑1 (MQ=255)
                                                                                                                                                                                               ccATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAAc                                                                                                                                      >  1:887783/1‑161 (MQ=255)
                                                                                                                                                                                                cATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAGACATCGAAAAAACCACCAATCACGACGTAAAAg                                           <  1:263877/251‑1 (MQ=255)
                                                                                                                                                                                                            ttATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAGACATCGAAAAAACCACCAATCACGACGTAAAAGCGGTGGAatat                                <  1:875187/250‑1 (MQ=255)
                                                                                                                                                                                                              aTCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAGACATCGAAAAAACCACCAATCACGACGTAAAAGCGGTGGAATATTTc                             <  1:479331/251‑1 (MQ=255)
                                                                                                                                                                                                                  ggATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGcac                                                                                        <  1:798857/188‑1 (MQ=255)
                                                                                                                                                                                                                   gATAAAACCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAGACATCGAAAAAACCACCAATCACGACGTAAAAGCGGTGGAATATTTCCTCa                         <  1:814688/250‑1 (MQ=255)
                                                                                                                                                                                                                      aaaaCCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTg                                                                                              >  1:532713/1‑178 (MQ=255)
                                                                                                                                                                                                                      aaaaCCGCCGAATTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAGACATCGAAAAAACCACCAATCACGACGTAAAAGCGGTGGAATATTTCCTCAAAGa                     >  1:479296/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aTTACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAg                                                                            <  1:744659/185‑1 (MQ=255)
                                                                                                                                                                                                                                   tACGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACAGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAGACATCGAAAAAACCACCAATCACGACGTAAAAGCGGTGGAATATTTCCTCAAAGAAAAAAGCGCAGCg        <  1:398636/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aCGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAGa                                                                           >  1:454923/1‑183 (MQ=255)
                                                                                                                                                                                                                                    aCGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAGACATCGAAAAAACCACCAATCACGACGTAAAAGCGGTGGAATATTTCCTCAAAGAAAAAAGCGCAGCGt       >  1:397002/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aCGCAGCATTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAGACATCGAAAAAACCACCAATCACGACGTAAAAGCGGTGGAATATTTCCTCAAAGAAAAAAGCGCAGCGt       <  1:243883/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gcaTTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAGACATCGAAAAAACCACCAATCACGACGTAAAAGCGGTGGAATATTTCCTCAAAGAAAAAAGCGCAGCGTTGccg  >  1:239129/1‑249 (MQ=255)
                                                                                                                                                                                                                                         gcaTTTTCAGTGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAGACATCGAAAAAACCACCAATCACGACGTAAAAGCGGTGGAATATTTCCTCAAAGAAAAAAGCGCAGCGTTGccg  >  1:561400/1‑249 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
ATCGCATTCGTGCCGCCTTGTTAGCCGGTATTCGCGCGGCGATTTTATGGCAACAGGTGGGCGGCACCAAATGGCAGTTGTTGTTCTTTCGAGGCAAAATCGCCGCCACCGCACAACAAATTCTCTCTTCTCTTTAATTTAACTCATGGAGCTTCACGCTATGCAACTCAGCGCATTAACCGCACTCTCCCCCATTGACGGACGTTATCAGGATAAAACCGCCGAATTACGCAGCATTTTCAGCGAATTCGGTTTATTAAAATTCCGTGTAACGGTAGAAGTTCGTTGGCTGCAAAAACTGGCGGCGACCGCACAAATCACGGAAGTGCCTTCTTTGTCGCAGGAAGCAAACGATTACCTTGACGGTATTGTAGAAAATTTCAGCATAATTGATGCACAACGCATTAAAGACATCGAAAAAACCACCAATCACGACGTAAAAGCGGTGGAATATTTCCTCAAAGAAAAAAGCGCAGCGTTGCTG  >  1/610026‑610509

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