Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 914,108 G→A Y124Y (TAC→TAT)  unknown ← Integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1914,1080GA100.0% 154.2 / NA 57Y124Y (TAC→TAT) unknownIntegrase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/30);  total (27/30)

CGGGGTTAGAACGCGTTCGCGCTTTTTGAACGTTGCAATAGATGAATTTGCGATTTCATCGGCTGGATTTGGCAGTTTATGTCCGCGCTGAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCGTAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCAGGAGAAATGCCTTCACTGACCTGCTTCCTTGCGATCATTAATTTTTCGCGCGCTTCTGCGAGATTTATTCCGTCGGCGCCGTATTTACCGATAGTTAATGTTTCGCGTCTACCGTTGATT  >  1/913862‑914354
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
cGGGGTTAGAACGCGTTCGCGCTTTTTGAACGTTGCAATAGATGAATTTGCGATTTCATCGGCTGGATTTGGCAGTTTATGTCCGCGCTGAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAg                                                                                                                                                                                                                                                    <  1:539611/251‑1 (MQ=255)
cGGGGTTAGAACGCGTTCGCGCTTTTTGAACGTTGCAATAGATGAATTTGCGATTTCATCGGCTGGATTTGGCAGTTTATGTCCGCGCTGAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAg                                                                                                                                                                                                                                                    <  1:578886/251‑1 (MQ=255)
cGGGGTTAGAACGCGTTCGCGCTTTTTGAACGTTGCAATAGATGAATTTGCGATTTCATCGGCTGGATTTGGCAGTTTATGTCCGCGCTGAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAg                                                                                                                                                                                                                                                    <  1:645151/251‑1 (MQ=255)
  gggTTAGAACGCGTTCGCGCTTTTTGAACGTTGCAATAGATGAATTTGCGATTTCATCGGCTGGATTTGGCAGTTTATGTCCGCGCTGAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTg                                                                                                                                                                                                                                                  >  1:697301/1‑251 (MQ=255)
  gggTTAGAACGCGTTCGCGCTTTTTGAACGTTGCAATAGATGAATTTGCGATTTCATCGGCTGGATTTGGCAGTTTATGTCCGCGCTGAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTg                                                                                                                                                                                                                                                  >  1:378174/1‑251 (MQ=255)
      tAGAACGCGTTCGCGCTTTTTGAACGTTGCAATAGATGAATTTGCGATTTCATCGGCTGGATTTGGCAGTTTATGTCCGCGCTGAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCa                                                                                                                                                                                                                                              >  1:828198/1‑251 (MQ=255)
                     cTTTTTGAACGTTGCAATAGATGAATTTGCGATTTCATCGGCTGGATTTGGCAGTTTATGTCCGCGCTGAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTg                                                                                                                                                                                                                               <  1:857433/251‑1 (MQ=255)
                        tttGAACGTTGCAATAGATGAATTTGCGATTTCATCGGCTGGATTTGGCAGTTTATGTCCGCGCTGAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTAc                                                                                                                                                                                                                            <  1:9928/251‑1 (MQ=255)
                        tttGAACGTTGCAATAGATGAATTTGCGATTTCATCGGCTGGATTTGGCAGTTTATGTCCGCGCTGAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTAc                                                                                                                                                                                                                            <  1:397955/251‑1 (MQ=255)
                        tttGAACGTTGCAATAGATGAATTTGCGATTTCATCGGCTGGATTTGGCAGTTTATGTCCGCGCTGAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTAc                                                                                                                                                                                                                            <  1:594267/251‑1 (MQ=255)
                                ttGCAATAGATGAATTTGCGATTTCATCGGCTGGATTTGGCAGTTTATGTCCGCGCTGAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGcc                                                                                                                                                                                                                    <  1:506311/251‑1 (MQ=255)
                                ttGCAATAGATGAATTTGCGATTTCATCGGCTGGATTTGGCAGTTTATGTCCGCGCTGAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGc                                                                                                                                                                                                                     >  1:75234/1‑250 (MQ=255)
                                   cAATAGATGAATTTGCGATTTCATCGGCTGGATTTGGCAGTTTATGTCCGCGCTGAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAg                                                                                                                                                                                                                                                    >  1:731675/1‑216 (MQ=255)
                                                          tCGGCTGGATTTGGCAGTTTATGTCCGCGCTGAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTc                                                                                                                                                                                          <  1:307631/251‑1 (MQ=255)
                                                          tCGGCTGGATTTGGCAGTTTATGTCCGCGCTGAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTc                                                                                                                                                                                          <  1:517811/251‑1 (MQ=255)
                                                           cGGCTGGATTTGGCAGTTTATGTCCGCGCTGAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAg                                                                                                                                                                                                                                                    <  1:799321/192‑1 (MQ=255)
                                                                          gTTTATGTCCGCGCTGAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAg                                                                                                                                                                          >  1:667450/1‑251 (MQ=255)
                                                                                         gAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAg                                                                                                                                                                                                                                                    >  1:590971/1‑162 (MQ=255)
                                                                                          aaTCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCgg                                                                                                                                                           <  1:655163/250‑1 (MQ=255)
                                                                                              gCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAgg                                                                                                                                                                                                  >  1:771401/1‑207 (MQ=255)
                                                                                               cATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAgg                                                                                                                                                     <  1:440842/251‑1 (MQ=255)
                                                                                                    cGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGCATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCttt                                                                                                                                                >  1:99927/1‑251 (MQ=255)
                                                                                                      gTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGc                                                                                                                                                             <  1:512179/236‑1 (MQ=255)
                                                                                                       tAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTg                                                                                                                                             <  1:45612/251‑1 (MQ=255)
                                                                                                        aGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAg                                                                                                                                                                                                                                                    >  1:137524/1‑147 (MQ=255)
                                                                                                                   aTTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATc                                                                                                                                                                     <  1:245183/215‑1 (MQ=255)
                                                                                                                                  aaaaTGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGCATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGaaaa                                                                                                                  >  1:41647/1‑250 (MQ=255)
                                                                                                                                   aaaTGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATg                                                                                                                 <  1:727993/251‑1 (MQ=255)
                                                                                                                                      tGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCt                                                                                                              <  1:394821/251‑1 (MQ=255)
                                                                                                                                            tGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGc                                                                                                                                                                  >  1:849756/1‑193 (MQ=255)
                                                                                                                                             ggAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGa                                                                                                       >  1:390824/1‑251 (MQ=255)
                                                                                                                                                aGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCCCTGAc                                                                                                      >  1:689406/1‑249 (MQ=255)
                                                                                                                                                     cACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAgg                                                                                                                                                                                                  >  1:447102/1‑152 (MQ=255)
                                                                                                                                                     cACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGt                                                                                                                                                                                                 <  1:857062/153‑1 (MQ=255)
                                                                                                                                                     cACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGACCTGCTTc                                                                                               <  1:34005/251‑1 (MQ=255)
                                                                                                                                                            aTCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAg                                                                                                                                                                                                                                                    >  1:700777/1‑95 (MQ=255)
                                                                                                                                                            aTCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGGGAAATGCCTTCACCGACCTGCTTTCtttct                                                                                         >  1:460723/1‑247 (MQ=255)
                                                                                                                                                            aTCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGACCTGCTTCCtttcgg                                                                                        >  1:116833/1‑247 (MQ=255)
                                                                                                                                                             tCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAg                                                                                                                                                                                                   <  1:275275/143‑1 (MQ=255)
                                                                                                                                                              cTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGAt                                                                                                                                           >  1:82385/1‑198 (MQ=255)
                                                                                                                                                                   aTTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGACCTGCTTcc                                                                                              >  1:478846/1‑238 (MQ=255)
                                                                                                                                                                       tATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGACCTGCTTCCTTGCGATCATTAAtttt                                                                             >  1:801969/1‑251 (MQ=255)
                                                                                                                                                                         ttGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAAAGCCTTCCCTGACCTGCTTCCTTGTGAATATTAAtttttt                                                                           >  1:396583/1‑250 (MQ=255)
                                                                                                                                                                                       gATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGACCTGCTTCCTTGCGATCATTAATTTTTCGCGCGCTTCTGCg                                                              <  1:620379/250‑1 (MQ=255)
                                                                                                                                                                                               ccGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGACCTGCTTCCTTGCGATCATTAATTTTTCGCGCGCTTCTGCGAGATTTAtt                                                     <  1:856667/251‑1 (MQ=255)
                                                                                                                                                                                                 gTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATg                                                                                                                 <  1:238452/189‑1 (MQ=255)
                                                                                                                                                                                                     tGATTTCGGACATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGACCTGCTTCCTTGCGATCATTAATTTTTCGCGCGCTTCTGCGAGATTTATTCCGTc                                                <  1:509733/250‑1 (MQ=255)
                                                                                                                                                                                                             gTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGACCTGCTTCCTTGCGATCATTAATTTTTCGCGCGCTTCTGCGAGATTTATTCCGTCGGCGCCGTa                                       >  1:538499/1‑251 (MQ=255)
                                                                                                                                                                                                                        aTTCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAg                                                                                                                                                                          >  1:298622/1‑109 (MQ=255)
                                                                                                                                                                                                                          tCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGcc                                                                                                               <  1:57942/166‑1 (MQ=255)
                                                                                                                                                                                                                          tCCCAAACGTGCCTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGACCTGCTTCCTTGCGATCATTAATTTTTCGCGCGCTTCTGCGAGATTTATTCCGTCGGCGCCGTATTTTCCGATAGtt                          >  1:492969/1‑251 (MQ=255)
                                                                                                                                                                                                                                      cTTTAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGACCTGCTTTCTTTCGATCATTAATTTTTCGCGCGCTTCTGCGAGATTTATTCCGTCGGCGCCGGATTTACCGATAg                            >  1:65576/1‑237 (MQ=255)
                                                                                                                                                                                                                                        ttAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGACCTGCTTCCTTGCGATCATTAATTTTTCGCGCGCTTCTGCgaga                                                           <  1:137041/204‑1 (MQ=255)
                                                                                                                                                                                                                                        ttAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGACCTGCTTCCTTGCGATCATTAATTTTTCGCGCGCTTCTGCGAGATTTATTCCGTCGGCGCCGTATTTACCGAt                              <  1:779187/233‑1 (MQ=255)
                                                                                                                                                                                                                                        ttAATATCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGACCTGCTTCCTTGCGATCATTAATTTTTCGCGCGCTTCTGCGAGATTTATTCCGTCGGCGCCGGATTTACCGATAGTTAATGTTTCGCGt              >  1:401494/1‑249 (MQ=255)
                                                                                                                                                                                                                                             atCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGACCTGCTTCCTTGCGATCATTAATTTTTCGCGCGCTTCTGCga                                                             <  1:599225/197‑1 (MQ=255)
                                                                                                                                                                                                                                             atCTCGTTCATAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGACCTGCTTCCTTGCGATCATTAATTTTTCGCGCGCTTCTGCGAGATTTATTCCGTCGGCGCCGTATTTACCGATAGTTAATGTTTCGCGTCTACCGt       <  1:281010/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 cgttcatAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGACCTGCTTCCTTGCGATCATTAATTTTTCGCGCGCTTCTGCGAGATTTATTCCGTCGGCGCCGGATTTACCGATAGTTAATGTTTCGCGTCTACCGTTGAt   >  1:701147/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  gttcatAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGACCTGCTTCCTTGCGATCATTAATTTTTCGCGCGCTTCTGCGAGATTTATTCCGTCGGCGCCGTATTTACCGATAGTTAATGTTTCGCGTCTACCGTTGAtt  >  1:829285/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  gttcatAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGACCTGCTTCCTTGCGATCATTAATTTTTCGCGCGCTTCTGCGAGATTTATTCCGTCGGCGCCGTATTTACCGATAGTTAATGTTTCGCGTCTACCGTTGAt   >  1:316508/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  gttcatAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAATGCCTTCACTGACCTGCTTCCTTGCGATCATTAATTTTTCGCGCGCTTCTGCGAGATTTATTCCGTCGGCGCCGTATTTACCGATAGTTAATGGTTTGCGTCTACCGTTGAtt  >  1:287486/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  gttcatAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCGGGAGAAAAGTCTTCACTGAACCGCTTCCTTGCGGTCATTAATTTTTCGCGCGCTTCTGCGAGAGTTATTTCGTCGGCGCCGTATTTACCGAGAGTTAATGTTTCGCGTCTACCGTTGAtt  >  1:487766/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
CGGGGTTAGAACGCGTTCGCGCTTTTTGAACGTTGCAATAGATGAATTTGCGATTTCATCGGCTGGATTTGGCAGTTTATGTCCGCGCTGAATCGCATATCGGTAGATATTCGCGATTAAATCCCGTACGAAAATGGCCGTGGAAGGCGCACCGCGATCTTTAATTTTATTGCAATGGCGGCGGATTTCATCCGTAGTGATTTCGGTCATTAAGCGATTCCCAAACGTGCCTTTAATATCTCGTTCGTAAGTGGCCACGCGGAGGGCTTTTGTACTTTCTGCCAGTTGCACCTCCTCAAGGTATTTTTCAGCGAAAGCGGAAAAGCGATCTGCATTGCGGATAAGGTCTTTTTGATTACGCTTTTCACCGGCAGGAGAAATGCCTTCACTGACCTGCTTCCTTGCGATCATTAATTTTTCGCGCGCTTCTGCGAGATTTATTCCGTCGGCGCCGTATTTACCGATAGTTAATGTTTCGCGTCTACCGTTGATT  >  1/913862‑914354

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