Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,216 C→A 3.7% P366P (CCC→ACA)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,2160C→A 3.7% 9.4 ‑0.8 ‑0.0189P366P (CCC→CCA) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (2/5):  ref base (92/90):  total (94/95)
Fisher's exact test strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.94e-01

Only 100 of 206 total aligned reads displayed.

GGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTCCCCCCGCGG  >  NC_014034/1397116‑1397316
                                                                                                     |                                                                                                    
gggCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGccc                                                                                                      <  2:1926881/101‑1 (MQ=255)
    gTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCa                                                                                                  >  1:4064221/1‑101 (MQ=255)
     tGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCAc                                                                                                 <  1:1032626/101‑1 (MQ=255)
      ggATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCAcc                                                                                                >  1:1418733/1‑101 (MQ=255)
       gATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCa                                                                                               <  1:2632400/101‑1 (MQ=255)
       gATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCa                                                                                               <  2:2496035/101‑1 (MQ=255)
         tCGAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCGCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGa                                                                                             <  2:1450949/101‑1 (MQ=255)
          cGAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGAt                                                                                            <  2:240093/101‑1 (MQ=255)
          cGAAACGTCCAGCTCGCGCCACATGTCTTCGACCGCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGAt                                                                                            >  1:1417829/1‑101 (MQ=255)
           gAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGTCA.CGACGCCCTGCACCACATc                                                                                           >  2:2995101/1‑101 (MQ=255)
            aaaCGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCCGCACCAGATcc                                                                                          >  1:2015707/1‑101 (MQ=255)
            aaaCGTCCAGCCCGCGCCACATGTCTTCGCCCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATcc                                                                                          <  2:1001680/101‑1 (MQ=255)
             aaCGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCcg                                                                                         <  1:2384942/101‑1 (MQ=255)
             aaCGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCcg                                                                                         >  2:3297422/1‑101 (MQ=255)
             aaCGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCcg                                                                                         >  1:863998/1‑101 (MQ=255)
             aaCGTCCAGCTCGCGCCACATGTCTTCCACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCcg                                                                                         >  2:1751101/1‑101 (MQ=255)
              aCGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCcgc                                                                                        >  1:3424039/1‑101 (MQ=255)
                gTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCgcgc                                                                                      >  1:1523018/1‑101 (MQ=255)
                 tCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACAAGATCCGCGcc                                                                                     >  1:1960526/1‑101 (MQ=255)
                   cAGCTCGCGCCACATGTCTTCTACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGt                                                                                   >  1:3289558/1‑101 (MQ=255)
                   cAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGt                                                                                   <  2:4413818/101‑1 (MQ=255)
                    aGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTc                                                                                  <  2:113481/101‑1 (MQ=255)
                     gCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCg                                                                                 >  2:3330053/1‑101 (MQ=255)
                       tCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCt                                                                               <  2:1964656/101‑1 (MQ=255)
                        cgcgCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTc                                                                              >  2:243160/1‑101 (MQ=255)
                          cgcCACATGTCTTCGACCGTGCGCCGTCGCTCCCACCGCGGGGGCGTGGCTCGAAACGTGTTTGCCA.AGAGGCACTGCGCCAGATCCACGACGACGCTccc                                                                            >  2:1985085/1‑101 (MQ=255)
                           gcCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTcccc                                                                           >  1:1293050/1‑101 (MQ=255)
                           gcCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCCGCACCAGATCCGCGCCGTCGCTcccc                                                                           >  1:3777710/1‑101 (MQ=255)
                             cacaTGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTcccccc                                                                         >  1:713689/1‑101 (MQ=255)
                             cacaTGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTcccccc                                                                         <  2:2674251/101‑1 (MQ=255)
                              acaTGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCcg                                                                        <  2:1775542/101‑1 (MQ=255)
                               caTGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCcgc                                                                       >  2:3781167/1‑101 (MQ=255)
                                  gTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCggg                                                                    >  2:3322467/1‑101 (MQ=255)
                                   tCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCgggg                                                                   <  2:133359/101‑1 (MQ=255)
                                   tCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCgggg                                                                   >  2:2291532/1‑101 (MQ=255)
                                   tCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCgggg                                                                   <  2:863913/101‑1 (MQ=255)
                                     ttCGACCGTGCGCCGGCGCCCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGc                                                                 <  1:2657932/101‑1 (MQ=255)
                                        gACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTg                                                              >  1:350903/1‑101 (MQ=255)
                                          ccGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGa                                                            >  1:630929/1‑101 (MQ=255)
                                           cGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                           <  1:987221/101‑1 (MQ=255)
                                            gTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                          <  1:2999405/101‑1 (MQ=255)
                                            gTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                          >  2:2789395/1‑101 (MQ=255)
                                             tGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                         >  2:1081026/1‑101 (MQ=255)
                                              gcggCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGCTGCTGACACCGCCAACAAATATGTGGTCAACACTGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                         <  2:3261086/101‑1 (MQ=2)
                                              gcggCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGCTGCTGACACCGCCAACAAATATGTGGTCAACACTGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                         <  2:2240384/101‑1 (MQ=2)
                                              gcggCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGCTGCTGACACCGCCAACAAATATGTGGTCAACACTGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                         <  2:2035313/101‑1 (MQ=2)
                                                gcCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                      <  2:1629972/101‑1 (MQ=255)
                                                 ccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                                     <  1:1467277/101‑1 (MQ=255)
                                                  cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaca                                                    >  2:2159683/1‑101 (MQ=255)
                                                  cGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaca                                                    <  2:4170766/101‑1 (MQ=255)
                                                   gtggCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGTCA.CGCCGACCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGCAGGCGTGGATCGCAacac                                                   >  2:743493/4‑101 (MQ=37)
                                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACcgc                                                <  1:3146109/101‑1 (MQ=255)
                                                       cTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACcgcg                                               >  1:4146097/1‑101 (MQ=255)
                                                        tCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCgcgc                                              <  2:4302671/101‑1 (MQ=255)
                                                         ccccccGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCa                                             <  1:785766/101‑1 (MQ=255)
                                                         ccccccGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCa                                             >  2:182559/1‑101 (MQ=255)
                                                          cccccGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCAc                                            <  2:3135337/101‑1 (MQ=255)
                                                           ccccGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCAcg                                           <  1:1081114/101‑1 (MQ=255)
                                                            cccGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCAcgc                                          <  2:1542881/101‑1 (MQ=255)
                                                               gcgGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCAcgcgcg                                       >  2:165486/1‑101 (MQ=255)
                                                               gcgGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCAcgcgcg                                       <  2:3082264/101‑1 (MQ=255)
                                                                cgGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGa                                      >  1:1200440/1‑101 (MQ=255)
                                                                cgGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGa                                      <  1:504343/101‑1 (MQ=255)
                                                                 gggggCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGaa                                     >  1:1729268/1‑101 (MQ=255)
                                                                  ggggCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCCCGCGCGacc                                    >  2:3965155/1‑99 (MQ=255)
                                                                  ggggCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAAc                                    >  2:3543041/1‑101 (MQ=255)
                                                                  ggggCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAAc                                    >  2:693068/1‑101 (MQ=255)
                                                                  ggggCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAAc                                    <  2:776835/101‑1 (MQ=255)
                                                                      cGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCg                                >  1:3990799/1‑101 (MQ=255)
                                                                       gTGGATTGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCgc                               >  2:635504/1‑100 (MQ=255)
                                                                       gTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGt                               >  2:3469219/1‑101 (MQ=255)
                                                                       gTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGt                               >  1:4182542/1‑101 (MQ=255)
                                                                       gTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGt                               >  1:2758706/1‑101 (MQ=255)
                                                                        tGGTTCGAAACGTGTCTTCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTc                              <  2:3777457/101‑1 (MQ=255)
                                                                        tGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTc                              >  1:1431893/1‑101 (MQ=255)
                                                                         ggATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCa                             <  2:4020187/101‑1 (MQ=255)
                                                                         ggATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCa                             <  2:3405031/101‑1 (MQ=255)
                                                                             cGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGccc                         <  1:3703354/101‑1 (MQ=255)
                                                                              gAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGcccc                        >  2:3043805/1‑101 (MQ=255)
                                                                              gAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCCTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCGCGCACGCGCGAACGGCGTCAGcccc                        <  2:3661014/101‑1 (MQ=255)
                                                                              gAAACGTGGCTGCCA.CGACCCCCTGCCCCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGcccc                        <  2:4154507/101‑1 (MQ=255)
                                                                                  cGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGa                    <  1:424805/101‑1 (MQ=255)
                                                                                  cGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGa                    <  1:1856971/101‑1 (MQ=255)
                                                                                   gtgtCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGaa                   >  1:24532/1‑101 (MQ=255)
                                                                                    tgtCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAAc                  >  1:2579994/1‑101 (MQ=255)
                                                                                     gtCTGCCA.CGACGCCCTGCACCAGATCCGCACCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACg                 >  1:3975233/1‑101 (MQ=255)
                                                                                            a.cgacgCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTcc          <  2:1200349/101‑1 (MQ=255)
                                                                                             ccgacgCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTcc          <  2:1325117/100‑1 (MQ=255)
                                                                                              cgacgCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTccc         <  1:3601195/101‑1 (MQ=255)
                                                                                              cgacgCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTccc         <  2:2890915/101‑1 (MQ=255)
                                                                                              cgacgCCCTGCACCAGATCCGCGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTccc         <  1:2813430/101‑1 (MQ=255)
                                                                                                acgCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTccccc       >  2:3951028/1‑101 (MQ=255)
                                                                                                acgCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTccccc       >  1:2334297/1‑101 (MQ=255)
                                                                                                acgCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTccccc       >  1:3877527/1‑101 (MQ=255)
                                                                                                acgCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTccccc       <  1:1678168/101‑1 (MQ=255)
                                                                                                acgCCCTGCACCAGATACGCGCCGCCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGCACGGCGTCAGCCCCGCGAACGGGGCTcgccc       >  2:3505672/1‑101 (MQ=37)
                                                                                                 cgCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTcccccc      <  1:3229235/101‑1 (MQ=255)
                                                                                                 cgCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTcccccc      >  1:60461/1‑101 (MQ=255)
                                                                                                    ccTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTCCCCCcgcg   >  1:4412079/1‑101 (MQ=255)
                                                                                                     cTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTCCCCCCGCgg  >  2:2854809/1‑101 (MQ=255)
                                                                                                     |                                                                                                    
GGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACCGTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCA.CGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGCGCACGCGCGAACGGCGTCAGCCCCGCGAACGTCGCTCCCCCCGCGG  >  NC_014034/1397116‑1397316

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG