Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,823,325:1 +G 17.5% intergenic (+62/+66) norW → / ← hypF NADH:flavorubredoxin oxidoreductase/carbamoyl phosphate phosphatase and maturation protein for [NiFe] hydrogenases

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,823,3251.G17.5% 175.0 / 45.8 57intergenic (+62/+66)norW/hypFNADH:flavorubredoxin oxidoreductase/carbamoyl phosphate phosphatase and maturation protein for [NiFe] hydrogenases
Reads supporting (aligned to +/- strand):  ref base . (20/27);  new base G (4/6);  total (24/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

GATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGT  >  Exported/1823257‑1823389
                                                                     |                                                                
gATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑Cg                                                                <  1:2565379/71‑1 (MQ=255)
gATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑Cg                                                                <  1:5339693/71‑1 (MQ=255)
 aTGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGc                                                               <  1:3050926/71‑1 (MQ=255)
 aTGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGc                                                               <  1:5402383/71‑1 (MQ=255)
   gTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCAt                                                             >  1:555788/1‑71 (MQ=255)
    tAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCAtg                                                            >  1:5614669/1‑71 (MQ=255)
      ggTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCAtgtg                                                          <  1:2626366/71‑1 (MQ=255)
       gTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATgtgt                                                         <  1:1008766/71‑1 (MQ=255)
       gTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATgtgt                                                         <  1:5467348/71‑1 (MQ=255)
          ggCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAgg                                                      >  1:3997628/1‑71 (MQ=255)
          ggCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAgg                                                      <  1:2687919/71‑1 (MQ=255)
               cTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGc                                                     >  1:3974633/1‑67 (MQ=255)
               cTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGa                                                 <  1:2086222/71‑1 (MQ=255)
                  tgCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGAt                                                <  1:60408/69‑1 (MQ=255)
                  tgCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATaa                                              <  1:2136262/71‑1 (MQ=255)
                   gCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAg                                             >  1:5442851/1‑71 (MQ=255)
                     cTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCAtgt                                                           <  1:4765372/55‑1 (MQ=255)
                     cTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGc                                           <  1:61392/71‑1 (MQ=255)
                     cTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGc                                           <  1:839346/71‑1 (MQ=255)
                       aaaaTGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATgtgt                                                         <  1:4665323/55‑1 (MQ=255)
                        gggtGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGc                                           <  1:3300383/66‑1 (MQ=18)
                        gggtGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCg                                          <  1:5224143/67‑1 (MQ=18)
                        aaaTGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGtt                                        <  1:2659720/71‑1 (MQ=255)
                        aaaTGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGtt                                        >  1:3545806/1‑71 (MQ=255)
                        aaaTGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGt                                         <  1:5583964/70‑1 (MQ=255)
                         ggtGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGtt                                        <  1:2093581/69‑1 (MQ=21)
                          aTGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTa                                      >  1:3678413/1‑71 (MQ=255)
                           tGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTAc                                     >  1:4195416/1‑71 (MQ=255)
                             tCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTAcg                                    <  1:226097/70‑1 (MQ=255)
                             tCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTAcg                                    <  1:1206975/70‑1 (MQ=255)
                              cGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTAcgcc                                  >  1:1009605/1‑71 (MQ=255)
                               ggATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTAcgcc                                  <  1:1886594/70‑1 (MQ=255)
                                gATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcgccg                                 <  1:3502571/71‑1 (MQ=18)
                                gATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcgccg                                 <  1:2220199/71‑1 (MQ=18)
                                 aTGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCg                                                   >  1:106473/1‑51 (MQ=255)
                                   gCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCAt                              <  1:342817/70‑1 (MQ=255)
                                    cGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATc                             <  1:5260542/70‑1 (MQ=255)
                                    cGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATc                             <  1:2446669/70‑1 (MQ=255)
                                     gACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATcc                            >  1:4586736/1‑71 (MQ=18)
                                                gTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTc                >  1:3734995/1‑71 (MQ=255)
                                                gTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTc                >  1:4468053/1‑71 (MQ=255)
                                                   ttATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAt                      >  1:3773046/1‑62 (MQ=255)
                                                   ttATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAt                      >  1:2486958/1‑62 (MQ=255)
                                                      tCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATg          >  1:1056081/1‑71 (MQ=255)
                                                       ccGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGgtgt                 >  1:381741/1‑63 (MQ=255)
                                                        cGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGtt                                        <  1:521248/40‑1 (MQ=17)
                                                        cGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGc         <  1:2103605/70‑1 (MQ=255)
                                                        cGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGc         <  1:1676696/70‑1 (MQ=255)
                                                        cGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGCGTCCAAATGc         <  1:150420/70‑1 (MQ=255)
                                                           ccTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTa                                      >  1:170489/1‑39 (MQ=14)
                                                            cTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAacac    >  1:3485102/1‑71 (MQ=255)
                                                              aCGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcgcc                                  >  1:2492745/1‑40 (MQ=11)
                                                              aCGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcgcc                                  >  1:617926/1‑40 (MQ=11)
                                                              aCGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGt  >  1:1425518/1‑71 (MQ=255)
                                                               cGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                         >  1:4126872/1‑47 (MQ=255)
                                                               cGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                         >  1:1035101/1‑47 (MQ=255)
                                                               cGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACg   <  1:4216816/69‑1 (MQ=255)
                                                                     |                                                                
GATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGA‑CGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGT  >  Exported/1823257‑1823389

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑