Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,965,269 +G 4.0% intergenic (+37/+55) bdhA → / ← RCAP_rcc01818 3‑hydroxybutyrate dehydrogenase/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,965,2691.→G 4.0% 6.2 ‑1.9 ‑0.2100intergenic (+37/+55)bdhA/RCAP_rcc018183‑hydroxybutyrate dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (57/39):  total (61/39)
Fisher's exact test strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.92e-01

Only 100 of 140 total aligned reads displayed.

CGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GCGGG.CT.C.C.G.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTC  >  NC_014034/1965169‑1965366
                                                                                                         |                                                                                                       
cGAGGCGGCGGCGCAGATCACCGGCACTACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGCGATTTCTCCCCGGGCCAACCCC.GAGC.CG.ACccg.cc                                                                                                          >  2:2265359/1‑101 (MQ=255)
 gAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCCGGGGAAGCCC.GAGG.CG.ACtcg.tc..g                                                                                                       >  2:823105/1‑101 (MQ=255)
 gAGGCGGCGACGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACccg.ccc                                                                                                         >  1:7685642/1‑100 (MQ=255)
  aGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GCGG.CG.ACccg.cc..gg                                                                                                      >  2:2145851/1‑100 (MQ=255)
      ggcggcgCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACccg.cccggcgg                                                                                                    >  2:7594507/1‑95 (MQ=255)
      ggcggcgCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGgggg                                                                                                  >  1:5496589/1‑100 (MQ=255)
      ggcggcgCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.c                                                                                                 >  2:1120463/1‑101 (MQ=255)
        cggcgCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GGGG.CG.ACCCG.CC..CGGGG.CT.c                                                                                              >  2:1114634/1‑101 (MQ=39)
         ggcgCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGgggggcc                                                                                               >  2:1965779/1‑97 (MQ=255)
         ggcgCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGgggggcc                                                                                               >  1:8620595/1‑97 (MQ=255)
           cgcAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GT.G.C.G.c                                                                                        >  2:8181747/1‑101 (MQ=38)
           cgcAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CG.C.C.G.c                                                                                        >  1:7583846/1‑101 (MQ=255)
             cAGGTCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CG.C.C.C.ccc                                                                                      >  1:7229250/1‑101 (MQ=38)
             cAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG..C..GCggggggccgccc                                                                                          >  2:2114038/1‑92 (MQ=37)
              aGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.G.cccc                                                                                     >  2:3442388/1‑101 (MQ=255)
               gATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGggggggccccccc                                                                                         >  1:3620908/1‑91 (MQ=38)
               gATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GT.C.C.G.ccccc                                                                                    >  2:7459174/1‑100 (MQ=255)
                  cACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GG.C.C.C.cccccccg                                                                                 >  2:410114/1‑97 (MQ=37)
                  cACCGGCACGACGATGTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GC.C.C.G.ccccccgc                                                                                 >  1:7239045/1‑97 (MQ=25)
                    ccGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GC.C.C.G.ccccccggcc                                                                               >  1:8240335/1‑95 (MQ=37)
                     cGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCGGAccc                                                                               >  2:5941075/1‑101 (MQ=37)
                          cgacgaTTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..CGGGG.CT.CGG.G.CCCCGGTCCCCCGAg                                                                          >  2:501321/1‑101 (MQ=37)
                           gacgaTTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..CGGGG.GT.C.C.GCCCCCGGACCCCCGAgg                                                                         >  2:5963352/1‑101 (MQ=37)
                               aTTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCGCCC..GGGGG.CC.C.C.C.CCCCCGGccccccgggggat                                                                     >  1:6586114/1‑93 (MQ=25)
                               aTTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..CGGGG.CT.C.C.GCCCCCCGACCCCCGAGGGAtt                                                                     >  1:2802488/1‑101 (MQ=25)
                                tttCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GTCC.C.C.CCCCCGGCCCCCCAGGGAttt                                                                    >  1:5936307/1‑101 (MQ=21)
                                 tggtggggggggcggcGGGCGGGCGTTGGGAGAGATTCCGCCGGGGGGCACCCCC.GGGG.CC.CCCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTc                                                                   <  2:7221878/91‑1 (MQ=9)
                                  tCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GT.C.C.GCCCCCCGGCCCCCCAGGGATTTCa                                                                  >  2:3288021/1‑101 (MQ=21)
                                    tgtgGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GT.T.C.GCCCCCCGGCCCCCCCGGTATTTTAgg                                                                >  2:6616805/1‑101 (MQ=18)
                                    tgtgGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGCGGGATTTCAgg                                                                >  2:1662911/1‑101 (MQ=25)
                                    tgggggtGGCGGCGGGCCGGCGCTGGGAGACTTCCCCCCGGGGGGCACGCCG.GGGC.CC.CCCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAgg                                                                <  1:2897386/95‑1 (MQ=2)
                                     gtgGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GTCC.C.C.CCCCCGGCCCCCCAGGGGTTTCAGGc                                                               >  1:5802370/1‑101 (MQ=21)
                                     gtgGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CTCC.C.C.CCCCCGACCCCCGCGGGATTTCCGGc                                                               >  1:2249429/1‑101 (MQ=21)
                                      tgGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GT.C.C.GCCCCCGGACCCCCCGGGGATTTCAGGCa                                                              >  2:5636575/1‑101 (MQ=21)
                                        gggggcggcGGGCGGGCGTGGTGAGGGTTTCCGCCGGGGGGCACCCCC.GGGC.CC.CCCCG.CC..GGGGG.CT.C.T.GCCCCCGGACCCCCGAGGTATTTCAGGCAAg                                                            <  1:6107488/98‑1 (MQ=2)
                                        gATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GT.C.C.GCCCCCGGGCCCCCGAGGGATTTCAGGCAAg                                                            >  2:7276820/1‑101 (MQ=25)
                                         gggggggcGGGGCGGGGCTGGGGGAGCTTCCGCCGGGGGGCACGCCG.GGGG.CC.CCCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGa                                                           <  2:5642761/96‑1 (MQ=9)
                                         ggggcggcGGGGCGGCGCTGTGAAAGTTTCCCCCGCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGa                                                           <  1:473817/99‑1 (MQ=18)
                                         aTGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCGGGCCCCCCGGGGATTTCCgggcagg                                                           >  2:5237564/1‑96 (MQ=18)
                                          tGGCGGCTGGACGGCGCTGTGAGAGATATCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GT.C.CGG.CCCCCGCCCCCCGCGGGAtttgcagccagaa                                                          >  2:2112683/1‑91 (MQ=16)
                                          gggcggcGGGCGGGGGTTGTGAGAGTTTCCCCCGGGGGGCACGCCG.GGGG.CC.CCCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGAt                                                          <  2:5539946/100‑1 (MQ=9)
                                           ggcggcTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGGATTTCAGGGCAGATg                                                         >  2:6675437/1‑101 (MQ=25)
                                              ggcTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGGATTTCAGGCAAGATGaaa                                                      >  1:1299692/1‑101 (MQ=21)
                                               gggggcgggCGGTGGGAGAGTTTCCGCCGGGGGCCAACCCG.GGGC.CC.CCCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAag                                                     <  2:4718566/94‑1 (MQ=2)
                                               gcgggccGGCGCTGTGAGAGATCCCGCCTCGGGGCACGCCC.GGGC.CC.CCCCG.CC..GGGTG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAag                                                     <  2:3196919/95‑1 (MQ=16)
                                               gcTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGGCCC.C.C.C.CCCCCGGCCCCCCCGGGATTTTCGGCAAGATGaaaa                                                     >  1:621864/1‑100 (MQ=12)
                                                cTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CCGGGGGGG.CC.C.C.C.CCCCGGCCCCCCCGGGGTTTTCAGCGAAGAGGaaaa                                                     >  1:2513049/1‑100 (MQ=12)
                                                         gcTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CCGGGGGGG.GC.C.C.C.CCCCGCGCCCCCGCGGGGTTTCAGCCAGGAGGAAAAAGGGgggcc                                            >  2:2069025/1‑100 (MQ=9)
                                                         gcTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GTCC.C.C.CCCCCGGCCCCCGCGGGATTTCAGGCAAGAAGGAAGAGGGgcgggc                                           >  2:3458240/1‑98 (MQ=16)
                                                         gcTGTGAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GT.C.C.GCCCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGGGCTc                                           >  2:4152718/1‑101 (MQ=21)
                                                           tgtgAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCGCCC..GGGGG.CC.C.C.C.CCCCCGGCCCCCCCGGGGTTTTAAGGGAGAAGGAAGAGGGGGGCTCAg                                         >  1:5986250/1‑101 (MQ=9)
                                                            gtgAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GTCC.C.C.CCCCCGGCCCCCCAGGGATTTTCGGGAAGATGAAAGAGGGgggccccgg                                        >  2:6929686/1‑96 (MQ=12)
                                                            gtgAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GT.C.C.GCCCCCGGGCCCCCCAGGGATTTTAGGCAAGATGAAAGAggggggggcagg                                        >  2:7787975/1‑93 (MQ=18)
                                                            gtgAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCAgg                                        >  2:8829520/1‑101 (MQ=21)
                                                             tgAGAGATTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCCggg                                       >  2:5277745/1‑100 (MQ=21)
                                                              gagagaTTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GT.C.C.GCCCCCGGGCCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCCGAggcg                                      >  2:3666811/1‑101 (MQ=21)
                                                               agagaTTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GG.C.C.GCCCCCCGCCCCCCGAGGAATTTCAGGCAAGATGAAAGAGGGGCGGGCAggggg                                     >  2:115514/1‑101 (MQ=18)
                                                               agagaTTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCAggcgg                                     >  2:1545964/1‑101 (MQ=21)
                                                               agagaTTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCCAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCCggcgg                                     >  2:8234648/1‑101 (MQ=21)
                                                                gagaTTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCGCACCCCCCCGGGATTTCAGGCAAGATGAAAGAGGGGGGCTCAggccgc                                    >  2:8400765/1‑101 (MQ=18)
                                                                  gaTTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GTCC.C.C.CCCCGGGCCCCCCCGGGATTTCTGGCAAGATGGAAGAGGGGGGCGCAGGGGGCct                                  >  1:8755106/1‑101 (MQ=12)
                                                                   aTTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCCAGGGATTTCAGGCAAGATGAAAGAGGGGCGCCCAGGCGGCctt                                 >  1:5798850/1‑101 (MQ=21)
                                                                   aTTTCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCCGGCCCCCCAGGGCTTTCGGGCAAGATGAAAAAGGGGGGCTCAGGCGGcccc                                 >  1:8385315/1‑99 (MQ=16)
                                                                    tctccgcgcGGGGCAACCCCGGCGG.CC.CCCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCctt                                 <  2:8669000/93‑1 (MQ=18)
                                                                     tcgcgcggggggCCGGCCC.GGGG.CGCCCCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCcttc                                <  2:5476407/93‑1 (MQ=18)
                                                                      tcgcCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GT.C.C.GCCCCCCGACCCCCCAGGCATTTCCGGCAAGATGAAAAAGGGGCGCGCAGGCGGCcttctt                              >  1:4035364/2‑101 (MQ=21)
                                                                       cgcGTCGGGGCAAGCGC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGAGCGCTCAGGCGGCCttcttc                             >  1:8825620/1‑101 (MQ=21)
                                                                        gcCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGACGGATTTCCGGCAAGAAGAAAGAGGGGCGCTCAGGCGGGCTTCTTcg                            >  1:8972582/1‑101 (MQ=21)
                                                                         gcccGGGGCACCCCCGGACG.CC.ACCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTcg                            <  1:5419641/98‑1 (MQ=18)
                                                                          ccgggggCCAGCCC.GGGGCCC.CCCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTcgc                           <  2:5360321/98‑1 (MQ=21)
                                                                           gggggccgcgCCC.GGCG.CG.CCCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTcgccg                         <  2:8070416/92‑1 (MQ=18)
                                                                            ggggcCAAGCCC.GGCG.CC.CCCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTcgccgc                        <  1:861366/100‑1 (MQ=21)
                                                                            cGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GT.C.C.GCCCCCCGACCCCCCAGGGATTTCCGGCAAGATGGAAGAGGGGGGCTCAGGCGGCCTTTCTcgccgc                        >  2:8142677/1‑101 (MQ=21)
                                                                            cGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGGGCTCAGGCGGCCTTCTTcgccgc                        >  2:9003112/1‑101 (MQ=21)
                                                                            cGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCCGGCCCCCCCGGGATTTCAGGCCAGATGAAAGAGGGGCGCTCAGGCGGCCCTCTTcgccgc                        >  1:1493935/1‑101 (MQ=21)
                                                                             ggggCAAGCCC.GAGG.CG.ACCCG.CCG.GGGGG.GC.C.C.C.CCCCCGACCCCCCCGGGATTTCAGGCAAAAAGAAAAAAGGGCGCTTCGGCGGCCTTCTTCgccccc                       >  2:2724938/1‑101 (MQ=12)
                                                                              gggCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CTCC.C.C.CCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCgccgccc                      >  2:8235987/1‑100 (MQ=25)
                                                                                  aaGCCC.GAGG.CG.ACCCG.CC..GGGGG.CTCC.C.C.CCCCCGGCCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCCGGCGGCCTTCTTCGCCGCCCTGcc                  >  2:8364642/1‑101 (MQ=21)
                                                                                  aaGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCCTCTTCGCCGCCTTGcc                  >  2:1371645/1‑101 (MQ=21)
                                                                                  aaGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCCGAGCCCCCCGGTATTTCAGGCAAGATGAAAGAGGGGCGGGCGGGCGGGCTTCCTCcccccccgtcc                  >  2:2834626/1‑96 (MQ=9)
                                                                                     ccc.GAGG.CG.ACCCG.CC..GGGGG.CTCC.C.C.CCCCCGACCCCCCAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCCTGccccg               >  2:7044733/1‑99 (MQ=21)
                                                                                       c.GAGG.CG.ACCCG.CC..GGGGG.GTCC.C.C.CCCCCGACCCCCCAGGGAATTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCCTGCCCgggg             >  2:5625296/1‑99 (MQ=21)
                                                                                       c.GAGC.CG.CCCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGgtg             <  2:3202607/101‑1 (MQ=25)
                                                                                        ccggggcccccccG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCgg               <  1:2298031/91‑1 (MQ=25)
                                                                                         gAGG.CG.ACCCG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGgtgt            >  2:7238006/1‑101 (MQ=25)
                                                                                         ccggccccccccG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGgt              <  1:1943294/92‑1 (MQ=25)
                                                                                         ccgcggccccccG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGgt              <  1:3687327/92‑1 (MQ=25)
                                                                                          aGG.CG.ACCCG.CC..GGGGG.GT.C.C.GCCCCCCGACCCCCGAGGGATTTCCGGCAAGATGAAAGAGGGGCGCCCAGGCGGCCCTCTTCGCCGCCCTGCCCCgtggc           >  2:2894662/1‑99 (MQ=18)
                                                                                            gggccccccG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtc           <  2:4123805/95‑1 (MQ=25)
                                                                                             gccgccccG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtct          <  2:8972161/96‑1 (MQ=25)
                                                                                             gcccccccG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtct          <  2:8754701/96‑1 (MQ=25)
                                                                                              cccccccG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtctc         <  1:131192/97‑1 (MQ=25)
                                                                                              cccccccG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtctc         <  1:5277929/97‑1 (MQ=25)
                                                                                              cccccccG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtctc         <  1:5406111/97‑1 (MQ=25)
                                                                                              cccccccG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtctc         <  1:7277124/97‑1 (MQ=25)
                                                                                               ccccccG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCg        <  1:8056243/98‑1 (MQ=37)
                                                                                                  cccG.CC..GGGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGc     <  1:1662958/101‑1 (MQ=255)
                                                                                                    cg.cc..ggGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGctt   <  2:3738524/101‑1 (MQ=255)
                                                                                                     g.cc..ggGGG.CT.C.C.GCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGcttc  <  2:7718799/101‑1 (MQ=255)
                                                                                                     g.cc..ggGGG.CC.C.C.GCCCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCCTCTTCGCCGCCTTGCCCGGGGTCTCGGGgttc  >  2:6491933/1‑101 (MQ=25)
                                                                                                         |                                                                                                       
CGAGGCGGCGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAGCCC.GAGG.CG.ACCCG.CC..GCGGG.CT.C.C.G.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTC  >  NC_014034/1965169‑1965366

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 25 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG