Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,965,279 +C 89.5% intergenic (+47/+45) bdhA → / ← RCAP_rcc01818 3‑hydroxybutyrate dehydrogenase/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,965,2791.→C89.5% 13.1 ‑0.9 ‑0.857intergenic (+47/+45)bdhA/RCAP_rcc018183‑hydroxybutyrate dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (30/21):  ref base (5/1):  total (35/22)
Fisher's exact test strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.76e-01

AGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GCGGGC.T.CCG..C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGC  >  NC_014034/1965183‑1965374
                                                                                                       |                                                                                                   
aGATCACCTGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGG.C.CCC..c.c.c.c                                                                                             >  1:3578392/1‑101 (MQ=38)
  aTCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.T.CCG..c.c.c.ccg                                                                                           >  2:3923592/1‑99 (MQ=255)
  aTCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.G.CCG..c.c.c.ccc                                                                                           >  1:7067367/1‑99 (MQ=255)
    cACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGG.G.CCC..c.c.c.cccgg                                                                                         >  1:749567/1‑97 (MQ=37)
        ggCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..CGGGGC.T.CCG..c.c.c.cccgacccc                                                                                     >  1:6959951/1‑93 (MQ=25)
          cACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGG.G.CCC..c.c.c.ccccccccccc                                                                                   >  1:538156/1‑91 (MQ=25)
           acgacgATTTCTGTGGATGGCGGCTGGTCGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGG.TACCGC.C.C.C.CGGAccccga                                                                                    >  2:3786489/1‑99 (MQ=25)
           acgacgATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGG.T.CCGC.C.C.C.CGGACCCCCgc                                                                                   >  2:3259254/1‑100 (MQ=37)
             gttgtttgcTGTGGGGGGGGGGTGGGCGGCGCTGTGAGGTTTTCCGCCGGGGGGCAGG.CCCGGGGCCCCCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAgg                                                                                 <  2:4813308/93‑1 (MQ=2)
              acgaTTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.C.CCGCCC.C.C.CGGAccccccggg                                                                                 >  1:4109494/1‑97 (MQ=21)
                 aTTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CCGGGGGGGC.C.CCC..C.C.C.CGGCcccccccgggtt                                                                              >  2:3678283/1‑94 (MQ=21)
                 aTTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.C.GCGC.C.C.C.CCGGCCCCCGAgcacgt                                                                             >  2:2282712/1‑96 (MQ=21)
                   ttCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.T.CCGC.C.C.C.CCGGCCCCCGAGGGATTTc                                                                           >  2:2108575/1‑101 (MQ=25)
                    tgtggggggggcggcGGGACGGCGGTGGGAGAGATTCCGCCGGGGGGCACC.CCGGGGGCGCCCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCa                                                                          <  2:1663383/92‑1 (MQ=9)
                     cTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCGCCC..GGGGGC.C.CCC..C.C.C.CCGGCCCCCCCGGGAATTCAg                                                                         >  2:2444425/1‑101 (MQ=21)
                       gtgGATGGCGGCTGTACGGCGCTGTGAGAGATTTCGCCTCGGGACAAG.CCCGAGGCGACCCGCCC..GGGGGG.C.CCGC.C.C.C.CGGCCCCCCGGGGAGTTTCAgc                                                                        >  1:6289584/1‑100 (MQ=16)
                       gtgGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGGAAG.CCCGAGGCGACCCG.CC..CGGGGCGC.CCC..C.C.C.CCGGCCCCCGAGGGGGTTCAGGc                                                                       >  2:6601945/1‑101 (MQ=19)
                       gtgGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGG.G.CCGC.C.C.C.CGGGCCCCCCCGGGATTTCAGGc                                                                       >  2:6753129/1‑101 (MQ=21)
                       gtgGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCaagg                                                                       >  2:3279957/1‑98 (MQ=25)
                        tgGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCGCCC..GGGGGT.C.CCC..C.C.C.CCGGCCCCCCAGGGGTTTCAgggc                                                                      >  1:7049404/1‑99 (MQ=16)
                        tgGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGGCCCCCGAGGGATTTCAGGCa                                                                      >  2:6582998/1‑101 (MQ=25)
                           aTGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.G.CCGC.C.C.C.CGGACCCCCCGGGGATTTTAGGCAAGa                                                                   >  2:3130774/1‑101 (MQ=21)
                            gggcggcGGGGCGGCGCGGGGAGAGTTTCCGCCGGGGGGCACC.CCCGGGGCCCCCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGAt                                                                  <  1:5512047/100‑1 (MQ=9)
                                ggcTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGGATTTCAGGCAAGATGaaa                                                              >  2:1148386/1‑101 (MQ=21)
                                ggcTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.T.CCGC.C.C.C.CCGGCCCCCGAGGGATTTCAGGCAAGATGaaa                                                              >  2:6997681/1‑101 (MQ=21)
                                  ggggcgggCGCTGGGGGGGCTTCCGCCCCGGGCCACC.CCGGGGGCCCCCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAaga                                                            <  1:5668433/95‑1 (MQ=2)
                                  cTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAaga                                                            >  1:3332131/1‑101 (MQ=21)
                                  cTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGGATTTCAGGCAAGATGAAaga                                                            >  1:1063692/1‑101 (MQ=21)
                                   ggggcGGCGCTGGGAGGGATTCCGCGCGGGGGCACG.CCGGGGGCCCCCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAagag                                                           <  2:5528325/97‑1 (MQ=9)
                                    ggACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CCGGGGGGGC.C.CCC..C.C.C.CGGCCCCCCCGGGGATTTTCGGCAGAATGaagagg                                                           >  2:6055765/1‑97 (MQ=9)
                                    ggACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGGCCCCCCAGGGATTTCAGGCAAGATGAAAGAgg                                                          >  2:1966842/1‑101 (MQ=21)
                                    ggACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.C.CCGC.C.C.C.CGGCCCCCCGGGGGATTTCAGGCAAGATGAAAGAgg                                                          >  2:1135338/1‑101 (MQ=21)
                                       gggCGCGGGGAGAGGTGCCGCCGGGGGGGAAG.CCCGGGGCGCCCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGgc                                                       <  1:5970606/100‑1 (MQ=12)
                                       cGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGG.C.CCGC.C.C.C.CCGACCCCCGAGGGATTTCAGGCAAGATGAAAGAggggg                                                       >  1:6749966/1‑100 (MQ=21)
                                       cGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGGATTTCAGGCAAGATGAAAAAGGGgc                                                       >  2:2052166/1‑101 (MQ=21)
                                            cTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGG.T.CCGC.C.C.C.CGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAgggggggcca                                                  >  1:5379402/1‑95 (MQ=21)
                                            cTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGGCCCCCGCGGGATTTCAGGCAAGATGAAAGAGGGGGGCTCa                                                  >  2:1104549/1‑101 (MQ=21)
                                             tgtgAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGG.T.CCGCCC.C.C.CGGACCCCC.AGGGATTTCAGGGAAGATGAAAGAGGGgggcgccg                                                 >  2:5590010/1‑97 (MQ=17)
                                             tgtgAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCCAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCAg                                                 >  1:2814987/1‑101 (MQ=21)
                                             tgtgAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.C.ACGC.C.C.C.CCGACCCCCGAGGTATTTCAGGCAAGATGCAAGAGGGgcgcggag                                                 >  2:4226006/1‑97 (MQ=21)
                                              gtgAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGG.T.CCGC.C.C.C.CCGGCCCCCCAGGGATTTCAGGCAAGATGAAAAAGGGGGGCTCAgg                                                >  1:6693854/1‑101 (MQ=18)
                                               tgAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGG.G.CCGC.C.C.C.CGGACCCCCCAGGGCATTCAGGAAAAAAGAAAAAAGGGGGGTCCggc                                               >  2:6207568/1‑101 (MQ=9)
                                               tgAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGGATTTCAGGCCAGATGAAAGAGGGGCGCTCCggg                                               >  2:4330509/1‑100 (MQ=21)
                                               tgAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..CGGGGC.T.CCG..CGC.C.CGGACCCCCGAGGGATTTCACGCTAAGTGAAAGATGGGGGCTctcga                                               >  1:1201134/1‑97 (MQ=16)
                                                gagagaTTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAggcg                                              >  1:1122943/1‑101 (MQ=21)
                                                 agagaTTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..CGCGGC.T.CCGC.C.C.C.CGGCCACCCGCGGTTTTTGCGGGAAGATGAAAGAGGGGCGCTCCggccg                                             >  1:3070346/1‑99 (MQ=16)
                                                   agaTTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCAggcgggc                                           >  2:2265389/1‑99 (MQ=21)
                                                       ttcgcCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGGCCCCCCAGGGATTTCCGGCAAGATGAAAGAGGGGCGCTCAGGGGGCcttct                                       >  1:3203575/3‑101 (MQ=21)
                                                        cccccggggggCAAGCCCCGGGCCCCCCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCcttct                                       <  2:5922423/94‑1 (MQ=18)
                                                          gcCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.TCCCC..C.C.C.CCGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCttcttcc                                    >  2:1645568/1‑100 (MQ=25)
                                                             gggggcCAGG.CCCGCGGCCCCCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTcgtcg                                 <  1:7001062/99‑1 (MQ=18)
                                                               ggggCAAG.CCCGAGGCGACCCG.CC..GGGGGG.TCCCC..C.C.C.CCGGCCCCCCAGGGATTTCAGGCAAAATGAAAGAAGGGGGCTCCGGCGGGCCTCTTCtccccc                               >  2:1956425/1‑101 (MQ=12)
                                                               ggggCAAG.CCCGAGGCGACCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCgccgcc                               >  2:3301308/1‑101 (MQ=255)
                                                                gggCAAG.CCCGAGGCGACCCG.CCG.GGGGGC.C.CCC..C.C.C.CCGGCCCCCCAGGGGTTTCCGGCAAGATGAAAGAGGGGCGCTCCGGCGGGCCTCTTCgccgccc                              >  2:5456036/1‑100 (MQ=19)
                                                                 ggCAAG.CCCGAGGCGACCCG.CC..GGGGGC.T.CCGCCC.C.C.CGGACCCCC.AGGGATTTCAGGCAAGATGCAAGAGGGGCGCTCAGGCGGCCCTCTTCgccgccct                             >  1:6737006/1‑99 (MQ=21)
                                                                   gccccccagAGGCGCCCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTg                            <  1:4226201/93‑1 (MQ=21)
                                                                    aaG.CCCGAGGCGACCCG.CC..GGGGGG.T.CCGC.C.C.C.CCGGCCCCCCCGGGATTTTCGGGCAGAAGGAAGAGGGGGGCGCCGGGGGGCCTCTTCGCCCCCCTGcc                          >  1:1672209/1‑101 (MQ=0)
                                                                       ccgcgcggCGCCCCG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCg                        <  2:4632253/95‑1 (MQ=25)
                                                                        ccggcgcccccccG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCgg                       <  2:7012487/91‑1 (MQ=25)
                                                                         cggggcccccccG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGgt                      <  1:4330731/92‑1 (MQ=25)
                                                                          cGAGGCGACCCG.CC..GGGGGC.TCCCC..C.C.C.CCGACCCCCGAGGGATTTCAGGCCAGAAGAAAGAAGGGGGCCCAGGCGGGCCTCTTCTCCGCCCTGccccggg                     >  2:2743228/1‑97 (MQ=12)
                                                                           gcgcccccccG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGgtgt                    <  2:2081807/94‑1 (MQ=25)
                                                                           gAGGCGACCCG.CC..GGGGGC.TCCCC..C.C.C.CCGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCCGGCGGCCCTCCTCCGCGCCCTGCCCggggc                    >  2:6781701/1‑98 (MQ=12)
                                                                              cccccccG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtctc                 <  1:7017370/97‑1 (MQ=25)
                                                                                 ccccG.CC..GGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCggg              <  1:4682194/100‑1 (MQ=255)
                                                                                 ccccG.CC..GGGGGC.T.CCG..C.C.CGCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCggg              <  1:5607254/100‑1 (MQ=255)
                                                                                   ccg.cc..gGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGct            <  1:2490659/101‑1 (MQ=255)
                                                                                   ccg.cc..gGGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGct            <  1:60196/101‑1 (MQ=255)
                                                                                    cg.cc..ggGGGC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGctt           <  2:4999147/101‑1 (MQ=255)
                                                                                          cgggggC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGcttctt        <  2:6866606/98‑1 (MQ=255)
                                                                                           gggggC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGGCCGGTGTCTCGGGCTTCTTg       >  1:4376058/3‑101 (MQ=255)
                                                                                           gggggC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTg       >  2:4448708/3‑101 (MQ=255)
                                                                                           gggggC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTg       >  1:5424354/3‑101 (MQ=255)
                                                                                           gggggC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTg       >  1:7056432/3‑101 (MQ=255)
                                                                                           gggggC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTg       >  2:5813367/3‑101 (MQ=255)
                                                                                           gggggC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTg       <  1:99470/99‑1 (MQ=255)
                                                                                           gggggC.T.CCG..C.CGC.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTg       >  1:4143678/3‑101 (MQ=255)
                                                                                            ggggC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt      >  1:1973075/2‑101 (MQ=255)
                                                                                            ggggC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt      >  1:7097575/2‑101 (MQ=255)
                                                                                            ggggC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt      >  1:5256579/2‑101 (MQ=255)
                                                                                             gggC.T.GCGC.C.C.C.CGGGCCCCCGAGGTATTTCAGGCAATATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTc     >  1:6856674/1‑101 (MQ=255)
                                                                                             gggC.T.CCGC.C.C.C.CGGACCCCCGAGGCATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTc     >  1:4110734/1‑101 (MQ=255)
                                                                                              ggC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCg    <  2:6605791/101‑1 (MQ=255)
                                                                                              ggC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCg    <  2:311260/101‑1 (MQ=255)
                                                                                              ggC.T.CCGC.C.C.C.CGGACCCCCGAGGCATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTCCGCCGCCTCGGCCGGTGGCTCTGCCTTCTTGTCg    >  2:1183511/1‑101 (MQ=25)
                                                                                               gC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCgg   >  1:6614943/1‑101 (MQ=255)
                                                                                                c.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGc  >  2:5636074/1‑101 (MQ=255)
                                                                                                       |                                                                                                   
AGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GCGGGC.T.CCG..C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGC  >  NC_014034/1965183‑1965374

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG