Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,961 G→C 5.2% T232T (ACC→ACG RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9610G→C 5.2% 50.0 ‑2.7 ‑0.7563T232T (ACC→ACGRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (21/10):  ref base (265/267):  total (286/277)
Fisher's exact test strand distribution p-value = 6.42e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.21e-01

Only 100 of 1259 total aligned reads displayed.

CGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCG  >  NC_014034/221861‑222061
                                                                                                    |                                                                                                    
ccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtgggctccgtgggatcggtcgggtccgtgggatcg                                                                                                      <  2:4500952‑M1/6‑1 (MQ=255)
   tGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtc                                                                                                   >  2:1217782/1‑101 (MQ=255)
     ggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgg                                                                                                 >  2:2264653‑M1/91‑101 (MQ=255)
        gACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtc                                                                                              <  1:746206/101‑1 (MQ=255)
         cggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtcc                                                                                             >  1:659573‑M1/87‑101 (MQ=255)
         ccgtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtcc                                                                                             <  1:1576407‑M1/15‑1 (MQ=255)
           gtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgt                                                                                           <  1:2813040‑M1/17‑1 (MQ=255)
             cggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggacggggccgtgg                                                                                         <  2:1032007‑M1/19‑1 (MQ=255)
              ggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtggggtcggtcgggtcggtgggatcggtcgggtcggtggg                                                                                        <  2:1110883‑M1/20‑1 (MQ=255)
                   cgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtcggtcgggtccgtgggatcggtcgggtccgtgggatcgg                                                                                   <  2:4076970‑M1/25‑1 (MQ=255)
                    gtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggt                                                                                  <  1:3493373‑M1/26‑1 (MQ=255)
                    gtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggccggtcgggtccgtggggtcggg                                                                                  <  1:1927817‑M1/26‑2 (MQ=255)
                     tgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtc                                                                                 >  2:2915067‑M1/75‑101 (MQ=255)
                     aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                                                 <  1:1784756/101‑1 (MQ=255)
                        gatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcgga                                                                              <  1:4547994‑M1/30‑2 (MQ=255)
                        gatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgga                                                                              >  2:427392‑M1/72‑100 (MQ=255)
                         atcggtcgggtccgtcggatcggtgggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggt                                                                             <  2:2009565‑M1/31‑1 (MQ=255)
                            ggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccg                                                                          <  2:3351831‑M1/34‑1 (MQ=255)
                            ggTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccg                                                                          >  2:1504134/1‑101 (MQ=255)
                            ggTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccg                                                                          >  1:2259851/1‑101 (MQ=255)
                             gtgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgt                                                                         >  1:2573025/1‑101 (MQ=255)
                               cgggtccgtctgatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgg                                                                       >  2:888606‑M1/65‑101 (MQ=255)
                                gggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggg                                                                      <  1:4089689‑M1/38‑1 (MQ=255)
                                ggattcgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtggg                                                                      >  1:937807‑M1/64‑101 (MQ=255)
                                 tggtggTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGGGTCGGATGCGTGGGCTCGGCCGGAtcgggcgggcccgcgggatcggtcgggtccgtggga                                                                     >  1:2020003/1‑101 (MQ=38)
                                 ggtccgtcggatcggtgggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtggga                                                                     >  1:4280082‑M1/63‑101 (MQ=255)
                                    tggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGgtcgggtccgtgggatcggtcgggtccgtgggatcg                                                                  <  1:820189/101‑1 (MQ=255)
                                    ccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcg                                                                  >  1:3041360‑M1/60‑101 (MQ=255)
                                     cgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcgggcgggtccgtgggatcgggcgggtccgtgggatcggtcgggtccgtgggattgg                                                                 <  2:3202744‑M1/43‑4 (MQ=255)
                                        gggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcg                                                              >  2:3652543‑M1/56‑101 (MQ=255)
                                          gatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgca                                                            <  1:1239578‑M1/48‑3 (MQ=255)
                                          gCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGa                                                            >  1:1699744/1‑101 (MQ=255)
                                            gTCGCCGGATCGCGCTCATCGACGGCGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                                                          <  2:718053/101‑1 (MQ=255)
                                             cggtgggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                         <  2:1382769‑M1/51‑1 (MQ=255)
                                             cggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                                         <  1:1447044‑M1/51‑1 (MQ=255)
                                                  gggtccgtgggatcggtcgggtccgtgggaacggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg                                                    <  1:3936208‑M1/56‑1 (MQ=255)
                                                  ggATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                                                    >  2:511185/1‑101 (MQ=255)
                                                    gtcggtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAt                                                  >  1:2748263‑M1/44‑101 (MQ=255)
                                                    aTCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt                                                  >  1:3138158/1‑101 (MQ=255)
                                                     tCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc                                                 <  1:3286763/101‑1 (MQ=255)
                                                       gcgcTCATCGACGGTGGGCGTCGGATGCGTTGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcgg                                               >  2:3829440/1‑101 (MQ=255)
                                                       gcgcTCATCGACGGTGGGCGTCGGATGCGTGGGGTCGCTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcgg                                               <  1:3993012/101‑1 (MQ=255)
                                                        gtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGatcggc                                              <  2:2768152‑M1/62‑2 (MQ=255)
                                                          gggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCGGTCGGGTCCGTGGGAtcggtcg                                            <  2:1871703‑M1/64‑1 (MQ=255)
                                                           tCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGAtcggtggg                                           >  2:2445464/1‑101 (MQ=255)
                                                             atcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGgatcggtcggat                                         >  1:538925‑M1/35‑99 (MQ=255)
                                                              tcggtcgggtccgtggggtcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                        >  2:1775029‑M1/34‑101 (MQ=255)
                                                              tCGGCGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                        <  1:1194/101‑1 (MQ=255)
                                                               cggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                                       >  1:2973884‑M1/33‑101 (MQ=255)
                                                               cggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                       >  1:2465749‑M1/33‑101 (MQ=255)
                                                                ggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTGGGGTCCg                                      <  2:1286208‑M1/70‑1 (MQ=255)
                                                                cgtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg                                      <  1:4407853‑M1/70‑1 (MQ=255)
                                                                 gtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                                     >  2:3726034‑M1/31‑101 (MQ=255)
                                                                    gggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTggg                                  >  1:4625878‑M1/28‑101 (MQ=255)
                                                                     ggggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGa                                 <  2:4357590/100‑1 (MQ=255)
                                                                         cgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcgg                             >  2:2558139‑M1/23‑101 (MQ=255)
                                                                           tgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCGGTCGGGGCCGTGGGgtcggtg                           <  2:1731227‑M1/81‑2 (MQ=255)
                                                                           tCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTATGTAGGATCGGTCGGGTCCGTCGGGtccgtg                           <  1:4150964/101‑2 (MQ=11)
                                                                            ggggtcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                          >  1:4354485‑M1/20‑101 (MQ=255)
                                                                            gggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                          <  1:1009700‑M1/82‑1 (MQ=255)
                                                                             ggatcggtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                         >  2:2143791‑M1/19‑101 (MQ=255)
                                                                             ggATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcgggcgg                         >  1:497183/1‑101 (MQ=255)
                                                                                tGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                      >  2:2626400/1‑101 (MQ=255)
                                                                                 cggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATcc                     >  1:384615‑M1/15‑101 (MQ=255)
                                                                                 cggtcggatccgtgGGATCGGTCGGGTCCGTGGGCTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAtcg                     <  1:3437475‑M1/87‑2 (MQ=255)
                                                                                  ggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCg                    <  1:2330757‑M1/88‑1 (MQ=255)
                                                                                  ggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCg                    <  1:2143941‑M1/88‑1 (MQ=255)
                                                                                  cGTGGGATCGGTCGGGTCGGTGGGATCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCg                    <  2:3087880/101‑1 (MQ=11)
                                                                                   gTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGt                   >  1:2903042/1‑101 (MQ=255)
                                                                                    tGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                  <  1:4466354/101‑1 (MQ=255)
                                                                                      ggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg                >  2:3991101/10‑101 (MQ=11)
                                                                                       tccgatctggatcggtgggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGa               <  1:2654211/93‑1 (MQ=11)
                                                                                       ggtccgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGATCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGGGCCGTGGGa               >  2:459549/6‑101 (MQ=11)
                                                                                        gtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAt              >  2:3013876‑M1/8‑101 (MQ=255)
                                                                                         tcggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc             >  1:3796690/1‑101 (MQ=255)
                                                                                         tcggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc             <  1:3347650/101‑1 (MQ=255)
                                                                                          cggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcg            >  2:356400/1‑101 (MQ=11)
                                                                                          ccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcg            >  1:240473‑M1/6‑101 (MQ=255)
                                                                                           cgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcgg           <  2:3937460/100‑1 (MQ=11)
                                                                                            gtgggatcggtcggggccgtgggatcggtcgggtccgtgGGATCTGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGtccgt          <  2:2356913/98‑1 (MQ=11)
                                                                                            gtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGCGGGATCGGTGGGGCCCGTGGGGTCGGTCGGGTGCCTCGGAtgggt          >  1:1371487/1‑101 (MQ=11)
                                                                                             tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc         <  1:1308284‑M1/99‑1 (MQ=255)
                                                                                             tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc         >  1:3608911‑M1/3‑101 (MQ=255)
                                                                                             tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc         <  1:1089336‑M1/99‑1 (MQ=255)
                                                                                             tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc         >  1:2187494‑M1/3‑101 (MQ=255)
                                                                                             tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGgatccgtggggtccgtg         >  2:2213516‑M1/3‑94 (MQ=255)
                                                                                             tagagtcggtcggatcgttgggatcggtcgtgttcgtggtgTCGGGCTGATCCTTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgttc         <  2:584171/96‑1 (MQ=2)
                                                                                             cgagctcggtcggctccgttcgatcggtcgggtccgtgggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc         <  2:4418090/96‑1 (MQ=11)
                                                                                              gGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg        >  1:2066170‑M1/2‑101 (MQ=255)
                                                                                              cggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGGCCGTGGGATCGGGCGGGTCCGTCGTAtcggtcg        >  2:801271/1‑101 (MQ=11)
                                                                                              cgattcggttggtgccgtgggatcggtcgggtccgtgggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg        <  2:328172/97‑1 (MQ=11)
                                                                                               ggatcggtcgggtccgtggggtccgtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAGCGGTCGGGTCCGTGGGAtcggtcgg       >  2:414947/1‑101 (MQ=11)
                                                                                               ggatcggtcgggtccgtgggatcggtcgggtccgtcGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtcgg       <  2:2690769/101‑1 (MQ=11)
                                                                                               ggatcggtcgggtccgtgggatcggtcgggtccgtcGGATCGGTGGGATCCGGGGGATCGGTCGGGTCCGTGGGAGCGGTCGGGTCCGTGGGAtcgggcgg       >  2:2066030/1‑101 (MQ=14)
                                                                                               ggatcggtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcggtcgg       <  2:1585429/101‑1 (MQ=11)
                                                                                                 gtgggtcgggtccgtgggatcggtcgggtccgtcggatcGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGt     <  1:3693485/98‑1 (MQ=11)
                                                                                                 atcggtcgggtcagtgggatgcgtcgggtcagtgggATCGGTCGTGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGt     <  2:3465081/101‑1 (MQ=11)
                                                                                                 atcggtcggatccgtgggatcggtcggatccgtgggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtcggat     >  1:1653854/1‑99 (MQ=11)
                                                                                                   cggtcgggtccgtggggtccgtcggatccgtgggatCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGcc   >  2:1657380/1‑101 (MQ=11)
                                                                                                   cggtcgggtccgtgggatcggtctggtccgtgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATcc   <  1:605364/101‑1 (MQ=11)
                                                                                                    |                                                                                                    
CGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCG  >  NC_014034/221861‑222061

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