Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 493,704 T→C intergenic (+2789/‑203004) 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
*1493,7040TC100.0% 233.1 / NA 69intergenic (+2789/‑203004)unknown/unknowntRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/39);  total (30/39)

TGGTGGCGCTTTTCTCCGCACCCAACAAATTGCGTTGCGCCCGATGCCAATCAATTAAAATTTTCACTTCCAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGTTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGCCCACCTTGGTAGAAAAAAGGTTTTTTTTCGTTATAATGCTCACACTTTTACACTAAAAATGAGCAAAATCCGTTCATTTTATCGCTATCCATATC  >  1/493466‑493950
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
tggtggCGCTTTTCTCCGCACCCAACAAATTGCGTTGCGCCCGATGCCAATCAATTAAAATTTTCACTTCCAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAAt                                                                                                                                                                                                                                            <  1:309393/251‑1 (MQ=255)
   tggCGCTTTTCTCCGCACCCAACAAATTGCGTTGCGCCCGATGCCAATCAATTAAAATTTTCCCTTCCAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATgg                                                                                                                                                                                                                                          <  1:65182/250‑1 (MQ=255)
            tctcCGCACCCAACAAATTGCGTTGCGCCCGATGCCAATCAATTAAAATTTTCACTTCCAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATgg                                                                                                                                                                                                                                 >  1:67569/1‑250 (MQ=255)
            tctcCGCACCCAACAAATTGCGTTGCGCCCGATGCCAATCAATTAAAATTTTCACTTCCAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGc                                                                                                                                                                                                                                <  1:118725/251‑1 (MQ=255)
             ctcCGCACCCAACAAATTGCGTTGCGCCCGATGCCAATCAATTAAAATTTTCACTTCCAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCa                                                                                                                                                                                                                               <  1:401278/251‑1 (MQ=255)
                    cccAACAAATTGCGTTGCGCCCGATGCCAATCAATTAAAATTTTCACTTCCAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCttt                                                                                                                                                                                                                         >  1:593368/1‑250 (MQ=255)
                       aacaaATTGCGTTGCGCCCGATGCCAATCAATTAAAATTTTCACTTCCAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAg                                                                                                                                                                                                                     <  1:659597/251‑1 (MQ=255)
                        acaaATTGCGTTGCGCCCGATGCCAATCAATTAAAATTTTCACTTCCAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGa                                                                                                                                                                                                                    <  1:51143/251‑1 (MQ=255)
                                   tgcgCCCGATGCCAATCAATAAAAATTTTCCCTTCCATTTCCGGTTTACCCTGCTTCGCAAAATAAATTTCATCCAGAATTCCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTgcgcgc                                                                                                                                                                                                         <  1:379910/251‑1 (MQ=255)
                                      gcCCGATGCCAATCAATTAAAATTTTCACTTCCAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACg                                                                                                                                                                                                      >  1:150159/1‑251 (MQ=255)
                                      gcCCGATGCCAATCAATTAAAATTTTCACTTCCAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACg                                                                                                                                                                                                      >  1:258674/1‑251 (MQ=255)
                                         cGATGCCAATCAATTAAAATTTTCACTTCCAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTgcgcgc                                                                                                                                                                                                         <  1:872530/245‑1 (MQ=255)
                                           aTGCCAATCAATTAAAATTTTCCCTTCAAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCT‑‑‑AATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGAtt                                                                                                                                                                                               <  1:840141/250‑1 (MQ=255)
                                           aTGCCAATCAATTAAAATTTTCACTTCCAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAACCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGa                                                                                                                                                                                                 >  1:445243/1‑251 (MQ=255)
                                           aTGCCAATCAATTAAAATTTTCACTTCCAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAACCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGa                                                                                                                                                                                                 >  1:429848/1‑251 (MQ=255)
                                              ccAATCAATTAAAATTTTCACTTCCAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGAtt                                                                                                                                                                                               <  1:264956/250‑1 (MQ=255)
                                                  tcaatTAAAATTTTCCCTTCCAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATc                                                                                                                                                                                          <  1:17883/251‑1 (MQ=255)
                                                           aTTTTCACTTCCATTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACAt                                                                                                                                                                                 <  1:800015/251‑1 (MQ=255)
                                                           aTTTTCACTTCCAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACAt                                                                                                                                                                                 >  1:584151/1‑251 (MQ=255)
                                                                ccctcccAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAAtt                                                                                                                                                                            <  1:261976/246‑1 (MQ=255)
                                                                cACTTCCATTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAAtt                                                                                                                                                                            <  1:732457/251‑1 (MQ=255)
                                                                          tCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCtat                                                                                                                                                                  <  1:659303/251‑1 (MQ=255)
                                                                          tCCGGTTTATCCTGCTTCGCACAATAAACTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCtat                                                                                                                                                                  <  1:428292/251‑1 (MQ=255)
                                                                             ggTTTATCCTGCTTCGCACAATAAATTTCATCCAGATTTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTtc                                                                                                                                                                       <  1:122133/243‑1 (MQ=255)
                                                                             ggTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCa                                                                                                                                                                                                                                        <  1:579911/178‑1 (MQ=255)
                                                                             ggTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTaa                                                                                                                                                                                     >  1:555020/1‑229 (MQ=255)
                                                                             ggTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAAtt                                                                                                                                                                            >  1:358345/1‑238 (MQ=255)
                                                                              ttttATCCTGCTTCGCACAATAAATTTCATCCAAAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGC‑TTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATaaaaa                                                                                                                                                             <  1:754745/250‑1 (MQ=255)
                                                                                  aTCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAAtt                                                                                                                                                                            <  1:366416/233‑1 (MQ=255)
                                                                                       gCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGa                                                                                                                                                     >  1:343502/1‑251 (MQ=255)
                                                                                       gCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGa                                                                                                                                                     >  1:734931/1‑251 (MQ=255)
                                                                                              acaataatttccacaaaGAATTTCTTTTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGaa                                                                                                                                              <  1:540515/235‑1 (MQ=255)
                                                                                              acaATAAATTCCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGaa                                                                                                                                              <  1:511629/251‑1 (MQ=255)
                                                                                                                 aTTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTa                                                                                                                           >  1:351157/1‑251 (MQ=255)
                                                                                                                   tttttGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAg                                                                                                                                                <  1:835187/225‑1 (MQ=255)
                                                                                                                   ttCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAg                                                                                                                                                >  1:360056/1‑228 (MQ=255)
                                                                                                                   ttCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAg                                                                                                                          >  1:563566/1‑250 (MQ=255)
                                                                                                                               cTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTAc                                                                                                             <  1:869463/251‑1 (MQ=255)
                                                                                                                                 cGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACaa                                                                                                           >  1:471109/1‑251 (MQ=255)
                                                                                                                                   tCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAat                                                                                                          <  1:641400/250‑1 (MQ=255)
                                                                                                                                            aaaTAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGaa                                                                                                                                              <  1:114314/204‑1 (MQ=255)
                                                                                                                                             aaTAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGAt                                                                                                                                                                                                <  1:240206/154‑1 (MQ=255)
                                                                                                                                              aTAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGCCCa                                                                                              >  1:380327/1‑251 (MQ=255)
                                                                                                                                                        gTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACACTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGc                                                                                                 >  1:441153/1‑238 (MQ=255)
                                                                                                                                                           acaATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAat                                                                                                          <  1:802765/226‑1 (MQ=255)
                                                                                                                                                           acaATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGCCCACCTTGGTAGaaaa                                                                                 >  1:874458/1‑251 (MQ=255)
                                                                                                                                                            caATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTtc                                                                                                                                                                       >  1:51540/1‑164 (MQ=255)
                                                                                                                                                            caATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCtat                                                                                                                                                                  >  1:552599/1‑169 (MQ=255)
                                                                                                                                                            caATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGaa                                                                                                                                              >  1:516524/1‑189 (MQ=255)
                                                                                                                                                                aatCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTAc                                                                                                             <  1:189945/218‑1 (MQ=255)
                                                                                                                                                                              cTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAAtt                                                                                                                                                                            >  1:756158/1‑141 (MQ=255)
                                                                                                                                                                                tGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGCCCACCTTGGTAGAAAAAAGGTTTTTTTTCGTTATAAt                                                            >  1:231709/1‑251 (MQ=255)
                                                                                                                                                                                     aGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGCCCACCTTGGTAGAAAAAAGGTTTTTTTTCGTTATAATGCTc                                                        <  1:834441/250‑1 (MQ=255)
                                                                                                                                                                                       tCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGCCCACCTTGGTAGAAAAAAGGTTTTTTTTCGTTATAATGCTcaca                                                     <  1:548208/251‑1 (MQ=255)
                                                                                                                                                                                             gAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGAttt                                                                                                                                                                                              <  1:646375/108‑1 (MQ=255)
                                                                                                                                                                                                taatTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGCCCACCTTGGTAGAAAAAAGGTTTTTTTTCGTTATAATGCTCACACTTTTacac                                            <  1:110109/251‑1 (MQ=255)
                                                                                                                                                                                                 gatTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGCCCACCTTGGTAGAAAAAAGGTTTTTTTTCGTTATAATGCTCACACTTTTaca                                             >  1:300117/2‑249 (MQ=255)
                                                                                                                                                                                                        cTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGCCCACCTTGGTAGAAAAAAGGTTTTTTTTCGTTATAATGCTCACACTTTTACACTAAAAATg                                    <  1:128822/251‑1 (MQ=255)
                                                                                                                                                                                                        cTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTAT‑AAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGCCCACCTTGGTAGAAAAAAGGTTTTTTTTCGTTATAATGCTCACACTTTTACACTAAAAATGa                                   >  1:542280/1‑251 (MQ=255)
                                                                                                                                                                                                          tttAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAAc                                                                                                                                                            <  1:595941/129‑1 (MQ=255)
                                                                                                                                                                                                          tttAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGCCCACCTTGGTAGAAAAAAGGTTTTTTTTCGTTATAATGCTCACACTTTTACACTAAAAATGAg                                  >  1:901555/1‑251 (MQ=255)
                                                                                                                                                                                                           ttAAAGGCAGCCGGGCTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGCCCACCTTGGTAGAAAAAAGGTTTTTTTTCGTTATAATGCTCACACTTTTACACTAAAAATGAGc                                 <  1:897907/251‑1 (MQ=255)
                                                                                                                                                                                                                          cTATAAATAAATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGCCCACCTTGGTAGAAAAAAGGTTTTTTTTCGTTATAATGCTCACACTTTTACACTAAAAATGAGCAAAATCCGTTCAttt                  >  1:292063/1‑251 (MQ=255)
                                                                                                                                                                                                                            ataaataaATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGccc                                                                                               >  1:378149/1‑172 (MQ=255)
                                                                                                                                                                                                                            ataaataaATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGCCCACCTTGGTAGAAAAAAGGTTTTTTTTCGTTATAATGCTCACACTTTTacac                                            >  1:777217/1‑223 (MQ=255)
                                                                                                                                                                                                                            ataaataaATTCCACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGCCCACCTTGGTAGAAAAAAGGTTTTTTTTCGTTATAATGCTCACACTTTTACACTAAAAATGAGCAAAATCCGTTCATTTTa                >  1:763156/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ccACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGc                                                                                                                                                          <  1:853533/102‑1 (MQ=255)
                                                                                                                                                                                                                                        cACTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGCCCACCTTGGTAGAAAAAAGGTTTTTTTTCGTTATAATGCTCACACTTTTACACTAAAAATGAGc                                 <  1:859503/222‑1 (MQ=255)
                                                                                                                                                                                                                                          cTTGCTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGCCCACCTTGGTAGAAAAAAGGTTTTTTTTCGTTATAATGCTCACACTTTTACACTAAAAATGAGCAAAATCCGTTCATTTTATCGCTATCCATATc  <  1:878733/251‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
TGGTGGCGCTTTTCTCCGCACCCAACAAATTGCGTTGCGCCCGATGCCAATCAATTAAAATTTTCACTTCCAGTTCCGGTTTATCCTGCTTCGCACAATAAATTTCATCCAGAATTTCTTGTCCCGCCTCGTCATTTTGCCAATAAAGTGCGGTCACAATAATCCGCGTTTTTGCTTGTCGAGTCAGCTGAATAATTTCGCTTTTAAAGGCAGCCGGGCTATAAATAAATTCCACTTGTTCCGGCTGTAATGGCAAACATGGCAGGCTTTCAAGACGTTTGCGCGCACGTTCGATTTTATCTGTAAACATAAATTAGTTCTCTATAAAAACGCTGTGAGTAAGAATGTAAATGAATTTTAATTAGCGCCGTAGTTTACAATATTTATCGCCCACCTTGGTAGAAAAAAGGTTTTTTTTCGTTATAATGCTCACACTTTTACACTAAAAATGAGCAAAATCCGTTCATTTTATCGCTATCCATATC  >  1/493466‑493950

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