Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,894,195 A→C R243S (AGA→AGC)  unknown → Hydrogenase‑4 component I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,894,1950AC98.5% 192.7 / NA 67R243S (AGA→AGC) unknownHydrogenase‑4 component I
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (36/29);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAAAAACTGAAAGGCAAACAGGAAAAAGCCGATCCGAATGCCGACGCCAAATTGCGTTTTCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGAGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCAGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCTGTG‑CCCG‑CCC‑GAATTTAC‑G‑TTT‑GGATAATGGGCGGGCACA‑GGACCCGCCCCTACAAAAGCAACTCTGTTGCTTAGATGTAAAGCAGTAAGGACGTATCATGACCACCGTAATCTTCTATCTTCTCAACGAACGCTTTGTCGAAAACGAAGAACAAGTGCAGGA  >  1/1893954‑1894442
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                
aTCAAAAACTGAAAGGCAAACAGGAAAAAGCCGATCCGAATGCCGACGCCAAATTGCGTTTTCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGt                                                                                                                                                                                                                                                         >  1:249783/1‑250 (MQ=255)
               gCAAACAGGAAAAAGCCGATCCGAATGCCGACGCCAAATTGCGTTTTCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGc                                                                                                                                                                                                                                          <  1:112517/250‑1 (MQ=255)
               gCAAACAGGAAAAAGCCGATCCGAATGCCGACGCCAAATTGCGTTTTCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGc                                                                                                                                                                                                                                          <  1:513489/250‑1 (MQ=255)
                      ggAAAAAGCCGATCCGAATGCCGACGCCAAATTGCGTTTTCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTa                                                                                                                                                                                                                                  >  1:736398/1‑251 (MQ=255)
                       gAAAAAGCCGATCCGAATGCCGACGCCAAATTGCGTTTTCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAg                                                                                                                                                                                                                                    <  1:867442/248‑1 (MQ=255)
                       gAAAAAGCCGATCCGAATGCCGACGCCAAATTGCGTTTTCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTa                                                                                                                                                                                                                                  <  1:36713/250‑1 (MQ=255)
                       gAAAAAGCCGATCCGAATGCCGACGCCAAATTGCGTTTTCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTa                                                                                                                                                                                                                                  <  1:148951/250‑1 (MQ=255)
                           aaGCCGATCCGAATGCCGACGCCAAATTGCGTTTTCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGccg                                                                                                                                                                                                                             >  1:697356/1‑249 (MQ=255)
                                  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‑CGTCAAc                                                                                                                                                                                                                                                     <  1:141316/183‑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)
                                                                                     tATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGccg                                                                                                                                                                                                                             >  1:495921/1‑191 (MQ=255)
                                                                                          gCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGaaa                                                                                                                                                                                                                                                              <  1:854750/156‑1 (MQ=255)
                                                                                          gCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGc                                                                                                                                                                                                                               >  1:268532/1‑186 (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)
                                                                                                 aCGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGt                                                                                                                                                                                                                                   >  1:104634/1‑175 (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)
                                                                                                                           aTCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGCGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGc                                                                                                                                                                                                                               >  1:844445/1‑153 (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‑CGTCAACCGTCTGCATGCCATTAGTa                                                                                                                                                                                                                                  <  1:779778/122‑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‑CGTCAACCGTCTGCATGCCATTAGTa                                                                                                                                                                                                                                  <  1:508865/77‑1 (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)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                
ATCAAAAACTGAAAGGCAAACAGGAAAAAGCCGATCCGAATGCCGACGCCAAATTGCGTTTTCCGGGCATTCCGCTGGATTTGCGTATCAGCTTAGAACGGGAGGCTCGTCGTCTTGCCGGTTATCGCCAGGGTGGCAATATCGCCAATCAATTTATGGCGATGATGTCGGAAGATGACAGCGTGCCGTTC‑GGTGTGCGTTTAGGGGAATACCTGGAACGGGAAGCGGATCCGCGCCTGAGAGAAAT‑CGTCAACCGTCTGCATGCCATTAGTATGCAGTTTTCGGGGTAAGTTATTGG‑CTTACATTCTCCCGTAGGGGCGGGTCCTGTG‑CCCG‑CCC‑GAATTTAC‑G‑TTT‑GGATAATGGGCGGGCACA‑GGACCCGCCCCTACAAAAGCAACTCTGTTGCTTAGATGTAAAGCAGTAAGGACGTATCATGACCACCGTAATCTTCTATCTTCTCAACGAACGCTTTGTCGAAAACGAAGAACAAGTGCAGGA  >  1/1893954‑1894442

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: ‑