Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,964 +A 14.7% coding (846/1578 nt) RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,9642.→A14.7% 170.5 ‑0.1 ‑0.0279coding (846/1578 nt)RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (19/22):  ref base (116/122):  total (135/144)
Fisher's exact test strand distribution p-value = 8.66e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.28e-01

Only 100 of 311 total aligned reads displayed.

AAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAA  >  NC_014034/1396864‑1397064
                                                                                                              |                                                                                                      
aaaCCGGTCCAGAGATAGACCGCGCGTACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTg                                                                                                         <  1:2547354/101‑1 (MQ=255)
 acccGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..c                                                                                                      <  2:3899835/99‑1 (MQ=255)
     ggTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CCG..CGTg.g                                                                                                 >  1:2891896/1‑101 (MQ=255)
       tCCAGAGATAGACAGCGCGGACCTCCTCAGACGAGGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTc                                                                                               <  1:3634689/101‑1 (MQ=255)
         cAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GGCGc                                                                                             >  2:4038710/1‑101 (MQ=255)
              aTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTcccc                                                                                        <  2:3502560/101‑1 (MQ=255)
              aTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GGCGCTcccc                                                                                        >  2:3204558/1‑101 (MQ=255)
               tAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTccccc                                                                                       >  2:4131773/1‑101 (MQ=255)
               tAGACCGCGCGGACCTCCTCAGACGCCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTccccc                                                                                       <  2:2766031/101‑1 (MQ=255)
                  aCCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCcgc                                                                                    >  1:4084411/1‑101 (MQ=255)
                         ggACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCg                                                                             <  1:1080464/101‑1 (MQ=255)
                         ggACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCg                                                                             >  2:950168/1‑101 (MQ=255)
                          gccctcctCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGt                                                                            <  1:1601367/99‑1 (MQ=255)
                            ctcgcggttgcGTCGCTCCCCCCGCTGGGGCGTGGATCGAAACATGG.CTAT.GGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTg                                                                           <  2:1854079/91‑1 (MQ=0)
                             ggtcaacactgTCGCTCCACCCGCGGGGGCGTGGATCGAAAC..GC.TG.G.CGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTgg                                                                          <  1:836338/91‑1 (MQ=9)
                             ctcctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                         <  2:281490/101‑1 (MQ=255)
                             attggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                         >  1:2052845/11‑101 (MQ=0)
                             attggccccggTCGCCCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                         <  1:1786308/91‑1 (MQ=0)
                              tcctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                        >  1:1018959/1‑101 (MQ=255)
                               tcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GC.TG.G.CGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                        >  2:3599946/9‑101 (MQ=9)
                               gcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGG.CTAT.GGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                        <  2:1243074/94‑1 (MQ=2)
                               gcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGG.CTAT.GGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                        >  1:1291637/8‑101 (MQ=2)
                               ccccgcccGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACT.GC.GT.G.CGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTgg                                                                          <  1:268214/94‑1 (MQ=0)
                                 tcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                     >  1:3723826/1‑101 (MQ=255)
                                 ccgcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GC.GT.G.CGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                        <  2:1557810/96‑1 (MQ=2)
                                 aagactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GC.TG.G.CGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                      >  1:2420254/7‑101 (MQ=9)
                                  gttgcggcGCTCCCCCCGCGGGGGCGTGGATCGAAACATGG.CTAT.GGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                     >  1:1921921/8‑101 (MQ=2)
                                  cgcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GC.GT.G.CGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                       <  2:764498/97‑1 (MQ=2)
                                  ccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                    <  2:2439711/96‑1 (MQ=2)
                                  ccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                    <  2:2243330/96‑1 (MQ=2)
                                  cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                    >  1:3035506/1‑101 (MQ=255)
                                   gaccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                   >  2:1579462/4‑101 (MQ=255)
                                   cccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                   <  2:2648408/97‑1 (MQ=2)
                                   aGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                   >  2:1681531/1‑101 (MQ=255)
                                     gcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGG.CTAT.GGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                  <  1:3900838/100‑1 (MQ=2)
                                     cgggcGCCCCCCCCGCGGGGGCGTGGATCGNAACA..A.CT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                                 <  1:2799515/97‑2 (MQ=2)
                                     cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                                 <  2:2016190/99‑2 (MQ=2)
                                     cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                                 >  1:2148916/3‑100 (MQ=2)
                                      ggTCGCTCCCCCCGCGGCGGCGTGGATCGAAACA..A.CT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActa                                                                >  1:549536/2‑99 (MQ=2)
                                      cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                                                >  1:1358714/1‑101 (MQ=255)
                                       gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActac                                                               <  2:2528041/101‑4 (MQ=2)
                                        tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGG.CTAT.GGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccg                                                               >  2:437598/1‑98 (MQ=2)
                                        tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActaca                                                              <  1:7030/101‑5 (MQ=2)
                                        tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActaca                                                              >  1:500958/1‑97 (MQ=2)
                                        tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActaca                                                              <  1:2235987/101‑5 (MQ=2)
                                         cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GC.GT.G.CGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACac                                                                <  2:96623/101‑2 (MQ=2)
                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACATGG.CTAT.GGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGGGGATCGAAAcccgcc                                                             >  2:2928340/1‑96 (MQ=2)
                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGc                                                            <  1:417249/101‑1 (MQ=255)
                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGGTCGAAACATCGGc                                                            >  2:3566385/1‑101 (MQ=255)
                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                                            >  2:1309375/1‑95 (MQ=2)
                                          gCTCCCCCCGCGGGGGCGTGGATCGAAAC..GC.TG.G.CGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggc                                                             >  2:3038203/1‑96 (MQ=9)
                                           cTCCCCCCGCGGGGGCGTGGATCGAAACATGG.CTAT.GGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgcca                                                            <  2:3087523/101‑7 (MQ=2)
                                           cTCCCCCCGCGGGGGCGTGGATCGAAACATGG.CTAT.GGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgcca                                                            >  1:466959/1‑95 (MQ=2)
                                           cTCCCCCCGCGGGGGCGTGGATCGAAACA.CCAAT.GCCGAA.GC.CAAG.GCCGCA.GCAGC.GGC..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGAACGAAACaaaaa                                                             >  1:1778079/1‑97 (MQ=2)
                                           cTCCCCCCGCGGGGGCGTGGATCGAAACA..A.CT.G.ACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccg                                                           <  2:1634861/101‑8 (MQ=2)
                                           cTCCCCCCGCGGGGGCGTGGATCGAAAC..GC.TG.G.CGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcc                                                            >  2:501510/1‑95 (MQ=9)
                                            tCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtc                                                          >  1:1045238/1‑101 (MQ=255)
                                            tCCCCCCGCGGGGGCGTGGATCGAAAC..GC.TG.G.CGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcct                                                           >  1:3920103/1‑94 (MQ=9)
                                            tCCCCCCGCGGGGGCGTGGATCGAAAC..GC.TG.G.CGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcct                                                           >  1:3196182/1‑94 (MQ=9)
                                            cccccccGCGGGGGCGTGGATCGAAACA.GC.CT.G.CACG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtc                                                          <  1:2167685/100‑1 (MQ=255)
                                             ccccccGCGGGGGCGTGGATCGAAACT.GC.GT.G.CGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaat                                                            >  1:1720171/1‑99 (MQ=2)
                                             ccccccGCGGGGGCGTGGATCGAAAC..GC.TG.G.CGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctc                                                          <  1:1551679/101‑9 (MQ=9)
                                              cccccGCGGGGGCGTGGATCGAAACT.GC.GT.G.CGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatg                                                           <  2:1183173/101‑4 (MQ=2)
                                              cccccGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcgt                                                        <  2:3391650/101‑1 (MQ=255)
                                                cccGCGGGGGCGTGGATCGAAACT.GC.GT.G.CGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgcc                                                         <  1:422083/101‑6 (MQ=0)
                                                  cgcgGGGGCGTGGATCGAAACT.GC.GT.G.CGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgccga                                                       >  1:910829/1‑94 (MQ=0)
                                                       gggCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCa                                               <  2:3165396/101‑1 (MQ=255)
                                                         gCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGc                                             <  2:1797123/101‑1 (MQ=255)
                                                            tGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCgg                                          <  1:3361384/101‑1 (MQ=255)
                                                             ggATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGa                                         >  1:641595/1‑101 (MQ=255)
                                                             ggATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGa                                         >  2:2612140/1‑101 (MQ=255)
                                                              gATCGAAACC.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGAc                                        <  1:4132035/101‑1 (MQ=255)
                                                              gATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGAc                                        <  2:2067419/101‑1 (MQ=255)
                                                              gATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGAc                                        <  2:3314361/101‑1 (MQ=255)
                                                               aTCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACg                                       <  2:1077494/101‑1 (MQ=255)
                                                               aTCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACg                                       <  2:2586443/101‑1 (MQ=255)
                                                                 cGAAGCA.GC.CG.G.CGCG.GC.CGCA.TCTGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGAt                                     <  2:2401862/101‑1 (MQ=255)
                                                                     aCA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTcc                                 <  1:3699062/101‑1 (MQ=255)
                                                                      cA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGGCGATGTCCg                                >  2:381142/1‑101 (MQ=255)
                                                                      cA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCg                                <  1:3945236/101‑1 (MQ=255)
                                                                         gC.CT.G.CGCG.GC.CGCA.TCGCCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTc                              <  2:1934989/101‑1 (MQ=255)
                                                                          c.cT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCg                             <  2:641518/101‑1 (MQ=255)
                                                                             t.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCt                           <  2:1108690/101‑1 (MQ=255)
                                                                             t.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCt                           <  2:3257704/101‑1 (MQ=255)
                                                                                  gcg.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTccc                        <  2:1358614/101‑1 (MQ=255)
                                                                                      gC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcccccc                     <  1:1410445/101‑1 (MQ=255)
                                                                                          gCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGGAACATAAGCTCGTCGTTGGTCAGCAGGACGATGTCCGTCGCTCCCCCcgcg                  >  1:3555133/1‑101 (MQ=255)
                                                                                               cGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCCCCCCCCGCggggg              >  2:334141/1‑101 (MQ=255)
                                                                                                ggCT.TCCGC.CTG..CGTG.GTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGc             <  2:3714353/101‑1 (MQ=255)
                                                                                                  cT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCgg           >  2:544704/1‑100 (MQ=255)
                                                                                                     tCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTgg         <  1:1959975/101‑1 (MQ=255)
                                                                                                       cGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGAt       <  1:1222848/101‑1 (MQ=255)
                                                                                                        gC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATc      <  1:1126987/101‑1 (MQ=255)
                                                                                                           cTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGa    <  2:1038798/101‑1 (MQ=255)
                                                                                                            tG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa   <  2:4180519/101‑1 (MQ=255)
                                                                                                            tG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa   <  2:3450813/101‑1 (MQ=255)
                                                                                                            tG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa   <  1:2830518/101‑1 (MQ=255)
                                                                                                            tG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa   <  2:801482/101‑1 (MQ=255)
                                                                                                            tG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa   <  1:3472301/101‑1 (MQ=255)
                                                                                                             g..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa  >  1:3544371/1‑101 (MQ=255)
                                                                                                              |                                                                                                      
AAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GC.CT.G.CGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAA  >  NC_014034/1396864‑1397064

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG