Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 682,714 T→C intergenic (+185161/‑16629) 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
*1682,7140TC98.5% 184.5 / ‑5.1 68intergenic (+185161/‑16629)unknown/unknowntRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (41/26);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

GGCTATTATGTAAGCACGGTAGGGCTAAATACAAAGGTAGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCCTTATAGAGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACTCTATCCATTATGGGTAGGTTCGTCCCGCAAGGGTGGTTTTTTAATTTAACGG  >  1/682520‑682960
                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
ggCTATTATGTAAGCACGGTAGGGCTAAATACAAAGGTAGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTAttttatgat                                                                                                                                                                                                >  1:203148/1‑246 (MQ=2)
ggCTATTATGTAAGCACGGTAGGGCTAAATACAAAGGTAGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTg                                                                                                                                                                                                >  1:719766/1‑251 (MQ=17)
ggCTATTATGTAAGCACGGTAGGGCTAAATACAAAGGTAGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTg                                                                                                                                                                                                >  1:262100/1‑251 (MQ=17)
                            aTACAAAGGTAGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAGGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAAc                                                                                                                                                                               >  1:705038/1‑240 (MQ=255)
                             tACAAAGGTAGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAAtt                                                                                                                                                                    >  1:807770/1‑250 (MQ=255)
                                    gTAGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACAt                                                                                                                                                                             <  1:704376/234‑1 (MQ=255)
                                     tAGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCa                                                                                                                                                            >  1:121018/1‑250 (MQ=255)
                                       gtgtAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACtt                                                                                                                                                         >  1:54500/1‑251 (MQ=255)
                                       gtgtAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTAGACCGCACtt                                                                                                                                                         >  1:432470/1‑251 (MQ=255)
                                               gTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTtgttgt                                                                                                                                                 >  1:114133/1‑251 (MQ=255)
                                               gTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTtgttgt                                                                                                                                                 >  1:297376/1‑251 (MQ=255)
                                               gTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTtgttgt                                                                                                                                                 >  1:79598/1‑251 (MQ=255)
                                               gTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTtgttgt                                                                                                                                                 >  1:807378/1‑251 (MQ=255)
                                               gTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATACGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTaaaa                                                                                                                                                                                  >  1:649152/1‑218 (MQ=255)
                                                       ggAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGt                                                                                                                                           >  1:7322/1‑249 (MQ=255)
                                                        gAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGtct                                                                                                                                         >  1:404344/1‑250 (MQ=255)
                                                        gAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGtct                                                                                                                                         <  1:793161/250‑1 (MQ=255)
                                                        gAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGtc                                                                                                                                          >  1:707668/1‑249 (MQ=255)
                                                              ggAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGc                                                                                                                                                                                       >  1:728366/1‑198 (MQ=255)
                                                                            tGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTGTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGaa                                                                                                                    <  1:433788/251‑1 (MQ=255)
                                                                              aTTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAg                                                                                                                  >  1:38189/1‑251 (MQ=255)
                                                                              aTTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAg                                                                                                                  >  1:156630/1‑251 (MQ=255)
                                                                                      tAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCTCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAg                                                                                                          <  1:448156/251‑1 (MQ=255)
                                                                                      tAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAg                                                                                                          <  1:189806/251‑1 (MQ=255)
                                                                                      tAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAg                                                                                                          <  1:765615/251‑1 (MQ=255)
                                                                                                        atatGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGAc                                                                                                                                                                >  1:567567/1‑179 (MQ=255)
                                                                                                         tatGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAattc                                                                                                                                  >  1:642302/1‑208 (MQ=255)
                                                                                                         tatGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAACAGCTTATTg                                                                                        >  1:738604/1‑250 (MQ=255)
                                                                                                         tatGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGc                                                                                       <  1:182577/251‑1 (MQ=255)
                                                                                                           tgtgGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAg                                                                                                          <  1:611023/230‑1 (MQ=255)
                                                                                                                   ccTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAattcattc                                                                                                                              <  1:461475/202‑1 (MQ=255)
                                                                                                                   ccTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATg                                                                             >  1:368567/1‑251 (MQ=255)
                                                                                                                    cTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAg                                                                                               >  1:883993/1‑232 (MQ=255)
                                                                                                                         aGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAg                                                                                                          <  1:791961/216‑1 (MQ=255)
                                                                                                                                 gTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTaaaaaa                                                                                                                                                                                >  1:163679/1‑138 (MQ=255)
                                                                                                                                 gTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAg                                                                                                                  <  1:667338/200‑1 (MQ=255)
                                                                                                                                 gTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGa                                                               >  1:770265/1‑251 (MQ=255)
                                                                                                                                  tCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCG‑TATGTTTTTTGTGCTAGaa                                                              <  1:161021/250‑1 (MQ=255)
                                                                                                                                       gagTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCg                                                          >  1:784670/1‑250 (MQ=255)
                                                                                                                                       gagTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCg                                                          >  1:378821/1‑250 (MQ=255)
                                                                                                                                           tATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCAccc                                                     <  1:604799/251‑3 (MQ=255)
                                                                                                                                           tATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCAccc                                                     <  1:60690/251‑3 (MQ=255)
                                                                                                                                                cACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTCTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATcc                                                <  1:442594/251‑1 (MQ=255)
                                                                                                                                                cACGGTTGTTGGTATTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGAGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATcc                                                <  1:592589/251‑1 (MQ=255)
                                                                                                                                                    gttgttGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTc                                                                                                                                                                                                   >  1:114875/1‑100 (MQ=255)
                                                                                                                                                     ttgttgGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATCCATTAt                                           <  1:596947/251‑1 (MQ=255)
                                                                                                                                                         tgGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATCCATTATGGGt                                       <  1:110779/251‑1 (MQ=255)
                                                                                                                                                              ttCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATCCATTATGGGTAGGtt                                  <  1:656965/251‑1 (MQ=255)
                                                                                                                                                              ttCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATCCATTATGGGTAGGtt                                  <  1:822286/251‑1 (MQ=255)
                                                                                                                                                              ttCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATCCATTATGGGTAGGtt                                  <  1:637833/251‑1 (MQ=255)
                                                                                                                                                                cATATGCCCTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATCCa                                               <  1:377388/236‑1 (MQ=255)
                                                                                                                                                                      cccTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATCCATTATg                                          >  1:619991/1‑235 (MQ=255)
                                                                                                                                                                       ccTTTGAAGGGCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATCCATTa                                            <  1:304165/232‑1 (MQ=255)
                                                                                                                                                                               gggCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATCCATTAt                                           >  1:603433/1‑225 (MQ=255)
                                                                                                                                                                                ggCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAatt                                                                                                                                   >  1:475590/1‑136 (MQ=255)
                                                                                                                                                                                ggCATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGc                                                                                       >  1:139948/1‑180 (MQ=255)
                                                                                                                                                                                  cATGTGAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATCCATTATGGGTAGGTTCGTCCCGCAAGGGTGGtttt              >  1:564492/1‑251 (MQ=255)
                                                                                                                                                                                     gtgAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATCCATTAt                                           >  1:179212/1‑219 (MQ=255)
                                                                                                                                                                                     gtgAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATCCATTATGGGTAGGTTCGTCCCGCAAggg                     >  1:323021/1‑241 (MQ=255)
                                                                                                                                                                                     gtgAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATCCATTATGGGTAGGTTCGTCCCGCAAGGGTGGTTTttta           >  1:647553/1‑251 (MQ=255)
                                                                                                                                                                                     gtgAAGTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATCCATTATGGGTAGGTTCGTCCCGCAAGGGTGGTTTttta           >  1:841238/1‑251 (MQ=255)
                                                                                                                                                                                          gTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAGCATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTg                                                                                        <  1:865443/169‑1 (MQ=255)
                                                                                                                                                                                          gTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCAccc                                                     <  1:718687/204‑3 (MQ=255)
                                                                                                                                                                                          gTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATCCATTATGGGTAGGTTCGTCCCGCAAGGGTGGTTTtttaattta      >  1:213221/1‑251 (MQ=255)
                                                                                                                                                                                          gTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATCCATTATGGGTAGGTTCGTCCCGCAAGGGTGGTTTtttaattta      >  1:61688/1‑251 (MQ=255)
                                                                                                                                                                                          gTACTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATCCATTATGGGTAGGTTCGTCCCGCAAGGGTGGTTTtttaattta      >  1:391296/1‑251 (MQ=255)
                                                                                                                                                                                            aCTAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATCCATTATGGGTAGGTTCGTCCCGCAAGGGTGGTTTtttaatttaa     >  1:626927/1‑250 (MQ=255)
                                                                                                                                                                                              tAGGCCCTTATAGGGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACCCTATCCATTATGGGTAGGTTCGTCCCGCAAGGGTGGTTTTTTAATTTAACgg  <  1:121738/251‑1 (MQ=255)
                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
GGCTATTATGTAAGCACGGTAGGGCTAAATACAAAGGTAGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCCTTATAGAGTAGCCTAAAAACCGCACGTTTTACGTGCGGATTGTTATTTTTCTTGATTCACAGCTAAAAAACATCTGGAAATTCGACCGCACTTTTTGTTGTTTACGTCTCAAATTCATTCTTTTCATGAAAGCGTATAAGTTATTTAAAAGCTTATTGCTTTCGTTATGTTTTTTGTGCTAGAATCCGCACTCTATCCATTATGGGTAGGTTCGTCCCGCAAGGGTGGTTTTTTAATTTAACGG  >  1/682520‑682960

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