Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,894,230 A→C Q255P (CAG→CCG)  unknown → Hydrogenase‑4 component I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,894,2300AC98.2% 150.0 / NA 55Q255P (CAG→CCG) unknownHydrogenase‑4 component I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (31/23);  minor base T (0/1);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGAATGCCGACGCCAAATTGCGTTTTCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGAGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCAGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCTGTG‑CCCG‑CCC‑GAATTTAC‑G‑TTT‑GGATAATGGGCGGGCACA‑GGACCCGCCCCTACAAAAGCAACTCTGTTGCTTAGATGTAAAGCAGTAAGGACGTATCATGACCACCGTAATCTTCTATCTTCTCAACGAACGCTTTGTCGAAAACGAAGAACAAGTGCAGGAGCAGGCGC  >  1/1893988‑1894450
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                    
tCCGAATGCCGACGCCAAATTGCGTTTTCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGCTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGCTTTCgg                                                                                                                                                                                                                              >  1:340072/1‑251 (MQ=255)
 ccGAATGCCGACGCCAAATTGCGTTTTCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCggg                                                                                                                                                                                                                             <  1:525139/251‑1 (MQ=255)
        ccGACGCCAAATTGCGTTTTCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGtt                                                                                                                                                                                                                      <  1:773771/251‑1 (MQ=255)
        ccGACGCCAAATTGCGTTTTCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGtt                                                                                                                                                                                                                      <  1:284370/251‑1 (MQ=255)
              ccAAATTGCGTTTTCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGTTGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑c                                                                                                                                                                                                               <  1:528512/251‑1 (MQ=255)
                    tGCGTTTTCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTTGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACAt                                                                                                                                                                                                         <  1:618368/251‑1 (MQ=255)
                      cGTTTTCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATtc                                                                                                                                                                                                       <  1:632478/251‑1 (MQ=255)
                         tttCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTcc                                                                                                                                                                                                    <  1:49479/251‑1 (MQ=255)
                                     ccGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCggg                                                                                                                                                                                        >  1:494258/1‑251 (MQ=255)
                                          ggATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAGT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTc                                                                                                                                                                                                                                >  1:440126/1‑207 (MQ=255)
                                          ggATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCggggc                                                                                                                                                                                      >  1:529437/1‑246 (MQ=255)
                                            aTTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCtgtg                                                                                                                                                                                 <  1:71214/251‑1 (MQ=255)
                                              ttGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCggggcccgtgg                                                                                                                                                                                >  1:501128/1‑242 (MQ=255)
                                              ttGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCTGTg‑cc                                                                                                                                                                              <  1:593532/251‑1 (MQ=255)
                                                gCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAgggg                                                                                                                                                                                            >  1:718415/1‑236 (MQ=255)
                                                        gCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCgggg                                                                                                                                                                                                                            <  1:603928/197‑1 (MQ=255)
                                                        gCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑c                                                                                                                                                                                                               >  1:532893/1‑209 (MQ=255)
                                                        gCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑c                                                                                                                                                                                                               >  1:197950/1‑209 (MQ=255)
                                                        gCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAgggg                                                                                                                                                                                            <  1:756692/228‑1 (MQ=255)
                                                        gCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCTGTg‑cccc‑cccccagt                                                                                                                                                                   >  1:504547/1‑246 (MQ=255)
                                                        gCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCTGTg‑cc                                                                                                                                                                              <  1:866842/241‑1 (MQ=255)
                                                        gCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCTGTG‑CCCG‑CCC‑GAAtt                                                                                                                                                                  >  1:787513/1‑251 (MQ=255)
                                                        gCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCCTTg‑cccg‑cccccatt                                                                                                                                                                   >  1:147568/1‑246 (MQ=255)
                                                        gCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCACCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTcc                                                                                                                                                                                                    >  1:3485/1‑220 (MQ=255)
                                                                       cTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCTGTg‑cccg‑cc                                                                                                                                                                         >  1:704483/1‑230 (MQ=255)
                                                                         cgtcgtCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCTGTG‑CCCG‑CCC‑GAATTTAC‑G‑TTT‑GGATAATggg                                                                                                                                              <  1:561533/251‑1 (MQ=255)
                                                                           tcgtcTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGGCCTGTG‑GCCG‑CCC‑GGAttttaggtttggtaaaggggg                                                                                                                                               >  1:530841/1‑233 (MQ=255)
                                                                             gtcTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCTGTG‑CCCG‑CCC‑GAATTTAC‑G‑TTT‑GGATAATggg                                                                                                                                              >  1:790578/1‑247 (MQ=255)
                                                                                      gTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAg                                                                                                                                                                                               <  1:709401/195‑1 (MQ=255)
                                                                                      gTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCTGTG‑CCCG‑CCC‑GAATTTAC‑G‑TTT‑GGATAATggg                                                                                                                                              >  1:691621/1‑238 (MQ=255)
                                                                                      gTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGGCCTGTG‑GCCG‑CCC‑GCATTTTC‑G‑TTT‑TGGTAATgggcgggggcaggag                                                                                                                                  >  1:44150/1‑242 (MQ=255)
                                                                                       ttATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTGAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCTGTG‑CCCG‑CCC‑GAATTTAC‑G‑TTT‑GGATAATGGGCGGGCACA‑GGAcc                                                                                                                                <  1:441037/250‑1 (MQ=255)
                                                                                          tCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCTGTG‑CCCG‑CCC‑GAATTTAC‑G‑TTT‑GGATAATGGGCGGGCACA‑GGACCCGcc                                                                                                                            <  1:259864/251‑1 (MQ=255)
                                                                                                          cccgcaaataattttagggggatgatGTGGGAAGATGA‑GGGGGGCCGTTGGGGTGTGCGTTTGGGGAAATCCCGGGAAGGGGAAGGGGATCCGCGCGTGAGAAAAAT‑CGACACCCGTCGGCATGCAATTAGTATGCCGTTTGGGGGGTAAGTTATGGG‑CTTACATTCCCCCGGGGGGGGGGGTCCTGTG‑CCCG‑CCC‑GAATTTCC‑G‑TTG‑GGATAAGGGGGGGGCTCA‑GGACCCGCCCCTACAAAAGCAActc                                                                                                            <  1:100806/231‑1 (MQ=255)
                                                                                                                     atTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGGCCTGTG‑GCCG‑CCCGGAATTTA‑‑G‑TTT‑TGGTAAAGGGGGGGGACA‑GGAcccccccccccaaaaaacaatcttgttgtttgg                                                                                                 >  1:28201/1‑226 (MQ=255)
                                                                                                                     atTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGGCCTGTG‑CCCG‑CCC‑CGATTTTC‑G‑TTT‑TGGTAAATGGCGGGGACA‑GGAcccccccccccaaaaaacaacccgtttctttga                                                                                                 >  1:734981/1‑226 (MQ=255)
                                                                                                                             cgagggtgTCGGAGGATGACGGCGTGCGGTTC‑GGTGTGCGTTGGGGGAAATCCGTGGGACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGACACCCGTCTGCATGCCTTTTGTAGGCCTTTTGCGGGGTAAGTTATTGG‑CTTACTTTCTCCGGTAGGGGGGGTTCCTGTG‑CCCG‑CCC‑GATTTTAC‑G‑TGT‑GGATAAGGGGGGGGCAGA‑GGACCCGCCCCTACAAAAGCAACTCTGTTGCTTAGATGTAAAGc                                                                                         <  1:901129/245‑1 (MQ=255)
                                                                                                                                  atgTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGCCCTGTG‑CCCGCCCC‑GAA‑TTAC‑GCTTT‑GG‑TAAAGGGGGGGCCGA‑GGGCCCCCCCCCCCAAAAACAAATCTGTTGCTTAGATGGaaaaaagtaa                                                                                    >  1:483677/1‑244 (MQ=255)
                                                                                                                                           gATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCTGTG‑CCCG‑CCC‑GAATTTAC‑G‑TTT‑GGAt                                                                                                                                                    <  1:487490/179‑1 (MQ=255)
                                                                                                                                                    gTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGGCCCGTG‑GCCGCCCC‑CAATTT‑C‑G‑TTT‑TGGTAAAAGGGGGGCGCA‑GGACCCCCCCCCCCAAAAAAAACACTGTTTCTTTGATGTAAAAAAGTAAGGACggagcctttctctt                                                                  >  1:861987/1‑238 (MQ=255)
                                                                                                                                                     gggcggttgggtgggctttggggggAATCCCTGGAGGGGGAGGCGGACCCGCGCCGGAG‑CAAATCCGCCACCCGTCTGCTTCCCTTTAGTTTGCTGTTTGGGGGTTAATTTTTTGG‑CTTCCTTTCCCCCGGGGGGGGGGGTCCTGCC‑CCCG‑CCC‑AATTTTCC‑T‑TTT‑GAATAAGGGGGGGGCACA‑GGCCCCGCCCCTACAAAAGCAACTCTGTTGCTTAGATGTAAAGCAGTAAGGACGTATCATGaccac                                                                   <  1:239164/229‑1 (MQ=255)
                                                                                                                                                       ccGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCCGTG‑CCCC‑CCC‑CCATTTTCGG‑TTT‑GGTAAATGGG‑GGGCACA‑AGACCCCCCCCCCCAAAAACAACTCTTTTGCTTTGATGTAAAGCAGGAAGGACGGATCATTaccacccca                                                               >  1:42518/1‑248 (MQ=255)
                                                                                                                                                       ccGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGGCCTGTG‑CCCG‑CCC‑GAATTTAC‑G‑TTT‑GGATAAAGGGCGGGCACA‑GGACCCCCCCCCCCAAAAACAAACCTGTTGCTTAGATGTAAAACAGTAAGGACGGATCATGaccacccta                                                               >  1:513918/1‑248 (MQ=255)
                                                                                                                                                               gtgtGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCTGTG‑CCCG‑CCC‑GAATTTAC‑G‑TTT‑GGATAATGGGCGGGCACA‑g                                                                                                                                    >  1:489523/1‑175 (MQ=255)
                                                                                                                                                               gtgtGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGGCCCGTGCCCCC‑CCC‑GAATTT‑C‑G‑TTT‑GGGTAATTGGCGGGGACC‑GGACCCCCCCCCCCAAAAACAACTCTTATTCTTAGATGTAAAGCAGGAAAGACGGATCATGAACACCCTAatctttt                                                        >  1:756094/1‑249 (MQ=255)
                                                                                                                                                                   gCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTGAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCCGTG‑CCCG‑CCC‑GCATTTTC‑G‑TTT‑TGGTAATGGGGGGGCACA‑AGAACCCCCCCCACAAAAAAAAACCTGTTTCTTTGATGTAAAACAGGAAAGAAGGATCATTAACACCCTAatcttttat                                                      >  1:536486/1‑249 (MQ=255)
                                                                                                                                                                   gCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTcc                                                                                                                                                                                                    >  1:140348/1‑114 (MQ=255)
                                                                                                                                                                   gCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGGCCCGTG‑CCCC‑CCC‑CGATTTTC‑G‑TTT‑TGTTTATGGGCGGGCACA‑GGACCCGCCCCCACCAAAACAAATCTGTTTCTTTGATGTAAAACAGGAAGGACGGATCacctgctctttttcacaactcc                                                     >  1:301716/1‑229 (MQ=255)
                                                                                                                                                                     gTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCTGTG‑CCCG‑CCC‑GAATTTAC‑G‑t                                                                                                                                                           <  1:259133/148‑1 (MQ=255)
                                                                                                                                                                     gTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGGCCTGTG‑CCCG‑CCC‑CGATTTAC‑G‑GTT‑TGATAATGGGGGGGCACA‑GGACCCCCCCCCCCCAAAACAACACTGTTGGTTTGATGTAAAAGAGGAAAGACGGAACATTAACACCGCAaacttctatcttt                                                  >  1:742703/1‑250 (MQ=255)
                                                                                                                                                                                tACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCTGTG‑CCCG‑CCC‑GAATTTAC‑G‑t                                                                                                                                                           <  1:292450/137‑1 (MQ=255)
                                                                                                                                                                                                              ggagaAAT‑CGTCACCCGTCTGCATGCAATTAGTAGGCCTTTTTGGGGGT‑AGTTTTTGGTCTTACTTTCCCCGGGGGGGGGGGGTCGTGCC‑CCCG‑‑CC‑AAATTTTG‑G‑TTTGGGAT‑AGGGGGGGGCAAAGGGCCCCGCCCCTACAAAAGCAACTCTGTTGCTTAGATGTAAAGCAGTAAGGACGTATCATGACCACCGTAATCTTCTATCTTCTCAACGAACGCTTTGTCGAAAACGAAGAACAAGTGCAGGa          <  1:210994/249‑1 (MQ=255)
                                                                                                                                                                                                                    aT‑CGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTcc                                                                                                                                                                                     >  1:226345/1‑80 (MQ=255)
                                                                                                                                                                                                                       cGTCAACCGTCTGCATGCCATTAGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCCGGG‑GCCC‑CCC‑GAAATTTC‑G‑TTT‑TGCTAAAAGGGGGGCGCA‑GGACCCCCCCCCCCAAAACCAACTCCGTTTGTTAGATGTTAAAAAGAAAGAACGTCTCATGACCCCCCCAAAATTTTATCTTCTCCACGAACGCCTTTTCGAAAACgaagaaaaaaggcaggggccagcgc  >  1:444691/1‑231 (MQ=255)
                                                                                                                                                                                                                                             aGTATGCCGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAg                                                                                                                                                                                               <  1:459973/46‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                    
TCCGAATGCCGACGCCAAATTGCGTTTTCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGAGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCAGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCTGTG‑CCCG‑CCC‑GAATTTAC‑G‑TTT‑GGATAATGGGCGGGCACA‑GGACCCGCCCCTACAAAAGCAACTCTGTTGCTTAGATGTAAAGCAGTAAGGACGTATCATGACCACCGTAATCTTCTATCTTCTCAACGAACGCTTTGTCGAAAACGAAGAACAAGTGCAGGAGCAGGCGC  >  1/1893988‑1894450

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