Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1630,4740TG49.2% 56.4 / 24.9 66L54V (TTG→GTG) unknownN‑acetylmannosamine‑6‑phosphate 2‑epimerase (EC 5.1.3.9)
Reads supporting (aligned to +/- strand):  ref base T (5/28);  new base G (32/0);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 1.39e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAATGTCTCCATTGACTAGAATGACTGCCAAATGCAGCCGCAAAAAAGAGGAGTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTACGTATCACCCCGTTTTTACA‑GGACATTGAAGACTTGGCGCACGCCGGAGCGGACATTATCGCCGTGGATGGCACTCGTCGCCCGCGCCCGGTAGAGATTGAAAGTGCGGTTAAAAAAATTCATG  >  1/630257‑630696
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                    
gAATGTCTCCACTGACTAGAATGACTGCCAAATGCAGCCGCAAAAAAGAGGAGTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑ATTGGAGGGGTGGGC‑AATTTaaaaagcg                                                                                                                                                                                                   >  1:265338/1‑247 (MQ=255)
   tGTCTCCACTGACTAGAATGACTGCCAAATGCAGCCGCAAAAAAGAGGAGTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAc                                                                                                                                                                                                                 <  1:134311/235‑1 (MQ=255)
   tGTCTCCACTGACTAGAATGACTGCCAAATGCAGCCGCAAAAAAGAGGAGTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGc                                                                                                                                                                                                     <  1:868139/246‑1 (MQ=255)
   tGTCTCCACTGACTAGAATGACTGCCAAATGCAGCCGCAAAAAAGAGGAGTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑CTTGGAGGGGGGGGC‑AATTTTAAAGACGTgg                                                                                                                                                                                                >  1:386460/1‑250 (MQ=255)
        ccactGACTAGAATGACTGCCAAATGCAGCCGCAAAAAAGAGGAGTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGGGGGGGTGCGC‑CTTGGAGGGGGGGGC‑AATTTAAAAACCGTG‑C‑Gccc                                                                                                                                                                                          >  1:408515/5‑250 (MQ=255)
             gACTAGAATGACTGCCAAATGCAGCCGCAAAAAAGAGGAGTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGG‑CTTGAAGGGGGGGGC‑AATTTTAAAAACGTG‑C‑GCCCGACCgg                                                                                                                                                                                    >  1:786831/1‑250 (MQ=255)
             gACTAGAATGACTGCCAAATGCAGCCGCAAAAAAGAGGAGTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑ATTGAAGGTGTGGGC‑AATTTTAAAGCCGTG‑C‑GCCCCACCg                                                                                                                                                                                     >  1:146642/1‑250 (MQ=255)
                   aaTGACTGCCAAATGCAGCCGCAAAAAAGAGGAGTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAAtg                                                                                                                                                                               <  1:25475/250‑1 (MQ=255)
                   aaTGACTGCCAAATGCAGCCGCAAAAAAGAGGAGTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAAt                                                                                                                                                                                <  1:23156/249‑1 (MQ=255)
                   aaTGACTGCCAAATGCAGCCGCAAAAAAGAGGAGTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTGCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAAt                                                                                                                                                                                <  1:428334/249‑1 (MQ=255)
                        cTGCCAAATGCAGCCGCAAAAAAGAGGAGTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑ATTGAAGGGGGGGGC‑AAATTAAAAGCCGGG‑C‑GGCCCACCGTAAATGTGccca                                                                                                                                                                         >  1:141963/1‑249 (MQ=255)
                          gCCAAATGCAGCCGCAAAAAAGAGGAGTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑ATTTAAgg                                                                                                                                                                                                                        >  1:409482/1‑205 (MQ=255)
                          gCCAAATGCAGCCGCAAAAAAGAGGAGTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑ATTGAAGGGGGGGGC‑AATTTAAAAACCGGG‑C‑GCCCGACCGTAAATGTGCCGat                                                                                                                                                                        >  1:146362/1‑250 (MQ=255)
                                  cAGCCGCAAAAAAGAGGAGTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTa                                                                                                                                                                <  1:247757/250‑1 (MQ=255)
                                        cAAAAAAGAGGAGTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑ATTGAAGGGGTGGAC‑AATTTAAAAACCGTG‑C‑GCCCGACCGCAAATGTGCCGATCATCGGGATCGTCaa                                                                                                                                                         >  1:8468/1‑251 (MQ=255)
                                                 ggagTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑ATTGAAGGGGGGGGC‑AATTTAAAAGCCGTG‑C‑GCCCGCCCGTAAATGGGGCGATCATCGGGATCGTCAAACggggt                                                                                                                                                  >  1:237467/1‑247 (MQ=255)
                                                  gagTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑ATTGAAGGTGTGGGC‑CATTTAAAAAACGTG‑C‑GCCCCACCGTAAATGTGCCGAGCATCGGGATCGTCAAACGGGATC‑tt                                                                                                                                              >  1:632455/1‑251 (MQ=255)
                                                    gTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGGGGCGGGGGTGCGC‑ATTGAAGGTGTGGGC‑AATTTAAAAACCGCG‑C‑GCCCGCCCGTAAAAGTGCCGATCATCGGGATCGTCAAACAGGATC‑TTc                                                                                                                                             >  1:790854/1‑250 (MQ=255)
                                                     tttCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑ATTGAAGGGGGGGAC‑AATTTTAAAACCGGG‑C‑GCCCGACCGGAAATGTGCCGatc                                                                                                                                                                       >  1:163451/1‑224 (MQ=255)
                                                        cTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGAc                                                                                                                                        <  1:383009/251‑1 (MQ=255)
                                                           tGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGAc                                                                                                                                        <  1:141668/248‑1 (MQ=255)
                                                                caGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGG‑CTTGAAGGGGGGGGC‑AAATTAAAAGCCGTG‑C‑GCCCGACCGTAAAt                                                                                                                                                                                >  1:318948/1‑204 (MQ=255)
                                                                caGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑ATTGAAGGGGTGGAC‑AATTTAAAAACCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAAAGGGGTC‑TTCCCCACAGccccgg                                                                                                                                >  1:695615/1‑248 (MQ=255)
                                                                   cTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTAc                                                                                                                              <  1:413553/250‑1 (MQ=255)
                                                                   cTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTCAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTAc                                                                                                                              <  1:271564/250‑1 (MQ=255)
                                                                       ctcAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGGGGGGGGGGTGGGG‑ATTGGAGGGGGGGGCAAATTTAAAAG‑CGGG‑G‑GCCCGACCGGAAATGTGGCGGTCATCGGGATCGTCAAAAGGGGAC‑TTTCCGACAGCCCGGTACggat                                                                                                                          >  1:897889/1‑247 (MQ=255)
                                                                       ctcAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTc                                                                                                                                                           <  1:683747/218‑1 (MQ=255)
                                                                       ctcAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTACGTATc                                                                                                                         <  1:474853/251‑1 (MQ=255)
                                                                       ctcAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTACGTATc                                                                                                                         <  1:18091/251‑1 (MQ=255)
                                                                       ctcAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑AATGTAAGTGGGGAC‑AAATTAAAAACCATG‑C‑GCCCGACCGCGAATGTTCCGATCATCGGTAATGTCAAACGGGATA‑TTCCCGACAGACCCGTACGTAtt                                                                                                                         >  1:35786/1‑250 (MQ=255)
                                                                                    tcatcaAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTACGTATCACCCCGTTTTTAc                                                                                                            <  1:239835/251‑1 (MQ=255)
                                                                                      atcaAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑ACTGAAGGTGTGGAC‑AATTTAAAAACCGCG‑C‑GCCCGCCCGTAAATGTGCCGATCATCCGTATCGTCAAACGGGATC‑TTCCCCACAGCCCCGGACGTATCACCCCGCTTTTACA‑g                                                                                                         >  1:91346/1‑251 (MQ=255)
                                                                                         aaGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTACGTATCACCCCGTTTTTACA‑GGAc                                                                                                      <  1:658673/251‑1 (MQ=255)
                                                                                              gTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑ATTGAAGGGGTGGAC‑AATTTAAAAACCGTG‑C‑GCCCGACCGGAAATGTGCCGATCATCGGGATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTACGTATc                                                                                                                         >  1:21594/1‑228 (MQ=255)
                                                                                               ttttCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTACGTATCACCCCGTTTTTACA‑GGACATTGa                                                                                                 <  1:451254/250‑1 (MQ=255)
                                                                                                   cGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTACGTATCACCCCGTTTTTACA‑GGACATTGAAGACt                                                                                            <  1:539883/251‑1 (MQ=255)
                                                                                                         aCTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGGGGGGGTGCGG‑ATTGAAggggggggc                                                                                                                                                                                                                 >  1:113638/1‑126 (MQ=255)
                                                                                                          cTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGCAGGGGGGGCGGGG‑TGCGCAATTTAAGGGGGGGGA‑AATTAAAAAGGC‑TG‑CGGCCCGACGGTAAAATTGCCG‑TCATCGGTATCCTCCAAAGGGGTT‑TTTCCGCCAGCCCCGGATGTATCAACCCGCTTTTTCATGAACATTG‑AGACTTGgcgg                                                                                      >  1:113418/1‑249 (MQ=255)
                                                                                                          cTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑CTTGAAGGGGTGGAC‑AATTTTAAAACCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGGATCGTCAAACGGGATC‑TTCCCGAc                                                                                                                                        >  1:690844/1‑201 (MQ=255)
                                                                                                               gAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGGGGGGGTGCGC‑ATTTGAGGGGGGGAC‑AATTTAAAAACCGTG‑C‑GCCCGACCGGAAATGTGGCGATCCTCCGTATCGGCAAACGGGATCTTTCCCG‑CAGCCCGGTGCGTATCCCCCCGTTTTTTCA‑GGACAATTAAGACTTGGcgcccgcc                                                                                 >  1:257526/1‑250 (MQ=255)
                                                                                                                      ttAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑CTTTAAGGGGTGGAC‑AAATTAAAAACCGTG‑C‑GCCCGGCCGTAAATGTGCCGATCATCGGGATCGTCAAACGGGGAC‑TTCCCCACAGCCCCGTAc                                                                                                                              >  1:492870/1‑199 (MQ=255)
                                                                                                                             ccTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGGGGCGGGGGTGGGC‑ATTGTAGGGGTGGGC‑AAATTAAAAACCGTG‑G‑GCCCGGCCGTAAATGTGCCGAACATCGGGGTTGTCAAAAGGGGTC‑TTTCCGACAGCCCGGTACGCATCACCCCGCTTTTTCA‑GGGAATTGAAGAATTTGGGCACGcccgggcgggcccgt                                                                    >  1:150409/1‑236 (MQ=255)
                                                                                                                             ccTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGGGGGGGTGCGC‑CTTGAAGGGGTGGAC‑AATTTAAAAGACGTG‑C‑GCCCGACCGTAAATGTGGCGATCATCGGTATCGTCAAAAGGGGTC‑TTCCCGACAGCCCCGGACGTATCACCCCGTTTTTTCA‑GGACATTGAAGACTTGGCGCACGCCGGAGCGGGCATTAt                                                                   >  1:548378/1‑250 (MQ=255)
                                                                                                                             ccTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTACGTATCACCCCGTTTTTACA‑GGACATTGAAGACTTGGCGCACGCCGGAGCGGACATTATc                                                                  <  1:576867/251‑1 (MQ=255)
                                                                                                                               tctTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTACGTATCACCCCGTTTTTACA‑GGACATTGAAGACTTGGCGCACGCCGGAGCGGACATTATCGc                                                                <  1:109771/251‑1 (MQ=255)
                                                                                                                               tctTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTACGTATCACCCCGTTTTTACA‑GGACATTGAAGACTTGGCGCACGCCGGAGCGGACATTATCGc                                                                <  1:508371/251‑1 (MQ=255)
                                                                                                                                ctTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGGGGGGTTGCGC‑ATTGGAGGGGGGGGC‑AATTTAAAAACCGTG‑C‑GCCCGACCGTCAATGTGCCGATCATCGGTATCGTCAAACGGGGTC‑TTCCCCACAGCCCGGTACGTATCACCCCGTTTTTTTC‑GGGCATTGAAGACTTGGCGCACGCCGGAGCGGGCATTATCGc                                                                >  1:152315/1‑250 (MQ=255)
                                                                                                                                    tCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTACGTATCACCCCGTTTTTACA‑GGACATTGAAGACTTGGCGCACGCCGGAGCGGACATTATCGCCGTg                                                            <  1:688812/250‑1 (MQ=255)
                                                                                                                                       gCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑ATTGAAGGGGTGGAC‑AATTTAAAAACCGGG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGGATCGTCAAACGGGGTC‑TTCCCGACAGCCCgggac                                                                                                                              >  1:638053/1‑179 (MQ=255)
                                                                                                                                           gTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGGGGCGGGGGTGCGC‑ATTGAAGGGGTGGGC‑AAATTAAAAGCCGGG‑C‑GCCCGACCGGAAATGTGGCGATCATTGGTATCGGCAAAAGGGGTA‑TTCCCCACAGCCCGGGAAGGATCAcccccgttt                                                                                                               >  1:675361/1‑188 (MQ=255)
                                                                                                                                           gTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑CTTGAAGGGGGGGAA‑AATTTAAAAACCGTG‑C‑GCCCCCCCGTAAATGTGCCGATCATCGGGATCGTCAAACGGGATC‑TTCCCGACAGCCCCGTACGTATCACCCCCCTTTTTCA‑GGACATTTAAGAc                                                                                             >  1:277529/1‑210 (MQ=255)
                                                                                                                                              gatgatGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑ATTGAAGGGGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGGACGGATCAc                                                                                                                       >  1:331252/1‑182 (MQ=255)
                                                                                                                                               atgatgGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑CTTGAAGGTGTGGGC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAAAGTGCCGATCATCGGGATCGTCAAAAGGGGTC‑TTTCCCACAGCCCGGGACGCATCACCCCCTTTTTTCA‑GGACATTGAAGACTTTGCGCCCGCCGGGGCCGACAATATCGCCGTGGGTGGCACtcgt                                                >  1:584297/1‑251 (MQ=255)
                                                                                                                                                  atgGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAc                                                                                                                                                                                                                 <  1:493989/92‑1 (MQ=255)
                                                                                                                                                    ggTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑CTTGAAGGGGTGGGC‑AATTTTAAAACCGCGAC‑CCCCCCCCG‑‑AATGTGCCGATCATCGGTTTCGGCCAAAGGGGTC‑TTCCCGCCAGCCCCGGACGGATCACCCCGGTTTTTCA‑GGACATTGAAGAATTGGCGCACGCCGGAGCGGGCAATATCGCCCTGGGTGGCGCTCGTcgccc                                           >  1:220468/1‑250 (MQ=255)
                                                                                                                                                               tAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTACGTATCACCCCGTTTTTACA‑GGACATTGAAGACTTGGCGCACGCCGGAGCGGACATTATCGCCGTGGATGGCACTCGTCGCCCGCGCCCGGTag                                <  1:168125/251‑1 (MQ=255)
                                                                                                                                                                  aCCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTACGTATCACCCCGTTTTTACA‑GGACATTGAAGACTTGGCGCACGCCGGAGCGGACATTATCGCCGTGGATGGCACTCGTCGCCCGCGCCCGGTAgaga                             <  1:650004/251‑1 (MQ=255)
                                                                                                                                                                   ccGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGc                                                                                                                                                                                                     <  1:109970/86‑1 (MQ=255)
                                                                                                                                                                   ccGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑CTTGTAGGGGTGGGC‑AATTTTAAAACCGTG‑C‑GCCTGACCGGAAAAGTGCCGATCATCGGGATCGTCAAACGGGATC‑TTCCCGACAGCCCGGGACGGATCACCCCGTTTTTACC‑GGACATTGAAGAc                                                                                             >  1:359044/1‑186 (MQ=255)
                                                                                                                                                                       aaaTTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTACGTATCACCCCGTTTTTACA‑GGACATTGAAGACTTGGCGCACGCCGGAGCGGACATTATCGCCGTGGATGGCACTCGTCGCCCGCGCCCGGTAGAGATTGaa                        <  1:632137/251‑1 (MQ=255)
                                                                                                                                                                           tGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTACGTATCACCCCGTTTTTACA‑GGACATTGAAGACTTGGCGCACGCCGGAGCGGACATTATCGCCGTGGATGGCACTCGTCGCCCGCGCCCGGTAGAGATTGAAAGTg                    <  1:214982/251‑1 (MQ=255)
                                                                                                                                                                                ctgctATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTACGTAt                                                                                                                          <  1:826715/145‑1 (MQ=255)
                                                                                                                                                                                ctgctATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGGGGGC‑AAATTAAAAACCCTG‑C‑GCCCCACCGGAAATGTGCCGATCATCGGTATCGTCAAACGGGGTC‑TTCCCCACCGCCCCGGACGGATCACCCCCGTTTTACA‑GGACATTGAAGACTTGGCGCACGCCCGAGCGGGCATTATCGCCCTGGGTGGCCCTCCTCGCCCCCCCCCCGGAGAGATTggaaaggcgggt               >  1:724703/1‑240 (MQ=255)
                                                                                                                                                                                             ggCTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑ATTGAAGGGGTGGAC‑AATTTAAAAACCGGG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGGTC‑TTTCCGACAGCCCGGGACGTATCACCCCGTTTTTTCA‑GGACATTGAAGACTTGGCGCACGCCGGAGCGGACATTATCGCCGTGGATGGCACTCGTCGCCCGCGCCCGGGAGAGATTGAAAGTGCGGGTAAAAAAATTCATg  >  1:467550/1‑251 (MQ=255)
                                                                                                                                                                                               cTTCCGTTATCGGC‑GGGGCGGCGGGGGTGCGC‑ATTGGAGGTGTGGGCAAATTTAAAA‑CCGCG‑C‑GCCCCACCGTAAAAGTGCCGATCATCGGTATTGTCAAACGGGATT‑TTCCCCACAGCCCGGTACGTATCACCCCGTTTTTACA‑GGACATTGAAGACTTGGCGCACGCCCGGGCGGACATTATCGCCGTGGATGGGCCTCGTCTCCCCCCCCCGGTAgagagtggaaa                      >  1:606526/1‑222 (MQ=255)
                                                                                                                                                                                                               gggcggcggggTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTACGTATCACCCCGTTTTTACA‑GGACATTGAAGACTTGGCGCACGCCGGAGCGGACATTATCGCCGTGGATGGCACTCGTCGCCCGCGCCCGGTAGAGATTGAAAGTGCGGTTAAAAAAATTCATg  >  1:849542/1‑234 (MQ=255)
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                    
GAATGTCTCCATTGACTAGAATGACTGCCAAATGCAGCCGCAAAAAAGAGGAGTTTCTATGTCACAGCTATCTCAAAAATTATCTCATCAAGCCGTTTTCGCCCAACTTCAGAACGGATTAATTGCCTCTTGTCAGCCGGTGGATGATGGTCCGATGGATAAACCGGAAATTGTCGCTGCTATGGCGCAGGCTTCCGTTATCGGC‑GGGGCGGCGGGGTTGCGC‑ATTGAAGGTGTGGAC‑AATTTAAAAGCCGTG‑C‑GCCCGACCGTAAATGTGCCGATCATCGGTATCGTCAAACGGGATC‑TTCCCGACAGCCCGGTACGTATCACCCCGTTTTTACA‑GGACATTGAAGACTTGGCGCACGCCGGAGCGGACATTATCGCCGTGGATGGCACTCGTCGCCCGCGCCCGGTAGAGATTGAAAGTGCGGTTAAAAAAATTCATG  >  1/630257‑630696

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