Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,551,157 A→G Y219Y (TAT→TAC)  unknown ← Phosphoenolpyruvate carboxykinase [ATP] (EC 4.1.1.49)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,551,1570AG98.5% 177.9 / NA 65Y219Y (TAT→TAC) unknownPhosphoenolpyruvate carboxykinase [ATP] (EC 4.1.1.49)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (37/27);  minor base T (0/1);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TGCGCGATAAATATCCGGTTCGTTTTCTTCGGAAAGGTGAATGGTTTTCGCGTAACAACCACCTTCAAAGTTGAAGATACCCACATCATCCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAATAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAACCGTTCATCACCACGAAACTTGGTTTGAAATCTTTTAATTCGGCGTCGGTCGGACGAATAAACATATT  >  1/1550911‑1551372
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                       
tGCGCGATAAATATCCGGTTCGTTTACTTCGGAAAGGTGAATGGTTTTCGCGTAACAACCACCTTCAAAGTTGAAGATACCCACATCATCCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAagtagt                                                                                                                                                                                                                     <  1:774769/251‑1 (MQ=255)
  cgcgATAAATATCCGGTTCGTTTTCTTCGGAAAGGTGAATGGTTTTCGCGTAACAACCACCTTCAAAGTTGAAGATACCCACATCATCCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTtc                                                                                                                                                                                                                   <  1:704471/251‑1 (MQ=255)
      aTAAATATCCGGTTCGTTTTCTTCGGAAAGGTGAATGGTTTTCGCGTAACAACCACCTTCAAAGTTGAAGATACCCACATCATCCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTtcatca                                                                                                                                                                                                               >  1:415783/1‑251 (MQ=255)
      aTAAATATCCGGTTCGTTTTCTTCGGAAAGGTGAATGGTTTTCGCGTAACAACCACCTTAAAACTTGACGATACCCACACCACCCCAACCGCGTTCATCACCACCAATTAGCTCACGTTCCGGATCGGAAGAAAGGGGGGTTTTACCGGTGCCGGACAAGCCGAAGAACCCAGCAACATCACCGCCTTTGCCTACGTTAGCGGAACAGTGCATTGGGCCCACGCCTTTTAGCGGCAGGAAGAAGTtcatca                                                                                                                                                                                                               >  1:57718/1‑251 (MQ=255)
              ccGGTTCGTTTTCTTCGGAAAGGTGAATGGTTTTCGCGTAACAACCACCTTCAAAGTTGAAGATACCCACATCATCCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAAc                                                                                                                                                                                                       <  1:249554/251‑1 (MQ=255)
                 gTTCGTTTTCTTCGGAAAGGTGAATGGTTTTCGCGTAACAACCACCTTCAAAGTTGAAGATACCCACATCATCCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATa                                                                                                                                                                                                    >  1:576573/1‑251 (MQ=255)
                  ttCGTTTTCTTCGGAAAGGTGAATGGTTTTCGCGTAACAACCACCTTCAAAGTTGAAGATACCCACATCATCCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGAC‑AGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATAcc                                                                                                                                                                                                  <  1:223408/251‑1 (MQ=255)
                   tCGTTTTCTTCGGAAAGGTGAATGGTTTTCGCGTAACAACCACCTTCAAAGTTGAAGATACCCACATCATCCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATAcc                                                                                                                                                                                                  <  1:408857/251‑1 (MQ=255)
                                   ggTGAATGGTTTTCGCGTAACAACCACCTTCAAAGTTGAAGATACCCACATCATCCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTtc                                                                                                                                                                                                                   <  1:754952/218‑1 (MQ=255)
                                   ggTGAATGGTTTTCGCGTAACAACCACCTTCAAAGTTGAAGATACCCACATCATCCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATAc                                                                                                                                                                                                   >  1:833271/1‑234 (MQ=255)
                                   ggTGAATGGTTTTCGCGTAACAACCACCTTCAAAGTTGAAGATACCCACATCATCCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTcgccg                                                                                                                                                                                  >  1:171792/1‑251 (MQ=255)
                                   ggTGAATGGTTTTCGCGTAACAACCACCTTCAAAGTTGAAGATACCCACATCATCCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGGAGTTCATCATAGAGAACATACCTTTCTTCATTTcgccg                                                                                                                                                                                  >  1:458870/1‑251 (MQ=255)
                                              ttCGCGTATCAACCACCTTCAAAGTTGAAGATACCCACATCATCCAAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCTCCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCGGGAATTAGTTCAGCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGt                                                                                                                                                                       <  1:896368/251‑1 (MQ=255)
                                                             ccTTCAAAGTTGAAGATACCCACATCATCCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTg                                                                                                                                                        >  1:487770/1‑251 (MQ=255)
                                                              cTTCAAAGTTGAAGATACCCACATCATCCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATAcc                                                                                                                                                                                                  >  1:533524/1‑208 (MQ=255)
                                                                   aaGTTGAAGATACCCACATCATCCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCg                                                                                                                                                  >  1:824829/1‑251 (MQ=255)
                                                                             tACCCACATCATCCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTtc                                                                                                                                                                                                                   <  1:246333/176‑1 (MQ=255)
                                                                              aCCCACATCATCCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTttcttc                                                                                                                                                                                           <  1:265735/199‑1 (MQ=255)
                                                                                      catcCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACg                                                                                                                                 >  1:647686/1‑249 (MQ=255)
                                                                                        tcCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGc                                                                                                                                <  1:71243/248‑1 (MQ=255)
                                                                                           cAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGaa                                                                                                                          >  1:811887/1‑251 (MQ=255)
                                                                                             aCCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGt                                                                                                                        <  1:696811/251‑1 (MQ=255)
                                                                                             aCCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGt                                                                                                                        <  1:459289/251‑1 (MQ=255)
                                                                                                gtgtTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGtttt                                                                                                                     >  1:284778/1‑251 (MQ=255)
                                                                                                gtgtTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGtttt                                                                                                                     >  1:858751/1‑251 (MQ=255)
                                                                                                 tgtTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTttcttc                                                                                                                                                                                           >  1:88607/1‑180 (MQ=255)
                                                                                                     catcatcaCCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGt                                                                                                                                                                       <  1:601405/196‑1 (MQ=255)
                                                                                                      atcatcaCCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAAt                                                                                                               >  1:720823/1‑251 (MQ=255)
                                                                                                             ccAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGc                                                                                                         >  1:537253/1‑250 (MQ=255)
                                                                                                                  tAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGgttaagt                                                                                                                                       >  1:577818/1‑215 (MQ=255)
                                                                                                                    gCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGtt                                                                                                   >  1:207072/1‑249 (MQ=255)
                                                                                                                     cTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTc                                                                                                  <  1:401046/249‑1 (MQ=255)
                                                                                                                                gATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAAc                                                                                                                                                                                                       >  1:650296/1‑137 (MQ=255)
                                                                                                                                     gTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGTGAAc                                                                                                                                                                                                       >  1:67459/1‑132 (MQ=255)
                                                                                                                                     gTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGt                                                                                                                                                                       <  1:322145/164‑1 (MQ=255)
                                                                                                                                     gTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCAAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGgtta                                                                                                                                          <  1:876605/193‑1 (MQ=255)
                                                                                                                                         aaaGGGTGGTTTTCCCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGc                                                                              <  1:289885/249‑1 (MQ=255)
                                                                                                                                         aaaGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGc                                                                              <  1:730225/249‑1 (MQ=255)
                                                                                                                                               tggtTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGaa                                                                      >  1:789386/1‑251 (MQ=255)
                                                                                                                                                 gtTTTACCGGTGCGGGACAAGCCGAAGAGCACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAAc                                                                     <  1:578155/250‑1 (MQ=255)
                                                                                                                                                 gtTTTACCGGTGCCGGACAAGCCGAAGAACGCAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAAc                                                                     <  1:133623/250‑1 (MQ=255)
                                                                                                                                                 gtTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAAc                                                                     <  1:92791/250‑1 (MQ=255)
                                                                                                                                                 gtTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAAc                                                                     <  1:153819/250‑1 (MQ=255)
                                                                                                                                                     tcccGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAACCGTtc                                                                <  1:379927/249‑1 (MQ=255)
                                                                                                                                                        cGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAACCGTtcatc                                                             <  1:68677/251‑1 (MQ=255)
                                                                                                                                                          gTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGc                                                                                                                                <  1:20511/182‑1 (MQ=255)
                                                                                                                                                                    aGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGACGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAACCGTTCATCACCACGAAACtt                                                 >  1:120367/1‑251 (MQ=255)
                                                                                                                                                                     gCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAACCGTTCATCACCACGAAACTTg                                                >  1:521597/1‑251 (MQ=255)
                                                                                                                                                                     gCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGCTCGGGTTAGTGCATTTAGAACCGTTCATCACCACGAAACTTg                                                >  1:811573/1‑251 (MQ=255)
                                                                                                                                                                                        tCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATAcc                                                                                                                                                                                                  >  1:397024/1‑86 (MQ=255)
                                                                                                                                                                                             gTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGt                                                                                                                                                                       >  1:760508/1‑108 (MQ=255)
                                                                                                                                                                                                   gCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAAc                                                                                                                                                                                                       <  1:42366/70‑1 (MQ=255)
                                                                                                                                                                                                        cGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATagag                                                                                                                                                                                                          <  1:876197/62‑1 (MQ=255)
                                                                                                                                                                                                           tAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTc                                                                                                  >  1:223825/1‑163 (MQ=255)
                                                                                                                                                                                                           tAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAACCGTTCATCACCACGAAACTTGGTTTGAAATCTTTTAATTCGGCGTCGGTCGGACGAATa          >  1:590405/1‑251 (MQ=255)
                                                                                                                                                                                                            aGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAACCGTTCATCACCACGAAACTTGGTTTGAAATCTTTTAATTCGGCGTCGGTCGGACGAATaa         >  1:200234/1‑251 (MQ=255)
                                                                                                                                                                                                                gAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAACCGTTCATCACCACGAAACTTGGTTTGAAATCTTTTAATTCGGCGTCGGTCGGACGAATAAACat     >  1:421553/1‑251 (MQ=255)
                                                                                                                                                                                                                 aaCAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAACCGTTCATCACCACGAAACTTGGTTTGAAATCTTTTAATTCGGCGTCGGTCGGACGAATAAACat     <  1:689222/250‑1 (MQ=255)
                                                                                                                                                                                                                  aCAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAACCGTTCATCACCACGAAACTTGGTTTGAAATCTTTTAATTCGGCGTCGGTCGGACGAATAAACatat   >  1:343065/1‑251 (MQ=255)
                                                                                                                                                                                                                  aCAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAACCGTTCATCACCACGAAACTTGGTTTGAAATCTTTTAATTCGGCGTCGGTCGGACGAATAAACatat   >  1:584850/1‑251 (MQ=255)
                                                                                                                                                                                                                  aCAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCCGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAACCGTTCATCACCACGAAACTTGGTTTGAAATCTTTTAATTCGGCGTCGGTCGGACGAATAAACatat   >  1:586756/1‑251 (MQ=255)
                                                                                                                                                                                                                   cAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAACCGTTCATCACCACGAAACTTGGTTTGAAATCTTTTAATTCGGCGTCGGTCGGACGAATAAACATAtt  >  1:423665/1‑251 (MQ=255)
                                                                                                                                                                                                                     gTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCtt                                                                                                                                                                                                >  1:276422/1‑59 (MQ=255)
                                                                                                                                                                                                                     gTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAACCGTTCATcac                                                           <  1:818029/192‑1 (MQ=255)
                                                                                                                                                                                                                       gCATTGCGCCCACGCCTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAACCGTTCATCACCACGAAACTTGGttt                                            >  1:891174/1‑205 (MQ=255)
                                                                                                                                                                                                                                     ccTTTTAGCGGCAGGAAGTAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAACCGTTCATcac                                                           >  1:221818/1‑176 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                       
TGCGCGATAAATATCCGGTTCGTTTTCTTCGGAAAGGTGAATGGTTTTCGCGTAACAACCACCTTCAAAGTTGAAGATACCCACATCATCCCAACCGTGTTCATCATCACCAATTAGCTCACGTTTCGGATCGGTAGAAAGGGTGGTTTTACCGGTGCCGGACAAGCCGAAGAACACAGCAACATCACCGTCTTTGCCTACGTTAGCGGAACAGTGCATTGCGCCCACGCCTTTTAGCGGCAGGAAATAGTTCATCATAGAGAACATACCTTTCTTCATTTCGCCGCCGTACCAGGTGCCGCCGATTAATTGAATGCGTTCGGTTAAGTTAAACGCTACGAAGTTTTCGGAATTCAAGCCTTGTTCTTTCCAGTTCGGGTTAGTGCATTTAGAACCGTTCATCACCACGAAACTTGGTTTGAAATCTTTTAATTCGGCGTCGGTCGGACGAATAAACATATT  >  1/1550911‑1551372

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