Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 854,914 T→G 17.7% intergenic (+103/‑45) hypD → / → hypE hydrogenase expression/formation protein HypD/hydrogenase expression/formation protein HypE

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034854,9140T→G17.7% 3.4 ‑2.4 0.072intergenic (+103/‑45)hypD/hypEhydrogenase expression/formation protein HypD/hydrogenase expression/formation protein HypE
Reads supporting (aligned to +/- strand):  new base (16/1):  ref base (40/15):  total (56/16)
Fisher's exact test strand distribution p-value = 9.50e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CCCAAGCACGGCATCGCCTCCCCCCGGGCCCCCCACCGACCGGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAGCGAT  >  NC_014034/854865‑855014
                                                 |                                                                                                    
ccggggggaggCATCGCCCCCCCCCGGGCCCCCCACCGACCGGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGt                                                   <  1:5256792/92‑1 (MQ=255)
    tctgaggCATCGCCTCCCCCCGGGCCCCCCACCGACCGGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTgag                                               <  1:2486929/96‑1 (MQ=255)
    ggggaggCATCGCCTCCCCCCGGGCCCCCCACCGACCGGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTgag                                               >  2:5967425/6‑101 (MQ=255)
     gggaggCATCGCCTCCCCCCGGGCCCCCCACCGACCCGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTgaga                                              >  2:3558431/5‑101 (MQ=255)
     ggcagccaccGCCCCCCCCCGGGCCCCCCACCGACCGGGGGAGGTGCGACCTCCCCCGGTCCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTgaga                                              <  1:7476422/92‑1 (MQ=255)
      ggaggCATCGCCTCCCCCCGGGCCCCCCACCGACCGGGGGAGGTGGGAACTCCCCCCGCCACCCGCGACGCATTTCGAGCAAGAGGAGATGGCGTTgcgag                                             >  2:8256476/4‑101 (MQ=255)
      ggaggCATCGCCTCCCCCCGGGCCCCCCACCGACCGGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTgagag                                             <  2:2802524/98‑1 (MQ=255)
       gaggCGTCGCCTCCCGCCGGGCCCCCCACCGACCGGGGGAGGGGCGACCTCCCCCGGACCCCCGCGACGGATTTGGAGCCAGATGAGATGGCGTTgggaga                                            >  1:8849860/3‑101 (MQ=255)
       gaggCATCGCCTCCCCCCGGGCCCCCCACCGACCGGGGGGGGGGCGACCTCCCCCGGACCCCCGCGACGGATTTGGAGCAAGATGAGATGGCGTTgggaga                                            >  2:2486856/3‑101 (MQ=255)
         ggCATCGCCTCCCCCCGGGCCCCCCACCGACCGGGGGGGGGGAGACCTCCCCCGGACCCCCACCTCCTAGTTGGAGCAAGATTGGATGGGGGTgagggacg                                          >  2:4740741/1‑99 (MQ=255)
         ggCATCGCCTCCCCCCGGGCCCCCCACCGACCGGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATg                                          <  1:446611/101‑1 (MQ=255)
          gCATCGCCTCCCCCCGGGCCCCCCACCGACCGGGGGAGGGGCGACCTCCCCCGGACCCCCGCGACGGATTTGGAGCAAGATGAGATGGCGTTGAGAGATGa                                         >  2:5330446/1‑101 (MQ=255)
            aTCGCCTCCCCCCGGGCCCCCCACCGACCGGGGGGGGGGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGGGTTGAGAGATGAGc                                       >  2:6479541/1‑101 (MQ=255)
             tCGCCTCCCCCCGGGCCCCCCACCGACCGGGGGGGGGGCGGGCCCCCCCCGACCCCCCCGCCCCATTTTGGGGAAGAAGAGAGGGGGGTGGGGGAGGAGCg                                      >  1:781024/1‑101 (MQ=255)
              cGCCTCCCCCCGGGCCCCCCACCGACCGGGGGGGGGGCGACCTCCCCCGGACCCCCCCGACGTATTTGGAGCAAGATGGGATGGCGTTGAGAGATGAGCgg                                     >  2:2489047/1‑101 (MQ=255)
              cGCCTCCCCCCGGGCCCCCCACCGACCCGGGGGGGGGGGACCTCCCCCCGACCCCCCCCCCGCATTTTGAGCAAGATGAGATGGCGTGGAGAGGAGAGCgg                                     >  1:1314738/1‑101 (MQ=255)
              cGCCTCCCCCCGGGCCCCCCACCGACCCGGGGAGGGGCGACCTCCCCCGGACCCCCGCGACGGATTTGGAGCAAGATGAGATGGCGGTGCGAGATGAGCgg                                     >  2:7107681/1‑101 (MQ=255)
               gCCTCCCCCCGGGCCCCCCACCGACCCGGGGGGGTGGCCCCCCCCCCCGGCCCCCCCGCCGCCTTTGGAGCCAGACAAGAAGGCGGTGAGAGAAGAGCggg                                    >  2:8985101/1‑101 (MQ=255)
                ccTCCCCCCGGGCCCCCCACCGACCGGGGGAGGGGCGACCTCCCCCGGACCCCCGCGACGCATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCgggg                                   >  1:4266671/1‑100 (MQ=255)
                 cTCCCCCCGGGCCCCCCACCGACCGGGGGGGGGGGGGCCCCCCCCCGCCCCCCCCGCCGCATTTGGAGCAAGATGAGATGGCGTTGGGAGATGGGGggggg                                  >  1:2209312/1‑99 (MQ=255)
                 cTCCCCCCGGGCCCCCCACCGACCGGGGGGGGGGCGAACTCCCCCCGGCCCCCGCGACGCATTTGGGGCAAGATGGGGTGGGGGTGGGGGATGAGCgggcg                                  >  2:5302371/1‑99 (MQ=255)
                 cTCCCCCCGGGCCCCCCACCGACCGGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTg                                  >  1:2985359/1‑101 (MQ=255)
                  tCCCCCCGGGCCCCCCACCGACCGGGGGGGGTGCGACCTCCCCCGGACCCCCCCGACGTATTTGGAGCAAGATGAGAAGGGGGTGAGGGATGGGCgggggg                                 >  2:5896011/1‑98 (MQ=255)
                   ccccccGGGCCCCCCACCGACCGGGGGAGGTGCGACCTCCCCCGGTCCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGAc                                <  1:5433183/101‑1 (MQ=255)
                     ccccGGGCCCCCCACCGACCCGGGGGGGGGCGACCTCCCCCCGACCCCCCCGCCGGATTTTGGGGAAGATGAGATGGGGGTGGGGGGATAGCGGGGGACGc                              >  2:5004193/1‑101 (MQ=255)
                       ccGGGCCCCCCACCGACCGGGGGAGGGGCGAGCTCCCCCCGGCCCCCACGGCCCACTTGGAGCAAGAAGAAATGGCGTGTAGAGATAGACGGGTGGCGCTg                            >  2:484905/1‑101 (MQ=255)
                         gggCCCCCCACCGACCGGGGGAGGGGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGc                          <  2:1071285/101‑1 (MQ=255)
                         gggCCCCCCACCGACCGGGGGAGGGGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGc                          >  2:4765552/1‑101 (MQ=255)
                         gggCCCCCCACCGACCGGGGGAGGGGCGACCTCCCCCGGACCCCCCCGACGTATTTGGGGCAAGATGAGANNNNNTTGAGAGATGAGCGGGGGGCGCCGGc                          >  1:5376191/1‑101 (MQ=255)
                         gggCCCCCCACCGACCCGGGGGGGGGGGGCCCCCCCCCGGCCCCCCCGGCGCATTTTGGGGAAGAAGAGATGGGGGTGAGAGATGGGCGGGGGGAGCTGGc                          >  2:2160372/1‑101 (MQ=255)
                          ggCCCCCCACCGACCGGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGca                         >  2:5643779/1‑101 (MQ=255)
                          ggCCCCCCACCGACCGGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGAGGGGGTTGAGAGGTGGGCGGGGGAGGCTGGca                         >  2:4273297/1‑101 (MQ=255)
                          ggCCCCCCACCGACCGGGGGAGGGGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGca                         >  2:7460780/1‑101 (MQ=255)
                          ggCCCCCCACCGACCGGGGGAGGGGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGGTGAGAGATGAGCGGGGGACGCTGGca                         >  2:6504577/1‑101 (MQ=255)
                           gCCCCCCACCGACC.GGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGcaca                       >  1:8736520/1‑101 (MQ=255)
                            ccccccACCGACCGGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGcaca                       <  1:5896246/101‑1 (MQ=255)
                               cccACCGACCGGGGGGGGGGGGGCCCCCCCCCGGCCCCCCCGCCGGATTTGGGGGCAGATGAGAAGGGGGTGAGAGAGGAGCGGGGGGCGCTggggcaagg                    >  2:4431453/1‑94 (MQ=255)
                               cccACCGACCGGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATgg                    >  2:6459748/1‑101 (MQ=255)
                               cccACCGACCGGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATgg                    <  1:4037519/101‑1 (MQ=255)
                                ccACCGACCGGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATggc                   >  2:3799100/1‑101 (MQ=255)
                                ccACCGACCGGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATggc                   >  2:2905178/1‑101 (MQ=255)
                                ccACCGACCGGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATggc                   >  2:2299768/1‑101 (MQ=255)
                                 cgccgaccgGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCGCATggcg                  <  2:55025/99‑1 (MQ=255)
                                 cACCGACCGGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATggcg                  >  2:5068465/1‑101 (MQ=255)
                                   ccgaccgGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATggcggc                >  2:2299337/1‑101 (MQ=255)
                                   ccgaccgGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATggcggc                >  2:2829540/1‑101 (MQ=255)
                                   ccgaccgGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATggcggc                <  1:5330630/101‑1 (MQ=255)
                                     gaccgGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGcggcgg              >  2:1275482/1‑101 (MQ=255)
                                     gaccgGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGGGACGCTGGCACATGGcggcgg              >  1:931571/1‑101 (MQ=255)
                                      ctcgGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGCGATGAGCGGGTGACGCTGGCACATGGCGGCGGt             <  1:6447425/99‑1 (MQ=255)
                                      accgGGGGGGGGGCGACCCCCCCCCGGCCCCCCCGGCGTATTTGGAGCAAGATGAGATGGCGGTGAGAGGTGAGCGGGGGACGCCGGCACATGGcgggggg             >  2:6177431/1‑100 (MQ=255)
                                      accgGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGcggcggg             >  2:1895949/1‑100 (MQ=255)
                                      accgGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGcggcggg             >  1:4991630/1‑100 (MQ=255)
                                        cgGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGGGACGCTGGCACATGGCGGCGGGgg           >  2:5798161/1‑101 (MQ=255)
                                        cgGGGGAGGGGCGACCTCCCCCGGACCCCCGCGCCGTATTTGGAGCAAAATGAGATGGCGTTGAGAGATGAGCGGGAGGCGCTGGCACATGGCGGCGGTgg           >  2:7433415/1‑101 (MQ=255)
                                         gggggAGGTGCGACCTCCCCCGGACCCCCGCGAGGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGc          >  1:6443475/1‑101 (MQ=255)
                                         gggggAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGc          >  1:6342440/1‑101 (MQ=255)
                                         gggggAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGc          >  1:5736486/1‑101 (MQ=255)
                                         gggggAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGc          >  2:6928411/1‑101 (MQ=255)
                                          ggggAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCa         >  1:2499930/1‑101 (MQ=255)
                                          ggggAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCa         >  1:6432922/1‑101 (MQ=255)
                                          ggggAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCa         >  1:3784305/1‑101 (MQ=255)
                                           gggAGTTGCTCCCTCCCCCGCCCCCCCGCGACGTATTTGGAGCAAGATTAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCaa        <  1:8932332/101‑1 (MQ=255)
                                           gggAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCaa        >  1:8688536/1‑101 (MQ=255)
                                            ggAGGTGCGACCTCCCCCGGACCGCCGCGGCGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCaaa       <  1:7671533/101‑1 (MQ=255)
                                            ggAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCaaa       >  2:5124500/1‑101 (MQ=255)
                                             gAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGTGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAg      <  2:3062528/101‑1 (MQ=255)
                                             gAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAg      <  1:2022460/101‑1 (MQ=255)
                                             gAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAg      <  2:8849454/101‑1 (MQ=255)
                                             gAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAg      >  1:2059218/1‑101 (MQ=255)
                                             gAGGTGCGACCTCCCACGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAg      <  1:8564736/101‑1 (MQ=255)
                                             gAGGTCCGACCTCCCCCGGGCCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAg      <  1:790338/101‑1 (MQ=255)
                                              aGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGTACATGGAGACGGAGGCaagac     >  1:5396348/1‑98 (MQ=255)
                                              aGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAGc     >  1:2827748/1‑101 (MQ=255)
                                              aGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAGc     >  2:4277828/1‑101 (MQ=255)
                                              aGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAGc     >  1:407389/1‑101 (MQ=255)
                                              aGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAGc     >  2:7726742/1‑101 (MQ=255)
                                               ggTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCGAAGCg    >  2:6218032/1‑101 (MQ=255)
                                               ggTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAGCg    >  2:3258235/1‑101 (MQ=255)
                                               ggTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAGCg    >  2:239193/1‑101 (MQ=255)
                                               ggTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAGCg    >  2:7725512/1‑101 (MQ=255)
                                                gTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAGCGa   >  1:3489801/1‑101 (MQ=255)
                                                gTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAGCGa   >  2:3243830/1‑101 (MQ=255)
                                                gTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAGCGa   >  1:7059547/1‑101 (MQ=255)
                                                 tGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAGCGAt  >  2:5955751/1‑101 (MQ=255)
                                                 tGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAGCGAt  >  2:8656707/1‑101 (MQ=255)
                                                 |                                                                                                    
CCCAAGCACGGCATCGCCTCCCCCCGGGCCCCCCACCGACCGGGGGAGGTGCGACCTCCCCCGGACCCCCGCGACGTATTTGGAGCAAGATGAGATGGCGTTGAGAGATGAGCGGGTGACGCTGGCACATGGCGGCGGTGGCAAAGCGAT  >  NC_014034/854865‑855014

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 23 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG