Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 136,941 +G 2.7% coding (608/2064 nt) RCAP_rcc00130 → radical SAM family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034136,9411.→G 2.7% 7.6 ‑4.1 ‑0.2188coding (608/2064 nt)RCAP_rcc00130radical SAM family protein
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (79/104):  total (84/104)
Fisher's exact test strand distribution p-value = 1.66e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.43e-01

Only 100 of 227 total aligned reads displayed.

TGGTCGGACAAGGTCCGCCGCTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCACCGAACTGCCCGCCGATGATCTG  >  NC_014034/136841‑137041
                                                                                                     |                                                                                                       
 ggTCGGGCAAGGTCCGCCGCTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.g                                                                                                        >  2:1995372/1‑101 (MQ=255)
 ggTCGGACAAGGTCCGCCGCTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.g                                                                                                        >  2:964415/1‑101 (MQ=255)
 ggTCGGACAAGGTCCGCCGCTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.g                                                                                                        >  1:4278927/1‑101 (MQ=255)
     ggACAAGGTCCGCCGCTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.ggc.gg                                                                                                   >  1:4977966/1‑101 (MQ=255)
     ggACAAGGTCCGCCGCTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.ggc.gg                                                                                                   <  1:604466/101‑1 (MQ=255)
       aCAAGGTCCGCCGCTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.ggg.gggg                                                                                                 >  2:1608233/1‑99 (MQ=255)
       aCAAGGTCCGCCGCTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.Ggc.ggcg                                                                                                 <  2:3175309/101‑1 (MQ=255)
        cAAGGTCCGCCGGTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGc.ggggg                                                                                                >  1:2737589/1‑101 (MQ=255)
        cAAGGTCCGCCGCTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGg.ggcgg                                                                                                >  2:3748486/1‑101 (MQ=255)
         aaGGTCCGCCGCTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.Ggc.ggcgtg                                                                                               >  2:5000526/1‑99 (MQ=255)
          aGGTCCGCCGCTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGG.ggggggg                                                                                              >  2:1967311/1‑100 (MQ=255)
           ggTCCGCCGCTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.ggcgggcg                                                                                             <  2:3943012/101‑1 (MQ=255)
              ccgccgCTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.ggcggggggaa                                                                                          >  1:2240628/1‑96 (MQ=255)
              ccgccgCTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGaaa                                                                                          <  2:4908528/101‑1 (MQ=255)
                gccgcTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGG.GGGGGGGGGAAGc                                                                                        >  1:4101735/1‑101 (MQ=255)
                   gcTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTCCGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.TGCGGGGGGGGGAAAGccc                                                                                      >  1:1997804/1‑101 (MQ=255)
                   gcTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGG.ggggggggggagaccc                                                                                     >  1:4401600/1‑91 (MQ=255)
                   gcTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGcccc                                                                                     <  1:3223220/101‑1 (MQ=255)
                     tCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGcccccc                                                                                   >  2:2208179/1‑100 (MQ=255)
                     tCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCcg                                                                                   <  2:13373/101‑1 (MQ=255)
                     tCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCcg                                                                                   <  2:4051391/101‑1 (MQ=255)
                     tCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCcg                                                                                   <  1:4655181/101‑1 (MQ=255)
                      cGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGCGGGGGGGGGGAAGcccccc                                                                                   >  1:100022/1‑100 (MQ=255)
                      cGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGccccccc                                                                                  >  1:3093902/1‑99 (MQ=255)
                        aTCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCTGGGG.GGGGGGGGGAAGCCCCcccg                                                                                 >  1:2936587/1‑101 (MQ=255)
                         tCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGAt                                                                               <  1:959831/101‑1 (MQ=255)
                           cTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGG.GGGCGGGGCAAGACCCcccgcgcg                                                                             >  1:2969640/1‑97 (MQ=255)
                            tGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCTGGGG.GGGGGGGGAAAGccccccccgccg                                                                             >  2:1059671/1‑94 (MQ=255)
                            tGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGG.GGCGGGGGGAAGccacccccgctcg                                                                            >  2:3067586/1‑93 (MQ=255)
                                cGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGGGAAAGCCCCCGCGATCTGGCCg                                                                        >  1:3124745/1‑101 (MQ=255)
                                  aTGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGat                                                                      <  1:3825724/101‑1 (MQ=255)
                                    gCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGata.t                                                                   >  2:1511861/1‑101 (MQ=255)
                                     ccAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGGGGGGGGAAGGCCCCGCGAACTGGGCGAGA.Tc                                                                  >  1:1656712/1‑101 (MQ=255)
                                     ccAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.Tc                                                                  <  1:2856419/101‑1 (MQ=255)
                                      gatctcGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.Tcc                                                                 <  2:260310/96‑1 (MQ=255)
                                      cAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.Tcc                                                                 <  2:898011/101‑1 (MQ=255)
                                       acggCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCcg                                                                <  1:1206966/99‑1 (MQ=255)
                                       aaGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCTAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCcg                                                                <  2:4373981/101‑1 (MQ=255)
                                       aaGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCcg                                                                <  2:2797971/101‑1 (MQ=255)
                                        aGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGGCCCCGCGATCTGGCCGATA.TCcgc                                                               >  2:312085/1‑101 (MQ=255)
                                        aGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCcgc                                                               <  2:4895669/101‑1 (MQ=255)
                                          gCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGGAAGGCCCCCCCGTCTGGCCGATA.TTCgcgg                                                             >  2:1345926/1‑101 (MQ=255)
                                           cGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCTGGGG.GGGGGGGGAAAGCCCCCCCGCCCTGCCCGAAA.CCCgggg                                                             >  2:3838001/1‑101 (MQ=255)
                                           cGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGGGAAAACCCCCGCGATCTGGCCGATA.TCCccggc                                                            >  1:123995/1‑101 (MQ=255)
                                           cGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.CGCGGGGGGAAGCCCCCGCGGTCTTGCCGagacccagatc                                                             >  1:260344/1‑91 (MQ=255)
                                           cGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCAGCT.GGC.GGCGGGCGAAAGCCCCCCCGACCTGGCCGATA.Tcccccgc                                                            >  2:1360962/1‑96 (MQ=255)
                                           cGGATCTTTTGCTTTACGGCAATGCCGAACACGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCgcggc                                                            <  1:4586334/101‑1 (MQ=255)
                                              aTCTTTTGCTTTACGGCAATGCCGAACGCGCCGCGGTCGAGGTGGCCAACCGGCG.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTg                                                        <  2:2270604/101‑1 (MQ=255)
                                               tCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTgg                                                       <  2:2175884/101‑1 (MQ=255)
                                                cTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGGGGAAACCCCCCCGCTCTTGCCGATA.TCCgcgggggggg                                                       >  1:4192782/1‑95 (MQ=255)
                                                  tttGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGccc                                                    <  1:5000856/101‑1 (MQ=255)
                                                    tGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGGAAGCCCCCGCCCTCTAGGCGCTA.TCCGCGGC.GGGGccccc                                                  >  2:5049175/1‑99 (MQ=255)
                                                      cTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGG.GGGGGGGGAAAGCCCCCCCCCGCTGGCCGAGA.TCCGCGGCGGTGGCCCCGt                                                 >  2:531877/1‑101 (MQ=255)
                                                      cTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTTtt                                                >  2:1466829/1‑101 (MQ=255)
                                                      cTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGtt                                                <  2:1777016/101‑1 (MQ=255)
                                                       tttACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGGGGAAGACCCCCCCGTCTGGCCCACA.TACGCCAG.GAGGCCCCGTTc                                               >  1:3383859/1‑101 (MQ=255)
                                                        ttACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGGGAAAAACCCCGCGATCTGGCCGGAA.TCCCCGGC.GGGGCCCCGGTcc                                              >  1:1311399/1‑101 (MQ=255)
                                                        ttACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTcc                                              <  2:2417381/101‑1 (MQ=255)
                                                        ttACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTcc                                              <  1:1062113/101‑1 (MQ=255)
                                                         tACGGCAATTCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCg                                             <  2:4698118/101‑1 (MQ=255)
                                                         tACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGGGAAAAGCCCCCCGCTCTGGGCGATA.TCCCCGGG.GGGGGCCCGTTccc                                             >  1:4321593/1‑100 (MQ=255)
                                                          aCGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAACCCCCCCGATCTCGCCCATA.TCCCCGGC.GGGGCCCCGTTCCGc                                            >  1:2132510/1‑101 (MQ=255)
                                                           cGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGG.GGCGGGGGAAAGACCCCGCCAGCGGGGCGAGA.TTGGCCGC.GGGGCCCTCTTCCCcc                                           >  1:4972361/1‑101 (MQ=255)
                                                            ggCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGGGGGGGGAAGACCCCCCGGTCTGGCCGGTA.TCCCCGGG.GGGGGCCTGTTCCGccc                                          >  1:1157830/1‑101 (MQ=255)
                                                            ggCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGGGGAAGCCCCCCCGATATGGCCGATA.TCCCCGGC.GGGGGCCCGTTCCGccc                                          >  1:2489598/1‑101 (MQ=255)
                                                             gCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGG.GGGGGGGGGAAACCCCCCCGACCCGGGCGGTA.TCCCCGGC.GGGGCCCCGTTCCCcccc                                         >  1:3445232/1‑100 (MQ=255)
                                                               aaTGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGt                                       <  2:1451960/101‑1 (MQ=255)
                                                                aTGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTg                                      >  2:340233/1‑101 (MQ=255)
                                                                  gCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGG.GGGGGGGGAAAAACCCCGCGATCTGGGCGATA.TCCGCGGC.GTGGGCCTGTTTCGCCCGGGGcc                                    >  2:4828880/1‑101 (MQ=255)
                                                                  gCCGAACGCGCCGTGGTCGAGGGGGCCAACCGGCT.GGC.GGGGGGCGGAAACCCCCCCGGCCTTGCCGGTATTCCGGGGG.GGGGCCCGTTTCCCCCCgcggc                                     >  2:2378108/1‑97 (MQ=255)
                                                                   ccGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGccc                                   <  1:3905526/101‑1 (MQ=255)
                                                                    cGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCg                                  <  2:4132246/101‑1 (MQ=255)
                                                                    cGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCg                                  <  1:4889403/101‑1 (MQ=255)
                                                                         gcgcCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGc                             <  2:795453/101‑1 (MQ=255)
                                                                           gcCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCtt                           >  1:1432430/1‑101 (MQ=255)
                                                                              gTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGGGCAAGGCCCCCCCGTCTGGCCGCTA.TCCCCCGC.GGGGGCCCCTTCCCCCCCGCGCCCCCGGGcgtgcc                        >  1:3540271/1‑96 (MQ=255)
                                                                              gTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCAc                        <  1:3975384/101‑1 (MQ=255)
                                                                               tGGTCGAGGTGGCCAACCGGCT.GGC.GGGGGGGGAAAACCCCCCCGCTCTGGCCGGTA.TCCCCGGC.GGGGGCCCGTTTCGCCCCGGGCCCCAGGGGTTCAcc                       >  2:3429794/1‑101 (MQ=255)
                                                                               tGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCAcc                       <  2:4581418/101‑1 (MQ=255)
                                                                               gggTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCAcc                       <  1:3614932/100‑1 (MQ=255)
                                                                                ggTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCACCg                      <  1:2340895/101‑1 (MQ=255)
                                                                                ggTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCACCg                      >  1:1891224/1‑101 (MQ=255)
                                                                                 gTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCACCGa                     <  2:1067588/101‑1 (MQ=255)
                                                                                   cGAGGTGGCCAACCGGCT.GGC.GGGGGGCGGAAGCCCCCCCGGACTGGGCCCTA.CCCCCGGC.GGGGCGCCGTTTCGCCCCGGGGCCGCGGGCGTCTCCGCAc                   >  2:1827936/1‑101 (MQ=255)
                                                                                   cGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCACCGAAc                   <  2:59583/101‑1 (MQ=255)
                                                                                   cGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCACCGAAc                   <  2:4591777/101‑1 (MQ=255)
                                                                                   cGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAACCCCCGCGATCTGGCCGATA.TCCGCGGC.GGGGCCCCGTTTCGCCCCGTGCCCCAGGGCTTCACCGaaa                   >  2:1036755/1‑100 (MQ=255)
                                                                                    gAGGTGGCCAACCGGCT.GGC.GGCGGGGGGAAGTCCCCCCGCAACGGCCCCTA.TCCGCGGG.CTGGGCCCTTTCCTCCCCCGGCCCCAGGGGGTCACCCAACt                  >  1:1878108/1‑101 (MQ=255)
                                                                                    gAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCACCGAACt                  <  1:71704/101‑1 (MQ=255)
                                                                                        tGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCACCGAACTGccc              <  1:2690154/101‑1 (MQ=255)
                                                                                            cAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCACCGAACTGCccgccg          <  2:3683146/101‑1 (MQ=255)
                                                                                            cAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCACCGAACTGCccgccg          <  2:1270438/101‑1 (MQ=255)
                                                                                            cAACCGGCT.GGC.GGCGGGCGAAAACCCCCGCGATCTGGCCGATA.TCCGCGGC.GGGGCCCCGTTCCGCCCGGGGCCCGAGGGCTTCACCGAACTGCccgccg          >  2:2884354/1‑101 (MQ=255)
                                                                                               ccGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCACCGAACTGCCCGCCgatg       <  1:4524698/101‑1 (MQ=255)
                                                                                                cGGCT.GGG.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCACCGAACTGCCCGCCgatga      <  2:1734283/101‑1 (MQ=255)
                                                                                                cGGCT.GGC.GGGGGGGGAAAGCCCCCCCGATCTGGCCGATA.TCCCCGGC.GGGGGCCTGTTCCGCCCCGGGCCCGAGGGCGTCCCCCAACCGCCCGCCcatga      >  2:4412517/1‑101 (MQ=255)
                                                                                                cGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCACCGAACTGCCCGCCgacga      <  1:2883686/101‑1 (MQ=255)
                                                                                                   cT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCACCGAACTGCCCGCCGATGATCt   <  1:531956/101‑1 (MQ=255)
                                                                                                   cT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCACCGAACTGCCCGCCGATGATCt   <  1:2560525/101‑1 (MQ=255)
                                                                                                    t.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCACCGAACTGCCCGCCGATGATCTg  <  2:4101418/101‑1 (MQ=255)
                                                                                                     |                                                                                                       
TGGTCGGACAAGGTCCGCCGCTCGATCCTGGCCGATGCCAAGGCGGATCTTTTGCTTTACGGCAATGCCGAACGCGCCGTGGTCGAGGTGGCCAACCGGCT.GGC.GGCGGGCGAAAGCCCCCGCGATCTGGCCGATA.TCCGCGGC.GTGGCCCTGTTCCGCCCGGTGCCCGAGGGCTTCACCGAACTGCCCGCCGATGATCTG  >  NC_014034/136841‑137041

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