Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,608,028 G→T 100% intergenic (+298/‑49) RCAP_rcc01484 → / → murD hypothetical protein/UDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,608,0280G→T100.0% 125.7 37intergenic (+298/‑49)RCAP_rcc01484/murDhypothetical protein/UDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase
Reads supporting (aligned to +/- strand):  new base (7/30):  ref base (0/0):  total (7/30)

ACCGGCGCCTGTCTTCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCCTTGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGGCCGCGCGGCGCGGGCCCAATTCACCCCCGCAGCTTACGGA  >  NC_014034/1607929‑1608068
                                                                                                   |                                        
gcggccgcCTGTC.TCCGCCCTTTTCCCAGCGGAGCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCcg                                       <  2:3655546/96‑1 (MQ=37)
   gccgcCTGTC.TCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCcgctg                                    <  2:3559641/99‑3 (MQ=38)
   gccgcCTGTC.TCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCcgctg                                    >  1:1398691/3‑99 (MQ=38)
   gccgcCTGTC.TCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCcgctg                                    <  1:329847/99‑3 (MQ=38)
          tgtctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCgc.tggcgcgg                             >  2:3297272/5‑101 (MQ=38)
           gtctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCCCTTGGTCCgc.tggcgcggc                            >  1:2448230/4‑100 (MQ=255)
            tctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCGCCCTCTTGGTCCgc.tggcgcggcc                           >  2:1392734/3‑99 (MQ=255)
             ctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcc                          <  2:2165764/100‑1 (MQ=38)
             ctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcc                          <  1:2377146/100‑1 (MQ=38)
             ctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcc                          >  2:181299/2‑101 (MQ=38)
             ctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAA.GGGGGCCGCGCCGCCCCCCCTTGGGTC.CGCTGGGCGGGccc                         >  2:211589/2‑101 (MQ=255)
              tCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcca                         <  2:266780/101‑2 (MQ=38)
              tCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcca                         <  2:1218787/101‑2 (MQ=38)
               ccGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCaa                       <  2:4017710/101‑1 (MQ=38)
                 gACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAAtt                     <  1:3330323/101‑1 (MQ=37)
                 gACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAAtt                     <  1:1044462/101‑1 (MQ=37)
                  aacTTTTCCCCGCCGCTCCCGCGGAGGGGGGCGCCGCC.TGCGAAGTCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGTCGC.GGCCCAATTc                    <  2:802030/99‑1 (MQ=12)
                  aCCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTc                    <  1:4333203/101‑1 (MQ=37)
                  aCCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTc                    >  2:4428933/1‑101 (MQ=38)
                  aCCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTc                    <  1:1583231/101‑1 (MQ=38)
                  aCCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTc                    <  2:1713126/101‑1 (MQ=37)
                   ccTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCCCCattca                    >  1:2391175/1‑97 (MQ=37)
                   ccTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCa                   <  1:1579310/101‑1 (MQ=38)
                    cTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAc                  <  2:17201/101‑1 (MQ=37)
                    cTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAc                  >  2:2046518/1‑101 (MQ=38)
                      tttCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccc                <  1:4429183/101‑1 (MQ=37)
                       ttCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAcccc               <  2:4616081/101‑1 (MQ=38)
                        tCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCCccccc              >  2:4523784/1‑101 (MQ=37)
                        tCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccccc              >  1:2689636/1‑101 (MQ=38)
                         ccgaGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCg             <  2:2195735/98‑1 (MQ=37)
                         cccAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCg             <  2:1673001/101‑1 (MQ=38)
                          ccAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGc            <  1:2600988/101‑1 (MQ=38)
                          ccAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.CTGGGC.GGCCCCATTCACCCCCGc            >  2:1739097/1‑101 (MQ=37)
                          ccAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCCCTTCGGCCGC.CGGGGC.GGCCCCATTCAccccccc            >  2:3393685/1‑99 (MQ=255)
                           cAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcg           <  1:4257637/101‑2 (MQ=38)
                            aGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgg          <  2:2351171/101‑3 (MQ=37)
                            aGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACACCCGcgg          >  2:1005174/1‑99 (MQ=37)
                              cGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcggga        <  2:1903421/101‑5 (MQ=37)
                               ggATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggat       <  1:4507243/101‑6 (MQ=37)
                               ggATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggat       <  1:3547005/101‑6 (MQ=37)
                                gATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggata      <  1:1056619/101‑7 (MQ=37)
                                 aTCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCCCCCCCGcgggatat     <  1:1270878/101‑8 (MQ=25)
                                  tCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatt    <  1:12134/101‑9 (MQ=37)
                                   cccGCGACGGGGTGCGCCGCC.CGCGGCGCCCCCCGCAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatattt   <  1:754508/101‑10 (MQ=19)
                                    ccGGGACGGGGTGCGACGCC.TGCGGCGCCGCCCTAAAGGGGGCCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatttc  <  1:2458652/101‑11 (MQ=255)
                                                                                                   |                                        
ACCGGCGCCTGTCTTCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCCTTGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGGCCGCGCGGCGCGGGCCCAATTCACCCCCGCAGCTTACGGA  >  NC_014034/1607929‑1608068

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG