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 (144/188):  ref base (122/71):  total (333/341)
Fisher's exact test strand distribution p-value = 1.32e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.50e-02

Only 100 of 861 total aligned reads displayed.

GTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGC  >  NC_014034/221962‑222163
                                                                                                     |                                                                                                    
 gcgtggccgtgggatcggtcggatccgggggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTa                                                                                                      <  2:665355/95‑1 (MQ=11)
   gggtccgtgggatcggtcgggtccgtgggatcggtcGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGgg                                                                                                    >  1:348084/1‑101 (MQ=11)
   gggtccgtgggatcggtcgggtccgtgggatcggtcGGGTCCGTGGGATCGGTCGGGTGCGAGGGATCGGTCGGGTCCGTGGGGGCCGTCGGGTGCGTGgg                                                                                                    >  2:3919840/1‑101 (MQ=11)
   gggtccgtgggatcggtcgggtccgtgggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGgg                                                                                                    <  1:3125997/101‑1 (MQ=11)
   gggtccgtgggatcggtcgggtccgtgggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGgg                                                                                                    <  1:1008170/101‑1 (MQ=11)
   gggtccgtgggatcggtcgggtccgtcggatctgtgGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGgg                                                                                                    <  2:4470923/101‑1 (MQ=11)
   ggatctgtgggatcggtcgggtgcgggggatcggtctgATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGgg                                                                                                    <  1:3340794/95‑1 (MQ=14)
     gtccgtgggatcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTTGGATCCGTGGGATCGGTCGGGTTCGTGGGAt                                                                                                  >  1:4182538/1‑101 (MQ=14)
     ctccgtgggtccggtcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAt                                                                                                  <  2:3741004/100‑1 (MQ=11)
     atccgtgggatcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGCCGGCTCCGTGGGAt                                                                                                  >  1:716710/2‑101 (MQ=11)
      gccgtggggtcgttgggttccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                                                 <  1:2456107/100‑1 (MQ=14)
       ccgtgggatcggtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGAtcc                                                                                                <  1:2042265/101‑2 (MQ=11)
       ccgtgggatcggtcggatccgggggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcg                                                                                                <  2:5054321/101‑1 (MQ=14)
        cgtgggatcggtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGAtccg                                                                                               <  1:1188791/101‑3 (MQ=11)
         gtggggtgggtcggacacgtgggatcggtcggCCCCGTGGGAACGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggt                                                                                              <  2:2014376/101‑1 (MQ=11)
         gtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGATACGGGGGAccggt                                                                                              >  2:2085050/1‑101 (MQ=14)
         gtgggatcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGAGGGATCGGTCGGATCCGTAGGATCGGTCGCGTCCGTAGGAtcggt                                                                                              >  1:4242424/1‑101 (MQ=11)
          tgggatgggtgtggtccgtggggtccgtcggATCGGTGGGAGCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                                             <  1:2145985/101‑1 (MQ=14)
          tgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTGGGGTCCGCGGGAtcggcc                                                                                             >  2:4980519/1‑99 (MQ=14)
          tgggatcggtcgggtccgtggggtccgtcggATCCGTGGCATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGTtcggtg                                                                                             >  1:2492561/1‑100 (MQ=11)
          tgggatcggtcgggtccgtgggatcggtcggATCCGTG.........GGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtc                                                                                    >  1:4779939/1‑101 (MQ=11)
          tgggatcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcccac                                                                                             <  2:4000632/101‑5 (MQ=11)
           ggggtccgtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGCCCGTGGGAtcgtgcg                                                                                            >  1:1480797/1‑97 (MQ=11)
           gggatcggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCAGGCGGATCGGTCGGGTCCGTGGGGtcgggcg                                                                                            >  2:4181146/1‑101 (MQ=11)
            ggtttgatcgtgtccgtggggtctttcggATCCGTGGGATTGGTCGGACCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                                           <  1:800239/94‑1 (MQ=11)
            ggatcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTGGGATCCGTAGGAtcgggcgg                                                                                           >  2:623000/1‑101 (MQ=11)
            ggaggggtcgggtccgtgggggccgtcggATCCGTGGGAGCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                                           <  1:1777129/96‑1 (MQ=11)
             gatcgggcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgga                                                                                          <  2:2175566/101‑2 (MQ=11)
              gtcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggt                                                                                         <  1:1534433/100‑1 (MQ=11)
              atcggtgggatccgtgggatcggccggGTCCGGGGGATCGGTCGGGTCCGTGGGACCGGTCAGGGCCGTGGTATCGGTCGGGGCCGTGGGTtgggtcggcc                                                                                         >  2:850949/1‑99 (MQ=2)
               tcggtagggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtc                                                                                        <  1:1437389/101‑1 (MQ=11)
                cggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtcc                                                                                       <  1:4435786/101‑1 (MQ=11)
                 ggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatccg                                                                                      >  2:1222591/1‑101 (MQ=11)
                 ggtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccg                                                                                      <  2:1075842/101‑1 (MQ=11)
                 cgtcgggtccggggggtccgtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGgtgggatccg                                                                                      <  2:4736538/100‑1 (MQ=11)
                  gtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtcggt                                                                                     >  2:3267742/1‑101 (MQ=11)
                   tcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCAGTGGGATCGGTCGGGCCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtg                                                                                    <  2:2981797/101‑2 (MQ=11)
                    cggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcg                                                                                   >  2:1282579/1‑101 (MQ=255)
                     ggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgg                                                                                  >  2:4646560/4‑101 (MQ=11)
                      ggtccgtcggatcggtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGCCGGGGCCgccgggtccggcggt                                                                                 >  1:1088586/1‑100 (MQ=14)
                          cgtaggacctgtcggGTCCGCAGGATCTGTCGGGTACGTCGGAGCCGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGtggggtccgtcgggtccg                                                                             <  2:2443729/101‑1 (MQ=9)
                           gtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTcgggtccgtcggatcggt                                                                            >  2:4311654/1‑101 (MQ=11)
                            taggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTg                                                                           <  1:4670624/99‑1 (MQ=11)
                            taggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGCTCCGTg                                                                           >  1:3270389/3‑101 (MQ=11)
                             ggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcgggtccgtggggtccgtcg                                                                          >  1:4809186/1‑99 (MQ=11)
                             ggggtccgtcggATCCGTGGGATCGGCCGGATCCGTGGGATCGGTGGGGTCCGTGGGGTCGGTGGGGTCCGTCGGCTCGgtcggctccgtggggcccccag                                                                          >  2:2334648/1‑97 (MQ=2)
                             gggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTcgggtccgtcggatcggtcg                                                                          >  1:5023177/1‑99 (MQ=11)
                             gggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTgg                                                                          >  2:2921417/1‑101 (MQ=11)
                             cgggtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTgg                                                                          >  2:170793/11‑101 (MQ=11)
                             cgggtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGGgg                                                                          >  1:2260910/11‑101 (MQ=11)
                              ggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTggg                                                                         <  1:4907258/101‑1 (MQ=14)
                              ggatcgggcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTcgg                                                                         <  2:3426194/101‑1 (MQ=11)
                              gctgccgtcgtATCCGTCGGATCGGTCGGCTCCGGCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcgg                                                                         <  2:2669122/95‑1 (MQ=11)
                               ggtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGa                                                                        <  1:165614/96‑1 (MQ=11)
                                gcgcggaagATCCGTGGGGTAGGTGGGGTGCGTGGGAGGGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAt                                                                       <  1:2608537/93‑1 (MQ=2)
                                atcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTggggc                                                                       >  2:4918604/1‑99 (MQ=11)
                                 tccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGAtc                                                                      <  1:2819288/98‑1 (MQ=11)
                                 gccgtgggATCGGTCGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtc                                                                      <  1:1934419/95‑1 (MQ=14)
                                  cggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcg                                                                     <  2:2124913/101‑1 (MQ=14)
                                  ccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtcc                                                                     <  2:1648481/99‑2 (MQ=11)
                                    gtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGATCCGTGGGGTCCGTGGGAtcggg                                                                   >  1:1271579/1‑100 (MQ=11)
                                    gtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGGtcggt                                                                   >  2:2835861/1‑101 (MQ=14)
                                     tcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtc                                                                  <  1:2700319/101‑1 (MQ=11)
                                     tcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtc                                                                  >  1:2991097/1‑101 (MQ=11)
                                     tcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtg                                                                  >  2:484697/1‑100 (MQ=14)
                                      cggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtccgtgg                                                                 >  1:3600135/1‑99 (MQ=11)
                                       ggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtccgtggg                                                                >  1:557655/1‑101 (MQ=11)
                                        ggtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCAGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcggg                                                               >  2:365174/3‑100 (MQ=14)
                                         aTCGGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGAt                                                              <  1:3708260/101‑1 (MQ=18)
                                          tCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATc                                                             <  2:2609418/101‑1 (MQ=11)
                                                    cggtcggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcg                                                   <  1:938009/101‑1 (MQ=14)
                                                    ccggggggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg                                                   <  1:3967990/96‑1 (MQ=255)
                                                      gtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGGCGTAGGAtcggc                                                 >  2:4057271/1‑100 (MQ=14)
                                                          ggTCCGTGGGATCGGCCGGGTCCGTGGGATCGGTCGGGTCCGCAGGATCGGTCCGGTCGCCCGGGCCCGGGGGATAGGTCGGCTGCGCGGGCtcggcctgg                                             >  2:335467/1‑96 (MQ=9)
                                                           gTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGGCCGTCGGGTCGGTCgggg                                            >  2:3532677/1‑100 (MQ=14)
                                                            tCCGTGGGGTCGGGCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTc                                           <  1:1330441/101‑1 (MQ=17)
                                                            tCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGACCGGTCGGGTCCGTCGGATCGGTCGGGTc                                           >  2:1438281/1‑101 (MQ=14)
                                                             ccGTGGGATCGGTCGGGCCCGTGGGATCGGTCGGATCCGTGGGCTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTcc                                          <  2:1382709/101‑1 (MQ=11)
                                                             ccGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                          >  1:4127559/1‑101 (MQ=17)
                                                             ccGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTcc                                          >  1:3879751/1‑101 (MQ=17)
                                                              cGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCg                                         >  2:4382459/1‑101 (MQ=11)
                                                                tcggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTg                                       >  2:694330/3‑101 (MQ=17)
                                                                 gggTTCGGTGGGATGCGTGGGATCGGTGGGGGCCGTCGGATCGGGGGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtag                                      <  1:634005/101‑3 (MQ=14)
                                                                 cggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgg                                      >  2:1534303/2‑101 (MQ=17)
                                                                  gggtcggtcggATCCGTGGGATCGGTGGGGTCCTTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTagg                                     <  1:4342425/98‑1 (MQ=17)
                                                                   gATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTagga                                    >  1:891784/1‑100 (MQ=17)
                                                                    aTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGCCCGTGGGATCGGTCGGATCCGTGGGATCGGGCGGGTCCGGGGGGc                                   >  2:4548588/1‑101 (MQ=18)
                                                                    aTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATAGGTCGGGTCCGTaggat                                   >  2:3836478/1‑99 (MQ=17)
                                                                    aTCCGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggat                                   <  2:276367/101‑3 (MQ=14)
                                                                      cggtcggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtg                                 >  2:3135755/1‑101 (MQ=255)
                                                                      cggtcggGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTcgggtcc                                 >  2:783493/1‑98 (MQ=14)
                                                                      ccgtcgcATCCGTGGTATCGGGCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGACCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcg                                 <  2:2519168/99‑5 (MQ=14)
                                                                        gtcggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtggg                               >  1:4042732/1‑100 (MQ=255)
                                                                          cggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtcg                             >  1:2381946/1‑93 (MQ=14)
                                                                                   gggATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggtcg                    >  2:1023297/1‑101 (MQ=255)
                                                                                    ggctcggttggCGCCGAAGGAGCGGGGGGCTCCGTCGGGGCCGTGGGACCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggtcgg                   <  1:4554014/98‑1 (MQ=18)
                                                                                      gtccgtggggTCCGTCGGATCTGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGAc                 <  2:913663/94‑1 (MQ=18)
                                                                                            cggGTCGGTAGGGTCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGt           <  2:186577/101‑1 (MQ=255)
                                                                                             gggTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTc          <  1:3796446/101‑1 (MQ=255)
                                                                                                 ccGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgccg      <  1:4598421/101‑1 (MQ=255)
                                                                                                     |                                                                                                    
GTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGC  >  NC_014034/221962‑222163

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 31 ≤  ATCG  < 40 ≤  ATCG