Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 566,288 T→C intergenic (+68735/‑133055) 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
*1566,2880TC100.0% 186.5 / NA 68intergenic (+68735/‑133055)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 (40/28);  total (40/28)

CGACAACACTTCCCCCAAGGGCGCGTGATAATACCCCGCCGCCCACTGCAAAAACTGCCAGCTTTGCGGAGAAAATAAGGAAACCTCATCAATGCAAGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAAGATATAATCAAAATAGTGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATAAGGGATAAAATACTGTCAAATTTAAGAAGCAAGAAATGCGTACAGGATGAATGAAAAAGCATCTTAGCTAAGAAGCAAATATGAAAATTTATATT  >  1/566041‑566533
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
cGACAACACTTCCCCCAAGGGCGCGTGATAATACCCCGCCGCCCACTGCAAAAACTGCCAGCTTTGCGGAGAAAATAAGGAAACCTCATCAATGCAAGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGAt                                                                                                                                                                                                                                                    >  1:59698/1‑251 (MQ=255)
     acacTTCCCCCAAGGGCGCGTGATAATACCCCGCCGCCCACTGCAAAAACTGCCAGCTTTGCGGAGAAAATAAGGAAACCTCATCAATGCAAGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAAc                                                                                                                                                                                                                                               <  1:531304/251‑1 (MQ=255)
           cccccAAGGGCGCGTGATAATACCCCGCCGCCCACTGCAAAAACTGCCAGCTTTGCGGAGAAAATAAGGAAACCTCATCAATGCAAGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGcaca                                                                                                                                                                                                                                         >  1:867413/1‑251 (MQ=255)
              ccAAGGGCGCGTGATAATACCCCGCCGCCCACTGCAAAAACTGCCAGCTTTGCGGAGAAAATAAGGAAACCTCATCAATGCAAGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGa                                                                                                                                                                                                                                      <  1:599089/251‑1 (MQ=255)
               cAAGGGCGCGTGATAATACCCCGCCGCCCACTGCAAAAACTGCCAGCTTTGCGGAGAAAATAAGGAAACCTCATCAATGCAAGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGaa                                                                                                                                                                                                                                     <  1:471265/251‑1 (MQ=255)
                         tGATAATACCCCGCCGCCCACTGCAAAAACTGCCAGCTTTGCGGAGAAAATAAGGAAACCTCATCAATGCAAGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCg                                                                                                                                                                                                                           >  1:892435/1‑251 (MQ=255)
                         tGATAATACCCCGCCGCCCACTGCAAAAACTGCCAGCTTTGCGGAGAAAATAAGGAAACCTCATCAATGCAAGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCg                                                                                                                                                                                                                           >  1:780210/1‑251 (MQ=255)
                         tGATAATACCCCGCCGCCCACTGCAAAAACTGCCAGCTTTGCGGAGAAAATAAGGAAACCTCATCAATGCAAGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCg                                                                                                                                                                                                                           >  1:245259/1‑251 (MQ=255)
                           ataataCCCCGCCGCCCACTGCAAAAACTGCCAGCTTTGCGGAGAAAATAAGGAAACCTCATCAATGCAAGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCggg                                                                                                                                                                                                                         >  1:311922/1‑251 (MQ=255)
                                         cccACTGCAAAAACTGCCAGCTTTGCGGAGAAAATAAGGAAACCTCATCAATGCAAGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGaa                                                                                                                                                                                                            <  1:265747/250‑1 (MQ=255)
                                                                tGCGGAGAAAATAAGGAAACCTCATCAATGCAAGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATgg                                                                                                                                                                                    <  1:250455/251‑1 (MQ=255)
                                                                            aaGGAAACCTCATCAATGCAAGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCGCCCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATaa                                                                                                                                                                                         >  1:285116/1‑235 (MQ=255)
                                                                               gAAACCTCATCAATGCAAGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGc                                                                                                                                                                      <  1:384687/250‑1 (MQ=255)
                                                                               gAAACCACATCAATGCAAGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGc                                                                                                                                                                      <  1:196446/250‑1 (MQ=255)
                                                                                   ccTCATCAATGCAAGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATggg                                                                                                                                                                                   <  1:283435/233‑1 (MQ=255)
                                                                                       atcaATGCAAGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGa                                                                                                                                                                                                                                      <  1:361616/178‑1 (MQ=255)
                                                                                                aGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGaa                                                                                                                                                                                                            >  1:565964/1‑195 (MQ=255)
                                                                                                aGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAat                                                                                                                                                                           >  1:151769/1‑228 (MQ=255)
                                                                                                aGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATtgttgt                                                                                                                                                    >  1:195368/1‑251 (MQ=255)
                                                                                                aGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATtgttgt                                                                                                                                                    >  1:24239/1‑251 (MQ=255)
                                                                                                aGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATtgttgt                                                                                                                                                    >  1:473478/1‑251 (MQ=255)
                                                                                                         tGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGTACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTggg                                                                                                                                           <  1:437130/251‑1 (MQ=255)
                                                                                                          gTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAg                                                                                                                                                                                                 >  1:726495/1‑196 (MQ=255)
                                                                                                          gTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTgggg                                                                                                                                          >  1:222231/1‑251 (MQ=255)
                                                                                                          gTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTgggg                                                                                                                                          >  1:393547/1‑251 (MQ=255)
                                                                                                          gTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTgggg                                                                                                                                          >  1:131291/1‑251 (MQ=255)
                                                                                                          gTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTgggg                                                                                                                                          >  1:879609/1‑251 (MQ=255)
                                                                                                            cTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGCTTGTTGTCTATTTGGGGga                                                                                                                                        >  1:494790/1‑251 (MQ=255)
                                                                                                               cAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGaa                                                                                                                                     >  1:667678/1‑251 (MQ=255)
                                                                                                                     cTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGctc                                                                                                                               <  1:510697/251‑1 (MQ=255)
                                                                                                                     cTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGctc                                                                                                                               <  1:614531/251‑1 (MQ=255)
                                                                                                                              cacaTCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCt                                                                                                                                                  <  1:806557/223‑1 (MQ=255)
                                                                                                                               acaTCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAAt                                                                                                                                    <  1:889957/236‑1 (MQ=255)
                                                                                                                                caTCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTggg                                                                                                                    >  1:723524/1‑251 (MQ=255)
                                                                                                                                  tCCGACTGCTCCGCCCCCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCGCGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCtgt                                                                                                                            <  1:381330/241‑1 (MQ=255)
                                                                                                                                  tCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCAt                                                                                                                                                                    >  1:865053/1‑201 (MQ=255)
                                                                                                                                                 aCCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGc                                                                                                                                                                      <  1:503640/184‑1 (MQ=255)
                                                                                                                                                 aCCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGctcc                                                                                                                              >  1:303116/1‑223 (MQ=255)
                                                                                                                                                  ccTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGGCGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGAt                                                                                                  <  1:44698/251‑1 (MQ=255)
                                                                                                                                                  ccTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGctc                                                                                                                               >  1:329096/1‑222 (MQ=255)
                                                                                                                                                   cTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACGCCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACAGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGc                                                                                                                                                                      <  1:537142/182‑1 (MQ=255)
                                                                                                                                                   cTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATa                                                                                                 >  1:786387/1‑251 (MQ=255)
                                                                                                                                                       accaccATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATAAggg                                                                                             <  1:106555/251‑1 (MQ=255)
                                                                                                                                                        ccaccaTGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATAAggg                                                                                             <  1:841895/250‑1 (MQ=255)
                                                                                                                                                          accaTGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATAAGGGATa                                                                                          >  1:169291/1‑251 (MQ=255)
                                                                                                                                                              tGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTAtt                                                                                                                                               <  1:530946/194‑1 (MQ=255)
                                                                                                                                                              tGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATAAGGGATAAAAt                                                                                      <  1:489436/251‑1 (MQ=255)
                                                                                                                                                               gCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGc                                                                                                                                                                      >  1:405886/1‑170 (MQ=255)
                                                                                                                                                                        cTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGt                                                                                                                  >  1:23511/1‑213 (MQ=255)
                                                                                                                                                                              cgcCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATAAGGGATAAAATACTGTCAAATTTaaga                                                                      <  1:382399/251‑1 (MQ=255)
                                                                                                                                                                                ccccAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTg                                                                                                                                             >  1:303591/1‑178 (MQ=255)
                                                                                                                                                                                      aCGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTg                                                                                                                                             >  1:375051/1‑172 (MQ=255)
                                                                                                                                                                                              cacaCCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATAAGGGATAAAATACTGTCAAATTTAAGAAGCAAGAAATGCGTAc                                                      <  1:416800/251‑1 (MQ=255)
                                                                                                                                                                                                 acCCTCG‑CCCCCACGACAGCATTCAGAGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATAAGGGATAAAATACTGTCAAATTTAAGAAGCAAGAAATGCGTACAgg                                                   <  1:804732/251‑1 (MQ=255)
                                                                                                                                                                                                  cccTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAAt                                                                                                             >  1:698429/1‑192 (MQ=255)
                                                                                                                                                                                                  cccTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATAAggg                                                                                             >  1:4908/1‑208 (MQ=255)
                                                                                                                                                                                                           cccACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATAAGGGATAAAATACTGTCAAATTTAAGAAGCAAGAAATGCGTACAGGAtgaatgaa                                          >  1:22450/1‑251 (MQ=255)
                                                                                                                                                                                                            ccACGACAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTggggg                                                                                                                                         >  1:330281/1‑155 (MQ=255)
                                                                                                                                                                                                                 aCAGCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTggg                                                                                                                                           >  1:392846/1‑148 (MQ=255)
                                                                                                                                                                                                                    gCATTCAGCGACGCCGGCAGGATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATAAGGGATAAAAt                                                                                      >  1:396367/1‑198 (MQ=255)
                                                                                                                                                                                                                                       ggATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGTATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGAt                                                                                                  >  1:696008/2‑167 (MQ=255)
                                                                                                                                                                                                                                       ggATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATAAGGGATAAAATAc                                                                                    >  1:657843/2‑181 (MQ=255)
                                                                                                                                                                                                                                       ggATATAATCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATAAGGGATAAAATACTGTCAAATTTAAGAAGCAAGAAATGCGTACAGGATGAATGAAAAAGCATCTTAGCTAAGAAGCAAATATg              >  1:728888/2‑251 (MQ=255)
                                                                                                                                                                                                                                             aaTCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATAAGGGATAAAATACTGTCAAATTTAAGAAGCAAGAAATGCGTACAGGATGAATGAAAAAGCATCTTAGCTaagaag                      <  1:591633/237‑1 (MQ=255)
                                                                                                                                                                                                                                             aaTCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATAAGGGATAAAATACTGTCAAATTTAAGAAGCAAGAAATGCGTACAGGATGAATGAAAAAGCATCTTAGCTAAGAAGCAAATATGAAAAtt        >  1:717303/1‑251 (MQ=255)
                                                                                                                                                                                                                                              aTCAAAATAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATAAGGGATAAAAt                                                                                      <  1:609174/172‑1 (MQ=255)
                                                                                                                                                                                                                                                  aaaTAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATAAGGGATAAAATACTGTCAAATTTAAGAAGCAAGAAATGCGTACAGGATGAATGAAAAAGCATCTTAGCTAAGAAGCAAATATGAAAATTTatat   <  1:434292/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   aaTAGCGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATAAGGGATAAAATACTGTCAAATTTAAGAAGCAAGAAATGCGTACAGGATGAATGAAAAAGCATCTTAGCTAAGAAGCAAATATGAAAATTTATAtt  >  1:220226/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
CGACAACACTTCCCCCAAGGGCGCGTGATAATACCCCGCCGCCCACTGCAAAAACTGCCAGCTTTGCGGAGAAAATAAGGAAACCTCATCAATGCAAGCCAGCACTGTCTTCAACTGCTCCGCAGGCACATCCGACTGCTCCGCCACCTCAACCACCATGCCAACACGCTTCTGCGCCCCAAACGGCACCCACACCCTCG‑CCCCCACGACAGCATTCAGCGACGCCGGCAAGATATAATCAAAATAGTGATGTAACGGCACAGGAAGCGCAATTCGGGCAAATGACATGAACGGAAATTCAGTGAAATAAAATGGGCTCGGAATATAGCATTTTGGGGGATTGTTGTCTATTTGGGGGAGAATAGCTCTGTTGTTTGGGGTTTAATTTCACTTGGATAAGGGATAAAATACTGTCAAATTTAAGAAGCAAGAAATGCGTACAGGATGAATGAAAAAGCATCTTAGCTAAGAAGCAAATATGAAAATTTATATT  >  1/566041‑566533

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