Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,036 C→G 4.2% P207P (CCG→CCC RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0360C→G 4.2% 106.4 ‑0.3 ‑0.01534P207P (CCG→CCCRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (31/37):  ref base (631/833):  total (663/871)
Fisher's exact test strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.71e-01

Only 100 of 1953 total aligned reads displayed.

TCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGG  >  NC_014034/221936‑222136
                                                                                                    |                                                                                                    
    atcggtcgggtccggtGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGgatcggtcgggt                                                                                                  >  1:3853090‑M1/17‑99 (MQ=255)
      cggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATcc                                                                                                >  1:3979980‑M1/15‑101 (MQ=255)
      cggtcggatccgtgGGATCGGTCGGGTCCGTGGGCTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAtcg                                                                                                <  1:3437475‑M1/87‑2 (MQ=255)
         tGGGTTAGGTCGGGTCTGTGGGGGAGGTGGGGTCGGGCGGGTCGGTGGGATCGGTCGGGTCCGTAGGGTCGGTCGGGTCCGTCGGGTCCGTGGGATCGGtc                                                                                             <  1:1353033/101‑2 (MQ=12)
          gggATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGCCGGGTCCGtcg                                                                                            >  2:506380/1‑99 (MQ=1)
            ggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGa                                                                                          >  2:3203830‑M1/9‑101 (MQ=255)
            gCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTcggg                                                                                          <  2:821425/101‑2 (MQ=11)
             gtccgggggatcggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAt                                                                                         <  1:2118693/94‑1 (MQ=11)
              ttcgccttatcggttgggtccgtgttattggtctggtccgtgcgATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                                        <  2:1217993/93‑1 (MQ=11)
              tcggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc                                                                                        >  2:1755049/1‑101 (MQ=255)
              tccggggggtcggtggggtccgtggggtccgtcgggtccgtgggaTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                                        <  1:3667658/92‑1 (MQ=11)
               ccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcg                                                                                       <  2:1415952‑M1/96‑1 (MQ=255)
                ggttggctgggtcgggtccgggggatccgtcgggtccgtgggatcGGTCGGCTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcgg                                                                                      <  2:2127137/92‑1 (MQ=14)
                  tggggtccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGAtcggtc                                                                                    >  1:4399298/9‑101 (MQ=14)
                  tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                                    >  2:389867‑M1/3‑101 (MQ=255)
                  gcggATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtc                                                                                    <  1:586608/100‑1 (MQ=11)
                  cgagctcggtcggctccgttcgatcggtcgggtccgtgggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                                    <  2:4418090/96‑1 (MQ=11)
                   gggatcggtcgggtccgtgggatcggtggggtcggtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                                                                   <  1:236250/100‑1 (MQ=11)
                   cggGTCGGTGGGATCCGTGGGATCCGTCGGGTCCGTCGGATCCGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                                                   <  2:2997153/101‑1 (MQ=11)
                    ggatcggtcgggtccgtgggatcggtcgggtccgtcGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtcgg                                                                                  <  2:2690769/101‑1 (MQ=11)
                    ggatcggtagggtcggtgggatcggtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                                  <  2:3550932/101‑1 (MQ=11)
                     gctcgggcggatgcgggggatcggtcgggtccgtgggaTCGGGCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgga                                                                                 <  1:4330564/99‑2 (MQ=11)
                       tcggtcgggtccgtgggatcggtcgggtccgtggggTCCGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATGGGTCGGATCCGTGGGGTCGGTCGGGTc                                                                               >  1:457593/1‑101 (MQ=14)
                       tcgggcgggtcggtgggatcggtcgggtccgtggggTCGGTCGGATCCGTGGGAGCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAGCGGTCGGATc                                                                               <  2:1522001/101‑1 (MQ=11)
                        cggtcgggtccgtgggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGGTCCGTCGGATcc                                                                              >  1:1827058/1‑101 (MQ=11)
                         tgcccggtccgtgtgatcggtcggggccgtggggtccgtcgGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCg                                                                             <  2:3048879/96‑1 (MQ=11)
                         ggtcggatccgtgggatcggtcgggtccgtgggatcTGTCGGATCCGTGGGATCTGTCGGGGCCGTGGGATCGGTCTGGTCCGTGGGATCAGGCGGGTctg                                                                             >  1:2266056/1‑99 (MQ=11)
                         gggcgggtccgggggagcggtcgggtccgggggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg                                                                             <  1:1736066/98‑1 (MQ=11)
                         gggcggatccgtgggatccgtggggtccgtgggatcggtCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCg                                                                             <  2:2050128/98‑1 (MQ=11)
                          gtcgggtccgtgggatcggtcgggtccgtggggccgGTCGGATCGATGGGATGAGCGGGGTCCGGGGAATGAGTCGGGTCCGCGGGGTCGGCCGGGTcatt                                                                            >  2:2827943/1‑98 (MQ=9)
                          ctggggcccgtgggatccgccgggtccgtgggatcggtcGGGGCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGt                                                                            <  2:4352146/98‑1 (MQ=11)
                           gcgggtgcgtgggatcggtcgggtccgtgggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGtg                                                                           <  1:1661349/100‑2 (MQ=11)
                            cgggtccgtgggatcggtcgggtccgtggggtccgtCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgg                                                                          <  1:2182037‑M2/101‑3 (MQ=255)
                            cgggtccgtgggatcggtcgggtccgtggggtccgtCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgg                                                                          <  2:3467986‑M2/101‑3 (MQ=255)
                            cggatccgtgggatcggtcggatccgtaggatcggtCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGtcg                                                                          >  1:2636330/1‑99 (MQ=11)
                             gggttcgttggaatggtcgggtccgtgggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGgg                                                                         <  2:1595632/101‑1 (MQ=11)
                             gggtccgtggggtccgtcggatccgtgggatcggtcGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGACCGTGgg                                                                         >  2:632964/1‑101 (MQ=11)
                             gggtccggggggtcggtcgggtccgtgggagcggtcGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGgg                                                                         <  1:2906625/101‑1 (MQ=11)
                              ggtccgtgggatcggtcgggtccgtggggtcggtcgGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggga                                                                        >  1:3333737‑M2/1‑97 (MQ=255)
                              ggtccgtgggatcggtcgggtccgtggggtccgtcgGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggga                                                                        >  1:476608‑M2/1‑97 (MQ=255)
                              ggtccgtgggatcggtcgggtccgtcggatcggtggGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGggg                                                                        >  1:2550950/1‑100 (MQ=11)
                                tccgtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatc                                                                      >  1:1260602‑M2/1‑95 (MQ=255)
                                tccgtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatc                                                                      >  2:130470‑M2/1‑95 (MQ=255)
                                tccgggggatcgggcgggtccgggggatcggtcggCTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                      <  1:4209479/101‑1 (MQ=14)
                                gctgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc                                                                      <  2:3828589/98‑1 (MQ=14)
                                 ccgtgggatttgtccgctccctgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcg                                                                     <  2:3726582/101‑1 (MQ=11)
                                 ccgtgggatcggtcgggtccgtggggtcggtcggGTCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGGTCCGTGGGCTCCGTCGGGTCCGTGGGAtcg                                                                     >  2:462690/1‑101 (MQ=11)
                                 ccgtgggatcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGTTCCGTGGGCtcg                                                                     >  2:4247832/1‑101 (MQ=11)
                                  cgtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcgg                                                                    <  1:3925535‑M2/101‑9 (MQ=255)
                                   gtgggatcggtcgggtccgtggggccggtgggGTCCGTGGGATCGGTCGCGCCCGTGGGCTCGGCAGGCTCCGTCGGATCGGCGGGGTCCgggcggaccgg                                                                   >  1:3596778/1‑91 (MQ=2)
                                   gtgggatcggtcgggtccgtgggatcggtcggATCCGTGGGGTCGGTCGGGTCCGTGGGGTCGGTCGGGTCGGTGGGATCGGGCGGGTCCGTGGGAtcggg                                                                   >  2:583252/1‑100 (MQ=14)
                                   gtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggt                                                                   <  2:743433‑M2/101‑10 (MQ=255)
                                    tgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                  >  1:924531/1‑101 (MQ=14)
                                    tgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtc                                                                  <  1:3238914‑M2/101‑11 (MQ=255)
                                     gggatcggccgggtccgtgggatcggtcggGTCCGGCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                                 <  1:377669/101‑1 (MQ=14)
                                          cggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGGCCGTGGGGTCGgtcggggcc                                                            >  2:3955425/1‑101 (MQ=11)
                                           ggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGGCGGGTCCGTGGGCTCGgtcgggtcgg                                                           >  2:908137/1‑99 (MQ=11)
                                              cgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGggcgggtccgtcg                                                        >  2:2539762/1‑101 (MQ=11)
                                                ggtccgtcgggtcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCAgtcgggtccgtcgga                                                      <  2:1989398/101‑2 (MQ=11)
                                                gatccgtaggatgggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcggg                                                      <  1:1003866/99‑1 (MQ=11)
                                                     gtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTcggggccgtcggatccgt                                                 >  1:4468212/1‑101 (MQ=11)
                                                     gtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGGCCGTCGGATCCGtggggtccgccgggtccgg                                                 >  2:4098854/4‑100 (MQ=11)
                                                      gggggtcggccggGTCCGTGGGGTCGGTCGGGGCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcgggtccgtggggtccgtc                                                <  2:2910577/100‑2 (MQ=11)
                                                       gggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTgggatcggtcgggtccgtcggatcggtcg                                               <  1:841825‑M2/101‑30 (MQ=255)
                                                        gggtccgtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgg                                              >  2:2212326‑M2/1‑71 (MQ=255)
                                                         ggtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTGCGGACCGTCGGAGCCGTGGGa                                             >  1:147304/6‑101 (MQ=11)
                                                         gatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGGAGGGTCCGTgggg                                             >  2:2366497/1‑100 (MQ=11)
                                                          atcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggt                                            >  1:2451517‑M2/1‑69 (MQ=255)
                                                           tcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtc                                           >  2:3137469‑M2/1‑68 (MQ=255)
                                                            cggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcg                                          >  1:3163767/1‑101 (MQ=11)
                                                            ccgtcggatCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtcc                                          <  1:2031803‑M2/99‑35 (MQ=255)
                                                             cgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcgg                                         >  1:15768/5‑101 (MQ=14)
                                                             cgtcggaTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccg                                         <  2:1129257‑M2/100‑36 (MQ=255)
                                                                cggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtcg                                      >  2:1673580/1‑101 (MQ=11)
                                                                 ggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTTGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccgtggg                                     <  1:754465/101‑1 (MQ=11)
                                                                    tCGGTGGGCTCGGTCTGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatc                                  >  2:2283486‑M2/1‑59 (MQ=255)
                                                                    tCCGGGGGATCGGTCGGGTCCGTGGGATCGGACGGATCCGTGGGATCGTTCGGGTCCGCCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTc                                  <  1:1185119/101‑1 (MQ=14)
                                                                     ccGTGGGATCGGTCGCCTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcg                                 >  2:2625692‑M2/1‑58 (MQ=255)
                                                                      ggtcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCg                                >  1:682716/5‑101 (MQ=17)
                                                                        tcggATCCGTGGGGTCCGTGTGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTg                              <  2:2755818/99‑1 (MQ=18)
                                                                         gggATGGGTCGAGCCGGCGGGATCGGTTGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTgg                             <  1:2834746/101‑1 (MQ=14)
                                                                         gggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcgggcgggtccgtgggatcggtcgggtccgtgggatcggtcg                             <  2:3305534‑M2/101‑48 (MQ=255)
                                                                         aggATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGtcg                             >  1:269421/2‑99 (MQ=14)
                                                                         aggATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGGCCGTCGGCTCCGTGGGGTCCGTCGGGGCCGTGGGATCGGTCGGGTCCGtcg                             >  2:3377988/2‑99 (MQ=11)
                                                                          ggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgg                            <  1:53311‑M2/101‑49 (MQ=255)
                                                                          ggATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtagg                            <  1:2993585‑M2/101‑49 (MQ=255)
                                                                            aTCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt                          >  1:1967427/1‑101 (MQ=255)
                                                                             tcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatc                         <  1:2756016‑M2/101‑52 (MQ=255)
                                                                                 tgggttcCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTcggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtg                     >  1:4064032‑M2/3‑46 (MQ=255)
                                                                                   gggTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGGCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGCCCGCCGGAtcggtccg                   >  1:1075131/1‑99 (MQ=14)
                                                                                     gTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggat                 >  2:1009617‑M2/1‑42 (MQ=255)
                                                                                     gTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggat                 >  1:3503168‑M2/1‑42 (MQ=255)
                                                                                     atcCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtggggt                 <  1:4117634‑M2/100‑60 (MQ=255)
                                                                                        cGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCg              <  2:4030125/101‑1 (MQ=17)
                                                                                          tcggatCGGTGGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtggggtcggtc            <  2:1909496‑M2/99‑65 (MQ=255)
                                                                                           gggATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcg           >  2:3652543‑M2/1‑36 (MQ=255)
                                                                                            ggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg          >  1:3043151/1‑101 (MQ=255)
                                                                                             gATCGGTCGGGTCCGTGGGATCGGTCGGCTCCGTgggagcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggg         >  2:4536240‑M2/1‑34 (MQ=255)
                                                                                               tcggtcggGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtc       <  2:2451351‑M2/101‑70 (MQ=255)
                                                                                                  gtcggGGCCGTGGGATCGGTCGGGTCCGTAGGATCGGGCGGGTCCGTCCGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGGCGGGTCCGCGGGGCtggg    >  2:4162520/1‑100 (MQ=21)
                                                                                                    |                                                                                                    
TCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGG  >  NC_014034/221936‑222136

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 40 ≤  ATCG