Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,731,801 Δ1 bp 2.4% coding (1526/2757 nt) RCAP_rcc02540 ← diguanylate cyclase/phosphodiesterase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,731,8010G→. 2.4% 7.3 ‑0.4 ‑0.3210coding (1526/2757 nt)RCAP_rcc02540diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  new base (2/3):  ref base (117/88):  total (119/91)
Fisher's exact test strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.83e-01

Only 100 of 219 total aligned reads displayed.

AGATCACCCGGTCGCCCAGCTGCGGGCCGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCGGGGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGCGTCGCGGCATGGCTGGC  >  NC_014034/2731706‑2731895
                                                                                               |                                                                                              
aGATCACCCGGTCGCCCAGCTGCGGGCCGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCgcgc                                                                                           >  2:422060/1‑101 (MQ=255)
aGATCACCCGGTCGCCCAGCTGCGGGCCGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCgcgc                                                                                           >  1:2456673/1‑101 (MQ=255)
 gATCACCCGGTCGCCCAGCTGCGGGCCGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCgcgcg                                                                                          >  2:2277062/1‑101 (MQ=255)
  aTCACCCGGTCGCCCAGCTGCGGGCCGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCgcgcgc                                                                                         >  1:1837291/1‑101 (MQ=255)
     aCCCGGTCGCCCAGCTGCGGGCCGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCgcgcgcgcg                                                                                      >  2:585159/1‑101 (MQ=255)
       ccGGTCGCCCAGCTGCGGGCCGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCa                                                                                    >  2:3201139/1‑101 (MQ=255)
        cGGTCGCCCAGCTGCGGGCCGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAg                                                                                   <  2:3736543/101‑1 (MQ=255)
          gTCGCCCAGCTGCGGGCCGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCt                                                                                 <  1:2082190/101‑1 (MQ=255)
          gTCGCCCAGCCGCGGGGCGAAACGGTCCTTGATCTGCTGGACCTGGTCGATCGACAGGAAGACCTGCGGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCt                                                                                 <  2:364873/101‑1 (MQ=255)
           tCGCCCAGCTGCGGGCCGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTg                                                                                <  1:4552266/101‑1 (MQ=255)
             gCCCAGCTGCGGGCCGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTgcg                                                                              <  2:1131572/101‑1 (MQ=255)
              cccAGCTGCGGGCCGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTgcgc                                                                             <  2:4442766/101‑1 (MQ=255)
              cccAGCTGCGGGCCGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGC..CGCGCGCGCGCAGCGGCGCCg                                                                           >  1:1363368/1‑101 (MQ=255)
                cggcggCGGGCCGAAACGGTCCTTGATGTGCTGGACCTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCg                                                                           <  2:3620121/96‑1 (MQ=255)
                cAGCTGCGGGCCGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCg                                                                           <  1:2189090/101‑1 (MQ=255)
                      cGGGCCGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGa                                                                     >  1:4161887/1‑101 (MQ=255)
                       gggCCGAAACGGTCCTTGATCTGCTGGAACTGGTGGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGAc                                                                    >  1:3856177/1‑101 (MQ=255)
                       gggCCGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGAc                                                                    <  1:1296242/101‑1 (MQ=255)
                       gggCCGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGAc                                                                    >  2:3836792/1‑101 (MQ=255)
                       gggCCGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGAc                                                                    <  2:388597/101‑1 (MQ=255)
                          ccGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTg                                                                 >  1:4497084/1‑101 (MQ=255)
                           cGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTgg                                                                <  1:2370382/101‑1 (MQ=255)
                           cGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTgg                                                                >  2:1419708/1‑101 (MQ=255)
                           cGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTgg                                                                >  2:4447135/1‑101 (MQ=255)
                            gAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTggg                                                               <  1:100425/101‑1 (MQ=255)
                             aaaCGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGgc                                                              >  1:1448421/1‑101 (MQ=255)
                             aaaCGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTCCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGgc                                                              <  2:3439484/101‑1 (MQ=255)
                              aaCGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGgcg                                                             <  2:1980067/101‑1 (MQ=255)
                               aCGGTCCTTGATCTGCTGGAACTGTGCGATCGACAGGAAGACCTTCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGgcgc                                                            <  2:3858379/101‑1 (MQ=255)
                               aCGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGgcgg                                                            >  1:4387931/1‑100 (MQ=255)
                                cGGTCCTTTATCTGCTGGAACGGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGcgcg                                                           <  2:3809526/101‑1 (MQ=255)
                                cGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGcgcg                                                           >  1:2746565/1‑101 (MQ=255)
                                   tCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGaa                                                        <  2:1837165/101‑1 (MQ=255)
                                   tCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGaa                                                        >  2:1492638/1‑101 (MQ=255)
                                   tCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGATGCGCCGGGCCGACCTGGGCGCGGaa                                                        >  1:3470353/1‑101 (MQ=255)
                                    ccTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAg                                                       <  2:2061908/101‑1 (MQ=255)
                                    ccTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAg                                                       >  2:535216/1‑101 (MQ=255)
                                       tGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGTAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGAt                                                    >  1:4017226/1‑101 (MQ=255)
                                        gATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGAtt                                                   <  2:1714934/101‑1 (MQ=255)
                                        gATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGAtt                                                   <  2:1776808/101‑1 (MQ=255)
                                         aTCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTc                                                  >  1:742887/1‑101 (MQ=255)
                                          tCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTcg                                                 <  1:3776726/101‑1 (MQ=255)
                                          tCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTcg                                                 >  1:2633669/1‑101 (MQ=255)
                                               tgGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGAtt                                            <  2:419360/101‑1 (MQ=255)
                                                ggAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTg                                           >  2:2651874/1‑101 (MQ=255)
                                                 gAACTGGTCGATTGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTgg                                          <  2:1530812/101‑1 (MQ=255)
                                                 gAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTgg                                          >  1:3039253/1‑101 (MQ=255)
                                                 gAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTgg                                          >  1:120240/1‑101 (MQ=255)
                                                    cTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCa                                       <  1:615553/101‑1 (MQ=255)
                                                       gTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGcc                                    >  1:670849/1‑101 (MQ=255)
                                                        tcgatcgaCAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGccc                                   <  2:4016941/101‑1 (MQ=255)
                                                        tcgatcgaCAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGccc                                   <  2:4486108/101‑1 (MQ=255)
                                                        tcgatcgaAAGGAAGACCGGCTGCAAGAGATCGCCCGCCGCTCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGccc                                   <  2:869416/101‑1 (MQ=255)
                                                         cgatcgaCAGGAAGGCCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCg                                  <  2:798810/101‑1 (MQ=255)
                                                         cgatcgaCAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCg                                  <  1:290719/101‑1 (MQ=255)
                                                         cgatcgaCAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCg                                  >  1:2864375/1‑101 (MQ=255)
                                                          gatcgaCAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCgg                                 >  2:2795951/1‑101 (MQ=255)
                                                          gatcgaCAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCgg                                 <  2:4356692/101‑1 (MQ=255)
                                                           atcgaCAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGt                                <  2:3307944/101‑1 (MQ=255)
                                                           atcgaCAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGt                                >  1:3323027/1‑101 (MQ=255)
                                                            tcgaCAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTc                               >  1:204307/1‑101 (MQ=255)
                                                            tcgaCAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTc                               <  1:4517385/101‑1 (MQ=255)
                                                              gaCAGGAAGACCTGCTGCAAGAGATCGCCCGC..CGCGCGCGCGCACATGCGCCCGGCCGAGCTGTGGGCGGAGGGGGTCGCGCGTTGGCAGccccgggacag                           >  2:2978716/1‑93 (MQ=255)
                                                               aCAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAgc                            >  2:2594037/1‑101 (MQ=255)
                                                               aCAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAgc                            >  1:3690151/1‑101 (MQ=255)
                                                               aCAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAgc                            >  1:4079238/1‑101 (MQ=255)
                                                                cAGGAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAgcg                           >  1:3564415/1‑101 (MQ=255)
                                                                   gAAGCCCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGcat                        <  1:1419803/101‑1 (MQ=255)
                                                                   gAAGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGcat                        <  2:2527174/101‑1 (MQ=255)
                                                                    aaGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGcatc                       >  2:887535/1‑101 (MQ=255)
                                                                     atgcttgctgcAAGAGACCGCCCACCGCGCGCGCGCGGCGCTGCGCCGGGCCGTCCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGcatca                      <  1:4453684/96‑1 (MQ=255)
                                                                     aGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGcatca                      >  1:1741715/1‑101 (MQ=255)
                                                                     aGACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGcatca                      >  2:3206887/1‑101 (MQ=255)
                                                                      gACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGcatcat                     >  2:4149318/1‑101 (MQ=255)
                                                                      gACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGcatcat                     >  2:1058352/1‑101 (MQ=255)
                                                                      gACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGcatcat                     <  1:2511983/101‑1 (MQ=255)
                                                                      gACCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGcatcat                     <  1:1542666/101‑1 (MQ=255)
                                                                       aCCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATg                    <  2:1448326/101‑1 (MQ=255)
                                                                       aCCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATg                    >  2:1381945/1‑101 (MQ=255)
                                                                       aCCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATg                    <  2:3391875/101‑1 (MQ=255)
                                                                       aCCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATg                    >  2:4101206/1‑101 (MQ=255)
                                                                       aCCTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATg                    <  1:3256183/101‑1 (MQ=255)
                                                                        ccTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGc                   <  1:1945499/101‑1 (MQ=255)
                                                                        ccTGCTGCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGc                   >  1:3315997/1‑101 (MQ=255)
                                                                         ctgctgCAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGCg                  >  2:488251/1‑101 (MQ=255)
                                                                          tgctgcAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGCGt                 >  1:3050380/1‑101 (MQ=255)
                                                                           gctgcAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGCGTc                <  1:2652059/101‑1 (MQ=255)
                                                                            ctgcAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGCGTcg               <  2:118945/101‑1 (MQ=255)
                                                                            ctgcAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGCGTcg               >  1:3000527/1‑101 (MQ=255)
                                                                               cAAGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGCGTCGCgg            >  1:4663461/1‑101 (MQ=255)
                                                                                aaGAGATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGCGTCGCGGc           >  2:1477793/1‑101 (MQ=255)
                                                                                 agagATCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGCGTCGCGGCa          <  2:265822/101‑1 (MQ=255)
                                                                                  gcgaTCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGCGTCGCGGCAt         <  1:1075145/99‑1 (MQ=255)
                                                                                   agaTCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGCGTCGCGGCAtg        <  1:3411293/101‑1 (MQ=255)
                                                                                   agaTCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGCGTCGCGGCAtg        >  1:4227203/1‑101 (MQ=255)
                                                                                     agcccccGCCGCGTGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGCGTCGCGGCAtggc      <  1:3119531/97‑1 (MQ=255)
                                                                                     aTCGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGCGTCGCGGCAtggc      >  2:1821328/1‑101 (MQ=255)
                                                                                       cGCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGCGTCGCGGCAtggctg    <  2:2434474/101‑1 (MQ=255)
                                                                                        gCCCGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGCGTCGCGGCAtggctgg   >  2:3463418/1‑101 (MQ=255)
                                                                                         cccGCCGCGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGCGTCGCGGCAtggctggc  >  2:4529468/1‑101 (MQ=255)
                                                                                               |                                                                                              
AGATCACCCGGTCGCCCAGCTGCGGGCCGAAACGGTCCTTGATCTGCTGGAACTGGTCGATCGACAGGAAGACCTGCTGCAAGAGATCGCCCGCGGGGCGCGCGCGCAGCTGCGCCGGGCCGACCTGGGCGCGGAAGGATTCGCGATTGGGCAGCCCGGTCAGCGCATCATGCGTCGCGGCATGGCTGGC  >  NC_014034/2731706‑2731895

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG