Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,255,868 +C 100% intergenic (+38/+204) RCAP_rcc01184 → / ← RCAP_rcc01185 CsbD family protein/methyl‑accepting chemotaxis sensory transducer

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,255,8681.→C100.0% 208.6 96intergenic (+38/+204)RCAP_rcc01184/RCAP_rcc01185CsbD family protein/methyl‑accepting chemotaxis sensory transducer
Reads supporting (aligned to +/- strand):  new base (44/52):  ref base (0/0):  total (44/52)

AAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGG....CTG.TCTCGGTCCGCCCGAAGAGGGGGGGCGCTG..TGGTAC..CCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCT  >  NC_014034/1255795‑1255963
                                                                               |                                                                                                   
aagaagGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAgggg                                                                               >  1:506884/1‑101 (MQ=255)
aagaagGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAgggg                                                                               >  2:2303470/1‑101 (MQ=255)
 agaagGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAggggg                                                                              >  2:1879727/1‑101 (MQ=255)
 agaagGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAggggg                                                                              >  1:2646011/1‑101 (MQ=255)
 agaagGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAggggg                                                                              <  2:3240971/101‑1 (MQ=255)
 agaagGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAggggg                                                                              >  2:848342/1‑101 (MQ=255)
  gaaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAgggggg                                                                             <  1:5232645/101‑1 (MQ=255)
  gaaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAgggggg                                                                             >  2:3828121/1‑101 (MQ=255)
  gaaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAgggggg                                                                             <  2:5499919/101‑1 (MQ=255)
  gaaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAgggggg                                                                             >  1:2664914/1‑101 (MQ=255)
   aaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGCAGAggggggg                                                                            >  2:2049770/1‑101 (MQ=255)
   aaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGCGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAggggggg                                                                            >  2:3683367/1‑101 (MQ=255)
    aggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAgggggggg                                                                           >  2:4166920/1‑100 (MQ=255)
    aggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGgc                                                                           <  2:663765/101‑1 (MQ=255)
    aggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGgc                                                                           >  1:758552/1‑101 (MQ=255)
    aggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGgc                                                                           <  2:1228894/101‑1 (MQ=255)
     ggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAggggggggg                                                                          >  1:5699906/1‑99 (MQ=255)
     ggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGgcg                                                                          <  1:4435851/101‑1 (MQ=255)
     ggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGgcg                                                                          <  1:1367121/101‑1 (MQ=255)
     ggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGgcg                                                                          >  2:2024387/1‑101 (MQ=255)
     ggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGgcg                                                                          <  2:407162/101‑1 (MQ=255)
      gaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGgcgc                                                                         >  2:1395929/1‑101 (MQ=255)
      gaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGgcgc                                                                         >  1:3600663/1‑101 (MQ=255)
       aagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAggggggggggg                                                                        >  2:4859956/1‑97 (MQ=255)
       aagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCt                                                                        <  1:445958/101‑1 (MQ=255)
       aagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCt                                                                        <  2:2964381/101‑1 (MQ=255)
       aagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGCCCGCCCGAAGAGGGGGGgcgcg                                                                        >  2:2693739/1‑100 (MQ=255)
        agTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtg                                                                       <  2:4119630/101‑1 (MQ=255)
          tCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgcc                                                                     <  1:3220583/101‑3 (MQ=255)
           cgacgaCTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccc                                                                    >  1:2476317/1‑98 (MQ=255)
           cgacgaCTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccc                                                                    <  2:2749667/101‑4 (MQ=255)
            gacgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgcccc                                                                   >  2:2715573/1‑97 (MQ=255)
            gacgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgcccc                                                                   <  2:506834/101‑5 (MQ=255)
            gacgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgcccc                                                                   >  1:4298038/1‑97 (MQ=255)
            gacgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGCGTGGCGCtgggcc                                                                   >  2:3596371/1‑97 (MQ=255)
             acgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAggggggggggggccccc                                                                  >  1:5602073/1‑91 (MQ=255)
             acgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccccc                                                                  <  1:4754261/101‑6 (MQ=255)
             acgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccccc                                                                  <  1:1587635/101‑6 (MQ=255)
             acgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccccc                                                                  <  1:2145570/101‑6 (MQ=255)
             acgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccccc                                                                  <  1:347492/101‑6 (MQ=255)
             accacTTCCTGTCGAACCTCTGCCCCGGGGGCCCCGGACCTGACCCGGGGGCCTT.CGGCCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccccc                                                                  <  1:2158050/98‑6 (MQ=255)
              cgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgcccccc                                                                 <  2:3196023/101‑7 (MQ=255)
              cgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgcccccc                                                                 <  2:3472167/101‑7 (MQ=255)
              cgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAACAGGGGGGGCGCtg..tccccc                                                               >  1:3149055/1‑96 (MQ=255)
               gacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccg                                                                <  1:2715722/101‑8 (MQ=255)
               gacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccg                                                                <  2:1974540/101‑8 (MQ=255)
                acTTTCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgt                                                               >  1:1189125/1‑93 (MQ=255)
                acTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgt                                                               <  1:1507166/101‑9 (MQ=255)
                  ttCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgtcc                                                             <  1:78467/101‑11 (MQ=255)
                  ttCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgtcc                                                             <  2:3505840/101‑11 (MQ=255)
                  ttCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgtcc                                                             >  1:470930/1‑91 (MQ=255)
                                 tGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGC..CCCTGCGGGACTccccc                                            >  2:5192652/1‑101 (MQ=9)
                                 tGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTccccc                                            >  2:3490795/1‑101 (MQ=9)
                                  gACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTcccccc                                           >  1:1865103/1‑101 (MQ=9)
                                  gACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTcccccc                                           >  2:4676503/1‑101 (MQ=9)
                                  gACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTcccccc                                           <  2:865800/101‑1 (MQ=9)
                                  gACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTcccccc                                           <  2:2446682/101‑1 (MQ=9)
                                   aCCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCg                                          <  2:192256/101‑1 (MQ=9)
                                   aCCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCg                                          >  2:4056481/1‑101 (MQ=9)
                                    ccTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCgg                                         <  2:5637525/101‑1 (MQ=2)
                                    ccTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCgg                                         >  1:1128540/1‑101 (MQ=2)
                                    accaGGGGTGCCCGGCCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCgg                                         <  1:2554583/98‑1 (MQ=0)
                                     cTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCggg                                        <  1:5186645/101‑1 (MQ=2)
                                     cTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCggg                                        <  1:4482837/101‑1 (MQ=2)
                                     cTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCggg                                        >  1:53408/1‑101 (MQ=2)
                                      tAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGGCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGa                                       >  1:2751710/1‑101 (MQ=2)
                                      tAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGa                                       <  1:2389759/101‑1 (MQ=2)
                                      tAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGa                                       >  2:612622/1‑101 (MQ=2)
                                      tAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGa                                       >  2:1695010/1‑101 (MQ=2)
                                      tAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGa                                       <  1:627467/101‑1 (MQ=2)
                                       aGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGat                                      <  2:5601781/101‑1 (MQ=2)
                                       aGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGat                                      >  2:4891922/1‑101 (MQ=2)
                                       aGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGat                                      >  2:4818232/1‑101 (MQ=2)
                                       aGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGat                                      <  2:5734155/101‑1 (MQ=2)
                                       aGGGGTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGat                                      >  1:2666213/1‑101 (MQ=2)
                                       aGGGGCGCCGGACCGGAGCGGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGat                                      <  2:5601459/101‑1 (MQ=0)
                                        ggggTGCCGGACCGGACCCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGata                                     <  2:2885611/101‑1 (MQ=0)
                                        ggggTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGata                                     <  1:3379839/101‑1 (MQ=2)
                                        ggggTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGata                                     >  2:4974777/1‑101 (MQ=2)
                                         gggTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGatat                                    <  1:1753513/101‑1 (MQ=2)
                                          tgTGCCGGTCCGGAACCCGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATAtt                                   <  2:5193154/100‑1 (MQ=0)
                                          ggggCCGGACCCGAACCGGGGGCCTTCCGGCCC.CCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATAtt                                   <  2:4242191/98‑1 (MQ=0)
                                          ggTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATAtt                                   <  2:1524391/101‑1 (MQ=2)
                                          ggTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATAtt                                   <  2:472069/101‑1 (MQ=2)
                                          ggTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCCTCCCTGCGGGACTCCCCCCGGCATCtt                                   >  2:1825287/1‑101 (MQ=2)
                                           gTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATAttt                                  >  1:994104/1‑101 (MQ=0)
                                           gTGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATAttt                                  >  1:2438331/1‑101 (MQ=0)
                                            tGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTc                                 <  2:4488980/101‑1 (MQ=2)
                                            tGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTc                                 >  1:3320902/1‑101 (MQ=2)
                                            tGCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTc                                 <  1:1039885/101‑1 (MQ=2)
                                             gCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTcc                                <  1:4565612/101‑1 (MQ=2)
                                             gCCGGACCGGAACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTcc                                <  2:868961/101‑1 (MQ=2)
                                              ccggaccggaCCCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCg                               <  2:4951220/101‑1 (MQ=2)
                                              ccggaccggaACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCg                               <  2:5361121/101‑1 (MQ=2)
                                               cggaccggaACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGc                              >  2:3483931/1‑101 (MQ=2)
                                               cggaccggaACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGc                              <  2:2970772/101‑1 (MQ=2)
                                                ggaccggaACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGcc                             >  1:1099252/1‑101 (MQ=2)
                                                 gaccggaACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCa                            <  1:4167170/101‑1 (MQ=2)
                                                  accggaACCGGGGGCCTT.CGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTACCTGCGGGACTCCCCCCGGGATATTTCCGCCaa                           >  1:4999388/1‑101 (MQ=2)
                                                                           ncNN.TNNCGGTCCGCCCGAAGAGGGGGGGCGCGG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCt  <  1:2067738/100‑1 (MQ=18)
                                                                               |                                                                                                   
AAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTT.CGG....CTG.TCTCGGTCCGCCCGAAGAGGGGGGGCGCTG..TGGTAC..CCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCT  >  NC_014034/1255795‑1255963

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