Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,256,189 +C 93.3% coding (1332/1734 nt) recJ → single‑stranded‑DNA‑specific exonuclease RecJ

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,256,1891.→C93.3% 15.3 ‑0.8 ‑0.1104coding (1332/1734 nt)recJsingle‑stranded‑DNA‑specific exonuclease RecJ
Reads supporting (aligned to +/- strand):  new base (53/44):  ref base (5/2):  total (58/46)
Fisher's exact test strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.70e-01

Only 100 of 126 total aligned reads displayed.

GGTGCAGCGCGTCGCGGCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCT.GCC.GAGG.GCAG.TGGAA.CCGCGATGG.AA.GGCT.GCCGAG.TGTTCGCCCGGCAGGGCG.C.GG..CAC.GGCGGCCCGGCCGATCTGCGGCTGGACGGGCTCCTGATGCCGGGCGCCGCGAC  >  NC_014034/3256093‑3256278
                                                                                                      |                                                                                                 
ggTGCAGCGCGTCGCGGCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATgg                                                                                                    >  2:2848773/1‑101 (MQ=255)
 gTGCAGCGCGTCGCGGCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATggc                                                                                                   <  1:2365595/101‑2 (MQ=255)
 gTGCAGCGCGTCGCGGCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATggc                                                                                                   <  1:1901283/101‑2 (MQ=255)
 gTGCAGCGCGTCGCGGCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATggc                                                                                                   <  1:1415792/101‑2 (MQ=255)
   gCAGCGCGTCGCGGGCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATggcac                                                                                                 <  2:448859/101‑4 (MQ=255)
   gCAGCGCGTCGCGGCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATggcac                                                                                                 >  2:2503809/1‑98 (MQ=255)
   gCAGCGCGTCGCGGCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATggcac                                                                                                 <  1:1666967/101‑4 (MQ=255)
     aGCGCGTCGCGGCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATggcacgg                                                                                               >  1:3818351/1‑96 (MQ=255)
     aGCGCGTCGCGGCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATggcacgg                                                                                               <  2:2639446/101‑6 (MQ=255)
       ggcgTCGCGGCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATggcacggct                                                                                             <  1:2975194/100‑8 (MQ=255)
       cgcgTCGCGGCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATggcacggct                                                                                             <  2:522620/101‑8 (MQ=255)
        gcgTCGCGGCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATggcacggctg                                                                                            <  2:3021663/101‑9 (MQ=255)
         cgTCGCGGCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATggcacggctgg                                                                                           >  1:3365151/1‑92 (MQ=255)
         cgTCGCGGCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATggcacggctgg                                                                                           >  1:3256455/1‑92 (MQ=255)
         cgTCGCGGCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATggcacggctgg                                                                                           <  1:2656973/101‑10 (MQ=255)
          gTCGCGGCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATggcacggctggc                                                                                          >  2:1588843/1‑91 (MQ=255)
          gTCGCGGCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATggcacggctggc                                                                                          >  1:1111058/1‑91 (MQ=255)
               ggCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGGC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAg.t                                                                                   <  1:1118370/101‑1 (MQ=255)
               ggCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAg.t                                                                                   >  1:3789868/1‑101 (MQ=255)
               ggCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAg.t                                                                                   >  2:10517/1‑101 (MQ=255)
                gCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAg.tt                                                                                  >  1:1343364/1‑100 (MQ=255)
                gCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAg.tt                                                                                  >  1:1238754/1‑100 (MQ=255)
                 ccGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAg.ttg                                                                                 >  1:6896/1‑99 (MQ=255)
                    aGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAg.ttgctc                                                                              >  2:3557312/1‑96 (MQ=255)
                    aGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAg.ttgctc                                                                              >  2:2781199/1‑96 (MQ=255)
                       ggCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAg.ttgctcgcc                                                                           >  2:3481095/1‑93 (MQ=255)
                         cTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCg                                                                          >  1:4071809/1‑101 (MQ=255)
                         cTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCg                                                                          <  2:2849962/101‑1 (MQ=255)
                         cTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCg                                                                          >  2:1062978/1‑101 (MQ=255)
                         cTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCg                                                                          <  2:1593111/101‑1 (MQ=255)
                         cTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCg                                                                          <  1:1544743/101‑1 (MQ=255)
                              gCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAgg                                                                     >  1:54871/1‑101 (MQ=255)
                              gCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAgg                                                                     >  2:886788/1‑101 (MQ=255)
                              gCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCCGCAgg                                                                     >  2:2886541/1‑101 (MQ=255)
                               cTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAggg                                                                    >  1:3155912/1‑101 (MQ=255)
                               cTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAggg                                                                    <  2:2578997/101‑1 (MQ=255)
                                 cAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAggggg                                                                  >  1:1729782/1‑99 (MQ=255)
                                 cAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAggggg                                                                  >  2:3984104/1‑99 (MQ=255)
                                 cAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGgcg                                                                  <  2:3256214/101‑1 (MQ=255)
                                 cAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCCGCAggggc                                                                  >  2:3239811/1‑99 (MQ=255)
                                 cAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCCCCCGGCAGGgcg                                                                  >  2:2242729/1‑101 (MQ=255)
                                 cAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGG.CACGGCTGGCCGAGTTGCTCGCCCGGCAGGgcg                                                                  <  2:3364906/101‑1 (MQ=255)
                                   aGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGG.CACGGCTGGCCGAGTTGCTCGCCCGGCAGGGcg.c.g                                                              <  1:612609/101‑1 (MQ=255)
                                    gggCGGCGGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.C.gg                                                             <  2:3236532/101‑1 (MQ=255)
                                      gcggcgGGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.C.gggc                                                           <  1:2315561/101‑3 (MQ=255)
                                        ggcggGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.C.gggcac                                                         <  1:354214/101‑5 (MQ=255)
                                         gcggGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGGG.C.ggggacg                                                        >  2:1642305/1‑96 (MQ=255)
                                         gcggGCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.C.gggcacg                                                        <  2:873799/101‑6 (MQ=255)
                                           gggCCCCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CAc.ggg                                                    <  1:3681722/101‑2 (MQ=255)
                                            ggCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CAc.gggc                                                   <  1:4126689/101‑3 (MQ=255)
                                            ggCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CAc.gggc                                                   <  2:1806084/101‑3 (MQ=255)
                                             gCACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGGG.C.ggggacggggg                                                    >  2:3437015/1‑92 (MQ=255)
                                              cACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGGG.CGGG..CAc.gggcgg                                                 >  2:1809786/1‑97 (MQ=255)
                                              cACCGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCCGCAGGGCC.CGGG..CAc.gggcgc                                                 >  1:2839833/1‑97 (MQ=255)
                                                ccGCATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGcggcc                                                <  2:2183027/101‑1 (MQ=255)
                                                   cATGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGcggcccgg                                             <  1:3147540/101‑1 (MQ=255)
                                                    gtGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGcggcccggc                                            <  1:3972471/100‑1 (MQ=255)
                                                    aTGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGcggcccggc                                            <  1:1588952/101‑1 (MQ=255)
                                                    aTGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGcggcccggc                                            >  1:1649584/1‑101 (MQ=255)
                                                     tGGCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGcgggccggcc                                           >  2:540962/1‑101 (MQ=255)
                                                       gCCGCGGGTCTGACCCTGGCC.GAGGCGCAGCGGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCGCC.AG..CACGGGCGGCCCGGCCGa                                         >  1:2394304/1‑101 (MQ=255)
                                                         cgcgGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGcggcccggcccgtc                                       >  1:3735473/1‑97 (MQ=255)
                                                         cgcgGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATc                                       <  1:256926/101‑1 (MQ=255)
                                                          gcgGGTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGGGGCCCGGGCGATcg                                      >  2:2744987/1‑100 (MQ=255)
                                                            gggTCTGACCCT.GGCAGAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..GACGGGCGGCCCGGCCGGTCttc                                    >  2:4043323/1‑99 (MQ=255)
                                                              gTCTGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCgg                                  <  1:1901779/101‑1 (MQ=255)
                                                                 tGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGGGGGCCGGCCGGTCTGCGGGTg                               >  2:2259703/1‑101 (MQ=255)
                                                                 tGACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTg                               >  1:3926522/1‑101 (MQ=255)
                                                                  gcccctggCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTgg                              <  1:876833/94‑1 (MQ=255)
                                                                  gACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGTCGATCGGCGGCGgg                              >  1:3656888/1‑101 (MQ=255)
                                                                  gACCCTGGCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCCGCCCGGCCGATCTGCGGCTgg                              >  1:2649460/1‑101 (MQ=255)
                                                                    ccctggCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGAc                            <  1:189797/96‑1 (MQ=255)
                                                                    ccctggCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGAc                            <  2:1238658/96‑1 (MQ=255)
                                                                     ccgggCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGACg                           <  2:4032150/97‑1 (MQ=255)
                                                                      ctggCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..GACGGGCGGCCCGGCCGGTCTGCGGCTgggagg                          >  1:3007456/4‑97 (MQ=255)
                                                                       tggCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGGTCTGCGGCTGGGCggg                         >  1:4073709/3‑101 (MQ=255)
                                                                       tggCC.GAGGCGCAGCTGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.C.GGGCCCC.GGCGGCCCGGCCGGTTCGCGGCttgcgtgg                         >  2:1287556/3‑94 (MQ=255)
                                                                             gaggcgcagctGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGACGGGCTCCt                    <  1:2314555/91‑1 (MQ=255)
                                                                               ggcgcagctGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGGGGCCCGGCCGGTCTGCGGCTGGGCGGGGGCCTgg                  >  2:1037200/9‑100 (MQ=255)
                                                                                gcgcagctGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGACGGGCTCCTGAt                 <  2:2735074/94‑1 (MQ=255)
                                                                                 cgcagctGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGGTCTGCGGCTGGGCGGGGGCCCGATg                >  2:1105499/7‑101 (MQ=255)
                                                                                  gcagctGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGCTCTGCGGCTGGACGGGCTCCTCATGc               >  2:650115/6‑101 (MQ=255)
                                                                                  gcagctGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGACGGGCTCCTGATGc               <  2:2338938/96‑1 (MQ=255)
                                                                                   cagctGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGACGGGCTCCTGATGcc              <  1:1096939/97‑1 (MQ=255)
                                                                                   cagctGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGAAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGACGGGCTCCTGATGcc              <  1:3557565/97‑1 (MQ=255)
                                                                                    agctGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGACGGGCTCCTGATGCCg             <  2:1649469/98‑1 (MQ=255)
                                                                                     gctGGAACCCGCGATGG.CACGGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGACGGGCTCCCGATGcccg            >  1:3850544/3‑99 (MQ=255)
                                                                                      ctGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CCCGGGCGGCCCGGGCGATCTGCGGCTGGGCGGGCTCCTGATGGCggg           >  1:201139/2‑101 (MQ=255)
                                                                                      ctGGAACCCGCGATGG.CACGGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCCGCGGCTGGACGGGCTCCCGATGCCggg           >  2:2663229/2‑101 (MQ=255)
                                                                                      ctGGAACCCGCGATGG.CACGGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCCGCCGATCTGCGGCTGGACCGGCTCCCGCTGccacg           >  2:1869062/2‑98 (MQ=255)
                                                                                       tGGAACCCGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGACGGGCTCCTGATGCCGGgc          <  2:3850285/101‑1 (MQ=255)
                                                                                        ggaacccGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGACGGGGTCCTGATTCCGGgcg         >  2:1718188/6‑101 (MQ=255)
                                                                                          aacccGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGACGGGCTCCTGATGCCGGGcgcc       <  1:500540/98‑1 (MQ=255)
                                                                                          aacccGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGACGGGCTCCTGATGCCGGGcgcc       <  2:811265/98‑1 (MQ=255)
                                                                                          aacccGCGATGG.CACGGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGACGGGCTCCTGATGCCGGGcgcc       >  2:3424115/4‑101 (MQ=255)
                                                                                           acccGCGATGG.CACGGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGACGGGCCCCTGATGCCGGGggccg      >  1:1146660/3‑101 (MQ=255)
                                                                                             ccGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGGCGGGCTCCTGATGCCCGGGGGcgcg    >  1:4119429/1‑101 (MQ=255)
                                                                                             ccGCGATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGACGGGCTCCCGATGCCGGGGGCcgcg    >  2:1364834/1‑101 (MQ=255)
                                                                                              cgcgATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGGCGGGGTCCTGCTGCCCGGCGCCGGGa   >  2:1027332/1‑101 (MQ=255)
                                                                                               gcgATGGCAC.GGCTGGCCGAGTTGCTCGCCCGGCAGGGCG.CGGG..CACGGGCGGCCCGGCCGATCTGCGGCTGGACGGGCTCCTGATGCCGGGCGCCGCGAc  <  1:2910544/101‑1 (MQ=255)
                                                                                                      |                                                                                                 
GGTGCAGCGCGTCGCGGCCGAGGGGCTTTTGCTCAAGGGCGGCGGGCACCGCATGGCCGCGGGTCTGACCCT.GCC.GAGG.GCAG.TGGAA.CCGCGATGG.AA.GGCT.GCCGAG.TGTTCGCCCGGCAGGGCG.C.GG..CAC.GGCGGCCCGGCCGATCTGCGGCTGGACGGGCTCCTGATGCCGGGCGCCGCGAC  >  NC_014034/3256093‑3256278

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG