Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,732,014 C→A R91L (CGC→CTC)  hypF carbamoyl phosphate phosphatase and maturation protein for [NiFe] hydrogenases

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,732,0140C→A100.0% 122.7 39R91L (CGC→CTC) hypFcarbamoyl phosphate phosphatase and maturation protein for [NiFe] hydrogenases
Reads supporting (aligned to +/- strand):  new base (17/22):  ref base (0/0):  total (17/22)

TTTGCGTATTCATGGTGCCGCCTGTGCTCTGGCGGATAGTGAACTCAGTGGGCAGTTGTGACCAG  >  REL606/2731982‑2732046
                                |                                
ttggCGTATTCATGGTGCCGCCTGTGCTCTGGAGGa                               <  1:898426/33‑1 (MQ=25)
 ttGCGTATTCGTGGTGCCGCCTGTGCTCTGGAGGAt                              <  1:5185498/36‑1 (MQ=25)
 ttGCGTATTCATGGTGCCGCCTGTGCTCTGGAGGAt                              <  1:439238/36‑1 (MQ=255)
 ttGCGTATTCATGGTGCCGCCTGTGCTCTGGAGGAt                              >  1:5314835/1‑36 (MQ=255)
  tGCGTATTCATGGTGCCGCCTGTGCTCTGGAGGATa                             <  1:5260053/36‑1 (MQ=255)
    cGTATTCATGGTGCCGCCTGTGCTCTGGAGGATAGt                           >  1:1831413/1‑36 (MQ=255)
      tATTCATGGTGCCGCCTGTGCTCTGGAGGATAGTGa                         <  1:4812545/36‑1 (MQ=255)
      tATTCATGGTGCCGCCTGTGCTCTGGAGGATAGTGa                         >  1:330075/1‑36 (MQ=255)
       aTTCATGGTGCCGCCTGTGCTCTGGAGGATAGTGaa                        >  1:1597807/1‑36 (MQ=255)
       aTTCATGGTGCCGCCTGTGCTCTGGAGGATAGTGaa                        >  1:2666628/1‑36 (MQ=255)
       aTTCATGGTGCCGCCTGTGCTCTGGAGGATAGTGaa                        <  1:354875/36‑1 (MQ=255)
          cATGGTGCCGCCTGTGCTCTGGAGGATAGTGAACTc                     <  1:7341578/36‑1 (MQ=255)
          cATGGTGCCGCCTGTGCTCTGGAGGATAGTGAACTc                     >  1:1737417/1‑36 (MQ=255)
          cATGGTGCCGCCTGTGCTCTGGAGGATAGTGAACTc                     >  1:1215684/1‑36 (MQ=255)
           aTGGTGCCGCCTGTGCTCTGGAGGATAGTGAACTCa                    >  1:4927001/1‑36 (MQ=255)
           aTGGTGCCGCCTGTGCTCTGGAGGATAGTGAACTCa                    <  1:5799695/36‑1 (MQ=255)
           aTGGTCCCGCCTGTGCTCTGGAGGATAGTGAACTCa                    <  1:5634224/36‑1 (MQ=25)
             ggTGCCGCCTGTGCTCTGGAGGATAGTGAACTCAGt                  <  1:1187813/36‑1 (MQ=255)
             aatGCCACCTATACTCTGAAGGATAGTGAACTCAGt                  <  1:956353/34‑1 (MQ=12)
               ggccgccTGTGCTCTGGAGGATAGTGAACTCAGTgg                <  1:507396/35‑1 (MQ=37)
                gccgccTGTGCTCTGGAGGATAGTGAACTCAGTggg               >  1:1096182/1‑36 (MQ=255)
                 ccggcTGTGCTCCGGAGGATAGTGAACTCAGTGGGc              <  1:4780856/36‑1 (MQ=21)
                   gccTGTGCTCTGGAGGATAGTGAACTCAGTGGGCAg            >  1:4746820/1‑36 (MQ=255)
                   gccTGTGCTCTGGAGGATAGTGAACTCAGTGGGCAg            >  1:7485611/1‑36 (MQ=255)
                     cTGTTCTCTGGAGTATAGTGAACTCAGTGGGCAGtt          <  1:1651585/36‑1 (MQ=21)
                     cTGTGCTCTGGAGGATAGTGAACTCAGTGGGCAGtt          >  1:2249605/1‑36 (MQ=255)
                      tgtgCTCTGGAGGATAGTGAACTCAGTGGGCAGTtg         >  1:2278198/1‑36 (MQ=255)
                       gtgCTCTGGAGGATAGTGAACTCAGTGGGCAGTtgt        <  1:948204/36‑1 (MQ=255)
                       gtgCTCTGGAGGATAGTGAACTCAGTGGGCAGTtgt        <  1:5947889/36‑1 (MQ=255)
                       gtgCTCTGGAGGATAGTGAACTCAGTGGGCAGTtgt        <  1:5324286/36‑1 (MQ=255)
                       gtgCTCTGGAGGATAGTGAACTCAGTGGGCAGTtgt        <  1:4891124/36‑1 (MQ=255)
                       gtgCTCTGGAGGATAGTGAACTCAGTGGGCAGCtgt        >  1:5642337/1‑36 (MQ=25)
                        tgCTCTGGAGGATAGTGAACTCAGTGGGCAGTtgtg       <  1:6485189/36‑1 (MQ=255)
                        tgCTCTGGAGGATAGTGAACTCAGTGGGCAGTtgtg       <  1:1806980/36‑1 (MQ=255)
                        tgCTCTGGAGGATAGTGAACTCAGTGGGCAGTtgtg       >  1:1088202/1‑36 (MQ=255)
                          ctctGGAGGATAGTGAACTCAGTGGGCAGTTGTGAc     >  1:20411/1‑36 (MQ=255)
                             tGGAGGATAGTGAACTCAGTGGGCAGTTGTGATCAg  >  1:5200831/1‑36 (MQ=25)
                             tGGAGGATAGTGAACTCAGTGGGCAGTTGTGACCAg  <  1:5438168/36‑1 (MQ=255)
                             tGGAGGATAGTGAACTCAGTGGGCAGTTGTGACCAg  <  1:865187/36‑1 (MQ=255)
                                |                                
TTTGCGTATTCATGGTGCCGCCTGTGCTCTGGCGGATAGTGAACTCAGTGGGCAGTTGTGACCAG  >  REL606/2731982‑2732046

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG