Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,293,915 G→A R167R (CGC→CGT)  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
*11,293,9150GA98.5% 189.7 / NA 65R167R (CGC→CGT) 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 G (0/1);  new base A (38/26);  total (38/27)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGGGTGGAATTGCACGCCGTAGAAACGGCGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGTGATCAGAGGTTTCCGTTAAGCCGCCAAGCTGCATTGCCATGGTTTGCATACCGT  >  1/1293680‑1294159
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    
cGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGAc                                                                                                                                                                                                                                       <  1:206468/251‑1 (MQ=255)
          aTGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATggc                                                                                                                                                                                                                                  >  1:623826/1‑246 (MQ=255)
          aTGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCaa                                                                                                                                                                                                                             >  1:347046/1‑251 (MQ=255)
          aTGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCaa                                                                                                                                                                                                                             >  1:435009/1‑251 (MQ=255)
          aTGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCaa                                                                                                                                                                                                                             >  1:593464/1‑251 (MQ=255)
          aTGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCaa                                                                                                                                                                                                                             >  1:265469/1‑251 (MQ=255)
            gTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAAt                                                                                                                                                                                                                            >  1:621315/1‑250 (MQ=255)
            gTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAAt                                                                                                                                                                                                                            >  1:41830/1‑250 (MQ=255)
             taataaCAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCg                                                                                                                                                                                                                          >  1:682178/1‑251 (MQ=255)
                                gaagaaTCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGc                                                                                                                                                                                                        <  1:268698/250‑1 (MQ=255)
                                     aTCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTgt                                                                                                                                                                                                  <  1:308434/251‑2 (MQ=255)
                                                       gCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCa                                                                                                                                                                                >  1:83888/1‑251 (MQ=255)
                                                       gCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCa                                                                                                                                                                                >  1:182173/1‑251 (MQ=255)
                                                        ccTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGAc                                                                                                                                                                                                                                       <  1:740402/195‑1 (MQ=255)
                                                                               tACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGCGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTgg                                                                                                                                                         >  1:880851/1‑250 (MQ=255)
                                                                                         ttAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCcacac                                                                                                                                              <  1:688875/251‑1 (MQ=255)
                                                                                         ttAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCcaca                                                                                                                                               >  1:814293/1‑250 (MQ=255)
                                                                                               cgcgcCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTaa                                                                                                                                        <  1:576735/251‑1 (MQ=255)
                                                                                                            cTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCg                                                                                                                           >  1:706301/1‑251 (MQ=255)
                                                                                                            cTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCg                                                                                                                           >  1:685194/1‑251 (MQ=255)
                                                                                                            cTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCg                                                                                                                           >  1:173880/1‑251 (MQ=255)
                                                                                                              tCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCTGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGt                                                                                                                          >  1:836052/1‑250 (MQ=255)
                                                                                                                   gatgTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGGAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCAtttt                                                                                                                    >  1:7607/1‑250 (MQ=255)
                                                                                                                      gTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGa                                                                                                                                                                                                  >  1:613204/1‑170 (MQ=255)
                                                                                                                      gTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCcacac                                                                                                                                              <  1:847943/222‑1 (MQ=255)
                                                                                                                       ttttCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAAt                                                                                                                 <  1:688759/250‑1 (MQ=255)
                                                                                                                                gTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGc                                                                                                                                                                                                                       <  1:334672/139‑1 (MQ=255)
                                                                                                                                           ttCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGAt                                                                                            <  1:99510/251‑1 (MQ=255)
                                                                                                                                                               ctacGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGGAACTTTATCGCCGTGGCTCATCCACACATCTAAtt                                                                                                                                      >  1:228081/1‑189 (MQ=255)
                                                                                                                                                               ctacGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGGCGGCGTGGTGCCGGCGATTTGGAAGTCCTTCGGCAGGCGGGTTACCTTATTTCCGGGGGCCCTCCACACATCTAATTTTCTGTTGCCCTCATTTAAAAGGGCAAAAAGTGCCGGCGGATTTTCTATTAAAACGGATGCa                                                                        >  1:561947/1‑251 (MQ=255)
                                                                                                                                                                         ggTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGc                                                                                                                                                        <  1:453522/161‑1 (MQ=255)
                                                                                                                                                                         ggTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCCAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGaa                                                                >  1:266514/1‑249 (MQ=255)
                                                                                                                                                                         ggTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGaa                                                                >  1:291204/1‑249 (MQ=255)
                                                                                                                                                                            taataaTTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTcg                                                            >  1:586957/1‑250 (MQ=255)
                                                                                                                                                                                  ttCCAAGCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCgg                                                         >  1:216018/1‑247 (MQ=255)
                                                                                                                                                                                        gCCTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGTGATCAGAgg                                               >  1:371714/1‑251 (MQ=255)
                                                                                                                                                                                         ccTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCCCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAAATCGCGGTGATCAGAggg                                              >  1:556578/1‑250 (MQ=255)
                                                                                                                                                                                          cTTTTTTGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGTGATCAGAGGtt                                             <  1:626197/251‑1 (MQ=255)
                                                                                                                                                                                              tttGGTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGTGATCAGAGGTTTCCg                                         <  1:723104/251‑1 (MQ=255)
                                                                                                                                                                                                 ggTGTGGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGTGATCAGATGTTTCCGTTa                                      >  1:836612/1‑251 (MQ=255)
                                                                                                                                                                                                    gtgGGTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGc                                                                                                                                                                                                        <  1:596461/86‑1 (MQ=255)
                                                                                                                                                                                                      gggTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGc                                                                                                                                                                                 <  1:351100/107‑1 (MQ=255)
                                                                                                                                                                                                      gggTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGTGATCAGAgg                                               <  1:11308/237‑1 (MQ=255)
                                                                                                                                                                                                       ggTTGCTTCAGTGTGGAATTGCACGCCGTAGAAACGGCGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCGTCCACACATCTAATTTGCTGTCGCCGTCATTTAAAt                                                                                                                 <  1:364332/170‑1 (MQ=255)
                                                                                                                                                                                                       ggTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGTGATCAGAGGTTTCCGTTAAgccg                                 >  1:43397/1‑250 (MQ=255)
                                                                                                                                                                                                       ggTTACTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGGGATCAGAGGGTTCCGTTAAgccg                                 >  1:183241/1‑250 (MQ=255)
                                                                                                                                                                                                            cTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAAcc                                                                   >  1:812825/1‑211 (MQ=255)
                                                                                                                                                                                                            cTTCAGTGTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGTGATCAGAgg                                               <  1:776421/231‑1 (MQ=255)
                                                                                                                                                                                                                 gtgTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCg                                                                  <  1:789612/205‑1 (MQ=255)
                                                                                                                                                                                                                   gTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGTGATCAGAgg                                               >  1:208470/1‑224 (MQ=255)
                                                                                                                                                                                                                   gTGGAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGTGATCAGAGGGTTCCGTTAAGCCGCCAAGCTGCATTg                    >  1:49741/1‑251 (MQ=255)
                                                                                                                                                                                                                      gAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATTCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGGGATCAGAGGTTTCCGTTAAGCCGCCAAGCTGCATTGCCa                 >  1:532790/1‑251 (MQ=255)
                                                                                                                                                                                                                      gAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGcc                                                                                                                            >  1:743212/1‑144 (MQ=255)
                                                                                                                                                                                                                      gAATTGCACGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGGAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGTGATCAGAGGTTTCCGTTAAGCCGCCAAGCTGCCTTGCCa                 >  1:694341/1‑251 (MQ=255)
                                                                                                                                                                                                                              cGCCGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGGAACTTTATCGCCGTGGCTCATCcacac                                                                                                                                              >  1:413210/1‑118 (MQ=255)
                                                                                                                                                                                                                                ccGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGc                                                                                                                                                              <  1:611462/100‑1 (MQ=255)
                                                                                                                                                                                                                                ccGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGTGATCAGAGGTTTCCGTTAAGCCGCCAAGCTGCATTGCCATGGtt            <  1:361021/246‑1 (MQ=255)
                                                                                                                                                                                                                                ccGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGTGATCAGAGGTTTCCGTTAAGCCGCCAAGCTGCATTGCCATGGTTTGCAt       <  1:324671/251‑1 (MQ=255)
                                                                                                                                                                                                                                ccGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGTGATCAGAGGTTTCCGTTAAGCCGCCAAGCTGCATTGCCATGGTTTGCAt       <  1:700219/251‑1 (MQ=255)
                                                                                                                                                                                                                                ccGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGTGATCAGAGGTTTCCGTTAAGCCGCCAAGCTGCATTGCCATGGTTTGCAt       <  1:660629/251‑1 (MQ=255)
                                                                                                                                                                                                                                ccGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGTGATCAGAGGTTTCCGTTAAGCCGCCAAGCTGCATTGCCATGGTTTGCAt       <  1:378369/251‑1 (MQ=255)
                                                                                                                                                                                                                                ccGTAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAAAGGAGGCATAACCGAACTCGCGGTGATCAGAGGTTTCCGTTAAGCCGCCAAGCTGCATTGCCATGGTTTGCAt       >  1:610640/1‑251 (MQ=255)
                                                                                                                                                                                                                                   tAGAAACGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGTGATCAGAGGTTTCCGTTAAGCCGCCAAGCTGCATTGCCATGGTTTGCATAc     <  1:361551/250‑1 (MQ=255)
                                                                                                                                                                                                                                      aaaCGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGTGATCAGAGGTTTCCGTTAAGCCGCCAAGCTGCATTGCCATGGTTTGCATACCGt  <  1:216163/250‑1 (MQ=255)
                                                                                                                                                                                                                                      aaaCGACGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGGAACTTTATTGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGGCGGATTTTCTATTAAAACGGAGGCATAAc                                                                    >  1:549417/1‑184 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    
CGATAGCGCGATGTAATAACAGCGCCACCACGGAAGAATCCACGCCGCCGGATAAGCCTAAAATCACTTCGTCATCACCTACTTGCGCTTTAATGCGCGCCACGGCATCTTCAATGATGTTTTCTGCCGTCCATTTGGTTTCGCAACCGCAAATGTTTACTACGAAATTGGTTAATAATTCCAAGCCTTTTTTGGTGTGGGTTACTTCAGGGTGGAATTGCACGCCGTAGAAACGGCGATTTTCATCAGACATGGCGGCAATCGGGCAGGTCGGCGTGGTGCCGGTGATTTGGAAGTCCTTCGGCAGGCGGGTAACTTTATCGCCGTGGCTCATCCACACATCTAATTTGCTGTCGCCGTCATTTAAATGGGCAAAAAGTGCGGTCGGATTTTCTATTAAAACGGAGGCATAACCGAACTCGCGGTGATCAGAGGTTTCCGTTAAGCCGCCAAGCTGCATTGCCATGGTTTGCATACCGT  >  1/1293680‑1294159

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