Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,478,786 C→T S26F (TCT→TTT)  unknown → Type IV prepilin peptidase TadV/CpaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,478,7860CT100.0% 256.2 / NA 88S26F (TCT→TTT) unknownType IV prepilin peptidase TadV/CpaA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (55/33);  total (55/33)

CTTCTGTTGAGTATGGGCTAATTGCTGTGGCTATTGCTGTTTTTGTAGTCGCCGTTTTAGTCGAGGATAACAGCTTTATTAAGGCGTTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCTCTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGATAAAAAAATGAACTGGCTTATTAATACTTTAGTAGTCATAACGATCTTACTGCTGACTACTTTATCTGTAACTGACATTCGTTCTCGGATTATTAGTAACCGAGTCGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGACTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGATTTTTACTGTTTTCTTTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTAATGCTAGCCGTTCCAAGTGACGAAATTTTTTCTTTTTTCTTCTTTACAACATTTGCCGGATTGGGAC  >  1/1478542‑1479028
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
cTTCTGTTGAGTATGGGTTAATTGCTGTGGCTATTGCTGTTTTTGTAGTCGCCGTTTTAGTCGGGGATAACAGCTTTATTAAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCggctt                                                                                                                                                                                                                                             >  1:674379/1‑248 (MQ=255)
  tCTGTTGAGTATGGGTTAATTGCTGTGGCTATTGCTGTTTTTGTAGTCGCCGTTTTAGTCGGGGATAACAGCTTTATTAAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGACCTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGcttat                                                                                                                                                                                                                                           >  1:430098/1‑251 (MQ=255)
          gTATGGGTTAATTGCTGTGGCTATTGCTGTTTTTGTAGTCGCCGTTTTAGTCGGGGATAACAGCTTTATTAAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAAc                                                                                                                                                                                                                                    >  1:359386/1‑250 (MQ=255)
           tATGGGTTAATTGCTGTGGCTATTGCTGTTTTTGTAGTCGCCGTTTTAGTCGGGGATAACAGCTTTATTAAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAAc                                                                                                                                                                                                                                    <  1:724868/249‑1 (MQ=255)
               ggttAATTGCTGTGGCTATTGCTGTTTTTGTAGTCGCCGTTTTAGTCGGGGATAACAGCTTTATTAAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGaatt                                                                                                                                                                                                                              >  1:779923/4‑248 (MQ=255)
                                      gTTTTTGTAGTCGCCGTTTTAGTCGGGGATAACAGCTTTATTAAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCat                                                                                                                                                                                                       >  1:85047/1‑251 (MQ=255)
                                      gTTTTTGTAGTCGCCGTTTTAGTCGGGGATAACAGCTTTATTAAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCat                                                                                                                                                                                                       >  1:261506/1‑251 (MQ=255)
                                                     gTTTTAGTCGGGGATAACAGCTTTATTAAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGttt                                                                                                                                                                                        >  1:725234/1‑251 (MQ=255)
                                                        ttAGTCGGGGATAACAGCTTTATTAAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTttatt                                                                                                                                                                                     <  1:843970/251‑1 (MQ=255)
                                                              ggggATAACAGCTTTATTAAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTaaa                                                                                                                                                                                  >  1:572688/3‑248 (MQ=255)
                                                               gggATAACAGCTTTATTAAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTTGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAatat                                                                                                                                                                              >  1:750316/2‑251 (MQ=255)
                                                               gggATAACAGCTTTATTAAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAatat                                                                                                                                                                              >  1:394607/2‑251 (MQ=255)
                                                               gggATAACAGCTTTATTAAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAatat                                                                                                                                                                              >  1:104107/2‑251 (MQ=255)
                                                                 gATAACAGCTTTATTAAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTc                                                                                                                                                                                                            >  1:367751/1‑219 (MQ=255)
                                                                    aaCAGCTTTATTAAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAggc                                                                                                                                                                         <  1:252038/251‑3 (MQ=255)
                                                                           ttattaAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTttatt                                                                                                                                                                                     >  1:520307/1‑232 (MQ=255)
                                                                            tattaAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAg                                                                                                                                                                                                                                        >  1:346725/1‑180 (MQ=255)
                                                                            tattaAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGaattg                                                                                                                                                                                                                             >  1:713371/1‑187 (MQ=255)
                                                                            tattaAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGt                                                                                                                                                                                                          >  1:312297/1‑210 (MQ=255)
                                                                            tattaAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTg                                                                                                                                                                 >  1:156350/1‑251 (MQ=255)
                                                                            tattaAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTg                                                                                                                                                                 >  1:866773/1‑251 (MQ=255)
                                                                            gattaAGGCAGTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTg                                                                                                                                                                 >  1:18966/2‑251 (MQ=255)
                                                                                      gtGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTccc                                                                                                                                                         <  1:463212/248‑1 (MQ=255)
                                                                                            ggTAAGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATg                                                                                                                                                 >  1:578120/1‑251 (MQ=255)
                                                                                               aaGTTCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTac                                                                                                                                              <  1:8210/251‑1 (MQ=255)
                                                                                                  ttCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATATCGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTacac                                                                                                                                            <  1:269642/250‑1 (MQ=255)
                                                                                                  ttCTCGGATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTacac                                                                                                                                            <  1:508286/250‑1 (MQ=255)
                                                                                                        gATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTttatta                                                                                                                                                                                    <  1:581467/204‑1 (MQ=255)
                                                                                                        gATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTc                                                                                                                                                           >  1:422692/1‑229 (MQ=255)
                                                                                                        gATTTAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTtattat                                                                                                                                     >  1:32452/1‑251 (MQ=255)
                                                                                                          tttAACGACCACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCg                                                                                                                                   >  1:287989/1‑251 (MQ=255)
                                                                                                                  ccACTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTc                                                                                                                                                                                             >  1:296583/1‑185 (MQ=255)
                                                                                                                     cTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTc                                                                                                                                                                                                            >  1:197941/1‑167 (MQ=255)
                                                                                                                     cTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAata                                                                                                                                                                               >  1:755739/1‑196 (MQ=255)
                                                                                                                     cTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATc                                                                                                                                    <  1:235106/239‑1 (MQ=255)
                                                                                                                     cTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTa                                                                                                                           >  1:23506/1‑248 (MQ=255)
                                                                                                                     cTGTCACTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTa                                                                                                                           >  1:124819/1‑248 (MQ=255)
                                                                                                                         cactGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTaa                                                                                                                                                                                   <  1:288495/186‑1 (MQ=255)
                                                                                                                          actGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCG‑TCAGTTTATTAAAATATCAGGCTATCTTTGTTATTc                                                                                                                                                           >  1:562427/2‑210 (MQ=255)
                                                                                                                               cGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTa                                                                                                                                                              <  1:144136/203‑1 (MQ=255)
                                                                                                                                      gTTAGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCAt                                                                                                       >  1:573738/1‑251 (MQ=255)
                                                                                                                                         aGCAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCAt                                                                                                       <  1:640233/248‑1 (MQ=255)
                                                                                                                                           cAAAAGTAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTg                                                                                                  <  1:81214/251‑1 (MQ=255)
                                                                                                                                                 tAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGgag                                                                                            <  1:879745/251‑1 (MQ=255)
                                                                                                                                                 tAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGgag                                                                                            <  1:319394/251‑1 (MQ=255)
                                                                                                                                                 tAGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGgag                                                                                            <  1:689030/251‑1 (MQ=255)
                                                                                                                                                  aGCTAAGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGgaga                                                                                           >  1:28355/1‑251 (MQ=255)
                                                                                                                                                      aaGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTc                                                                                                                                                           >  1:609532/1‑183 (MQ=255)
                                                                                                                                                      aaGAAAAATGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTacac                                                                                                                                            >  1:186427/1‑198 (MQ=255)
                                                                                                                                                           aaaTGAT‑AAAAAATGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAgct                                                                                  <  1:591419/251‑1 (MQ=255)
                                                                                                                                                                gataaaaaaTGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGc                                                                                                              >  1:413250/4‑219 (MQ=255)
                                                                                                                                                                     aaaaTGAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGt                                                                         >  1:901598/1‑251 (MQ=255)
                                                                                                                                                                          gAACTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTaa                                                                    >  1:279739/1‑251 (MQ=255)
                                                                                                                                                                            aCTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTAATg                                                                  <  1:872160/251‑1 (MQ=255)
                                                                                                                                                                            aCTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTAATg                                                                  <  1:153056/251‑1 (MQ=255)
                                                                                                                                                                            aCTGGCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTAATg                                                                  <  1:239420/251‑1 (MQ=255)
                                                                                                                                                                               ggCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCt                                                                                                                                                                    <  1:787907/150‑1 (MQ=255)
                                                                                                                                                                               ggCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTa                                                                                                                                                              <  1:261526/156‑1 (MQ=255)
                                                                                                                                                                               ggCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATc                                                                                                                                    >  1:804572/1‑182 (MQ=255)
                                                                                                                                                                               ggCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTAATgtta                                                               >  1:642253/1‑248 (MQ=255)
                                                                                                                                                                                gCTTATTAACGCTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTCCACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTAATGTTAg                                                              >  1:454836/1‑251 (MQ=255)
                                                                                                                                                                                           cTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTacac                                                                                                                                            <  1:749089/162‑1 (MQ=255)
                                                                                                                                                                                           cTTTAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTAATGTTAGCCGTTCCAAGt                                                   >  1:319621/1‑251 (MQ=255)
                                                                                                                                                                                             ttAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAgc                                                                                   <  1:156002/217‑1 (MQ=255)
                                                                                                                                                                                             ttAGTAATCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTAATGTTAGCCGTTCCAAGTGa                                                 >  1:838441/1‑251 (MQ=255)
                                                                                                                                                                                               agtaatCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAgg                                                                                                    <  1:592526/198‑1 (MQ=255)
                                                                                                                                                                                                  aatCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTacac                                                                                                                                            <  1:621817/153‑1 (MQ=255)
                                                                                                                                                                                                    tCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCt                                                                                                                                                                    >  1:160973/1‑129 (MQ=255)
                                                                                                                                                                                                    tCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTac                                                                                                                                              <  1:596611/151‑1 (MQ=255)
                                                                                                                                                                                                    tCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGgag                                                                                            <  1:491878/201‑1 (MQ=255)
                                                                                                                                                                                                    tCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTAATGTTAGCCGTTCCAAGTGAGGAAAtt                                          >  1:90960/1‑251 (MQ=255)
                                                                                                                                                                                                    tCATAACGCTCTCGCTGCTGGTTACTTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTAATGTTAGCCGTTCCAAGTGAGGAAAtt                                          >  1:288435/1‑251 (MQ=255)
                                                                                                                                                                                                                cgctgctggttactTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTAATGTTAGCCGTTCCAAGt                                                   <  1:678400/220‑1 (MQ=255)
                                                                                                                                                                                                                  ctgctggttactTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTaaaa                                                                                                                                                                                 >  1:881162/9‑102 (MQ=255)
                                                                                                                                                                                                                    gctggttactTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATAtt                                                                                       >  1:289766/7‑190 (MQ=255)
                                                                                                                                                                                                                    gctggttactTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTAATg                                                                  >  1:242086/7‑211 (MQ=255)
                                                                                                                                                                                                                     ctggttactTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGCTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTAATGTTAGCCGTTCCAAGTGAGGAAATTTTTTCTTTTTTcttctt                         >  1:845328/6‑251 (MQ=255)
                                                                                                                                                                                                                       ggttactTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTAATGTTAGCCGTTCCAAGTGAGGAAATTTTTTCTTTTTTCTTCTTTa                       >  1:525459/4‑251 (MQ=255)
                                                                                                                                                                                                                           actTTATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTac                                                                                                                                              >  1:720842/1‑128 (MQ=255)
                                                                                                                                                                                                                              ttATCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTAATGTTAGCCGTTCCAAGTGAGGAAATTTTTTCTTTTTTCTTCTTTACAACAtt                >  1:898579/1‑251 (MQ=255)
                                                                                                                                                                                                                                aTCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATTTCAGGCTATCTTTTTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGc                                                                                                              <  1:324953/155‑1 (MQ=255)
                                                                                                                                                                                                                                aTCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATc                                                                                                                                    <  1:146892/133‑1 (MQ=255)
                                                                                                                                                                                                                                aTCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTAATGTTAGc                                                             <  1:895121/204‑1 (MQ=255)
                                                                                                                                                                                                                                aTCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTAATGTTAGCCGTTCCAAGTGAGGAAATTTTTTCTTTTTTCTTCTTTACAACATTTg              >  1:743509/1‑251 (MQ=255)
                                                                                                                                                                                                                                aTCTATAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAAGAGATATTAAGCTTGCTTCGGTGCTAATGTTAGCCGTTCCAAGTGAGGAAATTTTTTCTTTTTTCTTCTTTACAACAt                 >  1:440929/1‑248 (MQ=255)
                                                                                                                                                                                                                                     tAACTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGATTATTAAAATATCAGGCTATCTTTGTTATTc                                                                                                                                                           >  1:351632/1‑105 (MQ=255)
                                                                                                                                                                                                                                       aCTGACATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAgg                                                                                                    >  1:342626/1‑158 (MQ=255)
                                                                                                                                                                                                                                            cATTCGTTTTCGGCTTATTAGTAACCGAATTGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGGCTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGGTTTTTACTGTTTTCTCTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTAATGTTAGCCGTTCCAAGTGAGGAAATTTTTTCTTTTTTCTTCTTTACAACATTTGCCGGATTGGGAc  <  1:291072/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
CTTCTGTTGAGTATGGGCTAATTGCTGTGGCTATTGCTGTTTTTGTAGTCGCCGTTTTAGTCGAGGATAACAGCTTTATTAAGGCGTTGTCAGGTAAGTTCTCGGATTTAACGACCACTGTCTCTGACGCCATAGTTAGCAAAAGTAGCTAAGAAAAATGATAAAAAAATGAACTGGCTTATTAATACTTTAGTAGTCATAACGATCTTACTGCTGACTACTTTATCTGTAACTGACATTCGTTCTCGGATTATTAGTAACCGAGTCGTTTTGTTTCTGTTATTCGTCATTATTCCGTTCAGTTTATTAAAATATCAGACTATCTTTGTTATTCCCGCCTTATGTACACTTATTATCGGATTTTTACTGTTTTCTTTGCACTTCATAGGTGCAGGAGATATTAAGCTTGCTTCGGTGCTAATGCTAGCCGTTCCAAGTGACGAAATTTTTTCTTTTTTCTTCTTTACAACATTTGCCGGATTGGGAC  >  1/1478542‑1479028

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