Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,185,841 +C 100% intergenic (+114/‑391) mtlK → / → RCAP_rcc02028 mannitol 2‑dehydrogenase/FeoA family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,185,8411.→C100.0% 17.9 129intergenic (+114/‑391)mtlK/RCAP_rcc02028mannitol 2‑dehydrogenase/FeoA family protein
Reads supporting (aligned to +/- strand):  new base (26/93):  ref base (7/2):  total (34/95)
Fisher's exact test strand distribution p-value = 1.07e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.80e-07

Only 100 of 158 total aligned reads displayed.

CCGCCCCGGAACAGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.G.CCCCGCACCCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTGTTTTTGATGGCACAAGAAACGCCGAACCGCCCCTTTC  >  NC_014034/2185741‑2185938
                                                                                                       |                                                                                                  
ccGCCCCGGAACAGGCGGCCCAGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.g                                                                                                     <  2:4583726/101‑1 (MQ=255)
 cGCCCCGGAACAGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.G.c                                                                                                   <  1:1162388/101‑1 (MQ=255)
   ccccGGAACAGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.G.ccc                                                                                                 <  1:5380448/101‑1 (MQ=255)
    cccGGAACAGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.G.cccc                                                                                                <  2:5586145/101‑1 (MQ=255)
     ccGGAACAGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.G.ccccc                                                                                               <  2:1167709/101‑2 (MQ=255)
     ccGGAACAGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.G.ccccc                                                                                               <  1:6701794/101‑2 (MQ=255)
      cGGAACAGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGGGCC.G.cccccc                                                                                              >  2:6678920/1‑99 (MQ=255)
       ggAACAGGCGGCCCCGCTCCCTCCCGGGTGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.G.cccccgc                                                                                             <  1:4966864/101‑4 (MQ=255)
         aaCAGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.G.cccccgcac                                                                                           <  2:5941277/101‑6 (MQ=255)
         aaCAGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.G.cccccgcac                                                                                           <  2:3928525/101‑6 (MQ=255)
          aCAGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCAcc                                                                                           >  2:885818/1‑101 (MQ=255)
          aCAGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCAcc                                                                                           <  1:2894949/101‑1 (MQ=255)
          aCAGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCAcc                                                                                           <  1:2004997/101‑1 (MQ=255)
          aCAGGCGGCCCCGCTCCCGCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCAcc                                                                                           <  2:4801910/101‑1 (MQ=255)
            aGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCAcccc                                                                                         <  2:494161/101‑1 (MQ=255)
             ggcggcCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCAccccc                                                                                        <  2:160742/101‑1 (MQ=255)
             ggcggcCCCCCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCATCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCAccccc                                                                                        <  1:4014894/101‑1 (MQ=38)
              gcggcCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGGCCC.G.cccccgccgccccc                                                                                      >  2:5189913/1‑91 (MQ=255)
               cggcCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCAcccccgg                                                                                      >  2:6207185/1‑99 (MQ=255)
                ggcCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCAcccccggc                                                                                     <  1:3559621/101‑4 (MQ=255)
                gatctcGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCAcccccggc                                                                                     <  2:2662228/96‑4 (MQ=37)
                    ccccccccccccGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTAt                                                                                <  2:6926760/96‑1 (MQ=25)
                    ccGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTAt                                                                                <  2:3392223/101‑1 (MQ=255)
                     ccctccccgcCGGGCGGGCCCTCTTGAGGTGCAAGGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTAtt                                                                               <  2:2296802/99‑1 (MQ=18)
                     cGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTAC.CGGGGGGGC.G.CCCCCCCcccccccccccgcc                                                                               >  1:2662324/1‑93 (MQ=255)
                      gCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGGGCC.GCCCCCCCACCCCCCgcgcgttt                                                                               >  1:3954812/1‑94 (MQ=255)
                      gCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTAttt                                                                              <  1:6190697/101‑1 (MQ=255)
                       ctccctccCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTa                                                                             <  2:2053227/101‑1 (MQ=255)
                       ctccctccCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTa                                                                             <  1:2940446/101‑1 (MQ=255)
                        tgcctcccGGGCGGGGCGGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTTCCGAGGCGTCCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTAt                                                                            <  1:1311992/99‑1 (MQ=25)
                        tccctcccGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTAt                                                                            >  2:5915591/1‑101 (MQ=255)
                           ctcccGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATccc                                                                         <  2:1772922/101‑1 (MQ=255)
                            tcccGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGGTCC.G.CCCCCCCCCCCCCCGCGGGTTTTTcccc                                                                        >  2:4234620/1‑100 (MQ=255)
                             cccGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCaa                                                                       <  2:6496531/101‑1 (MQ=255)
                              ccGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAg                                                                      <  2:6255344/101‑1 (MQ=255)
                               cgggcgggGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGa                                                                     <  1:6685308/101‑1 (MQ=255)
                                gggcgggGCCGCCTGATGTGTAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGaa                                                                    <  2:4174269/101‑1 (MQ=37)
                                gggcgggGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGaa                                                                    <  2:5747974/101‑1 (MQ=255)
                                gggcgggGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGaa                                                                    <  2:4722808/101‑1 (MQ=255)
                                   cgggGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGaaaaa                                                                 <  1:3820657/101‑1 (MQ=255)
                                      ggCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAagca                                                              <  2:5681770/101‑1 (MQ=255)
                                       gccgccTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAagcag                                                             <  2:4220492/101‑1 (MQ=255)
                                       gccgccTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.G.CCCCCGCCCCCCCGGCGTATTTATCCCAAGAAAAagcag                                                             <  1:4331316/101‑1 (MQ=255)
                                        ccgccTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAagcagc                                                            <  1:796400/101‑1 (MQ=255)
                                          gccTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.gg                                                         <  1:6598738/101‑1 (MQ=255)
                                            cTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCt                                                       <  2:6419391/101‑1 (MQ=255)
                                            cTGATGTGCAAAGCGGTTCCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGCATTTATCCCAAGGAAAAACAGC.GGcg                                                       >  2:4343516/1‑100 (MQ=25)
                                             tGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCtt                                                      <  2:2597716/101‑1 (MQ=255)
                                               aTGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGc                                                    <  1:4897975/101‑1 (MQ=255)
                                                tgtgCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGcc                                                   <  2:3705855/101‑1 (MQ=255)
                                                tgtgCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGcc                                                   <  1:4510249/101‑1 (MQ=255)
                                                 gtgCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGccc                                                  <  2:1445121/101‑1 (MQ=255)
                                                 gtgCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGccc                                                  <  2:1735096/101‑1 (MQ=255)
                                                 gtgCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.AGGGGGGCC.G.CCCCCCCCCCCCCCCCGCCTTTTTTCCACGAAAAAACAAC.GGCTGGccc                                                  >  1:6227615/1‑101 (MQ=255)
                                                  tgCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGGGTCC.C.CCCCCCCCCCCCCCCGGGCTTTTTCCCACCAAAAACAACC.GGCCTGccc                                                  >  1:2812034/1‑101 (MQ=255)
                                                  tgCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGcccc                                                 <  2:6536268/101‑1 (MQ=255)
                                                    cAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGa                                               <  1:1230441/101‑1 (MQ=255)
                                                     aaaGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGGATTTATCCCAAGAAAAAGCAGC.GGCTTGTCCCCaa                                              >  1:5638487/1‑101 (MQ=25)
                                                      aaGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAg                                             <  1:2699992/101‑1 (MQ=255)
                                                       aGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGGGCC.G.CCCCCCCCCCCCCCGCGGCTTTTTCCCCAAAAAAAAAAGCGGGGTTGcccccacg                                             >  1:4005322/1‑97 (MQ=255)
                                                         cGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCCGG.CCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGCAGGcc                                          >  1:3261327/1‑101 (MQ=37)
                                                           gTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCtt                                        <  2:3059095/101‑1 (MQ=255)
                                                           gTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCtt                                        <  2:597194/101‑1 (MQ=255)
                                                           gTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCCTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCtt                                        <  2:3454879/101‑1 (MQ=25)
                                                            tttCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTg                                       <  2:4005147/101‑1 (MQ=255)
                                                             ttCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTTCCCCCGCGGCCttcg                                      >  1:4176471/1‑99 (MQ=21)
                                                             ttCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTGt                                      <  1:5353638/101‑1 (MQ=255)
                                                              tCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTGtt                                     <  2:3623453/101‑1 (MQ=255)
                                                               cGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTGttt                                    <  1:3431267/101‑1 (MQ=255)
                                                                gCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCTTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTGtttt                                   >  1:832534/1‑101 (MQ=37)
                                                                gCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTGtttt                                   <  2:1787695/101‑1 (MQ=255)
                                                                   ccggcggTGCCGAGGCGACCCTGC.GGGGCCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTGttttgat                                <  2:1749330/100‑4 (MQ=37)
                                                                    cggcggTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCCAAGGGCTTGttttgatg                               >  1:4847851/1‑97 (MQ=25)
                                                                     ggcggTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATgg                             <  1:555829/101‑1 (MQ=37)
                                                                     ggcggTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGAAGC.GGCTTGGCCCGAAGGGCTTGTTGTTgtgt                              >  1:2283164/1‑98 (MQ=25)
                                                                      gcggTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGc                            <  2:6227296/101‑1 (MQ=37)
                                                                        ggTGCCGAGGCGGCCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGcac                          <  1:5190173/101‑1 (MQ=25)
                                                                         gTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCCCAAGGCTCG.TTTTGATGGcaca                         >  2:2304909/1‑101 (MQ=21)
                                                                         gTGCCGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCACGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGcaca                         <  1:885842/101‑1 (MQ=37)
                                                                            ccGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGa                      <  1:4408061/101‑1 (MQ=37)
                                                                             cGAGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGaa                     <  1:7060676/101‑1 (MQ=37)
                                                                              gAGGCGACCCTGC.GGGGGCGCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGaaa                    <  1:3566135/101‑1 (MQ=37)
                                                                               aGGCGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAAc                   <  1:4234815/101‑1 (MQ=37)
                                                                                  cGACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGcc                <  1:402922/101‑1 (MQ=37)
                                                                                   gACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCg               <  2:3566576/101‑1 (MQ=37)
                                                                                   gACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCCAAGGGCTTG.TTTTGATGGCACAAGAAAAGCCg               >  1:2046490/1‑101 (MQ=25)
                                                                                   gACCCTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCCAAGGCCTTG.TTTTGATGGCACAAGAAACGCCg               >  1:6535549/1‑101 (MQ=25)
                                                                                      ccTGC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAAc            <  1:2751019/101‑1 (MQ=37)
                                                                                         gC.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCCAAGGCCTTG.TTTTGATGGGACAAGAAAAGCCGAACCGc         >  1:31410/1‑101 (MQ=25)
                                                                                          c.GGGGGCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGGCTTG.TTTTGATCGCACAAAAAAAGCCGAAACggc        >  2:3291219/1‑99 (MQ=21)
                                                                                            gggggCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGccc       >  1:4855030/1‑101 (MQ=37)
                                                                                             tctgCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGcccc      <  2:2830474/98‑1 (MQ=37)
                                                                                             ggggCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGGACAAGAAACGCCGAACCGcccc      >  2:5186685/1‑101 (MQ=37)
                                                                                             ggggCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGcccc      <  1:1566177/101‑1 (MQ=37)
                                                                                              gggCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAAAGCCCAACCGccccc     >  2:73418/1‑100 (MQ=25)
                                                                                              gatcTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGccccc     <  2:1005038/98‑2 (MQ=37)
                                                                                                gCTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCctt   <  2:5545759/101‑1 (MQ=37)
                                                                                                 cTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCcttt  <  2:3954637/101‑2 (MQ=37)
                                                                                                 cTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCcttt  <  1:5806815/101‑2 (MQ=37)
                                                                                                 cTCC.GCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCCAACCGCCCCcttt  >  2:5491341/1‑100 (MQ=37)
                                                                                                       |                                                                                                  
CCGCCCCGGAACAGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGC.GGGGGCTCC.G.CCCCGCACCCCCCGGCGTATTTATCCCAAGAAAAAGCAGC.GGCTTGCCCCGAAGGCCTTGTTTTTGATGGCACAAGAAACGCCGAACCGCCCCTTTC  >  NC_014034/2185741‑2185938

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG