Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 533,647 G→T Q202K (CAA→AAA)  unknown ← Cytochrome c551 peroxidase (EC 1.11.1.5)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1533,6470GT100.0% 186.3 / NA 66Q202K (CAA→AAA) unknownCytochrome c551 peroxidase (EC 1.11.1.5)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (33/33);  total (33/33)

TCCTGCTCAAATCGTGCCACAATTTCTTCCCAGGAATGGGAACCCATTTCTACCGGATTTAATGGCGGTCCGCCGGCTTGATCGGCTAAATCCACGGCACGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTGAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTGCCGGTAAAAAGGCTTTCCGTTGTTCGCGAATCCAGTTTAGCAATGCGGTTTTTTCATCTTCATCAGGGCGGCTTCCCCAGTGAATATGAAGAAACTTG  >  1/533406‑533876
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     
tCCTGCTCAAATCGTGCCACAATTTCTTCCCAGGAATGGGAACCCATTTCTACCGGATTTAATGGCGGTCCGCCGGCTTGATCGGCTAAATCCACGGCACGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTAcc                                                                                                                                                                                                                              <  1:710042/251‑1 (MQ=255)
    gCTCAAATCGTGCCACAATTTCTTCCCAGGAATGGGAACCCATTTCTACCGGATTTAATGGCGGTCCGCCGGCTTGATCGGCTAAATCCACGGCACGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGGGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGGGTGACATTTAATCGTACCATCg                                                                                                                                                                                                                          >  1:42259/1‑251 (MQ=255)
                ccACAATTTCTTCCCAGGAATGGGAACCCATTTCTACCGGATTTAATGGCGGTCCGCCGGCTTGATCGGCTAAATCCACGGCACGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGa                                                                                                                                                                                                              <  1:543397/251‑1 (MQ=255)
                 cacaATTTCTTCCCAGGAATGGGAACCCATTTCTACCGGATTTAATGGCGGTCCGCCGGCTTGATCGGCTAAATCCACGGCACGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTTGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGa                                                                                                                                                                                                              >  1:651306/1‑250 (MQ=255)
                 cacaATTTCTTCCCAGGAATGGGAACCCATTTCTACCGGATTTAATGGCGGTCCGCCGGCTTGATCGGCTAAATCCACGGCACGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGAc                                                                                                                                                                                                             <  1:329163/251‑1 (MQ=255)
                            cccAGGAATGGGAACCCATTTCTACCGGATTTAATGGCGGTCCGCCGGCTTGATCGGCTAAATCCACGGCACGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAAc                                                                                                                                                                                                  <  1:648062/251‑1 (MQ=255)
                            cccAGGAATGGGAACCCATTTCTACCGGATTTAATGGCGGTCCGCCGGCTTGATCGGCTAAATCCACGGCACGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAAc                                                                                                                                                                                                  <  1:132399/251‑1 (MQ=255)
                                       gAACCCATTTCTACCGGATTTAATGGCGGTCCGCCGGCTTGATCGGCTAAATCCACGGCACGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAg                                                                                                                                                                                                <  1:552178/242‑1 (MQ=255)
                                         aCCCATTTCTACCGGATTTAATGGCGGTCCGCCGGCTTGATCGGCTAAATCCACGGCACGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGcc                                                                                                                                                                                        <  1:809437/248‑1 (MQ=255)
                                            cATTTCTACCGGATTTAATGGCGGTCCGCCGGCTTGATCGGCTAAATCCACGGCACGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAagag                                                                                                                                                                                  >  1:546815/1‑251 (MQ=255)
                                                                      cgccgGCTTGATCGGCTAAATCCACGGCACGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGAtt                                                                                                                                                        <  1:71228/251‑1 (MQ=255)
                                                                          ggCTTGATCGGCTAAATCCACGGCACGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATc                                                                                                                                                                                                                                   <  1:825891/172‑1 (MQ=255)
                                                                           gCTTGATCGGCTAAATCCACGGCACGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTtcat                                                                                                                                                    >  1:447781/1‑250 (MQ=255)
                                                                                         aTCCACGGCACGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGaa                                                                                                                                     >  1:285702/1‑251 (MQ=255)
                                                                                           ccACGGCACGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAAt                                                                                                                                             <  1:352103/241‑1 (MQ=255)
                                                                                               ggCACGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGAc                                                                                                                                                                                                             >  1:732197/1‑173 (MQ=255)
                                                                                                gCACGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAAc                                                                                                                               <  1:307381/250‑1 (MQ=255)
                                                                                                   cGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGc                                                                                                                                                                                         <  1:879858/189‑1 (MQ=255)
                                                                                                         tCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAAc                                                                                                                                                                                                  >  1:810586/1‑174 (MQ=255)
                                                                                                         tCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAAc                                                                                                                                                                                                  <  1:232020/174‑1 (MQ=255)
                                                                                                         tCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAAc                                                                                                                                                                                                  <  1:736916/174‑1 (MQ=255)
                                                                                                         tCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAg                                                                                                                                                                                                >  1:878020/1‑176 (MQ=255)
                                                                                                         tCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGAtt                                                                                                                                                        >  1:12316/1‑216 (MQ=255)
                                                                                                         tCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGAt                                                                                                                                                         <  1:77148/215‑1 (MQ=255)
                                                                                                         tCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAAc                                                                                                                               >  1:515248/1‑241 (MQ=255)
                                                                                                         tCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACgg                                                                                                                             <  1:59328/243‑1 (MQ=255)
                                                                                                         tCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGt                                                                                                                      >  1:354533/1‑250 (MQ=255)
                                                                                                         tCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGt                                                                                                                      >  1:459750/1‑250 (MQ=255)
                                                                                                         tCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGt                                                                                                                      >  1:717353/1‑250 (MQ=255)
                                                                                                         tCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGt                                                                                                                      >  1:43373/1‑250 (MQ=255)
                                                                                                         tCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGt                                                                                                                      >  1:399594/1‑250 (MQ=255)
                                                                                                         tCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTc                                                                                                                     <  1:474886/251‑1 (MQ=255)
                                                                                                            cAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGAtt                                                                                                                                                        >  1:156834/1‑213 (MQ=255)
                                                                                                              gAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTg                                                                                                                 >  1:221043/1‑250 (MQ=255)
                                                                                                                        aaGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAAc                                                                                                                                                                                                  <  1:235021/159‑1 (MQ=255)
                                                                                                                          gTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTAttta                                                                                                    >  1:277070/1‑250 (MQ=255)
                                                                                                                          gTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTAttta                                                                                                    >  1:403005/1‑250 (MQ=255)
                                                                                                                            tAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTAtttacc                                                                                                  <  1:132894/251‑4 (MQ=255)
                                                                                                                                    gcgTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAAc                                                                                                                                                                                                  <  1:122514/147‑1 (MQ=255)
                                                                                                                                    gcgTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTaaaaa                                                                                          >  1:702771/1‑251 (MQ=255)
                                                                                                                                      gTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTAAAAAgg                                                                                        >  1:526352/1‑251 (MQ=255)
                                                                                                                                                    gcgcGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTAAAAAGGCTTTCCGTTGTTcg                                                                          >  1:689399/1‑251 (MQ=255)
                                                                                                                                                      gcgTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTAAAAAGGCTTTCCGTTGTTcgcg                                                                        >  1:895770/1‑251 (MQ=255)
                                                                                                                                                                   cTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTAAAAAGGCTTTCCGTTGTTCGCGAATCCAGTTTAGc                                                           <  1:791587/251‑1 (MQ=255)
                                                                                                                                                                          gCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGCGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTAAAAAGGCTTTCCGTTGTTCGCGAATCCAGTTTAGCAATGCgg                                                    >  1:638047/1‑251 (MQ=255)
                                                                                                                                                                                    ccGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTc                                                                                                                     >  1:176775/1‑176 (MQ=255)
                                                                                                                                                                                          gATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTTCCATCGCCGGATAAACGACCATCCATAAAc                                                                                                                                                                                                  <  1:151450/93‑1 (MQ=255)
                                                                                                                                                                                          gATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTAAAAAGGCTTTCCgt                                                                                <  1:557091/207‑1 (MQ=255)
                                                                                                                                                                                              tATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTAAAAAGGCTTTCCGTTGTTCGCGAATCCAGTTTAGCAATGCGGTTTTTTCATCTTCATCAggg                                >  1:125571/1‑251 (MQ=255)
                                                                                                                                                                                              tATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTAAAAAGGCTTTCCGTTGTTCGCGAATCCAGTTTAGCAATGCGGTTTTTTCATCTTCATCAggg                                >  1:191268/1‑251 (MQ=255)
                                                                                                                                                                                              tATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTAAAAAGGCTTTCCGTTGTTCGCGAATCCAGTTTAGCAATGCGGTTTTTTCATCTTCATCAggg                                <  1:781303/251‑1 (MQ=255)
                                                                                                                                                                                                tCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACggtgg                                                                                                                          >  1:572134/1‑159 (MQ=255)
                                                                                                                                                                                                   aaaCGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTAAAAAGGCTTTCCGTTGTTCGCGAATCCAGTTTAGCAATGCGGTTTTTTCATCTTCATCAGggcgg                             >  1:207514/1‑249 (MQ=255)
                                                                                                                                                                                                   aaaCGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTAAAAAGGCTTTCCGTTGTTCGCGAATCCAGTTTAGCAATGCGGTTTTTTCATCTTCATCAGGGCGGCt                           <  1:665597/251‑1 (MQ=255)
                                                                                                                                                                                                       ggTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTAAAAAGGCTTTCCGTTGTTCGCGAATCCAGTTTAGCAATGCGGTTTTttcatcttcatc                                    >  1:575853/1‑238 (MQ=255)
                                                                                                                                                                                                          cAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTAAAAAGGCTTTCCGTTGTTCGCGAATCCAGTTTAGCAATGCGGTTTTTTCATCTTCATCAGGGCGGCTTCCCCAg                    <  1:486202/251‑1 (MQ=255)
                                                                                                                                                                                                          cAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTAAAAAGGCTTTCCGTTGTTCGCGAATCCAGTTTAGCAATGCGGTTTTTTCATCTTCATCAGGGCGGCTTCCCCAg                    <  1:680137/251‑1 (MQ=255)
                                                                                                                                                                                                            actccGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTc                                                                                                                     <  1:765827/149‑1 (MQ=255)
                                                                                                                                                                                                            acacCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTAAAAAGGCTTTCCGTTGTTCGCGAATCCAGTTTAGCAATGCGGTTTTTTCATCTTCATCAGGGCGGCTTCCCCAGTg                  <  1:497818/251‑1 (MQ=255)
                                                                                                                                                                                                                      tGTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTAAAAAGGCTTTCCGTTGTTCGCGAATCCAGTTTAGCAATGCGGTTTTTTCATCTTCATCAGGGCGGCTTCCCCAGTGAATATgaaga        <  1:382457/251‑1 (MQ=255)
                                                                                                                                                                                                                       gTCCAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTAtttacc                                                                                                  <  1:427996/160‑4 (MQ=255)
                                                                                                                                                                                                                         ccAGTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTAAAAAGGCTTTCCGTTGTTCGCGAATCCAGTTTAGCAATGCGGTTTTTTCATCTTCATCAGGGCGGCTTCCCCAGTGAATATGAAGAAAc     <  1:572538/251‑1 (MQ=255)
                                                                                                                                                                                                                            gTTGGTGGCAGGTGTGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTAAAAAGGCTTTCCGTTGTTCGCGAATCCAGTTTAGCAATGCGGTTTTTTCATCTTCATCAGGGCGGCTTCCCCAGTGAATATGAAGAAACTTg  >  1:69420/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gtgtGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAAc                                                                                                                               >  1:605650/1‑115 (MQ=255)
                                                                                                                                                                                                                                       gtgtGACATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTACCGGTAAAAAGGCTTTCCGTTGTTCGCGAATCCAGTTTAGCAATGCGGTTTTTTCATCTTCATCAggg                                <  1:486243/210‑1 (MQ=255)
                                                                                                                                                                                                                                            aCATTTAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGt                                                                                                                                         >  1:209328/1‑100 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     
TCCTGCTCAAATCGTGCCACAATTTCTTCCCAGGAATGGGAACCCATTTCTACCGGATTTAATGGCGGTCCGCCGGCTTGATCGGCTAAATCCACGGCACGCCCATCCCAGAATTGTGCGAAGTTAAATGCGGCGTTGAAGACAGTCGGCGCGTTAATGCCGCCTTTTTTGCCATCAATGCCGGTGGATGTATCTAAACGGTCAACACCGCCTTTGTCCAGTTGGTGGCAGGTGTGACATTGAATCGTACCATCGCCGGATAAACGACCATCCATAAACAGTTTATGCCCAAGAGCAACTTTTGCCGCATCGGTCATGATTTCATCAGGAATCGGTTGAACTAAACGGTGGTTGTCCGTGCCTTCGGTATTTGCCGGTAAAAAGGCTTTCCGTTGTTCGCGAATCCAGTTTAGCAATGCGGTTTTTTCATCTTCATCAGGGCGGCTTCCCCAGTGAATATGAAGAAACTTG  >  1/533406‑533876

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