Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,063 A→G 100% P198P (CCT→CCC RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0630A→G100.0% inf 674P198P (CCT→CCCRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (124/208):  ref base (111/92):  total (299/375)
Fisher's exact test strand distribution p-value = 1.10e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.42e-02

Only 100 of 897 total aligned reads displayed.

TCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCC  >  NC_014034/221963‑222161
                                                                                                    |                                                                                                  
  gggtccgtggggtccgtcggatccgtgggatcggtcGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGACCGTGgg                                                                                                  >  2:632964/1‑101 (MQ=11)
  gggaccgtgggatcggtcgggtccgtggggtccggcGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCCGGTCCGTGgg                                                                                                  >  2:128019/1‑101 (MQ=14)
  ggatcggggggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGgg                                                                                                  <  1:1048578/93‑1 (MQ=11)
  ggatccgtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGgg                                                                                                  >  2:531473/4‑101 (MQ=11)
  ggatccgtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGgg                                                                                                  <  2:2003864/98‑1 (MQ=11)
   gatccgtgggatctgtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGa                                                                                                 <  2:2108706/99‑1 (MQ=11)
    gtccgtgggatcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCGGTGgggt                                                                                                >  2:2275555/1‑99 (MQ=11)
     tccgtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGAGGGATCGGTCGGGTCAGTGGGAtc                                                                                               <  2:120643/101‑1 (MQ=14)
     tccgggggatcggtcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc                                                                                               <  2:836456/101‑1 (MQ=11)
     gccgtgggatcggtggggtgcgtcggttcggtcggGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                                               <  1:3174327/100‑1 (MQ=11)
      gcgtgggatcggtcgggtccgtcgggtcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGAtcg                                                                                              <  1:2275717/100‑1 (MQ=11)
      ccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCCGGGGGATCGGGCGGATCCGTGGGTtcg                                                                                              >  2:2657044/1‑101 (MQ=11)
      ccgtgggatcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGGTCCGTGGGATCGGTCTGGGCCGTGGGAgcg                                                                                              >  2:377632/1‑101 (MQ=11)
       cgtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcgg                                                                                             >  1:1367289/1‑101 (MQ=11)
        gtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                                                                            <  2:2409909/101‑1 (MQ=14)
        gtgggagcgggcggggccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                                                            <  1:4169110/101‑1 (MQ=11)
         ttggatcggtcgggtccctgggttcggtcggGTCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                                           <  2:1380258/99‑1 (MQ=14)
         gggggttgggcggggccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                                           <  2:2939002/91‑1 (MQ=11)
           ggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTGGGGGCCGTGGGAtcgggggg                                                                                         >  2:4031461/1‑97 (MQ=11)
             gtcgggcgggtccgtggggtcgggcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggt                                                                                       <  1:2724485/100‑1 (MQ=14)
             atcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGCGGGATCGGTCGGATCCGTGGGACCGCGCGGGACGGGGGGAtcgggccggc                                                                                       >  2:1690132/1‑95 (MQ=2)
              tcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatc                                                                                      <  2:2572114/101‑1 (MQ=11)
              tcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATGGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtc                                                                                      >  1:3919770/1‑101 (MQ=14)
                ggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatccg                                                                                    >  1:3425158/1‑101 (MQ=11)
                ggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGGCGGGTCCGTGGGCTCGgtcgggtcgg                                                                                    >  2:908137/1‑99 (MQ=11)
                  tcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtc                                                                                  <  1:125622/101‑1 (MQ=11)
                   ggggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcg                                                                                 <  2:341653/100‑1 (MQ=255)
                   gggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcg                                                                                 <  2:4435885/97‑1 (MQ=11)
                   cgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCGGTGGGATCGgtggggtccgcgg                                                                                 >  1:3068541/1‑98 (MQ=11)
                    gggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGGTAGGTCGGGCCCGTGGGATCGgtcgggtccgtcgg                                                                                >  1:3971322/1‑101 (MQ=11)
                    gggtccggggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatccgtggg                                                                                <  1:173160/101‑1 (MQ=11)
                    ggatccgtggggtccgtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGggggggtccgtcgg                                                                                >  2:586071/4‑101 (MQ=11)
                    ggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcgg                                                                                >  1:1830145/4‑101 (MQ=255)
                    ggatccgtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGgtcgggtccgtcgg                                                                                >  1:3551190/4‑101 (MQ=11)
                     ggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAAGAgtggggtccggggga                                                                               >  2:3906500/1‑96 (MQ=11)
                     ggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGGCCGTGGGTTCGgcgcgggccgtcgga                                                                               >  2:2699788/1‑100 (MQ=11)
                     gatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcgga                                                                               >  2:4238997/3‑100 (MQ=255)
                     gatccgtaggatgggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcggg                                                                               <  1:1003866/99‑1 (MQ=11)
                     gatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcggg                                                                               >  2:3017946/3‑101 (MQ=11)
                     gatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgccggg                                                                               >  2:4587036/3‑101 (MQ=11)
                        ccgtaggatcagtcggGTGTGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGTTCAGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgggtcc                                                                            <  1:4643572/101‑1 (MQ=11)
                         cgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGAACGGTCGGATCCGTGGGGTCGGTCGGGTCCGTGGGATCGgtggaggccgtggaatcgg                                                                           >  2:2834530/1‑93 (MQ=11)
                         cgtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtcggatcgg                                                                           >  1:2671544/1‑99 (MQ=11)
                          gtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcgggtccgtggggtccgt                                                                          <  1:978990/101‑1 (MQ=14)
                          gtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCCGTggggtccgtcggggccgt                                                                          >  1:2048653/4‑101 (MQ=11)
                           ttggagcggtccgGTCCGTGGGGTCGGTCCGGTCCGTGGGATCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcgggtccgtggggtccgtc                                                                         <  2:3666910/99‑2 (MQ=11)
                           tgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGGTCCGTg                                                                         >  2:481655/1‑101 (MQ=14)
                            ggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGGGGGCTCGGTCGGGTCCGTGGCATCGgttgcatccgcggggtcgtgcg                                                                        >  1:2166368/1‑96 (MQ=0)
                            cgggtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTgg                                                                        >  2:3157847/11‑101 (MQ=11)
                            cgggtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTgg                                                                        >  1:2406381/11‑101 (MQ=11)
                              ggtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGa                                                                      <  2:4103016/96‑1 (MQ=11)
                              gatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGa                                                                      <  1:4239285/101‑1 (MQ=14)
                               gtcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAt                                                                     <  1:1492107/100‑1 (MQ=11)
                               gtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATAGGTCGGGTCCGTCGGGCTCGTGGGAt                                                                     >  2:722339/8‑101 (MQ=11)
                               gtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAt                                                                     >  2:3327019/5‑101 (MQ=11)
                               atcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTggggt                                                                     >  1:2823377/1‑99 (MQ=11)
                                tccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtc                                                                    >  1:3026450/4‑101 (MQ=11)
                                tccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtc                                                                    <  2:2181887/98‑1 (MQ=11)
                                 ccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAGCGGTCGGGTCCGAGGGGTCCGTCGGAtcc                                                                   >  1:609971/3‑100 (MQ=11)
                                  ggtcggATCCGTGGGATCGGTGGGATCCGTGGGATAGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccg                                                                  <  1:1167358/101‑3 (MQ=11)
                                  cgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcgg                                                                  >  1:521701/5‑101 (MQ=14)
                                  cgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGCAGGATCGGTCGGGCCCGTCGGGGGCGTGGCAccgg                                                                  >  2:3613420/5‑101 (MQ=11)
                                   gccggGTCCGTGGGAGCTGTCGGCTCCGTGCGTTCGCGCGGCTCCGTCGAATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggt                                                                 <  2:1933650/99‑1 (MQ=9)
                                    gcggGTCCGTAGGAATGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtc                                                                <  2:2892977/100‑1 (MQ=14)
                                      gggtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcgg                                                              >  1:4346157/4‑101 (MQ=14)
                                      ggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcggtcgg                                                              <  2:4319772/101‑1 (MQ=17)
                                      ggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTTGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccgtggg                                                              <  1:754465/101‑1 (MQ=11)
                                       gATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGa                                                             <  1:2589870/101‑1 (MQ=11)
                                        aTCCGTGGGATCGGCCGGAGCCGTGGGATCGGTCGGATCCGGGGGATCGGTCGGGTCCGTGGGATGGGTCGGGGCCGTGGGGTCCGCCGGATCCGTGGGAt                                                            >  2:386572/1‑101 (MQ=11)
                                           ggtcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCg                                                         >  1:682716/5‑101 (MQ=17)
                                           cGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGGGGGGGAGGAGGGGCGCGTCGGGTCGGTGGGGtcagtggggtggg                                                         >  2:3234369/1‑96 (MQ=2)
                                           cGTAGGATCGGTCGGGTCCGTGGGATCCGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCCTGGGGTCGCTGGGGTCCg                                                         >  1:3636704/1‑101 (MQ=14)
                                            gtcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGt                                                        >  1:1560786/4‑101 (MQ=17)
                                            gtcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGGCCGTGGGATCGGTCGGATCCGt                                                        >  1:890763/4‑101 (MQ=17)
                                             tcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCCTCGGGCCCGTGGGGtccggcgggcacgcg                                                       >  2:1782207/3‑94 (MQ=14)
                                             nggngtcggtcggGGCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTg                                                       <  1:4103683/96‑1 (MQ=14)
                                              gggATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGCTCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTgg                                                      >  2:4164332/1‑101 (MQ=14)
                                                gagaggtcggGTCCGTCGGATCGGTGGGGTCCGTGGGATTGTTCGGGTTCGTGGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGa                                                    <  1:3845479/97‑1 (MQ=11)
                                                  tccttggggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc                                                  <  1:328885/95‑1 (MQ=255)
                                                   ccgtggggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg                                                 <  2:2362492/96‑1 (MQ=255)
                                                      tcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                              <  1:3236324/101‑1 (MQ=17)
                                                      tcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                              >  1:85166/1‑101 (MQ=17)
                                                       ggggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                             <  2:1525304/100‑1 (MQ=255)
                                                        gggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgg                                            >  1:3719130/1‑101 (MQ=255)
                                                         ggTCCGTGGGATGGGTCGGATCCGGGGGATCGGACGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCggg                                           <  1:4031746/101‑1 (MQ=17)
                                                               tGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTg                                     >  2:1619786/1‑101 (MQ=255)
                                                               tGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTg                                     >  1:2934503/1‑101 (MQ=255)
                                                                    tcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTGGGGGGCGGgggggc                                >  2:3472029/1‑99 (MQ=11)
                                                                    tcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatc                                <  1:4559407/101‑4 (MQ=14)
                                                                    tcggtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGGCGGATCGGGCGGGTCCGTgggatg                                >  2:2495668/1‑98 (MQ=17)
                                                                         cggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtggtgg                           >  1:4421138/1‑101 (MQ=255)
                                                                         cggATCGGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtcg                           <  2:3956416/101‑9 (MQ=14)
                                                                         cggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtcg                           >  2:3070588/1‑93 (MQ=14)
                                                                               cgggggAGCGGGGGGGTCCGTAGGATCGGGCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcgg                     <  1:1979261/98‑1 (MQ=21)
                                                                                gttggctcggtccgGTCCGTTGGGCTGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggt                    <  2:2768828/95‑1 (MQ=18)
                                                                                gTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggt                    >  1:2861731/1‑101 (MQ=255)
                                                                                        ggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGc            >  2:4498975/1‑101 (MQ=255)
                                                                                               tCCGTGGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgcc     <  1:2515347/101‑1 (MQ=255)
                                                                                                 cGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgccgc   >  2:1139186/1‑101 (MQ=255)
                                                                                                  gTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgccgcc  >  1:1649984/1‑101 (MQ=255)
                                                                                                    |                                                                                                  
TCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCC  >  NC_014034/221963‑222161

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