Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,949 C→A 8.8% E236D (GAG→GAT RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9490C→A 8.8% 47.1 ‑2.1 ‑0.3250E236D (GAG→GATRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (14/8):  ref base (105/123):  total (119/131)
Fisher's exact test strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.55e-01

Only 100 of 288 total aligned reads displayed.

CCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAT  >  NC_014034/221851‑222049
                                                                                                   |                                                                                                    
 cGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCtcg                                                                                                   >  2:1202657/1‑101 (MQ=255)
  ggTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGNNCncgn                                                                                                  >  1:1984196/1‑100 (MQ=255)
     ttCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCtcggtcg                                                                                               >  2:3212172/1‑101 (MQ=255)
      tCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCtcggtcgg                                                                                              <  1:3907545/101‑1 (MQ=255)
       cGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCgga                                                                                             <  1:4079656/101‑1 (MQ=255)
       cGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCgga                                                                                             <  1:2372898/101‑1 (MQ=255)
        gCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggat                                                                                            <  2:794683/101‑1 (MQ=255)
        gCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggat                                                                                            >  2:1519357/1‑101 (MQ=255)
         ccGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatc                                                                                           <  1:3655995/101‑1 (MQ=255)
         ccGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatc                                                                                           <  2:1970394/101‑1 (MQ=255)
         ccGGTGCCCGACACCGTCGCCGAGGCTTCGGTGATGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatc                                                                                           <  1:4974300/101‑1 (MQ=255)
          cGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcg                                                                                          <  1:2025273/101‑1 (MQ=255)
             tGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtc                                                                                       <  1:1076952/101‑1 (MQ=255)
              gCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcg                                                                                      >  1:4290882/1‑101 (MQ=255)
              gCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcg                                                                                      <  2:1304390/101‑1 (MQ=255)
               cccGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgg                                                                                     <  1:3732157/101‑1 (MQ=255)
                  gACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtc                                                                                  >  2:2404455/1‑101 (MQ=255)
                   acacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtcc                                                                                 <  1:4610673/101‑1 (MQ=255)
                   acacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtcc                                                                                 >  1:1165916/1‑101 (MQ=255)
                    ccccGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccg                                                                                <  2:1650911/99‑1 (MQ=255)
                     acCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccgt                                                                               <  2:2782169/101‑1 (MQ=255)
                           gCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccgtgggatc                                                                         <  2:142847/101‑1 (MQ=255)
                           gCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccgtgggatc                                                                         <  1:1649268/101‑1 (MQ=255)
                           gCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccgtgggatc                                                                         <  1:5667235/101‑1 (MQ=255)
                            ccGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccgtgggatcg                                                                        <  2:3580631/101‑1 (MQ=255)
                             cGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccgtgggatcgg                                                                       <  1:3302589/101‑1 (MQ=255)
                               gggCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                                     <  2:5340971/100‑1 (MQ=255)
                               aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                                     <  1:5730391/101‑1 (MQ=255)
                               aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                                     <  2:2974997/101‑1 (MQ=255)
                                ggCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcg                                                                    <  1:1124732/101‑1 (MQ=255)
                                ggCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccgggggatcgggcg                                                                    >  2:95196/1‑101 (MQ=255)
                                 gCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgg                                                                   >  2:5512189/1‑101 (MQ=255)
                                 gCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgg                                                                   >  1:5300326/1‑101 (MQ=255)
                                  cTTCGGTGTTGGTGGTCAGCGTCGCGGGTCCGCCCTCATCGCGGGGGGGGGCGGGAGGC.GGGGCCTCGGTGggatcggtggggcccgtggggcggggcggg                                                                  >  2:4997613/1‑91 (MQ=255)
                                   ttCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcggggccgtgggatcggtcgggt                                                                 >  2:4884190/1‑101 (MQ=255)
                                    tCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggccgggtccgtgggatcggtcgggtc                                                                >  2:1554566/1‑101 (MQ=255)
                                     cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                               <  2:1797298/101‑1 (MQ=255)
                                     cGGTGTTGGTGGTCAGCGGCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                               <  1:920978/101‑1 (MQ=255)
                                      ggTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccgtggggtcggtcgggtccg                                                              >  2:3988828/1‑101 (MQ=255)
                                       gtgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgg                                                             >  2:892314/1‑100 (MQ=255)
                                       gtgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GGGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgt                                                             >  2:5383788/1‑101 (MQ=255)
                                        tgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgtg                                                            <  2:2151345/101‑1 (MQ=255)
                                        tgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgtg                                                            >  2:100369/1‑101 (MQ=255)
                                         gtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGgatcggtcgggtccgtgggatcggtcgggtccgtgg                                                           >  2:3153284/1‑101 (MQ=255)
                                           tggtggTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGatcggtcgggtccgtgggatcggtcgggtccgtgggg                                                         >  1:4014891/1‑100 (MQ=255)
                                            ggtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATcggtcgggtccgtgggatcggtcgggtccgtgggat                                                        >  1:4568801/1‑101 (MQ=255)
                                            ggtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATcggtcgggtccgtgggatcggtcgggcccgtgggat                                                        >  1:4800205/1‑101 (MQ=255)
                                             gtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCggtcgggtccgtgggatcggtcgggtccgtgggatc                                                       <  1:5042632/101‑1 (MQ=255)
                                                 tCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCgggtccgtgggatcggtcgggtccgtgggatcggtc                                                   <  1:4129137/101‑1 (MQ=255)
                                                 tCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCgggtccgtgggatcggtcgggtccgtgggatcggtc                                                   <  1:2207169/101‑1 (MQ=255)
                                                   aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgg                                                 >  2:4449078/1‑101 (MQ=255)
                                                   aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgg                                                 >  1:2579486/1‑101 (MQ=255)
                                                     cGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt                                               <  2:4799924/101‑1 (MQ=255)
                                                      gTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCAGTGGGATCGGTCGGATc                                              >  2:3004491/1‑101 (MQ=255)
                                                      gTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCCGGTCCGTGGGGTCGGGCGGacc                                              >  2:5237227/1‑99 (MQ=255)
                                                       tCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                                             >  2:5750087/1‑101 (MQ=255)
                                                        cggcGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCg                                            <  1:4449350/98‑1 (MQ=255)
                                                           cGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                         <  1:4869751/101‑1 (MQ=255)
                                                            ggATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                                        >  2:2900764/1‑101 (MQ=255)
                                                              aTCGCGCTGATGGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt                                      <  1:3382555/101‑1 (MQ=255)
                                                              aTCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt                                      <  2:4586381/101‑1 (MQ=255)
                                                                cgcgCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg                                    <  1:3593976/101‑1 (MQ=255)
                                                                 gcgcTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcgg                                   >  1:4060305/1‑101 (MQ=255)
                                                                  cgcTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                  >  1:2457832/1‑101 (MQ=255)
                                                                   gcTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGTTCGGGCGGGTCCGGGGGATCGGTCGGATCCGTGGGAttggtc                                 >  2:2276944/1‑101 (MQ=25)
                                                                   gcTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                 <  1:1575881/101‑1 (MQ=255)
                                                                   gcTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                 <  2:306366/101‑1 (MQ=255)
                                                                     tCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgg                               <  1:648204/101‑1 (MQ=255)
                                                                       aTCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCgggg                             >  1:539305/1‑100 (MQ=255)
                                                                       aTCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGt                             >  2:4353329/1‑101 (MQ=255)
                                                                          gACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCg                          >  2:4076153/1‑101 (MQ=255)
                                                                           gcGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt                         <  1:1386912/100‑1 (MQ=255)
                                                                             gggatcggTCGGGTCC.GTGGGATCGGTCGGATCCGTGGGGCCCGTGGGATCGGTCGGGTCCGTGGGATCGGCCGGGTCCGTAGGAGCGGCTGGGTCCGtcg                       >  2:2599995/8‑99 (MQ=2)
                                                                             cgggtgCGTCGGATTCGGCGGGTTCGG.CGGGTCCGGCGGGTCCGGGGGATCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgg                       <  2:2687867/96‑1 (MQ=2)
                                                                              gTGGGCGTCGGATGG.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg                      <  1:2277070/101‑1 (MQ=255)
                                                                              gTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg                      <  1:3083631/101‑1 (MQ=255)
                                                                               tGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTGGGGTCGGTgggg                     >  2:2029514/1‑100 (MQ=18)
                                                                               tGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGa                     <  1:2797046/101‑1 (MQ=255)
                                                                                gggCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTAGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                    >  2:133832/1‑101 (MQ=255)
                                                                                gccgtgcGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                    <  2:4708678/95‑1 (MQ=255)
                                                                                    gTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                <  1:5088824/101‑1 (MQ=255)
                                                                                     tCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCGGGCGGGTCCGTggggcggggc               >  2:422487/1‑94 (MQ=18)
                                                                                       ggATGG.GTGGGATCCGTGGGGTCGGTCGGGTGTGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg             <  1:5496234/101‑1 (MQ=11)
                                                                                       ggATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg             >  2:34922/1‑101 (MQ=255)
                                                                                       ggATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg             >  1:583036/1‑101 (MQ=255)
                                                                                        cggtctgTGGGATCCGTGGGTTCCGGCGGTTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg             <  1:3239207/95‑1 (MQ=1)
                                                                                         gtGC.GTGGGATCGGGCGGATCCGTGGGGTCCGTGGGATCGGGCGGGGCCGTGGAATCGGTGGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAt           <  1:897333/100‑1 (MQ=11)
                                                                                         atcc.GTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGGTCCGTGGGAtcggtcgggg           >  2:764761/4‑99 (MQ=11)
                                                                                         atcc.GTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATAGGTCGGGTCGGTGGGAtcgggcgggc           >  1:1521754/4‑99 (MQ=11)
                                                                                         aTGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcggcc           >  1:5185634/1‑99 (MQ=255)
                                                                                         aTGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt           >  1:2877781/1‑101 (MQ=255)
                                                                                          tcc.GTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGGTCCGTGGGATCGCTCGGGTCCGTGGTATCGGTCGGGTCCGTAGGAtcggccgtttc          >  1:344889/3‑97 (MQ=14)
                                                                                          tGC.GTGGGCTCGGTCGGATCGCTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc          <  1:125234/101‑1 (MQ=255)
                                                                                            c.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCg        >  2:23262/1‑101 (MQ=255)
                                                                                            c.GTGGGATCGGTTGGTTTCGTGGGTTCGGTGGGATCCTTGGGGGCCGTGGGATCGGGCGGGTCCGGGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCg        <  1:853087/101‑1 (MQ=2)
                                                                                            c.GTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg        <  1:145030/101‑1 (MQ=11)
                                                                                              gTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt       <  1:194000/101‑1 (MQ=11)
                                                                                              gTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGGGGGAtcggtcggggccgg       >  2:841565/1‑95 (MQ=11)
                                                                                                gggCTCGGTCGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg     >  1:5551567/1‑101 (MQ=255)
                                                                                                 ggCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg    >  2:3366721/1‑101 (MQ=255)
                                                                                                   |                                                                                                    
CCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGC.GTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAT  >  NC_014034/221851‑222049

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