Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,955 C→G 1.6% P234P (CCG→CCC)  RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9550C→G 1.6% 2.6 0.0 ‑0.0373P234P (CCG→CCC) RCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (3/3):  ref base (198/168):  total (201/172)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.98e-01

Only 100 of 603 total aligned reads displayed.

TTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCG  >  NC_014034/221855‑222055
                                                                                                    |                                                                                                    
  tCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggtcgg                                                                                                    <  1:5548821/101‑1 (MQ=255)
      cGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcg                                                                                                >  2:4910320/1‑101 (MQ=255)
           cccGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggccgg                                                                                           >  2:8871877/1‑101 (MQ=255)
              gACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtc                                                                                        >  2:4551279/1‑101 (MQ=255)
              gACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtc                                                                                        <  1:2671083/101‑1 (MQ=255)
               acacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGATCATCGACGGGGGGCGTCGGATGCGTGGGCTCGGTCggatgggtcgggtcc                                                                                       >  2:6860260/1‑101 (MQ=255)
                 acCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgt                                                                                     >  2:1563949/1‑101 (MQ=255)
                    gTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtggg                                                                                  >  2:2351493/1‑101 (MQ=255)
                     tCGCCGCGGCTGCGGCGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGTTCGGTGggatccgttgggtccgtgggg                                                                                 >  2:382755/1‑100 (MQ=255)
                      cGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggat                                                                                <  2:3963544/101‑1 (MQ=255)
                      cGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggat                                                                                >  2:146536/1‑101 (MQ=255)
                      cGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggac                                                                                >  1:7043324/1‑100 (MQ=255)
                       gCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatc                                                                               >  1:212052/1‑101 (MQ=255)
                        ccGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcg                                                                              <  1:5145520/101‑1 (MQ=255)
                          gAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcggggccgtgggatcggt                                                                            >  2:6633572/1‑101 (MQ=255)
                           aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                                           >  1:5279297/1‑101 (MQ=255)
                             gCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgg                                                                         <  1:5963135/101‑1 (MQ=255)
                              cTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcggg                                                                        <  1:8516079/101‑1 (MQ=255)
                              cTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcggg                                                                        <  1:1183556/101‑1 (MQ=255)
                               ttCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggt                                                                       <  1:8339104/101‑1 (MQ=255)
                                 cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatgggtcgggtcc                                                                     >  1:575199/1‑101 (MQ=255)
                                 cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcgggcggggccgtgggatcggtcggggcc                                                                     >  2:1955226/1‑101 (MQ=255)
                                  ggTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccg                                                                    >  1:7988753/1‑101 (MQ=255)
                                  ggTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccg                                                                    >  1:4848402/1‑101 (MQ=255)
                                  ggTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccg                                                                    <  2:5359993/101‑1 (MQ=255)
                                  ggTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccg                                                                    <  1:3183506/101‑1 (MQ=255)
                                   gtgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcgggcgggtcccg                                                                   >  2:2199871/1‑99 (MQ=255)
                                         gtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCggtcgggtccgtgggatcggtcgggtccgtgggatc                                                             <  1:1673597/101‑1 (MQ=255)
                                          tggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTTGGATCGgtcgggtccgtgggatcggtcgggtccgtgggatcg                                                            <  1:729517/101‑1 (MQ=255)
                                           ggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGtcgggtccgtgggatcggtcgggtccgtgggatcgg                                                           <  1:768105/101‑1 (MQ=255)
                                           ggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGtcgggtccgtgggatcggtcgggtccgtgggatcgg                                                           <  1:1527402/101‑1 (MQ=255)
                                             tCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGGCGGATGCGTGGGCTCGGTCGGATCGGGCgggtccctgggatcggtcgggtccctgggatcggtc                                                         >  2:808754/1‑101 (MQ=38)
                                               aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgg                                                       <  2:4339297/101‑1 (MQ=255)
                                                  gTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                                                    <  1:540539/101‑1 (MQ=255)
                                                   tCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                                                   <  1:8646401/101‑1 (MQ=255)
                                                    cGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGATCCg                                                  <  2:1317708/101‑1 (MQ=255)
                                                       cGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                               >  1:4101454/1‑101 (MQ=255)
                                                         gATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGa                                             >  2:6138192/1‑101 (MQ=255)
                                                         gATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGa                                             <  2:7500611/101‑1 (MQ=255)
                                                           tCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtg                                           >  2:3404842/1‑100 (MQ=255)
                                                              cgcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                        <  2:3600413/101‑1 (MQ=255)
                                                                    tCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                  <  1:1857397/101‑1 (MQ=255)
                                                                     cGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                 >  2:1175358/1‑101 (MQ=255)
                                                                     cGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                 <  1:8162246/101‑1 (MQ=255)
                                                                      ggcGGTGGGCGTCGGATGCGTGGGGTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCg                                <  2:7077052/99‑1 (MQ=18)
                                                                       aCGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAGC.GTCggggtcggg                              >  2:2969506/1‑95 (MQ=21)
                                                                           tGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggt                           >  2:9187951/1‑100 (MQ=255)
                                                                             ggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                         <  2:2579393/101‑1 (MQ=255)
                                                                              gCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg                        <  2:8371105/101‑1 (MQ=255)
                                                                                gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                      >  1:4575618/1‑101 (MQ=255)
                                                                                gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggg                      >  2:3199881/1‑100 (MQ=255)
                                                                                 tCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                     <  1:370906/101‑1 (MQ=255)
                                                                                  cGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                    >  2:8550603/1‑101 (MQ=255)
                                                                                   ggATGCTTGGGCTCGGTCGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGTTCTGTCGGGTCCGTGGGAtcggtcag                   >  2:5959744/1‑99 (MQ=17)
                                                                                     gtGCGTGGGGTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt                 <  1:1346800/100‑1 (MQ=255)
                                                                                     aTGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt                 <  1:7844329/101‑1 (MQ=255)
                                                                                      tGCGTGGGCTCGGTCGGATCGGTCGGCTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                <  2:2752154/101‑1 (MQ=255)
                                                                                       tcGTGGGAGCGGTCGGGTCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc               <  2:3665732/100‑1 (MQ=11)
                                                                                       gCGTGGGATCGGTCGGATCGGTCGGGTCCGTCGGGTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc               <  1:4232956/101‑1 (MQ=255)
                                                                                       cggTGGGATCCGTCGGATCGGTGGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTcc               <  2:3123418/99‑1 (MQ=11)
                                                                                        ggtgccagcggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCg              <  1:2826597/93‑1 (MQ=255)
                                                                                        gggcggCTCCGTCCGAGCGCTCCGCTCCTTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCg              <  2:1743538/97‑1 (MQ=11)
                                                                                        ggTGGGCTCGGGGGGATCGGTCGGGGCCGTGGGAGCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCg              <  1:9050942/100‑1 (MQ=14)
                                                                                        cGTGGGGTCGGGCGGAGCGGCGGGGGCCGGGGGATCGGTGGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg              <  2:5579765/101‑1 (MQ=14)
                                                                                        cGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTTGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCg              >  1:965321/1‑101 (MQ=255)
                                                                                         gTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt             <  2:5877212/101‑1 (MQ=11)
                                                                                         gTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt             <  2:6704780/101‑1 (MQ=11)
                                                                                          ggggGTCGGTCGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg            <  2:7256576/100‑1 (MQ=255)
                                                                                          ggggCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg            <  2:8990829/100‑1 (MQ=255)
                                                                                           gggCTCCGTGGGTCCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCTGTCGGATCCGTTGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgg           <  1:6900158/101‑1 (MQ=11)
                                                                                            gggtcgggcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTcgg          <  1:1000386/98‑1 (MQ=11)
                                                                                              cTCGGTCGGTTCCGTGGGGTCGGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt        <  2:2111731/101‑1 (MQ=11)
                                                                                              cTCGGTCGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt        <  1:6259543/101‑1 (MQ=255)
                                                                                              atcggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAt        <  1:3101968/100‑1 (MQ=11)
                                                                                              atcggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTCGGAt        >  2:6313434/2‑101 (MQ=11)
                                                                                               tcggtgggAGCGGTCGGGTCCGTGGGGTCGGTCGGGGCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc       <  1:3885949/101‑1 (MQ=11)
                                                                                               tcgggcggGTCCGTCGGATCCGTGGGGTCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGCCGGGTCCGTTGGAtc       <  1:8848352/101‑1 (MQ=11)
                                                                                               tccgtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc       <  1:4973682/98‑1 (MQ=11)
                                                                                                cggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg      >  2:5425708/1‑101 (MQ=255)
                                                                                                cggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg      <  1:6779170/101‑1 (MQ=255)
                                                                                                cggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg      <  2:6545620/101‑1 (MQ=255)
                                                                                                cggtcggATCGGGCGGGTCCGTGGGATCCGTCGGGTCCGTGGGATCGGTTGGATCCGTGGGATCCGTCGGGGCCGTGGGATCGGTCGGATCCGCGGCAtcg      >  1:6518372/1‑101 (MQ=17)
                                                                                                ccgtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTGCGTGGGATCGGTCGGGTACGTGGGAtcg      >  2:7891360/3‑101 (MQ=11)
                                                                                                ccgtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg      <  1:4412978/99‑1 (MQ=11)
                                                                                                ccgtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg      >  1:4723852/3‑101 (MQ=11)
                                                                                                ccggcggATCGGTCGGGTCCGTGGGGGCCGTCGGCTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg      <  2:4748617/97‑1 (MQ=14)
                                                                                                 ggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcgg     >  2:4126883/1‑101 (MQ=255)
                                                                                                 ggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcgg     >  2:5098531/1‑101 (MQ=255)
                                                                                                 cgtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg     >  1:3625964/2‑101 (MQ=11)
                                                                                                 cggcggATCGGTGGGGGCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg     <  1:64569/98‑1 (MQ=11)
                                                                                                  gtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggt    <  2:7726360/101‑1 (MQ=11)
                                                                                                  gtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggg    >  2:4480485/1‑100 (MQ=11)
                                                                                                  gtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAGCGGTCGGGGCGGTGGGAgcggt    >  1:2210074/1‑101 (MQ=11)
                                                                                                  gtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt    >  1:5616286/1‑101 (MQ=255)
                                                                                                  gtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt    >  1:8270154/1‑101 (MQ=255)
                                                                                                  gtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTGGGAtcggg    >  2:8678673/1‑100 (MQ=14)
                                                                                                   tcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtg   >  2:8251737/1‑100 (MQ=255)
                                                                                                   gcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc   <  1:3404946/100‑1 (MQ=255)
                                                                                                    cggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg  <  1:6154336/101‑1 (MQ=255)
                                                                                                    cggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtcg  <  1:2391989/101‑1 (MQ=11)
                                                                                                    |                                                                                                    
TTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCG  >  NC_014034/221855‑222055

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG