Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,493,952 +G 97.8% coding (252/252 nt) RCAP_rcc02319 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,493,9521.→G97.8% 3.9 ‑0.6 ‑1.8183coding (252/252 nt)RCAP_rcc02319hypothetical protein
Reads supporting (aligned to +/- strand):  new base (50/129):  ref base (0/4):  total (50/133)
Fisher's exact test strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.61e-02

Only 100 of 226 total aligned reads displayed.

TGAGCCCCGAAGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AA.GGGGGG.G.CGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCG  >  NC_014034/2493852‑2494047
                                                                                                      |                                                                                                  
tGAGCCCCGAAGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.aa                                                                                                     <  2:3253500/101‑1 (MQ=255)
tGAGCCCCGAAGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.aa                                                                                                     <  2:3429073/101‑1 (MQ=255)
 gAGCCCCGAAGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AA.g                                                                                                   <  2:3339724/101‑1 (MQ=255)
 gAGCCCCGAAGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AA.g                                                                                                   >  1:1919621/1‑101 (MQ=255)
     cccGAAGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AA.ggggg                                                                                               <  1:379028/101‑1 (MQ=255)
      ccGAAGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AA.gggggg                                                                                              <  1:142699/101‑1 (MQ=255)
      ccGAAGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AA.gggggg                                                                                              >  2:1850893/1‑101 (MQ=255)
        gAAGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AA.gggggg.gg                                                                                           >  2:5555168/1‑100 (MQ=255)
         aaGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AA.gggggg.ggg                                                                                          >  1:5335410/1‑99 (MQ=255)
         aaGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AA.gggggg.ggc                                                                                          <  1:5071306/101‑3 (MQ=255)
          agagCGCCGCGCGGTGATGGCCGAGCTGCCCCCCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AA.gggggg.ggcg                                                                                         <  2:2569433/101‑4 (MQ=255)
          agagCGCCGCGCGGTGATGGCCGAGCTGCCCCCCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AA.gggggg.ggcg                                                                                         <  1:2289197/101‑4 (MQ=255)
          agagCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AA.gggggg.ggcg                                                                                         >  2:5758294/1‑98 (MQ=255)
            agCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGGGG.AA.gggggg.gggggg                                                                                       >  2:4922100/1‑96 (MQ=255)
            agCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AA.gggggg.ggcgct                                                                                       <  2:3454417/101‑6 (MQ=255)
             gcgcCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCTCGGCCCGACAGGCCGAGG.AA.aggggg.ggggccg                                                                                      >  1:4936426/1‑95 (MQ=255)
                ccgcGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGccc                                                                                   <  1:4757816/101‑1 (MQ=255)
                 cgcgcgGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGcccc                                                                                  <  1:634624/101‑1 (MQ=255)
                   cgcgGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGcccccc                                                                                <  1:4800813/101‑2 (MQ=255)
                    gcgGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.GGCTGccccccg                                                                               >  2:2298434/1‑99 (MQ=255)
                      ggTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGccccccgtc                                                                             <  2:2649277/101‑5 (MQ=255)
                      ggTGATGGCCGAGCTGCCCACCCGCGCGCAGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGccccccgtc                                                                             <  1:3538395/101‑5 (MQ=255)
                       gTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGccccccgtcc                                                                            <  1:2805853/101‑6 (MQ=255)
                        tGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCg                                                                           <  1:151217/101‑1 (MQ=39)
                        tGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGG.AAGGGGGGG.G.GGCTGccccccgcgcc                                                                           >  2:155700/1‑95 (MQ=255)
                          aTGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCa                                                                         <  1:3104331/101‑1 (MQ=39)
                           tGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCaa                                                                        <  2:3592959/101‑1 (MQ=38)
                            ggCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAg                                                                       <  1:1760994/101‑1 (MQ=38)
                             gCCGAGCTGCCCCCCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGc                                                                      <  2:4012335/101‑1 (MQ=38)
                               cGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGC.AAGGGGGGG.G.AGCTGCCCCCGCTCCGCAAGccg                                                                    >  1:5704098/1‑99 (MQ=255)
                                  gCTGCCCCCCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACt                                                                 <  2:4847232/101‑1 (MQ=37)
                                  gCTGCCCCCCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACt                                                                 <  1:4262109/101‑1 (MQ=37)
                                   cTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGAACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTc                                                                <  1:2134887/101‑1 (MQ=37)
                                    tGCCCCCCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTcc                                                               <  2:587110/101‑1 (MQ=37)
                                    tGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGGAGCCCCCCGTCCGCAACCGGACCcc                                                               >  2:1725876/1‑101 (MQ=255)
                                    tGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTcc                                                               <  1:1581136/101‑1 (MQ=37)
                                       ccccccGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTccccc                                                            <  1:3010057/98‑1 (MQ=25)
                                       ccACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGGCCCCCGGCCCCAAGCGGGCTccccc                                                            >  2:1763812/1‑101 (MQ=37)
                                        cACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTcccccc                                                           <  1:3665497/101‑1 (MQ=37)
                                        cACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTcccccc                                                           >  2:4430323/1‑101 (MQ=37)
                                        cACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTcccccc                                                           >  1:5246225/1‑101 (MQ=37)
                                        cACCCGCCCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTcccccc                                                           <  2:2145222/101‑1 (MQ=25)
                                           ccGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCggg                                                        >  1:3910379/1‑101 (MQ=25)
                                            cgcgcgCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGa                                                       <  2:1585252/101‑1 (MQ=25)
                                              cgcgcCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGata                                                     <  2:1020821/101‑1 (MQ=25)
                                              cgcgcCGCGTCTGACCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGCCCGCCACCGGACTCCCCCCAGGatt                                                     >  1:5727190/1‑100 (MQ=25)
                                               gcgcCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGatat                                                    <  1:753929/101‑1 (MQ=25)
                                                  ccgcGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTg                                                 <  2:3028916/101‑1 (MQ=25)
                                                   cgcgTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgc                                                >  2:3152893/1‑101 (MQ=25)
                                                    gcgTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgca                                               <  2:2392124/101‑1 (MQ=25)
                                                     ggTCTGGCCCAGCGCGGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcag                                              <  2:4669248/100‑1 (MQ=21)
                                                       tCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcagca                                            <  1:1097514/101‑1 (MQ=25)
                                                         tGGCCCGGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAg                                          <  1:4722034/101‑1 (MQ=21)
                                                         tGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAg                                          <  2:716772/101‑1 (MQ=21)
                                                          ggCCCCGCGCCGCGGCGGGCGCGCGGCACGACCGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGa                                         <  1:792750/101‑1 (MQ=18)
                                                          ggCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGa                                         <  1:2897694/101‑1 (MQ=21)
                                                           gCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGAt                                        <  2:4055277/101‑1 (MQ=21)
                                                            cccCGCGCCGCGGCGGGCGGGCGGCAGGACAGGCCGAGT.AA.GGGGGGCG.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATg                                       <  1:2682086/101‑1 (MQ=18)
                                                             ccAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGa                                      <  2:4612123/101‑1 (MQ=21)
                                                             ccAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGa                                      <  2:2800768/101‑1 (MQ=21)
                                                               aGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaaa                                    >  1:543753/1‑101 (MQ=21)
                                                                 cgccgcGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                  <  1:553261/101‑1 (MQ=21)
                                                                 cgccgcGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                  <  2:4360332/101‑1 (MQ=21)
                                                                 cgccgcGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                  <  1:2709425/101‑1 (MQ=21)
                                                                  gccgcGGCGGGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCg                                 >  2:5668893/1‑101 (MQ=21)
                                                                     gcggcgGGCGCGCGGCACGCCAGGCCGAGTCAG.GGGGGG.G.CGCGGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.At                             <  1:369795/101‑1 (MQ=18)
                                                                      gggggggcgCCGGGCCCGACAGGCCGAGT.AA.GGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGc                           <  1:2269614/97‑1 (MQ=18)
                                                                      cggcggGCGCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATg                            <  2:3747437/101‑1 (MQ=21)
                                                                       ggcgggcgCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGc                           <  1:1229830/101‑1 (MQ=21)
                                                                          gggcgCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCtt                        <  2:5331538/101‑1 (MQ=21)
                                                                          gggcgCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCtt                        <  1:195317/101‑1 (MQ=21)
                                                                          gggcgCGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCtt                        <  2:4272241/101‑1 (MQ=21)
                                                                           gggggggggCCCGCCAGGCCGAGT.AAGGGGGGG.G.CGCTCCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCttt                       <  1:3210721/94‑1 (MQ=17)
                                                                           ggcgTGCGGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGGTATTTGCAGCAAGATGAAACCGG.ATGCCttt                       >  1:5239086/1‑101 (MQ=18)
                                                                              gcgcgGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCAt                    <  2:5376790/101‑1 (MQ=21)
                                                                                 cgGCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCtt                 <  2:4126816/101‑1 (MQ=21)
                                                                                  ggCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTg                <  1:4043740/101‑1 (MQ=21)
                                                                                  ggCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTg                >  2:2517607/1‑101 (MQ=21)
                                                                                  ggCACGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTg                <  1:166654/101‑1 (MQ=21)
                                                                                    cgcGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGcc              <  1:334478/99‑1 (MQ=18)
                                                                                     aCGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGccc             <  2:4743496/101‑1 (MQ=21)
                                                                                     aCGACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGGCCGCCAGCGGACTCCCCCCCGGATATTTGCCGCAAGAAGAAAACGGAAGGCCTTTCTTCtgccc              >  1:1041465/1‑97 (MQ=9)
                                                                                       gACAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCCG.ATGCCTTTCATCTTGgccca           >  1:2620391/1‑99 (MQ=18)
                                                                                        aCAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGaa          <  2:5258250/101‑1 (MQ=21)
                                                                                        aCAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGaa          <  1:1126193/101‑1 (MQ=21)
                                                                                         cAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGaaa         <  1:2331066/101‑1 (MQ=21)
                                                                                         cAGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGaaa         >  2:1955983/1‑101 (MQ=21)
                                                                                          aGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAat        <  2:4539795/101‑1 (MQ=21)
                                                                                          aGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAat        <  2:4629544/101‑1 (MQ=21)
                                                                                          aGGCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAat        <  2:3911735/101‑1 (MQ=21)
                                                                                           ggCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAata       >  1:4318997/1‑101 (MQ=21)
                                                                                            gCCGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAatat      <  1:3596599/101‑1 (MQ=21)
                                                                                            gCCGAGT.AA.GGGGGG.GCCGCTGCCCCCCGTCCGCAAGCGGCCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAatat      <  2:5644917/101‑1 (MQ=21)
                                                                                             ccGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATc     <  2:5480580/101‑1 (MQ=21)
                                                                                              cGAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATcc    <  1:1768512/101‑1 (MQ=21)
                                                                                               gAGT.AAGGGGGGG.G.CGGTGCCCCCCGTCCGCAAGCGGTCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATccc   <  2:4730446/101‑1 (MQ=18)
                                                                                               gAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATccc   >  1:5668494/1‑101 (MQ=18)
                                                                                               gAGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATccc   >  2:5357442/1‑101 (MQ=18)
                                                                                                aGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCCAGCGGACTCCCACCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCg  <  2:2268504/101‑1 (MQ=18)
                                                                                                aGT.AAGGGGGGG.G.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCg  <  2:1405575/101‑1 (MQ=18)
                                                                                                      |                                                                                                  
TGAGCCCCGAAGAGCGCCGCGCGGTGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGT.AA.GGGGGG.G.CGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCG  >  NC_014034/2493852‑2494047

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 32 ≤  ATCG  < 40 ≤  ATCG