Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 519,267 A→G A88A (GCT→GCC)  unknown ← Acyl‑CoA thioesterase II (EC 3.1.2.‑)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1519,2670AG98.6% 177.0 / NA 69A88A (GCT→GCC) unknownAcyl‑CoA thioesterase II (EC 3.1.2.‑)
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (36/32);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACGGGTTATTCAGGTATTTGGTGCGAATATCAAAGGGGCGTTCCGCCGTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATAGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCTGAATGGCGGCAGAAAGGGCTTGCGCCACGACTTGTCCGCCGAAAACCTGACGGAAGCCTAAATCCTGACTTTCGCCGCGAAAGAGTAAATCGTCGATTTTTTCCAGTTTTAACAGGT  >  1/519027‑519505
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 
gAACGGGTTATTCAGGTATTTGGTGCGAATATCAAAGGGGCGTTCCGCCGTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGGAAt                                                                                                                                                                                                                                        >  1:176218/1‑251 (MQ=255)
  aCGGGTTATTCAGGTATTTGGTGCGAATATCAAAGGGGCGTTCCGCCGTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAAttt                                                                                                                                                                                                                                      <  1:853044/251‑1 (MQ=255)
       ttATTCAGGTATTTGGTGCGAATATCAAAGGGGCGTTCCGCCGTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGc                                                                                                                                                                                                                                  >  1:468732/1‑250 (MQ=255)
              ggTATTTGGTGCGAATATCAAAGGGGCGTTCCGCCGTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGGGAAATTTTTGCCTTcgggg                                                                                                                                                                                                                          >  1:710956/1‑248 (MQ=255)
              ggTATTTGGTGCGAATATCAAAGGGGCGTTCCGCCGTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGGGAAATTTTTGCCTTcgcgg                                                                                                                                                                                                                          >  1:384568/1‑250 (MQ=255)
                 aTTTGGTGCGAATATCAAAGGGGCGTTCCGCCGTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGCGGGGGAAATTTTTG‑CTTCCCG‑TAGa                                                                                                                                                                                                                      >  1:18861/1‑251 (MQ=255)
                         cGAATATCAAAGGGGCGTTCCGCCGTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCcac                                                                                                                                                                                                              <  1:217117/251‑1 (MQ=255)
                             tatCAAAGGGGCGTTCCGCCGTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATTTGGTAAATCGGTGATTTGTGCTGAATAGCTTTGGAGCGCATGGC‑GGTGTGATTTTTGCCGTCGCGATGGAATTTcacgtca                                                                                                                                                                                                           >  1:487081/1‑245 (MQ=255)
                             tatCAAAGGGGCGTTCCGCCGTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCAcatca                                                                                                                                                                                                          >  1:14380/1‑251 (MQ=255)
                                 aaaGGGGCGTTCCGCCGTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATaaa                                                                                                                                                                                                      <  1:41355/251‑1 (MQ=255)
                                    ggggCGTTCCGCCGTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAAATTCTGCCTTCGCG‑TAGAATTTTCACATCATaaa                                                                                                                                                                                                      >  1:852758/1‑248 (MQ=255)
                                    ggggCGTTCCGCCGTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGGGGAATTTCTGCCTTCGCG‑GAGAATTTCCACATCATAAATGa                                                                                                                                                                                                   >  1:528234/1‑251 (MQ=255)
                                     gggCGTTCCGCCGTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATg                                                                                                                                                                                                    <  1:791930/249‑1 (MQ=255)
                                      ggCGTTCCGCCGTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCTTTCGCG‑TAGAATTTCCACATCATAAATGA‑t                                                                                                                                                                                                 <  1:161377/250‑1 (MQ=255)
                                      ggCGTTCCGCCGTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑t                                                                                                                                                                                                 <  1:137274/250‑1 (MQ=255)
                                           tCCGCCGTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGata                                                                                                                                                                                           <  1:610651/251‑1 (MQ=255)
                                                cgTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGc                                                                                                                                                                                      <  1:826511/251‑1 (MQ=255)
                                                cgTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGc                                                                                                                                                                                      <  1:665183/251‑1 (MQ=255)
                                                               cACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGCAAAACGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGGGAAAATTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATAA‑ACGGATATTTGCTGTCcccccggcgt                                                                                                                                                                        >  1:126602/1‑244 (MQ=255)
                                                                  cGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTaaaa                                                                                                                                                                    <  1:747642/251‑1 (MQ=255)
                                                                              cATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGGGAAATTTTTGCCTTCGCG‑TAGAATTTTCACATCATAAATGG‑TCGGGTATTGGCTGGCGCCCGGGGGTaaaa                                                                                                                                                                    >  1:353775/1‑239 (MQ=255)
                                                                                    cgcCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGaa                                                                                                                                                   <  1:113600/250‑1 (MQ=255)
                                                                                      ccATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTTCACATCATAAATTA‑TCGGGTATTg                                                                                                                                                                                        >  1:538276/1‑211 (MQ=255)
                                                                                        aTTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTaa                                                                                                                                              <  1:473138/251‑1 (MQ=255)
                                                                                             cTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTTCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAg                                                                                                                                          <  1:427685/250‑1 (MQ=255)
                                                                                               ggATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTTTGCCTTCGCG‑CAGAATTTCCACATCATAAATTA‑TCGGATATTGGGTGTCGCCCGGGGCTAAAAAATACGCCTGGCAAGAATGTAAAAAGCGa                                                                                                                                       >  1:751162/1‑251 (MQ=255)
                                                                                                           tttCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGcc                                                                                                                           <  1:867510/251‑1 (MQ=255)
                                                                                                            ttCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTTCTGAATGGTCTTGAACCGCATAGC‑CGGGGAATTTTTCCCTTCGCG‑GAGCACTTTCACCTCATTACTGATTCGAATTTTGCCT‑TCTCCCGGGGCTAAAAAATACGCCTGGGAAGAAAGTAAAAAGGGCTCTGTCCGCGcc                                                                                                                           >  1:545007/1‑250 (MQ=255)
                                                                                                             tCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATg                                                                                                                                                                                                    >  1:680734/1‑177 (MQ=255)
                                                                                                              ccGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTaaaa                                                                                                                                                                    <  1:462636/207‑1 (MQ=255)
                                                                                                                gACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGGGAAATTTTTGCCCTCGCG‑GAGAAttttc                                                                                                                                                                                                                >  1:436719/1‑160 (MQ=255)
                                                                                                                gACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGGGAAATTTCTGCCTTCGCG‑TAGAATTTTCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGGTAAAAAAAACGCATGGCAAGAATGTAAAAAGCGAACTGTCGGCGCCAc                                                                                                                         >  1:7291/1‑248 (MQ=255)
                                                                                                                 aCATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTTTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTCAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCTGa                                                                                                                     >  1:3143/1‑251 (MQ=255)
                                                                                                                    tGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCcac                                                                                                                                                                                                              >  1:72768/1‑160 (MQ=255)
                                                                                                                     gAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGatat                                                                                                                                                                                          <  1:658045/178‑1 (MQ=255)
                                                                                                                                  gCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCTGAATGGCGGCAGAAAggg                                                                                                     <  1:531386/250‑1 (MQ=255)
                                                                                                                                   ccTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCTGAATGGCGGCAGAAAGGGCt                                                                                                   <  1:7819/251‑1 (MQ=255)
                                                                                                                                    cTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTTTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCTGAATGGCGGGAGAAAGGGCtt                                                                                                  >  1:300021/1‑251 (MQ=255)
                                                                                                                                         cTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTaaa                                                                                                                                             <  1:15865/203‑1 (MQ=255)
                                                                                                                                         cTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGGGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCCGGGCTAAAAAAAAAGCAAGGCaaaaaattaaaaaa                                                                                                                                          >  1:176908/1‑194 (MQ=255)
                                                                                                                                                      aCTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATg                                                                                                                                                                                                    <  1:878909/136‑1 (MQ=255)
                                                                                                                                                       cTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGGGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCT‑AAAAATACGCATGGCAAGAATGTAAAAAACGATCTGTCGGCGCCACCTGAATGGCGGCAGAAAGGGCTTGCGCCACGACTTGTccgccg                                                                              >  1:467037/1‑251 (MQ=255)
                                                                                                                                                         ttGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCTGAAt                                                                                                                   <  1:33538/213‑1 (MQ=255)
                                                                                                                                                              gCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGGAATTTTTGCCTTCGCG‑TAGAATTTTCACAACATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCTGAATGGGGGCAGAAAGGGGTTGCGCCACGAATTTTCCGCCGAAAAc                                                                         >  1:768762/1‑250 (MQ=255)
                                                                                                                                                              gCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCCGAAAGGCGGCAGAAAGGGCTTGCGCCACGACTTGTCCGCCGAAAAc                                                                         >  1:9154/1‑250 (MQ=255)
                                                                                                                                                              gCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGGAAATCGGGGATTTGTGCTGAATGGCCTTGACGCGCATGGC‑GGTGGAACTTCTTCCTTCCCG‑TAAAATTTCCACCTCAAAAATGA‑CCGGGCATTGGGTCGCGCCCCGGGGTAAAAAATATGCATGGCGAGAAAGTAAAAAACGGTCTGTCGGCGCCCCCCTAAAGGGGGGAAAAAGGGCTTGTGCCAcccgttgtcccgcggaaac                                                                         >  1:816576/1‑232 (MQ=255)
                                                                                                                                                               cTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATg                                                                                                                                                                                                    <  1:276438/127‑1 (MQ=255)
                                                                                                                                                               cTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGc                                                                                                                                                             <  1:477646/165‑1 (MQ=255)
                                                                                                                                                                                 ggTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTTTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGGCGCCCGGGGCTAAAAAATACGCAAGGCAAGAATGTAAAAAGCGAACTGGCGGCGCCACCTGAAAGGCGGCAGAAAGGGCTTGCGCCCCGACTTGTCCGCCGAAAACCTGACGGAAGCCTAAATcc                                                      >  1:271104/1‑250 (MQ=255)
                                                                                                                                                                                 ggTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCcac                                                                                                                                                                                                              <  1:62859/99‑1 (MQ=255)
                                                                                                                                                                                  gTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGatat                                                                                                                                                                                          >  1:220340/1‑117 (MQ=255)
                                                                                                                                                                                  gTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCTGAATGGCGGCAGAAAGGGCTTGCGCCACGACTTGTc                                                                                   <  1:735460/220‑1 (MQ=255)
                                                                                                                                                                                     tGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCTGAATGGCGGCAGAAAGGGCTTGCGCCACGACTTGTCCGCCGAAAACCTGACGGAAGCCTAAATCCTGACt                                                 <  1:211085/251‑1 (MQ=255)
                                                                                                                                                                                       gAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGc                                                                                                                                                             <  1:872600/141‑1 (MQ=255)
                                                                                                                                                                                            agGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGaa                                                                                                                                                   <  1:766815/146‑1 (MQ=255)
                                                                                                                                                                                              gCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGaa                                                                                                                                                   >  1:141681/1‑144 (MQ=255)
                                                                                                                                                                                                   gATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATg                                                                                                                                                                                                    >  1:536521/1‑91 (MQ=255)
                                                                                                                                                                                                   gATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCTGAATGGCGGCAGAAAGGGCTTGCGCCACGACTTGTCCGCCGAAAACCTGACGGAAGCCTAAATCCTGACTTTCGCCGCGAAaga                                   >  1:879962/1‑251 (MQ=255)
                                                                                                                                                                                                      gtgGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGaa                                                                                                                                                   <  1:4868/136‑1 (MQ=255)
                                                                                                                                                                                                      gtgGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCTGAATGGCGGCAGAAAGGGCTTGCGCCACGACTTGTCCGCCGAAAACCTGACGGAAGCCTAAATCCTGACTTTCGCCGCGAAAGAGTa                                >  1:467503/1‑251 (MQ=255)
                                                                                                                                                                                                           aaaTCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCTGAATGGCGGCAGAAAGGGCTTGCGCCACGACTTGTCCGCCGAAAACCTGACGGAAGCCTAAATCCTGACTTTCGCCGCGAAAGAGTAAAtcg                           <  1:779698/251‑1 (MQ=255)
                                                                                                                                                                                                           aaaTCGGTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCTGAATGGCGGCAGAAAGGGCTTGCGCCACGACTTGTCCGCCGAAAACCTGACGGAAGCCTAAATCCTGACTTTCGCCGCGAAAGAGTAAAtcg                           <  1:553815/251‑1 (MQ=255)
                                                                                                                                                                                                                ggTGATTTGTGCTGAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCTGAATGGCGGCAGAAAGGGCTTGCGCCACGACTTGTCCGCCGAAAACCTGACGGAAGCCTAAATCCTGACTTTCGCCGCGAAAGAGTAAATCGTCGAt                      >  1:837177/1‑251 (MQ=255)
                                                                                                                                                                                                                             gAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAg                                                                                                                                          <  1:432374/122‑1 (MQ=255)
                                                                                                                                                                                                                             gAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCTGAATGGCGGCAGAAAGGGCTTGCGCCACGACTTGTCCGCCGAAAACCTGACGGAAGCCTAAATCCTGACTTTCGCCGCGAAAGAGTAAATCGTCGATTTTTTCCAGtttt         >  1:662702/1‑251 (MQ=255)
                                                                                                                                                                                                                             gAATGGCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCTGAATGGCGGCAGAAAGGGCTTGCGCCACGACTTGTCCGCCGAAAACCTGACGGAAGCCTAAATCCTGACTTTCGCCGCGAAAGAGTAAATCGTCGATTTTTTCCAGtttt         >  1:865451/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gCTTTGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCTGAATGGCGGCAGAAAGGGCTTGCGCCACGACTTGTCCGCCGAAAACCTGACGGAAGCCTAAATCCTGACTTTCGCCGCGAAAGAGTAAATCGTCGATTTTTTCCAGTTTTAACAg    >  1:91279/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ttGACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCTGAATGGCGGCAGAAAGGGCTTGCGCCACGACTTGTCCGCCGAAAACCTGACGGAAGCCTAAATCCTGACTTTCGCCGCGAAAGAGTAAATCGTCGATTTTTTCCAGTTTTAACAGGt  >  1:571909/1‑250 (MQ=255)
                                                                                                                                                                                                                                       gACGCGCATGGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGc                                                                                                                                                                                      >  1:530413/1‑68 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 
GAACGGGTTATTCAGGTATTTGGTGCGAATATCAAAGGGGCGTTCCGCCGTAAATTTATCGCGCACCGCTTCCGGAATCATTTGCGCCATTTGCTGGATCACCGAGCTTTCCGACATGAGTGCATCGGGGGCCTCCACTTTTGGCATCTCACTTTGATGCTCAAAGCCTTGCTCCAAGGTTTGGAAAGAGGCGGTGATGTGGTAAATCGGTGATTTGTGCTGAATGGCTTTGACGCGCATAGC‑GGTGAAATTTCTGCCTTCGCG‑TAGAATTTCCACATCATAAATGA‑TCGGATATTGGCTGTCGCCCGGGGCTAAAAAATACGCATGGCAAGAATGTAAAAAGCGATCTGTCGGCGCCACCTGAATGGCGGCAGAAAGGGCTTGCGCCACGACTTGTCCGCCGAAAACCTGACGGAAGCCTAAATCCTGACTTTCGCCGCGAAAGAGTAAATCGTCGATTTTTTCCAGTTTTAACAGGT  >  1/519027‑519505

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