Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 243,633 G→T 1.3% intergenic (+62/+58) mogA → / ← ubiG molybdopterin biosynthesis mog protein/3‑demethylubiquinone‑9 3‑O‑methyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034243,6330G→T 1.3% 2.1 ‑2.0 ‑0.1298intergenic (+62/+58)mogA/ubiGmolybdopterin biosynthesis mog protein/3‑demethylubiquinone‑9 3‑O‑methyltransferase
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (169/125):  total (173/125)
Fisher's exact test strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.96e-01

Only 100 of 325 total aligned reads displayed.

CCGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATCG  >  NC_014034/243533‑243733
                                                                                                              |                                                                                                    
   gCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAg                                                                                                   >  1:3895310/1‑101 (MQ=255)
    cGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGc                                                                                                  <  2:3745695/101‑1 (MQ=255)
     gAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGcc                                                                                                 <  1:3578158/101‑1 (MQ=255)
      aaaaTCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCt                                                                                                >  1:171164/1‑101 (MQ=255)
        aaTCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGt                                                                                              <  1:5053053/101‑1 (MQ=255)
         aTCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGtt                                                                                             <  1:1221650/101‑1 (MQ=255)
           cgcgGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.GCTC.TGG.C..TTC.GCC.A..TTCACACCCCGGATATTTGAGGCAAGAT.AAAGGAGCCTGtttt                                                                                           >  1:4322195/1‑101 (MQ=255)
           cgcgGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGtttt                                                                                           <  1:1057115/101‑1 (MQ=255)
           cgcgGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGtttt                                                                                           <  1:312415/101‑1 (MQ=255)
           cgcgGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGtttt                                                                                           <  1:2025351/101‑1 (MQ=255)
             cggcgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTAt                                                                                         <  2:1132647/101‑1 (MQ=255)
             cggcgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTAt                                                                                         <  2:2374963/101‑1 (MQ=255)
             cggcgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTAt                                                                                         <  1:4935869/101‑1 (MQ=255)
              ggcgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATc                                                                                        <  2:4270120/101‑1 (MQ=255)
              ggcgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATc                                                                                        <  1:2942261/101‑1 (MQ=255)
               gcgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCt                                                                                       >  2:3956992/1‑101 (MQ=255)
               gcgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCt                                                                                       <  1:4144229/101‑1 (MQ=255)
                cgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCtt                                                                                      <  2:3245382/101‑1 (MQ=255)
                cgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCtt                                                                                      <  2:2675130/101‑1 (MQ=255)
                 gTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTg                                                                                     <  2:674911/101‑1 (MQ=255)
                  tttCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGc                                                                                    <  2:3653525/101‑1 (MQ=255)
                    tCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCt                                                                                  <  2:472282/101‑1 (MQ=255)
                     cGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.AGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTc                                                                                 <  1:3511076/101‑1 (MQ=255)
                      gCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCa                                                                                <  1:3295137/101‑1 (MQ=255)
                      gCCCGAAGACGAAGTAAAAGGGGGCGCTCCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCa                                                                                <  1:3847225/101‑1 (MQ=255)
                       cccGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCaa                                                                               <  2:3206193/101‑1 (MQ=255)
                       cccGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCaa                                                                               <  1:2432011/101‑1 (MQ=255)
                        ccGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCaaa                                                                              <  2:1024842/101‑1 (MQ=255)
                        ccGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCaaa                                                                              >  1:4012101/1‑101 (MQ=255)
                         cGAAGACGAAGTAAAAGGGGGCGCTTCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAAt                                                                             <  1:4590877/101‑1 (MQ=255)
                         cGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAAt                                                                             <  2:3426508/101‑1 (MQ=255)
                          gAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATg                                                                            >  2:3575222/1‑101 (MQ=255)
                           aaGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGt                                                                           >  1:2284258/1‑101 (MQ=255)
                             gACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTcc                                                                         <  2:487527/101‑1 (MQ=255)
                               gccccgaaaaGGGGGCGCTGCCCACCTCTTGGCC..TTCGGCC.AATTTAAACCCCGGGATATTTGCGGCAAGGTGAAAGGA.CGGGGTTTATCTTGCCCCAAAt                                                                             >  2:3892540/7‑101 (MQ=9)
                                gAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGCGATATTTGAGGCAAGAT.AAAGGATCCTGTTTTATCTTGCCTCAAATGTCCCgg                                                                      >  2:853031/1‑101 (MQ=255)
                                 aaGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCggg                                                                     >  1:3810552/1‑101 (MQ=255)
                                  aGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCgggg                                                                    <  1:306068/101‑1 (MQ=255)
                                    tAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCGAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGc                                                                  >  2:2892693/1‑101 (MQ=255)
                                    gggcaGGGGGCGCTGCCC.CCTCTTGG.CCTTTG.GCCAA..TTCACCCCTGGGATATTTCGGGCAAGATGAAATCATGCCCCTTCATCTTGCCGCAAATATCCCg                                                                       >  2:2059290/5‑101 (MQ=0)
                                    gggcaGGGGGCGCTGCCC.CCTCTTGG.CCTTTG.GCCAA..TTCACCCCTGGGATATTTCGGGCAAGATGAAATCATGCCCCTTCATCTTGCCGCAAATATCCCg                                                                       >  2:284437/5‑101 (MQ=0)
                                     aaaaGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGcc                                                                 >  2:2735687/1‑101 (MQ=255)
                                     aaaaGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGcc                                                                 <  2:4977770/101‑1 (MQ=255)
                                     aaaaGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGcc                                                                 >  2:1407155/1‑101 (MQ=255)
                                     aaaaGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGcc                                                                 <  2:4377347/101‑1 (MQ=255)
                                     aaaaGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGcc                                                                 <  1:3398034/101‑1 (MQ=255)
                                       aaGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTg                                                               >  1:2470075/1‑101 (MQ=255)
                                         gggggCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTggg                                                             >  1:1192276/1‑101 (MQ=255)
                                          ggggCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTgggg                                                            >  1:1844010/1‑101 (MQ=255)
                                          ggggCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTgggg                                                            >  1:1901732/1‑101 (MQ=255)
                                          ggggCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTgggg                                                            <  1:4739628/101‑1 (MQ=255)
                                          ggggCAATGCCC.CGTC.TGG.C..TTC.GCT.A..TTCTTCCCGGGGGTATTTGAGGCAGGAT.AAAGGAGGCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTgggg                                                            <  1:419654/101‑1 (MQ=21)
                                           gggCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTggggg                                                           <  2:3374492/101‑1 (MQ=255)
                                            ggCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATTTCCCGGGGGCCTGGGGGc                                                          <  2:1761154/101‑1 (MQ=255)
                                               gcTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAgg                                                       >  1:2568390/1‑100 (MQ=255)
                                               gcTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAgc                                                       <  1:436361/101‑1 (MQ=255)
                                                  gCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGcc                                                    >  2:2974169/1‑101 (MQ=255)
                                                   ccc.ccTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGccc                                                   >  2:3842010/1‑101 (MQ=255)
                                                     c.ccTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGccccc                                                 >  2:1577741/1‑101 (MQ=255)
                                                       ccTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCa                                                >  1:4054247/1‑101 (MQ=255)
                                                          c.tGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGcc                                             >  1:889293/1‑101 (MQ=255)
                                                            tGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGGAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCt                                            >  2:4181296/1‑101 (MQ=255)
                                                            cgg.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCt                                            <  2:4953482/100‑1 (MQ=255)
                                                             gg.C..TTC.GCC.A..TTCCCCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCtt                                           <  1:1570754/101‑1 (MQ=255)
                                                             gg.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCtt                                           <  2:2435210/101‑1 (MQ=255)
                                                             gg.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCtt                                           >  1:2096599/1‑101 (MQ=255)
                                                             gg.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCtt                                           >  2:737650/1‑101 (MQ=255)
                                                              g.C..TTC.GCC.A..TTCCCCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTa                                          <  2:1261246/101‑1 (MQ=255)
                                                                c..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAg                                         >  2:2175012/1‑101 (MQ=255)
                                                                c..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAg                                         >  2:4349370/1‑101 (MQ=255)
                                                                         c.A..TTCCCCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTc                                   <  2:1475335/101‑1 (MQ=255)
                                                                         c.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTc                                   >  2:4588565/1‑101 (MQ=255)
                                                                           a..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCa                                  >  1:3333716/1‑101 (MQ=255)
                                                                           a..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCa                                  >  2:3297943/1‑101 (MQ=255)
                                                                              ttgccccccGGGATATTTGAGGCAAGAT.AAAGGAGCGTGTTTAATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCAc                                 <  2:3618823/97‑1 (MQ=255)
                                                                              ttCCCCCCCGGGATATTTGAGGCGAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCAc                                 <  1:2892906/101‑1 (MQ=255)
                                                                                ccccccgGGGATATTTAGGGAAAGAA.AAGGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGc                               <  1:4410877/99‑1 (MQ=25)
                                                                                cACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGc                               >  1:4262442/1‑101 (MQ=255)
                                                                                     ccGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAggcg                          >  2:3106548/1‑101 (MQ=255)
                                                                                        ggATATTTGAGGGAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAggcggcg                       >  1:2745826/1‑101 (MQ=255)
                                                                                         gATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAggcggcgg                      >  1:5006945/1‑101 (MQ=255)
                                                                                            atTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGt                   >  1:63614/1‑101 (MQ=255)
                                                                                               tGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCCGGCGGCGGCGAAGt                >  1:2226869/1‑101 (MQ=255)
                                                                                               tGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGt                >  2:4757719/1‑101 (MQ=255)
                                                                                               tGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGt                >  2:1681203/1‑101 (MQ=255)
                                                                                                 aGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCCTAGCGTTTCACGCagg.ggcggccggagttc             >  1:1023533/1‑93 (MQ=37)
                                                                                                   gCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCAc            >  1:1465173/1‑101 (MQ=255)
                                                                                                   gCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCAc            <  1:4852730/101‑1 (MQ=255)
                                                                                                    cAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACg           >  2:4908677/1‑101 (MQ=255)
                                                                                                    cAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACg           >  1:365907/1‑101 (MQ=255)
                                                                                                        aT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTga       >  1:3567301/1‑101 (MQ=255)
                                                                                                        aT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTga       >  2:2312638/1‑101 (MQ=255)
                                                                                                         t.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTgag      >  2:4024868/1‑101 (MQ=255)
                                                                                                         t.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTgag      >  1:1269810/1‑101 (MQ=255)
                                                                                                         t.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGCTTCACGCAGGCGGCGGCGTAGCTCACGCTgag      >  2:1639011/1‑101 (MQ=255)
                                                                                                           aaaGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTgaga     >  1:542576/1‑101 (MQ=255)
                                                                                                           aaaGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTgaga     >  2:620804/1‑101 (MQ=255)
                                                                                                            aaGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGTCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGAt    <  2:2578628/101‑1 (MQ=255)
                                                                                                             aGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATc   >  2:4784213/1‑101 (MQ=255)
                                                                                                             aGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATc   <  2:2526476/101‑1 (MQ=255)
                                                                                                              |                                                                                                    
CCGGCGAAAATCGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCC.CCTC.TGG.C..TTC.GCC.A..TTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATCG  >  NC_014034/243533‑243733

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG