Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,948,851 +G 85.8% coding (1069/2472 nt) quiP ← acyl‑homoserine‑lactone acylase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,948,8511.→G85.8% 131.5 ‑0.7 ‑1.2288coding (1069/2472 nt)quiPacyl‑homoserine‑lactone acylase
Reads supporting (aligned to +/- strand):  new base (120/127):  ref base (17/24):  total (137/151)
Fisher's exact test strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.30e-02

Only 100 of 299 total aligned reads displayed.

TGATCCGACGGCAGCACCGGGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTC.G..CCC.A.CCCTCGGGGCTGCGGTA.CGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCGGCGAAGGCGGCGGTGATGC  >  NC_014034/1948751‑1948943
                                                                                                     |                                                                                                 
  aTCCGACGGCAGCACCGGGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTC.gg                                                                                                 >  2:30403/1‑100 (MQ=255)
    ccGACGGCAGCACCGGGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTC.G..gcc                                                                                             <  2:1982514/101‑1 (MQ=255)
     cGACGGCAGCACCGGGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTC.G..gccc                                                                                            <  2:3353065/101‑2 (MQ=255)
      gACGGCAGCACCGGGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTC.G..gccca                                                                                           >  1:1308876/1‑99 (MQ=255)
      gACGGCAGCACCGGGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTC.G..gccca                                                                                           <  2:5683616/101‑3 (MQ=255)
        cGGCAGCACCGGGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTC.G..gcccaac                                                                                         <  1:1917384/101‑5 (MQ=255)
        cGGCAGCACCGGGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTC.G..gcccaac                                                                                         <  2:5099741/101‑5 (MQ=255)
          gcagcaCCGTGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTC.G..gcccaaccc                                                                                       >  1:2403513/1‑95 (MQ=255)
          gcagcaCCGGGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTC.G..gcccaaccc                                                                                       >  2:867659/1‑95 (MQ=255)
          gcagcaCCGGGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTC.G..gcccaaccc                                                                                       >  1:1544621/1‑95 (MQ=255)
           cagcaCCGGGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTC.G..gcccaaccct                                                                                      <  2:3232119/101‑8 (MQ=255)
            agcaCCGGGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTC.G..gcccaaccctc                                                                                     >  1:3174209/1‑93 (MQ=255)
            agcaCCGGGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTC.G..gcccaaccctc                                                                                     <  1:303636/101‑9 (MQ=255)
            agcaCCGGGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTC.G..gcccaaccctc                                                                                     >  1:920223/1‑93 (MQ=255)
            agcaCCGGGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTC.G..gcccaaccctc                                                                                     <  1:3563148/101‑9 (MQ=255)
               aCCGGGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCggg                                                                                  >  2:2812775/1‑101 (MQ=255)
                 cGGGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGc                                                                                >  1:618434/1‑101 (MQ=255)
                 cGGGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGc                                                                                <  1:433847/101‑1 (MQ=255)
                  gggCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCt                                                                               >  2:1743544/1‑101 (MQ=255)
                  gggCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCt                                                                               >  2:5506611/1‑101 (MQ=255)
                    gCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGc                                                                             >  1:5224735/1‑101 (MQ=255)
                       gTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGt                                                                          >  2:1627546/1‑101 (MQ=255)
                        ttttCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTa                                                                         >  2:467417/1‑101 (MQ=255)
                         tttCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTag                                                                        >  1:5766663/1‑100 (MQ=255)
                          ttCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTagc                                                                       <  2:1587636/101‑3 (MQ=255)
                          ttCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTagc                                                                       >  1:442691/1‑99 (MQ=255)
                             ggTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGTGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGtggcgct                                                                    >  2:2293366/1‑95 (MQ=255)
                              gTCCAGCGCAGGATGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTc                                                                   >  1:4133983/1‑101 (MQ=255)
                               tCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCg                                                                  >  2:3567116/1‑101 (MQ=255)
                                 cAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCggg                                                                <  2:1276119/101‑1 (MQ=255)
                                  aGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGc                                                               >  1:3803656/1‑101 (MQ=255)
                                  aGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGc                                                               >  1:556197/1‑101 (MQ=255)
                                   gcgcAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCg                                                              >  1:2619221/1‑101 (MQ=255)
                                    cgcAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGCTCACCTTGCGGGTGGTGTACTCGG..CCCAA.CCCTCGGGGCGGCGGTAGCGCTCGGGCGt                                                             >  1:1157113/1‑101 (MQ=255)
                                     gcAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCCCCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGtt                                                            <  2:618376/101‑1 (MQ=255)
                                     gcAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGtt                                                            <  1:136716/101‑1 (MQ=255)
                                      cAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGttt                                                           >  2:43958/1‑101 (MQ=255)
                                      cAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGttt                                                           <  2:4435854/101‑1 (MQ=255)
                                       agggtgagggtgATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGtttt                                                          >  2:1479176/1‑101 (MQ=255)
                                       agggtgagggtgATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGtttt                                                          >  1:1609562/1‑101 (MQ=255)
                                       agggtgagggtgATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGtttt                                                          >  1:3881507/1‑101 (MQ=255)
                                        gggtgagggtgaTCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTc                                                         >  1:2423892/1‑101 (MQ=255)
                                        gggtgagggtgaTCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTc                                                         <  1:2544389/101‑1 (MQ=255)
                                            gagggtgaTCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGAt                                                     >  1:1198591/1‑101 (MQ=255)
                                               ggtgaTCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCaa                                                  >  2:3561557/1‑101 (MQ=255)
                                                gtgaTCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAg                                                 >  2:2997377/1‑101 (MQ=255)
                                                 tgaTCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGt                                                >  1:3591528/1‑101 (MQ=255)
                                                 tgaTCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGt                                                <  1:5348780/101‑1 (MQ=255)
                                                  gaTCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGtt                                               <  2:2901544/101‑1 (MQ=255)
                                                   aTCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTc                                              >  1:3666108/1‑101 (MQ=255)
                                                   aTCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTc                                              >  2:5562785/1‑101 (MQ=255)
                                                      ggCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTc                                           <  1:3331156/101‑1 (MQ=255)
                                                      ggCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTc                                           <  2:527441/101‑1 (MQ=255)
                                                      ggCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTc                                           >  2:5366133/1‑101 (MQ=255)
                                                       gcgcGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTcg                                          <  1:4734569/101‑1 (MQ=255)
                                                         gcgGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTcgac                                        >  1:4046991/1‑101 (MQ=255)
                                                             cATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAGCTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAAc                                    >  2:1083208/1‑101 (MQ=255)
                                                             cATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTTGGATCAAGTTCCTCGACGAAc                                    >  2:1332889/1‑101 (MQ=255)
                                                             cATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAAc                                    >  2:5664881/1‑101 (MQ=255)
                                                             cATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAAc                                    >  1:5648221/1‑101 (MQ=255)
                                                              aTCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACa                                   >  2:4294433/1‑101 (MQ=255)
                                                              aTCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACa                                   >  2:1012701/1‑101 (MQ=255)
                                                               tCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAg                                  >  2:1387491/1‑101 (MQ=255)
                                                               tCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAg                                  <  1:1895138/101‑1 (MQ=255)
                                                               tCCGTGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAg                                  <  2:4070852/101‑1 (MQ=255)
                                                                  ttGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTc                               >  2:5006761/1‑101 (MQ=255)
                                                                  ttGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGGAGCGCTCGGGGGTTTTCGGATCAAGTTCCTCGACGAACAGGTc                               >  1:4906957/1‑101 (MQ=255)
                                                                   tGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTcc                              >  2:2661977/1‑101 (MQ=255)
                                                                   tGACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCGCGGGGCTGCGGTAGGGCTCTGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTcc                              <  2:3817278/101‑1 (MQ=255)
                                                                    gACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCt                             >  2:2612706/1‑101 (MQ=255)
                                                                    gACCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCt                             <  2:1308788/101‑1 (MQ=255)
                                                                     aCCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTc                            <  1:1762316/101‑1 (MQ=255)
                                                                     aCCGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTc                            <  2:815078/101‑1 (MQ=255)
                                                                     aCCGTGACGATCAACTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTc                            <  2:4135203/101‑1 (MQ=255)
                                                                       cGTGACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGt                          >  2:2421802/1‑101 (MQ=255)
                                                                          gACGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCAt                       >  1:3358305/1‑101 (MQ=255)
                                                                           aCGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATc                      <  1:2743607/101‑1 (MQ=255)
                                                                           aCGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATc                      <  2:883633/101‑1 (MQ=255)
                                                                           aCGATCACCTTGCGGGTGGTGAACTC.GGACCA.A.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATc                      >  2:2099050/1‑101 (MQ=255)
                                                                            cGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCg                     <  2:3967863/101‑1 (MQ=255)
                                                                            cGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCg                     <  1:561259/101‑1 (MQ=255)
                                                                            cGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCg                     <  1:2812936/101‑1 (MQ=255)
                                                                            cGATCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCg                     <  2:4847649/101‑1 (MQ=255)
                                                                              aTCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCGGc                   >  2:3730638/1‑101 (MQ=255)
                                                                              aTCACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCGGc                   >  1:3457410/1‑101 (MQ=255)
                                                                              aTCACCTTGCGGGTGGGGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCGGc                   >  2:1508938/1‑101 (MQ=255)
                                                                                cACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCGGCGa                 <  1:2830660/101‑1 (MQ=255)
                                                                                cACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCGGCGa                 >  1:2176022/1‑101 (MQ=255)
                                                                                cACCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCGGCGa                 <  1:2064573/101‑1 (MQ=255)
                                                                                 aCCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCGGCGaa                <  1:4823430/101‑1 (MQ=255)
                                                                                 aCCTTGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCGGCGaa                <  2:1903031/101‑1 (MQ=255)
                                                                                     tGCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCGGCGAAggcg            >  2:4253225/1‑101 (MQ=255)
                                                                                      gCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCGGCGAAggcgg           <  2:5766368/101‑1 (MQ=255)
                                                                                      gCGGGTGGTGAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCGGCGAAggcgg           >  2:653713/1‑101 (MQ=255)
                                                                                           tggtgAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCGGCGAAGGCGGCGGTg      >  1:1441979/1‑101 (MQ=255)
                                                                                            ggtgAACTCGG..CCC.ACCCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCGGCGAAGGCGGCGGTGa     <  1:78404/101‑1 (MQ=255)
                                                                                             gtgAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCGGCGAAGGCGGCGGTGag    >  1:4916382/1‑100 (MQ=255)
                                                                                             gtgAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCGGCGAAGGCGGCGGTGAt    <  1:5007042/101‑1 (MQ=255)
                                                                                              tgAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCGGCGAAGGCGGCGGTGATg   <  1:2973916/101‑1 (MQ=255)
                                                                                               gAACTCGG..CCCAA.CCCTCGGGGCTGCGGTAGCGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCGGCGAAGGCGGCGGTGATGc  >  1:3084005/1‑101 (MQ=255)
                                                                                                     |                                                                                                 
TGATCCGACGGCAGCACCGGGCCGTTTTCGGTCCAGCGCAGGGTGAGGGTGATCGGCGCGGCATCCTTGACCGTGACGATCACCTTGCGGGTGGTGAACTC.G..CCC.A.CCCTCGGGGCTGCGGTA.CGCTCGGGCGTTTTCGGATCAAGTTCCTCGACGAACAGGTCCTCGTCATCGGCGAAGGCGGCGGTGATGC  >  NC_014034/1948751‑1948943

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 28 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG