Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,367 G→A 7.5% intergenic (+28/+790) sseA → / ← sseB 3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,653,3670GA 7.5% 198.8 / 3.4 68intergenic (+28/+790)sseA/sseB3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate
Reads supporting (aligned to +/- strand):  ref base G (26/35);  new base A (2/3);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

TCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCAC  >  NC_000913/2653069‑2653665
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                          
gggcgGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGaa                                                                                                                                                                                                                                                                                                          <  1:91240/299‑1 (MQ=255)
      tcATATTCCCGGAGCACTGAAGTTTGTGTGGACGTGACGGGTGCGCGAAGGCGATCTAAAAACGACCGATTATCTGGATGCTATATTTTTTGGTCGCGGCTGCAGCTCCGACACACCTATTATCGTCAGCTGGTGCTCTGGTGTATCGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGTTGTGCCAAACGTGAAACTGTCCGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGc                                                                                                                                                                                                                                                                                                    <  1:467874/300‑1 (MQ=255)
        agtgtCCGGGGGCACTGAAGGGTCCGTGGGCGGAACGGGTGCGCGAAGGCGTACTAAAAACGACCGATGAACTGGATGCTATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCgg                                                                                                                                                                                                                                                                                                  <  1:157595/297‑1 (MQ=255)
           ttCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTTAAATAAGTATTTTACAGGCAATAAAAAACAGCCGAATTTGGCGGttt                                                                                                                                                                                                                                                                                               >  2:585658/1‑301 (MQ=255)
            ccccggggcctttaaTGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATGTTTTGGTCGCGGCGTCAGCTACGACAAACCAATAATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGtttt                                                                                                                                                                                                                                                                                              <  1:76103/288‑1 (MQ=255)
              ccGGGGCACTGAATGTTCCGTGGCAGGAATTGTGGCGCGAAGGCGAACGAAAAACGACCGATGAACTGGAGTCGATATTTGTTGGTCGCGGCGTCAGCTACGACAACCCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACTGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGtttttt                                                                                                                                                                                                                                                                                            <  1:82266/301‑1 (MQ=255)
               cGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTa                                                                                                                                                                                                                                                                                           <  2:551066/301‑1 (MQ=255)
                ggAGTACTGAAGGTTCCGTGTACGTAACTGGTGCGCGAAGGCGAACTAAAAATGACCGATGAACTGGATGCGATATTTTTTGGTCGCGTCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGTCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTAt                                                                                                                                                                                                                                                                                          <  2:312674/301‑1 (MQ=255)
                  aGCTCTGAATGTTCCGTTGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTg                                                                                                                                                                                                                                                                                        <  2:333323/301‑1 (MQ=255)
                     cgcaccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTa                                                                                                                                                                                                                                                                                     >  2:325207/8‑301 (MQ=255)
                           agTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGGGGTTTTGTTCGCACTCGCGCCGCTGGATGTGCCAACCGGGAAACTGTACGCCGGCCCAGGGCGTGAATGGGGCGCGCGGGCAGCTTTACCGGTTGAGCCAGTTAAATAAGTCTTTTACAGGCAATAAAAAACCGCCAAATTAGGCAGTTTTTTagtccaactccgg                                                                                                                                                                                                                                                                               >  1:612077/2‑289 (MQ=255)
                                      cGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCtt                                                                                                                                                                                                                                                                    >  2:241817/1‑301 (MQ=255)
                                       ggAACTGGTGGGGGAAGGCGAACTAAAATGGACCGATGAACTGGATGCGTTATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCGGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCttt                                                                                                                                                                                                                                                                   <  1:57971/301‑1 (MQ=255)
                                        gAACTGGTGCGCGAATGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGACGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTc                                                                                                                                                                                                                                                                  >  2:154076/1‑301 (MQ=255)
                                        gAACTGGTGCGCGAAGGCGAACTAAAAACGACCGCTGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGGTGTGCCAAACGTGAAACTGTACGACGGCGCATGGCGTGACTGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTCAAATAAGTATTTTACAGGCAATAAAAACACGCACAAATTGGCGGATTTTTATTTCTAGTTGGGTTCCGCGCCtttt                                                                                                                                                                                                                                                                  >  1:347902/1‑300 (MQ=255)
                                        gAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTc                                                                                                                                                                                                                                                                  >  2:22997/1‑301 (MQ=255)
                                                 cgcgAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGtt                                                                                                                                                                                                                                                         >  2:222814/1‑301 (MQ=255)
                                                         ggagctgaaaagtacgtaTGAACTGGATGTGATAGTTTTTGGTGGCGGGGTCAGCTAGGACAACTCAATTATCGTCAGCTGCGGCTCTGGTGGAACGGCAGGCGTGGTTTTGTTAGCACTCGCGACGCTGGATGGGCCAAACGGGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgtg                                                                                                                                                                                                                                                 <  1:584934/284‑1 (MQ=255)
                                                            attAAAAACGACCGATGAATGGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTa                                                                                                                                                                                                                                              <  1:222571/299‑1 (MQ=255)
                                                                         gATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGACATAAGTATTTTACAGGCAATAAAACACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGAATTGTGTTACATGAGACACGCa                                                                                                                                                                                                                                 >  1:222276/1‑301 (MQ=255)
                                                                                 ggATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTc                                                                                                                                                                                                                         >  2:379021/1‑301 (MQ=255)
                                                                                    gttggtatGTTTTGGTTGGGGCGTCAGCGACGACAAACCAATTATCGTCAGCTGCGGTGCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGGGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTcaca                                                                                                                                                                                                                      <  1:532823/296‑1 (MQ=255)
                                                                                       tgttgttttttGTTGCTGCGTCAGCTACGACAAAACAATTTTCGTTAGCTGCGGCTTTGGTGTTACGGCAGCCGTGTTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCACCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCa                                                                                                                                                                                                                   <  2:268324/296‑1 (MQ=255)
                                                                                              tttGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGGGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACaaa                                                                                                                                                                                                            >  2:517827/1‑300 (MQ=255)
                                                                                               ttGGTTGCGGCTTTAGTTACGACAAACCAATTATCGTCAGTTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGGGCCAAACGTGAAACTGTACGACGGCGAATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTGGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATa                                                                                                                                                                                                           <  1:378525/301‑1 (MQ=255)
                                                                                                  gTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCCTGGAGTGAATGCGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCCATAAAAAACCACCGACTTTGGCGGATTTTTATTTCTAGTCTGGTTCGCGGCCTTTCCAGCCGGTTGACTTGTCTTACCTTCGCAACGCACGTCCTTCACCGAAACACAATACTc                                                                                                                                                                                                        >  1:87747/1‑301 (MQ=255)
                                                                                                    cgcgGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAAAGCAGGTGCTTAAAAGCAAAACAATACTCaa                                                                                                                                                                                                      >  2:642178/1‑300 (MQ=255)
                                                                                                     gcggcgTCTGCTACGACAAACCAATTATCGTCAGCTGCGTCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGGGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCAcc                                                                                                                                                                                                     <  2:6430/301‑1 (MQ=255)
                                                                                                                  cGATAAACCTGTTATCGTTAGCTGCGGCTTTGGTTTAACGGCATCCGTGGTTTTGTTAGCACTCGCGACGCTGTATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGAGTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCttt                                                                                                                                                                                        <  2:338220/301‑1 (MQ=255)
                                                                                                                  cGACAAGCCTATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCTTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCttt                                                                                                                                                                                        <  1:641301/301‑1 (MQ=255)
                                                                                                                   gACAAACCATTTATGGTCAGGTGCGGGTGTTGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCtttt                                                                                                                                                                                       <  1:517155/301‑1 (MQ=255)
                                                                                                                       aaCCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGAGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTACTTCACAGCAACACAATAATCACCAGAAACTCTCTTTTTGtt                                                                                                                                                                                   >  1:361882/1‑300 (MQ=255)
                                                                                                                        acgtattttggttggttgcggtttgtgtgtAAGCGCAGCCGTGGTTTTGTTAGCACGCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCa                                                                                                                                                                                  <  1:37360/275‑1 (MQ=255)
                                                                                                                          cAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAg                                                                                                                                                                                <  2:222519/301‑1 (MQ=255)
                                                                                                                                  gTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAagag                                                                                                                                                                        <  2:612935/301‑1 (MQ=255)
                                                                                                                                      ggtggggttttGGTGTAACGCCAGCCTTTGTTTTGTTAGCACTCGGGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTGCTTGACAAAAAGAGAGtt                                                                                                                                                                    <  1:140494/291‑2 (MQ=255)
                                                                                                                                              tttGGTGTAACGGCAGCCGTGGTTTTGTGAGCACTCGGGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGt                                                                                                                                                            <  1:241550/299‑1 (MQ=255)
                                                                                                                                                  ttttAACGGCAGCCTTTGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGtttat                                                                                                                                                        <  2:565797/298‑1 (MQ=255)
                                                                                                                                                    gtAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGAAAAACAATACTCAACAATAACTCTCTTTTTGTAAAGCAAAAGAAAGTAATTAATGtttaaag                                                                                                                                                      >  1:565083/1‑298 (MQ=255)
                                                                                                                                                        cGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCg                                                                                                                                                  <  2:409457/301‑1 (MQ=255)
                                                                                                                                                                ttGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCg                                                                                                                                          <  2:38988/300‑1 (MQ=255)
                                                                                                                                                                          tAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCAGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGccc                                                                                                                                <  2:616101/301‑1 (MQ=255)
                                                                                                                                                                          aggCACTGGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGccc                                                                                                                                <  1:22946/299‑1 (MQ=255)
                                                                                                                                                                            gCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATACAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCTGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAAACCTCTTTTTGTCAAGCAAAAGAGCGTAATTAATGTTTATTTAACGTATTATCAACACCGGACCtt                                                                                                                              >  2:87040/1‑301 (MQ=255)
                                                                                                                                                                                 ttgctcttccgATCTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTccgc                                                                                                                         <  2:303592/291‑1 (MQ=255)
                                                                                                                                                                                        gttgATGTGCCAATTGTGTAACTGTACGACGGCTTATGGAGTGAATGGGGCTCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCTTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTc                                                                                                                  <  2:402711/298‑1 (MQ=255)
                                                                                                                                                                                         gtgctgtgCCAATGGTGAAGCTGTACTACGGGGCATGGAGTGAATGGGGCGCGCGGGAAGATTGACCGGTTGAGCAAGTGAAATAAGTATTTTAAAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTGATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCg                                                                                                                 <  1:153925/297‑1 (MQ=39)
                                                                                                                                                                                               tgCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGATGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTAAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGtaata                                                                                                           >  2:559354/1‑301 (MQ=255)
                                                                                                                                                                                                gCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGtaataa                                                                                                          >  2:124020/1‑301 (MQ=255)
                                                                                                                                                                                                gCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGGCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTccgccaatccgaaacacct                                                                                                          >  1:303188/1‑287 (MQ=255)
                                                                                                                                                                                                             ttGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCtt                                                                                             <  2:362366/300‑1 (MQ=255)
                                                                                                                                                                                                                  cgacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCttattat                                                                                        >  1:119836/1‑300 (MQ=255)
                                                                                                                                                                                                                      ggttacggggttGAATGTGGGGGGTGGGCAGTGTGTCCGGTTGTGCTAGTGAAATAAGTATTTTGCAGGCAATAAAAAACCGCCTAATTTGTCCGTTTTTGATGGCTAGTCTGGTTCGCGGCCTGTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACGCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCg                                                                                    <  1:149460/290‑1 (MQ=255)
                                                                                                                                                                                                                                ttGAATGGGGCGCGTGGTCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACa                                                                          <  2:348365/300‑1 (MQ=255)
                                                                                                                                                                                                                                         gcgggcggGCAGATTTACCGGTTGAGCCAGTGAAATAATTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATataaat                                                                 <  1:123896/301‑1 (MQ=255)
                                                                                                                                                                                                                                              cgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATAAAAATAAATc                                                            >  2:342162/1‑301 (MQ=255)
                                                                                                                                                                                                                                                  cAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCCAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACCCCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGc                                                        >  1:476717/1‑301 (MQ=255)
                                                                                                                                                                                                                                                     aTTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTAATTAGCGTATTATCGACACCAGCCCTTTACGCCGTGTTCGGTAATAAAATAAACACACTTATTACTCCTAATTAAGACAACAAAAAATAAATcataatca                                                     >  1:148278/1‑294 (MQ=38)
                                                                                                                                                                                                                                                                gAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGGTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTAGTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAAt                                          <  2:87871/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                  gCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCCCAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCTCCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACACACATAAATAAATCCTACTCAAAATTAAAAAttt                                        >  2:188254/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                     aGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTaa                                     >  2:133997/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                          tcagacgtgtgctcttccgatctaaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGt                                <  2:149998/278‑1 (MQ=255)
                                                                                                                                                                                                                                                                              aGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAACATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAAAAAATCCTGCTCAAAATTAAAAATTCTAACCGGTaaaa                            >  1:613068/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                  ttttAAAGGCAATAAAAATCCGCCGAAGTTGGCGGGTTTTTAGTGCTAGTCGGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGAAACGCAGGTGCTTCAAAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGata                        <  1:133872/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         ggCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATAATACTTa                 >  2:194489/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                           cAATAAATAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTaaa               <  2:119960/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             aTAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCACCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAACATTAAAAATTCAAACCGGAAAAAGATATTACTAAAACa             >  1:248303/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                 aaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACAGGTagatcggaagagcgacgtgtagg         >  1:149847/1‑279 (MQ=255)
                                                                                                                                                                                                                                                                                                     ccgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTAACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCATTTCCGCCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAAATCAGACAAATAAAAATAAATCCTGCACAAAATTAAAAATTCAAACAGGAAAAAGTAATAACTTaacactataaatt     >  1:3446/1‑290 (MQ=38)
                                                                                                                                                                                                                                                                                                        ccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGAATTGTGTTACATGGGCAACGCAGGTGCTTCACAGCAAAAAAATACTCACCAGTAACGCTCTTTTTGTAAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTACGAATACAGACAAATATAAATAAATACTGCTCAAAATTAAAAATTCAAACCGGAAAAAGATATTACTTAAACATGTAAATTCAc  >  1:112939/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                          
TCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCAC  >  NC_000913/2653069‑2653665

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: