Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,412,228 C→G 7.1% C534S (TGT→TCT)  RCAP_rcc01318 ← outer membrane autotransporter barrel domain‑containing protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,412,2280C→G 7.1% 27.9 ‑3.5 ‑0.0210C534S (TGT→TCT) RCAP_rcc01318outer membrane autotransporter barrel domain‑containing protein
Reads supporting (aligned to +/- strand):  new base (4/11):  ref base (111/84):  total (115/95)
Fisher's exact test strand distribution p-value = 3.06e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.88e-01

Only 100 of 241 total aligned reads displayed.

CCCGCGATCGGCGCGGAGAAGGTCTTTACGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCGCCACGGTCGGTGCCGT  >  NC_014034/1412129‑1412328
                                                                                                   |                                                                                                    
  cgcgATCGGCGCGGAGAAGTTCTTTCCGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGAcaca                                                                                                   <  2:2934080/101‑1 (MQ=255)
     gATCGGCGCGGAGAAGGTCTTTACGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGc                                                                                                >  2:2637376/1‑101 (MQ=255)
     gATCGGCGCGGAGAAGGTCTTTACGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGc                                                                                                >  2:2629801/1‑101 (MQ=255)
     gATCGGCGCGGAGAAGGTCTTTACGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGc                                                                                                >  1:3214250/1‑101 (MQ=255)
      aTCGGCGCGGAGAAGGTCTTTACGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGcc                                                                                               >  2:4102667/1‑101 (MQ=255)
      aTCGGCGCGGAGAAGGTCTTTACGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGcc                                                                                               <  1:737007/101‑1 (MQ=255)
        gggCGGGGAGAAGGTCTTTCCGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGt                                                                                             <  2:2715008/100‑1 (MQ=255)
        gggCGCGGAAAAGGTCTTTCCGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGt                                                                                             <  1:1532568/100‑1 (MQ=255)
        cGGCGCGGAGAAGGTCTTTACGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGt                                                                                             >  2:2745141/1‑101 (MQ=255)
        cGGCGCGGAGAAGGTCTTTACGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGt                                                                                             >  2:1166621/1‑101 (MQ=255)
        cGGCGCGGAGAAGGTCTTTACGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGt                                                                                             <  1:1875131/101‑1 (MQ=255)
         ggCGCGGAGAAGGTCTTTACGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGCGCCGTCAACTGCGAAGGCGGTGCCCGCAAGACACAGACCGtc                                                                                            >  2:2168021/1‑101 (MQ=255)
          gcggGGAGAAGGTCTTTCCGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGtct                                                                                           <  1:337775/101‑1 (MQ=255)
            gcgGAGAAGGTCTTTACGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTcc                                                                                         >  1:2670619/1‑101 (MQ=255)
             gggAGAAGGTCTTTCCGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTccc                                                                                        <  1:2247316/100‑1 (MQ=255)
               gagaAGGTCTTTACGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCcgc                                                                                      >  2:1954855/1‑101 (MQ=255)
                agaAGGTCTTTACGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCcgcg                                                                                     >  1:67023/1‑101 (MQ=255)
                  aaGGCCTTTACGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGc                                                                                   >  2:449439/1‑101 (MQ=255)
                     gTCTTTCCGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATa                                                                                <  1:3373386/101‑1 (MQ=255)
                     gTCTTTACCCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATa                                                                                <  2:2393691/101‑1 (MQ=255)
                      tCTTTACGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTACCCGCAAGACACAGGCCGTCTCCCGCGGCATAg                                                                               >  1:3643496/1‑101 (MQ=255)
                       cTTTACGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAgg                                                                              >  1:63491/1‑101 (MQ=255)
                       cTTTACGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAgg                                                                              >  1:3097066/1‑101 (MQ=255)
                        tttACGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAggt                                                                             <  1:3411962/101‑1 (MQ=255)
                         ttgcgcccgccAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAggtg                                                                            <  1:4064210/98‑1 (MQ=255)
                           accccccccAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAggtggt                                                                          <  1:481545/98‑1 (MQ=255)
                             gcccgccAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGc                                                                        >  2:921603/1‑101 (MQ=255)
                             gcccgccAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGc                                                                        >  2:3790635/1‑101 (MQ=255)
                             ccccgccAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGc                                                                        <  2:3021215/100‑1 (MQ=255)
                                 gccAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCAc                                                                    <  2:3438513/101‑1 (MQ=255)
                                  ccAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCAcc                                                                   <  2:3554711/101‑1 (MQ=255)
                                   cAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACcg                                                                  <  2:1034321/101‑1 (MQ=255)
                                   cAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACcg                                                                  >  1:2383850/1‑101 (MQ=255)
                                   cAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACcg                                                                  >  2:2882221/1‑101 (MQ=255)
                                   cAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACcg                                                                  >  1:3488686/1‑101 (MQ=255)
                                      tgtCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCgcgc                                                               >  2:3980962/1‑101 (MQ=255)
                                       gtCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGcc                                                              >  2:3841905/1‑101 (MQ=255)
                                        tCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGccc                                                             <  2:526182/101‑1 (MQ=255)
                                        tCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGccc                                                             <  1:1195992/101‑1 (MQ=255)
                                          tGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCcgc                                                           >  1:1156274/1‑101 (MQ=255)
                                            aTCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCcgcgc                                                         >  1:3858490/1‑100 (MQ=255)
                                              cGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACt                                                       <  1:375746/101‑1 (MQ=255)
                                                  cGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTggg                                                   <  1:2705301/101‑1 (MQ=255)
                                                   gAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGc                                                  >  1:2284691/1‑101 (MQ=255)
                                                   gAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGc                                                  >  1:2929455/1‑101 (MQ=255)
                                                   gAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGc                                                  >  2:1799806/1‑101 (MQ=255)
                                                   gAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGc                                                  >  2:1588231/1‑101 (MQ=255)
                                                   gAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGc                                                  <  1:1249438/101‑1 (MQ=255)
                                                    aGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCg                                                 >  1:197555/1‑101 (MQ=255)
                                                      gTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTc                                               <  1:22438/101‑1 (MQ=255)
                                                      gTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTc                                               >  1:958435/1‑101 (MQ=255)
                                                        gcgcACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCaa                                             <  1:1415768/101‑1 (MQ=255)
                                                          gcACCACGGTGCGGTTAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCCTCAACa                                           <  2:592239/101‑1 (MQ=255)
                                                           caccacGGTGCCGTCAANTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACAt                                          <  2:1074407/101‑1 (MQ=255)
                                                               acGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGc                                      >  1:841989/1‑101 (MQ=255)
                                                               acGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGc                                      <  2:3926633/101‑1 (MQ=255)
                                                                 ggTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCt                                    <  2:204134/101‑1 (MQ=255)
                                                                  gTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCtt                                   <  1:1391236/101‑1 (MQ=255)
                                                                   ggCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTc                                  <  2:3107906/100‑1 (MQ=255)
                                                                    gCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGNCTTcn                                 >  1:917568/1‑100 (MQ=255)
                                                                    gCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCg                                 >  1:1035143/1‑101 (MQ=255)
                                                                      cGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAg                               >  2:1786099/1‑101 (MQ=255)
                                                                       gTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAgg                              >  1:212149/1‑101 (MQ=255)
                                                                       gTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAgg                              >  1:2915367/1‑101 (MQ=255)
                                                                        tCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGc                             <  1:3790892/101‑1 (MQ=255)
                                                                        tCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGc                             >  1:2759255/1‑101 (MQ=255)
                                                                          aaCTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCa                           >  2:1867180/1‑101 (MQ=255)
                                                                           aCTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCAc                          <  1:3260733/101‑1 (MQ=255)
                                                                           aCTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCAc                          >  1:2315466/1‑101 (MQ=255)
                                                                            cTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCAcc                         >  1:1953061/1‑101 (MQ=255)
                                                                             tGCGAAGGTGGTCCCCCCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCt                        <  1:2625217/101‑1 (MQ=255)
                                                                              gggAAAGTGGTGCCCGCAAGAGACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTg                       <  1:1091542/99‑1 (MQ=9)
                                                                              gCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTg                       >  2:4172160/1‑101 (MQ=255)
                                                                               cGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTgg                      >  1:584314/1‑101 (MQ=255)
                                                                               cGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTgg                      <  1:4102928/101‑1 (MQ=255)
                                                                               cGAAAGTGGTGCCCGCAAGAGACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTgg                      <  2:1179008/101‑1 (MQ=2)
                                                                                gAAAGTGGTGCCCGCAAGAGACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGc                     <  1:2514917/101‑1 (MQ=9)
                                                                                  aGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGccc                   >  1:2982517/1‑101 (MQ=255)
                                                                                   ggtggtGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCg                  >  2:3175679/1‑101 (MQ=40)
                                                                                   agtggtgCCCGCAAGAGACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCg                  <  2:2977851/100‑1 (MQ=18)
                                                                                    gtggtgCCCGCAAGAGACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCGc                 >  1:1813343/1‑101 (MQ=14)
                                                                                     tggtgCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCGcc                >  2:2311328/1‑101 (MQ=40)
                                                                                      ggtgCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCGCCa               <  1:101631/101‑1 (MQ=255)
                                                                                         gCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCGCCAcgg            <  2:1252831/101‑1 (MQ=255)
                                                                                         ccccGCAAGCCACAGCCCGTCTCCCGCGGCATGGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCGCCAcgg            <  2:94832/100‑1 (MQ=255)
                                                                                          cccGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCGCCAcggt           >  2:1428715/1‑101 (MQ=255)
                                                                                            cggaaGACACAGGCCGTTTCCCGCGGTATAGGTGGTGCCCACTGCACCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCCCCTGGCCCGCCAcggtcg         <  2:696159/98‑1 (MQ=17)
                                                                                            cGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCGCCAcggtcg         <  2:1482691/101‑1 (MQ=255)
                                                                                              cAAGACCCAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCGCCAcggtcggt       <  2:3690586/101‑1 (MQ=255)
                                                                                              cAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCGCCAcggtcggt       <  1:3044015/101‑1 (MQ=255)
                                                                                              cAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCCCCAcggtcggc       >  2:1091448/1‑100 (MQ=255)
                                                                                              cAAGACACAGGCCGCCCCCCGCGGCATAGGTGGTGCCCCCCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCGCCAcggtcggt       <  2:773089/101‑1 (MQ=255)
                                                                                               aaGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCGCCAcggtcgggg      >  1:796130/1‑99 (MQ=255)
                                                                                               aaGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTAAACATAGGCCTTCGAGGCCACCTGGCCCGCCACGGTCGGTg      >  2:169115/1‑101 (MQ=255)
                                                                                                aGACACAGGCCGCCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCGCCACGGTCGGTGc     <  2:2279982/101‑1 (MQ=255)
                                                                                                 gACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCCCCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCGCCACGGTCGGTGcc    <  1:1608981/101‑1 (MQ=255)
                                                                                                  acacAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCACGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCGGCACGGTCGGTGCCg   >  2:2806092/1‑101 (MQ=255)
                                                                                                   cacaGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCGCCACGGTCGGTGCCgg  >  2:3330343/1‑100 (MQ=255)
                                                                                                   cacaGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCGCCACGGTCGGTGCCGt  >  2:3647770/1‑101 (MQ=255)
                                                                                                   cacaGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCGCCACGGTCGGTGCCGt  <  2:212111/101‑1 (MQ=255)
                                                                                                   |                                                                                                    
CCCGCGATCGGCGCGGAGAAGGTCTTTACGCCCGCCAGTGTCTGATCGCCCGAGGTGCGCACCACGGTGCCGTCAACTGCGAAGGTGGTGCCCGCAAGACACAGGCCGTCTCCCGCGGCATAGGTGGTGCCCACCGCGCCCGCGACTTGGGCGTCAACATAGGCCTTCGAGGCCACCTGGCCCGCCACGGTCGGTGCCGT  >  NC_014034/1412129‑1412328

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG