Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,066 A→G 3.1% D197D (GAT→GAC RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0660A→G 3.1% 20.8 ‑1.7 0.0596D197D (GAT→GACRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (13/8):  ref base (255/315):  total (271/325)
Fisher's exact test strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 769 total aligned reads displayed.

GGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGG  >  NC_014034/221966‑222166
                                                                                                    |                                                                                                    
gctcggtgggctcgggctgccccgggggctccggcggATCCGTGCGCTCGCTCGCATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGa                                                                                                      <  2:1321930/96‑1 (MQ=2)
 ttccgtgggatcgggcggatccgtgggatctggcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAt                                                                                                     <  1:1767942/100‑1 (MQ=11)
 tcttccgatctcggtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGAt                                                                                                     <  2:1018492/91‑1 (MQ=11)
 gtccgtgggtttcgtccggtgcgtcggacgggtcgcCTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt                                                                                                     <  2:3162210/101‑1 (MQ=11)
  tctgtgggatgggtcgggtccgggggatcggccggGTCCGGGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                                                    <  1:243160/98‑1 (MQ=11)
   ccgtgggatcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcg                                                                                                   >  2:4104773/1‑101 (MQ=11)
    cgtgggctcggcggggtccgtgggctcggtcggGTCCGTGGGAGCCGTCGGCTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                                                                                  <  1:1294809/101‑1 (MQ=11)
      tgggctcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCCGCTCCGTGGGATCGGTCGGATCCCTCGTTTCTGTTGGATCCCTGGGTtcggtc                                                                                                >  2:2815210/1‑101 (MQ=2)
      tgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAGCGGGCGGGGCCGGGGGAtcggtc                                                                                                >  2:846358/1‑101 (MQ=11)
      tgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                                                <  2:2340455/101‑1 (MQ=14)
       gggatggggcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                                                               <  1:1841940/101‑1 (MQ=11)
       gggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTACGTGGGAtcggtcg                                                                                               >  2:487622/1‑101 (MQ=11)
       gggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGCtcggtcg                                                                                               >  1:2764250/1‑101 (MQ=14)
       gggagcggtcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                                                               <  1:819981/101‑1 (MQ=14)
        ggattggtcgggtcgcggggatccctcggATCCGTGGGATCGGTCGGTTCCGTGGGATCGGTCGGCTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                                              <  2:3289274/101‑1 (MQ=11)
        ggatcggtcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGGGGGAtcggtcgg                                                                                              >  2:551159/1‑101 (MQ=14)
         ggtcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgga                                                                                             <  1:1071050/99‑2 (MQ=11)
         gatcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGGGGGATCGgtcggg                                                                                             >  1:1529639/1‑101 (MQ=11)
          atcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGGTCCGGGGGATCGgtggggc                                                                                            >  2:3215973/1‑100 (MQ=11)
           tcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtc                                                                                           <  2:3535779/101‑1 (MQ=255)
            cggtcggttccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatcc                                                                                          <  2:2275065/101‑1 (MQ=11)
            cggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtcc                                                                                          <  2:3983283/101‑1 (MQ=11)
             ggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGAGCGGTCGGGTCCGTCGGATCGGTCGGGGCCGTGGGATCGggcggttccg                                                                                         >  2:2153333/1‑101 (MQ=11)
                cgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGTAGGGGCCGGGGGATCGgtcgggcccgtgg                                                                                      >  2:2539121/1‑99 (MQ=11)
                 ggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcgg                                                                                     <  1:1177392/98‑1 (MQ=255)
                  tatccgtggggtcggtccgATCGGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtggga                                                                                    <  2:1571360/99‑2 (MQ=11)
                  ggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCTGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtgggg                                                                                    >  1:3703656/1‑101 (MQ=11)
                  gatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcggg                                                                                    >  1:2991610/3‑101 (MQ=11)
                   atccgttgtatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcggat                                                                                   <  2:2871380/100‑3 (MQ=14)
                   atccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcggat                                                                                   <  1:1001929/100‑3 (MQ=255)
                   atccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCggggggtccgtcgggt                                                                                   >  2:3927383/2‑101 (MQ=11)
                    tgcgggggatcggtcggTGACGTGGGCTCGGTCGGGTCCGTGGGCTCGGTCGGATCCTTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtgggatc                                                                                  <  2:2012941/96‑1 (MQ=11)
                    tccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcggatc                                                                                  <  1:576126/101‑1 (MQ=255)
                      cgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcggatcgg                                                                                <  2:1530672/101‑3 (MQ=255)
                       gtgggatcgggcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcgggtccgtggggtccgt                                                                               <  1:457253/101‑1 (MQ=11)
                       gggggctcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcgggtccgtggggtccgt                                                                               <  2:4040744/99‑1 (MQ=11)
                        tgggttccgtcggTTCCGTGGGATCGGTCGGCTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcgggtccgtggggtccgtc                                                                              <  1:4250921/101‑2 (MQ=11)
                        tggggtcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcgggtccgtggggtccgtc                                                                              <  2:1916290/101‑2 (MQ=14)
                        tgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGtcgggtccgtcgggtccggg                                                                              >  2:1258929/1‑99 (MQ=11)
                          ggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTcgg                                                                            >  1:3044360/1‑101 (MQ=11)
                           ggtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCCGTCGGATCGGGGGGa                                                                           >  2:2027708/6‑101 (MQ=11)
                           gatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGGCCGTCCGGGCCGTgggg                                                                           >  2:3106460/1‑100 (MQ=11)
                            gtcggtcggGTCCGTAGGCTCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAt                                                                          <  2:1124629/100‑1 (MQ=11)
                             tccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtc                                                                         <  2:2751487/98‑1 (MQ=11)
                              ccgtcggTTCCGTGGGTTCGGTCGGATGCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtcc                                                                        <  2:3395885/99‑2 (MQ=11)
                               cgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtccg                                                                       <  2:1805137/100‑3 (MQ=11)
                                 tcggGTCGGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGATCCGTGGGGTCCGTGGGAtcggtc                                                                     <  1:3170470/101‑1 (MQ=11)
                                 tcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtc                                                                     >  2:661257/1‑101 (MQ=14)
                                 tcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtc                                                                     <  2:967360/101‑1 (MQ=14)
                                  cggGTCCGTTGGATCTGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcg                                                                    <  1:442132/101‑1 (MQ=14)
                                  cggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtccgtgg                                                                    >  1:1344907/1‑99 (MQ=11)
                                  cggATCCGACCGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcg                                                                    <  1:3951593/101‑1 (MQ=17)
                                   gggtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcgg                                                                   >  1:4103805/4‑101 (MQ=14)
                                   ggATCCGTGGGATCGGTCGGATCCGGGGGATCGGGCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtccgtggg                                                                   <  1:2641645/101‑1 (MQ=11)
                                   gaaTCCCTAGGTTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGCCCGTGGGATCGGTCGGAtccgtggg                                                                   <  2:384599/99‑1 (MQ=14)
                                    ggtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtcggg                                                                  >  1:1611775/3‑100 (MQ=11)
                                      tcggTGGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATc                                                                <  1:3440835/98‑1 (MQ=18)
                                      tCCGGGGGATCGGTCGGGTCCGTGGGATCGGACGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTc                                                                <  1:1461891/101‑1 (MQ=14)
                                        ggtcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCg                                                              >  2:2960407/5‑101 (MQ=17)
                                        ggtcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCg                                                              <  1:1239796/97‑1 (MQ=17)
                                        ggTGGGTTCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCg                                                              <  1:3253409/100‑1 (MQ=18)
                                        cGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTGGGATCGGTCGGGTCCg                                                              >  1:2337750/1‑101 (MQ=255)
                                          taggATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGCACGTCGACACCGTGCGGACCGTCGGAGCCGTg                                                            >  2:1481528/3‑101 (MQ=2)
                                           aggATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGtcg                                                           <  2:1725775/100‑3 (MQ=14)
                                            ggATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGGCCGTGGGATCGGTCGGGTCCGGggg                                                          >  2:3728798/1‑101 (MQ=14)
                                                 ggtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGCCCCTGGGAacgg                                                     >  1:876502/1‑101 (MQ=14)
                                                    ggggTCCGTGGGATCGGGCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                                  <  1:270629/100‑1 (MQ=255)
                                                     gggTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGGtcggttgg                                                 >  2:1387931/1‑101 (MQ=14)
                                                     gggTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtcgg                                                 <  1:3918430/101‑1 (MQ=17)
                                                      gggggggtggATCTGTCGGATCCGTGGGCGTGGGGGGGTCCGTGGGCTCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCggg                                                <  1:2851417/93‑1 (MQ=2)
                                                      ggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCggg                                                <  2:2412477/101‑1 (MQ=255)
                                                      gatCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATGGGCCGGGTCCGTGGGATCGGTCGGGTCCGCAGGGTGGGTCggg                                                >  2:2610413/3‑101 (MQ=14)
                                                       gTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGCTCCGTCGGTGCCGTGGGCGCCGTCGGGGCGGTGGGGtgggtggagg                                               >  2:4153104/1‑96 (MQ=2)
                                                         ccGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                             <  1:1388028/101‑1 (MQ=255)
                                                         ccGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGATCCGTGGGGTGCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGAcc                                             >  2:4403865/1‑101 (MQ=14)
                                                           gtcgggtaggtcgcATACGTAGGATCGGGTGGGTCAGTCGCATCCGTTGGATATGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGACCGGTCGGGTCCGt                                           <  2:756741/93‑1 (MQ=2)
                                                           gtcggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGGCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCCGTCGGGTCCGt                                           >  1:349569/4‑101 (MQ=17)
                                                           gTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCgg                                           >  1:4048617/1‑100 (MQ=255)
                                                            tcggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGCCCGACGGGGCCGCGGGGTCGGGCGGGGGCGGCGGATTGGTGGGGTCCGTg                                          >  2:1677357/3‑101 (MQ=9)
                                                             gggATCGGTCGGGTCCGTGGAATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAAAGGTCGGGTCCGTgg                                         >  2:29273/1‑101 (MQ=18)
                                                             gggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtag                                         <  1:2138237/101‑3 (MQ=17)
                                                              gggtcggtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGGTCGGTCGGGTCCGCCGGATCGGTCGGGCCCGTggg                                        >  1:2560366/4‑101 (MQ=17)
                                                              ggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTggg                                        >  1:694404/1‑101 (MQ=18)
                                                               gATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggg                                       <  1:4002599/101‑1 (MQ=255)
                                                                gtcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggat                                      <  2:42591/100‑3 (MQ=14)
                                                                 tcggtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatc                                     <  1:1989688/101‑4 (MQ=14)
                                                                 acggtcggAGCCGTGGGACCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatc                                     <  1:2251199/100‑4 (MQ=14)
                                                                  cggtcggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtg                                    >  1:4351640/1‑101 (MQ=255)
                                                                           ccGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGGCCGTGGGGCTGGTGGGTtcg                           >  1:3773651/1‑101 (MQ=255)
                                                                                  ttcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGAc                    <  2:3020501/100‑1 (MQ=255)
                                                                                  aTCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGAc                    >  1:759646/1‑101 (MQ=255)
                                                                                  aTCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGGCCGTGGGGCTGGTGGGATCGGTCGGAc                    >  1:2532582/1‑101 (MQ=255)
                                                                                   tcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGAcc                   <  1:1701155/101‑1 (MQ=255)
                                                                                      ggcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCt                <  1:3873868/99‑1 (MQ=255)
                                                                                         gggTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTc             >  1:3588800/1‑101 (MQ=255)
                                                                                            tCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgcc          <  1:3360620/101‑1 (MQ=255)
                                                                                             ccGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgccg         >  1:60896/1‑101 (MQ=255)
                                                                                              cGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgccgc        >  2:3459858/1‑101 (MQ=255)
                                                                                                  ggATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCgccgccgcc    <  1:968569/101‑1 (MQ=255)
                                                                                                   gATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGccgccgccg   <  1:203251/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
GGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGG  >  NC_014034/221966‑222166

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 22 ≤  ATCG  < 31 ≤  ATCG  < 40 ≤  ATCG