Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,647,331 G→A D224D (GAC→GAT)  unknown ← 2‑dehydro‑3‑deoxygluconokinase (EC 2.7.1.45)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,647,3310GA100.0% 163.7 / NA 58D224D (GAC→GAT) unknown2‑dehydro‑3‑deoxygluconokinase (EC 2.7.1.45)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/31);  total (27/31)

CGACATAATTGCCTTGTCGGCAGCAAATGTCTAAAGGTTTATTTCGTAAATAACCATTCAAAAAGCCGGCATTAAAAGAATCTCCTGCTGATGTGGTGTCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTGTCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTTATAAGCTGTTCAAGCAAAAAAGTGCGGTCATTTTTAGGTAAGATCGCAAGTGAAATGCCACTCAAGTAAATCATATCCACTTGT  >  1/1647084‑1647574
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
cgACATAATTGCCTTGTCGGCAGCAAATGTCTAAAGGTTTATTTCGTAAATAACCATTCAAAAAGCCGGCATTAAAAGAACCTCCTGCTGATGTGGTGTCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCg                                                                                                                                                                                                                                                  <  1:244728/251‑1 (MQ=255)
                                         tctCGTAAATAACCATTCCAAAAGCCGGCATTAACAGAATCTCCTGCTGATGTGGTGTCAACAACATTTGAAGTACGTTTCGGCACCACTCCTCCATATTTATAATCTTGAGAATCCGCACAAATCGCACCATTTCCACACTACTTTACAACTACTTTTGGTACGCCACGTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCTGTCCACaa                                                                                                                                                                                                                                          >  1:706871/3‑218 (MQ=255)
                                             gTAAATAACCATTCAAAAAGCCGGCATTAAAAGAATCTCCTGCTGATGTGGTGTCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACt                                                                                                                                                                                                      >  1:315541/1‑250 (MQ=255)
                                                aaTAACCATTCAAAAAGCCGGCATTAAAAGAATCTCCTGCTGATGTGGTGTCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTcaca                                                                                                                                                                                                  <  1:585408/251‑1 (MQ=255)
                                                     ccATTCAAAAAGCCGGCATTAAAAGAATCTCCTGCTGATGTGGTGTCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCtc                                                                                                                                                                                                                                 <  1:137611/215‑1 (MQ=255)
                                                     ccATTCAAAAAGCCGGCATTAAAAGAATCTCCTGCTGATGTGGTGTCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAgg                                                                                                                                                                                             >  1:498140/1‑251 (MQ=255)
                                                           aaaaaGCCGGCATTAAAAGAATCTCCTGCTGATGTGGTGTCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGtaataa                                                                                                                                                                                       <  1:803655/251‑1 (MQ=255)
                                                            aaaaGCCGGCATTAAAAGAATCTCCTGCTGATGTGGTGTCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTaataat                                                                                                                                                                                      <  1:323165/251‑1 (MQ=255)
                                                                 ccGGCATTAAAAGAATCTCCTGCTGATGTGGTGTCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCaa                                                                                                                                                                                 <  1:9449/251‑1 (MQ=255)
                                                                      aTTAAAAGAATCTCCTGCTGATGTGGTGTCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAAc                                                                                                                                                                            <  1:723323/251‑1 (MQ=255)
                                                                               aTCTCCTGCTGATGTGGTGTCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAg                                                                                                                                                                   >  1:610911/1‑251 (MQ=255)
                                                                                tctcCTGCTGATGTGGTGTCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGc                                                                                                                                                                  <  1:461611/251‑1 (MQ=255)
                                                                                 ctcCTGCTGATGTGGTGTCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGc                                                                                                                                                                  >  1:133348/1‑250 (MQ=255)
                                                                                      gctgATGTGGTGTCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTg                                                                                                                                                            >  1:162880/1‑251 (MQ=255)
                                                                                         gATGTGGTGTCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGtaataa                                                                                                                                                                                       >  1:691973/1‑221 (MQ=255)
                                                                                         gATGTGGTGTCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCaa                                                                                                                                                                                 <  1:553388/227‑1 (MQ=255)
                                                                                          aTGTGGTGTCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGgaa                                                                                                                                                         <  1:510405/250‑1 (MQ=255)
                                                                                          aTGTGGTGTCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGgaa                                                                                                                                                         <  1:47942/250‑1 (MQ=255)
                                                                                               gtgtCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGgaatgaat                                                                                                                                                    <  1:406175/250‑1 (MQ=255)
                                                                                                 gtCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTg                                                                                                                                                                                    >  1:650921/1‑216 (MQ=255)
                                                                                                 gtCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGcc                                                                                                                                                                 <  1:210069/235‑1 (MQ=255)
                                                                                                 gtCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGcc                                                                                                                                                                 <  1:424104/235‑1 (MQ=255)
                                                                                                         caTTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCt                                                                                                                                                                       <  1:100053/221‑1 (MQ=255)
                                                                                                              gAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCagag                                                                                                                                             >  1:457791/1‑242 (MQ=255)
                                                                                                              gAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTAACGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGt                                                                                                                                      >  1:589105/1‑249 (MQ=255)
                                                                                                                  tAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTaaaa                                                                                                                                <  1:402283/251‑1 (MQ=255)
                                                                                                                                 cTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCtc                                                                                                                                                                                                                                 <  1:608200/139‑1 (MQ=255)
                                                                                                                                 cTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCt                                                                                                                                                                       <  1:870567/197‑1 (MQ=255)
                                                                                                                                   tGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTaataat                                                                                                                                                                                      <  1:411045/180‑1 (MQ=255)
                                                                                                                                     tCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGAtt                                                                                                              <  1:599609/250‑1 (MQ=255)
                                                                                                                                        atatTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCa                                                                                                          >  1:833582/1‑251 (MQ=255)
                                                                                                                                                  gTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTg                                                                                                                                                                     >  1:397715/1‑182 (MQ=255)
                                                                                                                                                              gAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTtat                                                                                    >  1:483564/1‑251 (MQ=255)
                                                                                                                                                                          ccATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCtttttt                                                                                                >  1:111263/1‑227 (MQ=255)
                                                                                                                                                                            aTTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTTATAAGCTGTTCAAGCa                                                                      <  1:96313/251‑1 (MQ=255)
                                                                                                                                                                                cGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGAtt                                                                                                              <  1:49324/207‑1 (MQ=255)
                                                                                                                                                                                     gCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTTATAAGCTGTTCAAGCAAAAAAGTg                                                              >  1:654554/1‑250 (MQ=255)
                                                                                                                                                                                        cTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTTATAAGCTGTTCAAGCAAAAAAGTGCGGt                                                          >  1:659058/1‑251 (MQ=255)
                                                                                                                                                                                          ttACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTTATAAGCTGTTCAAGCAAAAAAGTGCGGTCa                                                        >  1:470005/1‑251 (MQ=255)
                                                                                                                                                                                                  actTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGtaataa                                                                                                                                                                                       <  1:365967/116‑1 (MQ=255)
                                                                                                                                                                                                  actTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTTATAAGCTGTTCAAGCAAAAAAGTGCGGTCATTTTTAgg                                                >  1:798953/1‑251 (MQ=255)
                                                                                                                                                                                                                aaTTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTTATAAGCTGTTCAAGCAAAAAAGTGCGGTCATTTTTAGGTAAGATCGCAAGTg                                  >  1:531660/1‑251 (MQ=255)
                                                                                                                                                                                                                aaTTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTTATAAGCTGTTCAAGCAAAAAAGTGCGGTCATTTTTAGGTAAGATCGCAAGTg                                  >  1:506657/1‑251 (MQ=255)
                                                                                                                                                                                                                      aTGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAAc                                                                                                                                                                            <  1:691395/107‑1 (MQ=255)
                                                                                                                                                                                                                      aTGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCATGCCTTTTTTAGCTAATTTTATAAGCTGTTCAAGCAAAAAAGTGCGGTCATTTTTAGGTAAGATCGCAAGTGAAATg                             >  1:492221/1‑250 (MQ=255)
                                                                                                                                                                                                                      aTGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTTATAAGCTGTTCAAGCAAAAAAGTGCGGTCATTTTTAGGTAAGATCGCAAGTGAAATGc                            <  1:491569/251‑1 (MQ=255)
                                                                                                                                                                                                                               cgTTAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTTATAAGCTGTTCAAGCAAAAAAGTGCGGTCATTTTTAGGTAAGATCGCAAGTGAAATGCCACTCAAGt                   >  1:453639/1‑251 (MQ=255)
                                                                                                                                                                                                                                 ttAGCGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTTATAAGCTGTTCAAGCAAAAAAGTGCGGTCATTTTTAGGTAAGATCGCAAGTGAAATGCCACTCAAGTaa                 <  1:243154/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTTATAAGCTGTTCAAGCAAAAAAGTGCGGTCATTTTTAGGTAAGATCGCAAGTGAAATGCCACTCAAGTAAAt               <  1:258905/249‑1 (MQ=255)
                                                                                                                                                                                                                                     cGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTTATAAGCTGTTCAAGCAAAAAAGTGCGGTCATTTTTAGGTAAGATCGCAAGTGAAATGCCACTCAAGTAAATc              <  1:56921/250‑1 (MQ=255)
                                                                                                                                                                                                                                     cGTGGCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTTATAAGCTGTTCAAGCAAAAAAGTGCGGTCATTTTTAGGTAAGATCGCAAGTGAAATGCCACTCAAGTAAATCa             <  1:636911/251‑1 (MQ=255)
                                                                                                                                                                                                                                        ggCCGGTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTTATAAGCTGTTCAAGCAAAAAAGTGCGGTCATTTTTAGGTAAGATCGCAAGTGAAATGCCACTCAAGTAAATCatat          >  1:102979/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAAt                                                                                                                               <  1:84092/129‑1 (MQ=255)
                                                                                                                                                                                                                                             gTTCATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTTATAAGCTGTTCAAGCAAAAAAGTGCGGTCATTTTTAGGTAAGATc                                         >  1:594244/1‑215 (MQ=255)
                                                                                                                                                                                                                                                cATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTTATGAGCTGTTCAAGCAAAAAAGTGCGGTCATTTTTAGGTAAGATCGCAAGTGAAATGCCACTCAAGTAAATCATATCCACTTGt  >  1:439165/1‑251 (MQ=255)
                                                                                                                                                                                                                                                cATTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTTATAAGCTGTTCAAGCAAAAAAGTGCGGTCATTTTTAGGTAAGATCGCAAGTGAAATGCCACTCAAGTAAATCATATCCACTTGt  >  1:477669/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 aTTTGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAAc                                                                                                                                                                                                                  <  1:200444/42‑1 (MQ=255)
                                                                                                                                                                                                                                                  tttGTATCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCAc                                                                                                         >  1:434906/1‑146 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
CGACATAATTGCCTTGTCGGCAGCAAATGTCTAAAGGTTTATTTCGTAAATAACCATTCAAAAAGCCGGCATTAAAAGAATCTCCTGCTGATGTGGTGTCAACAACATTTGAAATAGGTTTCGGCACAACTTGTCCATATTGAGAAGTTTGAGAATCTGAAAAAATCGCACCATTTCGACCGCACTTTACAACTACTTTTGGTATGCCAATTTGATGTAAGCGCGTTAGCGTGGCCGGTTCATTTGTGTCGTTCCACAATAATTGCTCATCATCAAAGGTAACAAGCGCGAGACTCACACAAGGTAATAATTGCAAATAACAATCTTGAGCCTCTTGGAATGAATCCCAGAGTTTAGGTCTAAAATTACTATCAAAAGAGATTTGCACGCCTTTTTTAGCTAATTTTATAAGCTGTTCAAGCAAAAAAGTGCGGTCATTTTTAGGTAAGATCGCAAGTGAAATGCCACTCAAGTAAATCATATCCACTTGT  >  1/1647084‑1647574

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