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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,965,2761.→C 5.1% 3.7 ‑1.3 ‑0.259intergenic (+44/+48)bdhA/RCAP_rcc018183‑hydroxybutyrate dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (34/22):  total (37/22)
Fisher's exact test strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.92e-01

CGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GCGGGC.T.CCG..C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGC  >  NC_014034/1965177‑1965374
                                                                                                         |                                                                                                       
cggcgCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.T.c                                                                                                        >  1:7106635/1‑101 (MQ=255)
 ggcgCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCCAG.CCCGAGGCGACCCG.CC..CGGGGC.G.cc                                                                                                       >  1:6606131/1‑101 (MQ=255)
 ggcgCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GGGGGC.G.cc                                                                                                       >  1:3241168/1‑101 (MQ=255)
      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:5813367/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:99470/99‑1 (MQ=255)
                                                                                                 gggggC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTg       >  1:7056432/3‑101 (MQ=255)
                                                                                                 gggggC.T.CCG..C.CGC.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTg       >  1:4143678/3‑101 (MQ=255)
                                                                                                  ggggC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt      >  1:5256579/2‑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)
                                                                                                   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:311260/101‑1 (MQ=255)
                                                                                                    ggC.T.CCGC.C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCg    <  2:6605791/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)
                                                                                                         |                                                                                                       
CGGCGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTCGGGGCAAG.CCCGAGGCGACCCG.CC..GCGGGC.T.CCG..C.C.C.CGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGC  >  NC_014034/1965177‑1965374

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