Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,218 G→A 3.1% C367Y (TGC→AAT)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,2180G→A 3.1% 7.5 ‑2.3 ‑0.6191C367Y (TGC→TAC) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (1/5):  ref base (99/86):  total (100/91)
Fisher's exact test strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.71e-01

Only 100 of 198 total aligned reads displayed.

GCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTCCCCCCGCGGGA  >  NC_014034/1397118‑1397318
                                                                                                     |                                                                                                    
 cGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGc                                                                                                     <  2:642108/101‑1 (MQ=255)
  gTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCa                                                                                                    >  1:2207567/1‑101 (MQ=255)
   tGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCGGCAc                                                                                                   <  2:3581947/101‑1 (MQ=255)
    ggATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCAcc                                                                                                  >  1:2814914/1‑101 (MQ=255)
    ggATCGAAACGTCCAGCTCGCCCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCAcc                                                                                                  <  1:586659/101‑1 (MQ=255)
        cGAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGAt                                                                                              >  1:3264612/1‑101 (MQ=255)
         gAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATc                                                                                             <  1:4208288/101‑1 (MQ=255)
         gAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATc                                                                                             <  1:3528004/101‑1 (MQ=255)
         gAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATc                                                                                             <  2:2053825/101‑1 (MQ=255)
         gAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATc                                                                                             <  1:1016275/101‑1 (MQ=255)
          aaaCGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATcc                                                                                            >  1:841627/1‑101 (MQ=255)
           aaCGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCcg                                                                                           <  2:1525554/101‑1 (MQ=255)
            aCGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCcgc                                                                                          <  1:1213022/101‑1 (MQ=255)
            aCGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCcgc                                                                                          >  2:1138894/1‑101 (MQ=255)
             cGTCCAGCTCGCTCCACATGTCGTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCcgcg                                                                                         <  2:2972746/101‑1 (MQ=255)
             cGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCcgcg                                                                                         >  2:962205/1‑101 (MQ=255)
             cGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCcgcg                                                                                         >  1:2091556/1‑101 (MQ=255)
              gTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCgcgc                                                                                        >  2:867784/1‑101 (MQ=255)
               tCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGCGTCTGCCA.CGACGCCCTGCACCAGATCCGCGcc                                                                                       >  1:3626133/1‑101 (MQ=255)
                ccAGCTCGCGCCACATGTCTTCGACCGTTCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCg                                                                                      <  2:367840/101‑1 (MQ=255)
                 cAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGt                                                                                     >  2:4201465/1‑101 (MQ=255)
                  aGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTc                                                                                    <  1:3357130/101‑1 (MQ=255)
                    cTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGc                                                                                  >  1:3892997/1‑101 (MQ=255)
                     tCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCt                                                                                 <  2:3779821/101‑1 (MQ=255)
                     tCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCt                                                                                 >  1:881533/1‑101 (MQ=255)
                      cgcgCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTc                                                                                <  1:297699/101‑1 (MQ=255)
                        cgcCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTccc                                                                              <  2:822904/101‑1 (MQ=255)
                        cgcCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTccc                                                                              >  2:3646698/1‑101 (MQ=255)
                         gcCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTcccc                                                                             <  1:4114419/101‑1 (MQ=255)
                          ccACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCCccccc                                                                            >  2:1022265/1‑101 (MQ=255)
                           cacaTGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTcccccc                                                                           >  2:645756/1‑101 (MQ=255)
                           cacaTGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTcccccc                                                                           >  1:2768632/1‑101 (MQ=255)
                             caTGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCcgc                                                                         >  2:1729174/1‑101 (MQ=255)
                               tGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCgg                                                                       <  2:3974602/101‑1 (MQ=255)
                                gTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCggg                                                                      <  1:1474182/101‑1 (MQ=255)
                                gTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCggg                                                                      >  1:2912488/1‑101 (MQ=255)
                                gTCTTCGACCGTGCGCCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCggg                                                                      <  1:2588503/101‑1 (MQ=255)
                                 tCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCANCAGATCCGCGCCGTCGCTCCCCCCGCgggg                                                                     >  1:3554254/1‑101 (MQ=255)
                                   ttCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGc                                                                   >  1:1392241/1‑101 (MQ=255)
                                   ttCGACCGTGCGCCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGc                                                                   <  1:3910400/101‑1 (MQ=255)
                                     cgcccGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGt                                                                 <  2:2207408/98‑1 (MQ=255)
                                      gACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTg                                                                <  2:3448630/101‑1 (MQ=255)
                                         cGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                             >  1:1009997/1‑101 (MQ=255)
                                          gTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                            >  2:1369687/1‑101 (MQ=255)
                                            gcggCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGCTGCTGACACCGCCAACAAATATGTGGTCAACACTGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                           <  1:1891777/101‑1 (MQ=2)
                                              ggcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGCTGCTGACACCGCCAACAAATATGTGGTCAACACTGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                         >  1:909279/3‑101 (MQ=2)
                                              gcCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                        >  1:2140038/1‑101 (MQ=255)
                                              gcCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                        >  2:487236/1‑101 (MQ=255)
                                              gcCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                        <  2:2814695/101‑1 (MQ=255)
                                                cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGCTGCTGACACCGCCAACAAATATGTGGTCAACACTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                                       <  1:2944038/101‑1 (MQ=2)
                                                cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGCTGCTGACACCGCCAACAAATATGTGGTCAACACTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                                       <  1:3083776/101‑1 (MQ=2)
                                                 gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacac                                                     <  1:2213481/101‑1 (MQ=255)
                                                    gCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACcgc                                                  <  2:3308181/101‑1 (MQ=255)
                                                    gCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACcgc                                                  >  2:3213392/1‑101 (MQ=255)
                                                    gCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACcgc                                                  >  1:703082/1‑101 (MQ=255)
                                                     cTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACcgcg                                                 >  2:4120235/1‑101 (MQ=255)
                                                      tCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCgcgc                                                >  2:2446535/1‑101 (MQ=255)
                                                       ccccccGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCa                                               <  1:3780327/101‑1 (MQ=255)
                                                        cccccGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGTGTGGATCGAAACACCGCGCAc                                              <  1:2407784/101‑1 (MQ=255)
                                                          cccGCGGGGGCGTGGATGGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCAcgc                                            <  2:714966/101‑1 (MQ=255)
                                                           ccGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCAcgcg                                           <  2:3873027/101‑1 (MQ=255)
                                                           ccGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCAcgcg                                           <  2:568093/101‑1 (MQ=255)
                                                             gcgGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCAcgcgcg                                         >  1:1678005/1‑101 (MQ=255)
                                                               gggggCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGaa                                       <  2:1607414/101‑1 (MQ=255)
                                                                ggggCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGacc                                      >  1:2737433/1‑99 (MQ=255)
                                                                ggggCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAAc                                      <  2:3534794/101‑1 (MQ=255)
                                                                  ggCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACgg                                    >  1:995983/1‑101 (MQ=255)
                                                                  ggCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACgg                                    <  1:3507646/101‑1 (MQ=255)
                                                                    cGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCg                                  >  2:2764776/1‑101 (MQ=255)
                                                                      tGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTc                                >  2:1663271/1‑101 (MQ=255)
                                                                      tGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTc                                >  1:1794138/1‑101 (MQ=255)
                                                                      tGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTc                                <  2:3625880/101‑1 (MQ=255)
                                                                         aTCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGCCAGc                             >  1:3170031/1‑101 (MQ=255)
                                                                           cGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGccc                           <  2:3101629/101‑1 (MQ=255)
                                                                           cGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGccc                           >  2:2663651/1‑101 (MQ=255)
                                                                           cGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGccc                           <  1:1754107/101‑1 (MQ=255)
                                                                             aaaCGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCcg                         >  2:3895157/1‑101 (MQ=255)
                                                                             aaaCGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCcg                         >  2:3445895/1‑101 (MQ=255)
                                                                             aaaCGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCcg                         >  1:2351815/1‑101 (MQ=255)
                                                                                 gtgtCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGaa                     <  1:645833/101‑1 (MQ=255)
                                                                                  tgtCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAAc                    <  1:110535/101‑1 (MQ=255)
                                                                                    tCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGt                  <  2:3355945/101‑1 (MQ=255)
                                                                                    tCTCCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGt                  <  2:1757326/101‑1 (MQ=255)
                                                                                     cTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTc                 <  2:1948196/101‑1 (MQ=255)
                                                                                      tGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCg                >  2:354782/1‑101 (MQ=255)
                                                                                      tGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCg                >  1:1628187/1‑101 (MQ=255)
                                                                                       gCCA.CGTCGCCCTGCACCAGATCCTCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGc               >  2:130167/1‑101 (MQ=255)
                                                                                       gCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGCCAGCCCCGCGAACGTCGc               >  2:1349242/1‑101 (MQ=255)
                                                                                         cA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTc             >  1:1274605/1‑101 (MQ=255)
                                                                                            cgacgCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTccc           >  2:1546243/1‑101 (MQ=255)
                                                                                             gacgCCCTGCCCCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTcccc          <  1:767941/101‑1 (MQ=255)
                                                                                             gacgCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTcccc          <  1:506061/101‑1 (MQ=255)
                                                                                              acgCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTccccc         >  2:3703755/1‑101 (MQ=255)
                                                                                              acgCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTccccc         <  1:1544525/101‑1 (MQ=255)
                                                                                               cgCCCTGCACCAGATCCGCGCCGTCGCTCCCCCGGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTcccccc        <  1:3642058/101‑1 (MQ=255)
                                                                                               cgCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTcccccc        <  2:3264371/101‑1 (MQ=255)
                                                                                                gCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTCCCCCcg       >  2:2589408/1‑101 (MQ=255)
                                                                                                   cTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTCCCCCCGCgg    >  1:3959681/1‑101 (MQ=255)
                                                                                                    tGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTCCCCCCGCggg   >  1:1551899/1‑101 (MQ=255)
                                                                                                     gCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTCCCCCCGCGGGa  <  2:2497818/101‑1 (MQ=255)
                                                                                                     |                                                                                                    
GCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTCCCCCCGCGGGA  >  NC_014034/1397118‑1397318

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