Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 915,604 T→C H214R (CAC→CGC)  unknown ← GMP synthase [glutamine‑hydrolyzing], amidotransferase subunit (EC 6.3.5.2) / GMP synthase [glutamine‑hydrolyzing], ATP pyrophosphatase subunit (EC 6.3.5.2)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1915,6040TC100.0% 185.7 / NA 67H214R (CAC→CGC) unknownGMP synthase [glutamine‑hydrolyzing], amidotransferase subunit (EC 6.3.5.2) / GMP synthase [glutamine‑hydrolyzing], ATP pyrophosphatase subunit (EC 6.3.5.2)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/31);  total (36/31)

TCGCTCAAGAAACGGCTTTCTGCATCAACACGGGTAATGTTCAAACCGAATTTATCGCCGAACATTTCCATCACCTGAACGCCTTCGTTTAAACGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGTGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGGGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACAC  >  1/915373‑915848
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
tCGCTCAAGAAACGGCTTTCTGCATCAACACGGGTAATGTTCAAACCGAATTTATCGCCGAACATTTCCATCACCTGAACGCCTTCGTTTAAACGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAat                                                                                                                                                                                                                                   <  1:231837/251‑1 (MQ=255)
              gCTTTCTGCATCAACACGGGTAATGTTCAAACCGAATTTATCGCCGAACATTTCCATCACCTGAACGCCTTCGTTTAAACGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCg                                                                                                                                                                                                                     >  1:740079/1‑251 (MQ=255)
                     gCATCAACACGGGTAATGTTCAAACCGAATTTATCGCCGAACATTTCCATCACCTGAACGCCTTCGTTTAAACGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCg                                                                                                                                                                                                                     >  1:469563/1‑244 (MQ=255)
                      cATCAACACGGGTAATGTTCAAACCGAATTTATCGCCGAACATTTCCATCACCTGAACGCCTTCGTTTAAACGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTg                                                                                                                                                                                                             <  1:677315/251‑1 (MQ=255)
                               gggTAATGTTCAAACCGAATTTATCGCCGAACATTTCCATCACCTGAACGCCTTCGTTTAAACGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCaa                                                                                                                                                                                                    >  1:233962/1‑251 (MQ=255)
                                 gTAATGTTCAAACCGAATTTATCGCCGAACATTTCCATCACCTGAACGCCTTCGTTTAAACGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAAc                                                                                                                                                                                                   >  1:859186/1‑250 (MQ=255)
                                    atcttCAAACCGAATTTATCGCCGAACATTTCCATCACCTGAACGCCTTCGTTTAAACGCAATAAACCGTTATCCACGCACACGCACTGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCCCCACGGAAGAATCCACGCCGCCGCCTCACCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTCATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTATCGCAACTGCa                                                                                                                                                                                               >  1:706901/4‑251 (MQ=255)
                                         cAAACCGAATTTATCGCCGAACATTTCCATCACCTGAACGCCTTCGTTTAAACGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGt                                                                                                                                                                                          <  1:646077/251‑1 (MQ=255)
                                                         ccGAACATTTCCATCACCTGAACGCCTTCGTTTAAACGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTgg                                                                                                                                                                           <  1:69306/250‑1 (MQ=255)
                                                             aCATTTCCATCACCTGAACGCCTTCGTTTAAACGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAatga                                                                                                                                                                                                                                 >  1:418134/1‑192 (MQ=255)
                                                               aTTTCCATCACCTGAACGCCTTCGTTTAAACGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTtaataa                                                                                                                                                                    <  1:670918/251‑1 (MQ=255)
                                                                tttCCATCACCTGAACGCCTTCGTTTAAACGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTaataat                                                                                                                                                                   >  1:69434/1‑251 (MQ=255)
                                                                  tCCATCACCTGAACGCCTTCGTTTAAACGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAAtt                                                                                                                                                                  >  1:748550/1‑250 (MQ=255)
                                                                       cACCTGAACGCCTTCGTTTAAACGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGTCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTtacta                                                                                                                                                                                    >  1:290369/1‑227 (MQ=255)
                                                                         ccTGAACGCCTTCGTTTAAACGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGt                                                                                                                                                                                                           <  1:633249/202‑1 (MQ=255)
                                                                         ccTGAACGCCTTCGTTTAAACGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGc                                                                                                                                                            >  1:340114/1‑249 (MQ=255)
                                                                              cattacccgtgttgatgcagaaagccGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCtttttt                                                                                                                                                     >  1:822410/25‑251 (MQ=255)
                                                                                  cTTCGTTTAAACGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGgtg                                                                                                                                                 <  1:291183/251‑1 (MQ=255)
                                                                                             cGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTc                                                                                                                                      <  1:737348/251‑1 (MQ=255)
                                                                                             cGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTc                                                                                                                                      <  1:763043/251‑1 (MQ=255)
                                                                                             cGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTc                                                                                                                                      <  1:812922/251‑1 (MQ=255)
                                                                                               cAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGgtg                                                                                                                                                 <  1:620263/238‑1 (MQ=255)
                                                                                               cAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGgtg                                                                                                                                                 <  1:738736/238‑1 (MQ=255)
                                                                                                   aaaCCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGtgtg                                                                                                                                >  1:864970/1‑251 (MQ=255)
                                                                                                      ccGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGaa                                                                                                                             >  1:278191/1‑251 (MQ=255)
                                                                                                         gtATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTg                                                                                                                          >  1:571874/2‑251 (MQ=255)
                                                                                                           aTCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTc                                                                                                                                      <  1:407820/237‑1 (MQ=255)
                                                                                                            tCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCAc                                                                                                                       <  1:182844/251‑1 (MQ=255)
                                                                                                                     acGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGt                                                                                                                                                                                                                              >  1:899089/1‑139 (MQ=255)
                                                                                                                      cGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTgg                                                                                                                                                   >  1:788649/1‑213 (MQ=255)
                                                                                                                      cGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAAGCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGaaa                                                                                                             <  1:443256/251‑1 (MQ=255)
                                                                                                                        cAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGaaa                                                                                                             <  1:511496/249‑1 (MQ=255)
                                                                                                                          aTGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTgg                                                                                                                               >  1:722906/1‑229 (MQ=255)
                                                                                                                          aTGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGaaa                                                                                                             <  1:41123/247‑1 (MQ=255)
                                                                                                                                      gCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTtactac                                                                                                                                                                                   <  1:405488/165‑1 (MQ=255)
                                                                                                                                        cGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGAc                                                                                           <  1:206468/251‑1 (MQ=255)
                                                                                                                                                  aTGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTg                                                                                                                                                                                                             >  1:291110/1‑127 (MQ=255)
                                                                                                                                                  aTGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTgg                                                                                                                                                   >  1:690732/1‑185 (MQ=255)
                                                                                                                                                  aTGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTc                                                                                                                                      <  1:168877/198‑1 (MQ=255)
                                                                                                                                                  aTGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTc                                                                                                                                      >  1:169399/1‑198 (MQ=255)
                                                                                                                                                  aTGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTc                                                                                                                                      >  1:335249/1‑198 (MQ=255)
                                                                                                                                                  aTGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTg                                                                                                                          <  1:311158/210‑1 (MQ=255)
                                                                                                                                                  aTGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGaaa                                                                                                             >  1:10968/1‑223 (MQ=255)
                                                                                                                                                  aTGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGaaa                                                                                                             <  1:164289/223‑1 (MQ=255)
                                                                                                                                                  aTGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATggc                                                                                      >  1:623826/1‑246 (MQ=255)
                                                                                                                                                  aTGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCaa                                                                                 >  1:435009/1‑251 (MQ=255)
                                                                                                                                                  aTGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCaa                                                                                 >  1:593464/1‑251 (MQ=255)
                                                                                                                                                  aTGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCaa                                                                                 >  1:347046/1‑251 (MQ=255)
                                                                                                                                                  aTGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCaa                                                                                 >  1:265469/1‑251 (MQ=255)
                                                                                                                                                  aTGTAATAACAGCGCCACCACAGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCAc                                                                                                                       <  1:232687/213‑1 (MQ=255)
                                                                                                                                                   tGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGaaa                                                                                                             <  1:280911/222‑1 (MQ=255)
                                                                                                                                                    gTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAAt                                                                                >  1:41830/1‑250 (MQ=255)
                                                                                                                                                    gTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAAt                                                                                >  1:621315/1‑250 (MQ=255)
                                                                                                                                                     taataaCAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCg                                                                              >  1:682178/1‑251 (MQ=255)
                                                                                                                                                           cAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGaaa                                                                                                             >  1:698819/1‑214 (MQ=255)
                                                                                                                                                                        gaagaaTCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGc                                                            <  1:268698/250‑1 (MQ=255)
                                                                                                                                                                             aTCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTgt                                                      <  1:308434/251‑2 (MQ=255)
                                                                                                                                                                                    ccgccgGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTtactac                                                                                                                                                                                   <  1:16392/119‑1 (MQ=255)
                                                                                                                                                                                         ggATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTtaataa                                                                                                                                                                    <  1:146926/129‑1 (MQ=255)
                                                                                                                                                                                               gCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCa                                    >  1:182173/1‑251 (MQ=255)
                                                                                                                                                                                               gCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCa                                    >  1:83888/1‑251 (MQ=255)
                                                                                                                                                                                                ccTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGAc                                                                                           <  1:740402/195‑1 (MQ=255)
                                                                                                                                                                                                          cTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGaa                                                                                                                             >  1:621839/1‑151 (MQ=255)
                                                                                                                                                                                                                 atcaCCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGaaa                                                                                                             <  1:691940/160‑1 (MQ=255)
                                                                                                                                                                                                                       tACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGCGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTgg             >  1:880851/1‑250 (MQ=255)
                                                                                                                                                                                                                                 ttAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCcacac  <  1:688875/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ttAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCcaca   >  1:814293/1‑250 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
TCGCTCAAGAAACGGCTTTCTGCATCAACACGGGTAATGTTCAAACCGAATTTATCGCCGAACATTTCCATCACCTGAACGCCTTCGTTTAAACGCAATAAACCGTTATCCACGAACACGCAATGTAAGTTTTTGCCGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGTGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGGGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACAC  >  1/915373‑915848

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: ‑