Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,824,288 G→A 54.5% D452D (GAC→GAT hypF ← 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,824,2880GA54.5% 20.8 / 120.2 88D452D (GAC→GAThypFcarbamoyl phosphate phosphatase and maturation protein for [NiFe] hydrogenases
Reads supporting (aligned to +/- strand):  ref base G (27/13);  new base A (18/30);  total (45/43)
Fisher's exact test for biased strand distribution p-value = 5.94e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

AGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCG  >  Exported/1824223‑1824348
                                                                 |                                                            
aGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCTCTGGCTGGAGACATAGCCCGGATGCGCATCATg                                                         <  1:939830/71‑1 (MQ=255)
aGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATg                                                         <  1:1306062/71‑1 (MQ=255)
aGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATg                                                         <  1:1352859/71‑1 (MQ=255)
        ttGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGca                                                       >  1:2544455/1‑65 (MQ=255)
         tGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTa                                                >  1:1222878/1‑71 (MQ=255)
         tGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTa                                                >  1:1649974/1‑71 (MQ=255)
          gCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTa                                                <  1:3291666/70‑1 (MQ=255)
             tCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCTTGCACAAc                                                   <  1:654133/64‑1 (MQ=255)
               ggCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCgg                                          >  1:2481099/1‑71 (MQ=255)
                 cAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCgggg                                        >  1:743865/1‑71 (MQ=255)
                 cAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCACGGATGCGCATCATGCACAACGTAt                                               >  1:2724064/1‑64 (MQ=255)
                   gATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGGACAACGTATTGCGGGGTa                                      <  1:2595544/71‑1 (MQ=255)
                    aTTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGt                                       <  1:3083935/69‑1 (MQ=255)
                       cATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTAt                                               <  1:2112227/58‑1 (MQ=255)
                        aTTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATc                                 <  1:2357117/71‑1 (MQ=255)
                        aTTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGAGGGGTAAAATc                                 <  1:526034/71‑1 (MQ=255)
                         tttCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCa                                                           <  1:297810/44‑1 (MQ=255)
                         tttCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTg                                             >  1:1024323/1‑58 (MQ=255)
                         tttCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTg                                             >  1:764486/1‑58 (MQ=255)
                         tttCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTg                                             >  1:1952831/1‑58 (MQ=255)
                         tttCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCATCATCCACAACGTATTg                                             >  1:1194438/1‑58 (MQ=255)
                           tCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCg                                >  1:3145744/1‑69 (MQ=255)
                             aCGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCg                                >  1:509301/1‑67 (MQ=255)
                             aCGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCg                                >  1:851598/1‑67 (MQ=255)
                             aCGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTa                              <  1:2763131/69‑1 (MQ=255)
                                cgcCCACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCg                                <  1:909990/64‑1 (MQ=255)
                                 gcCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGtt                        >  1:343642/1‑71 (MQ=255)
                                  cccACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGttt                       <  1:2408197/71‑1 (MQ=255)
                                   ccACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGtttt                      >  1:2838259/1‑71 (MQ=255)
                                   ccACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGtttt                      <  1:939539/71‑1 (MQ=255)
                                   ccACTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGtttt                      <  1:2916342/71‑1 (MQ=255)
                                     aCTGGCTGGAGACATAGCCCGTATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTg                     >  1:3296548/1‑70 (MQ=255)
                                     aCTGGCTGGAGACATAGCCCGGATGCGCGTCATGc                                                        <  1:538659/35‑1 (MQ=255)
                                     aCTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTa                                                >  1:2348441/1‑43 (MQ=255)
                                     aCTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTg                     >  1:54973/1‑70 (MQ=255)
                                     aCTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGtttt                      >  1:3624743/1‑69 (MQ=255)
                                     aCTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGtttt                      >  1:1004317/1‑69 (MQ=255)
                                     aCTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTg                     >  1:3031115/1‑70 (MQ=255)
                                     aCTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTg                     >  1:550794/1‑70 (MQ=255)
                                     aCTGGCTGGAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTg                     >  1:2161356/1‑70 (MQ=255)
                                      ctggctggAGACATAGCCCGGATGCGCATCATGCACAACg                                                  <  1:2861382/40‑1 (MQ=255)
                                      ctggctggAGACATAGCCCGGATGCGCATCATGCACAACg                                                  <  1:605682/40‑1 (MQ=255)
                                       tggctggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGtttt                      <  1:407830/67‑1 (MQ=255)
                                        ggctggAGACATAGCCCGGATGCGCGTCATGCACAAc                                                   >  1:2917736/1‑37 (MQ=255)
                                        ggctggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGc                                            >  1:30468/1‑44 (MQ=255)
                                        ggctggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGtttt                      <  1:3653981/66‑1 (MQ=255)
                                        ggctggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCAt                  >  1:1713364/1‑70 (MQ=255)
                                        ggctggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCAt                  >  1:3130176/1‑70 (MQ=255)
                                        ggctggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCAt                  >  1:837470/1‑70 (MQ=255)
                                         gctggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGc                    >  1:2813750/1‑67 (MQ=255)
                                         gctggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGc                    >  1:2396969/1‑67 (MQ=255)
                                         gctggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGc                    >  1:3629819/1‑67 (MQ=255)
                                           tggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGttt                       >  1:647456/1‑62 (MQ=255)
                                           tggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGttt                       >  1:1240828/1‑62 (MQ=255)
                                           tggAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGttt                       >  1:144015/1‑62 (MQ=255)
                                           tggAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGt                               <  1:2339049/54‑1 (MQ=255)
                                           tggAGACATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGt                               <  1:3284054/54‑1 (MQ=255)
                                             gagaCATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCAt                  <  1:2630556/65‑1 (MQ=255)
                                               gaCATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGt           >  1:3599406/1‑70 (MQ=255)
                                               gaCATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGt           >  1:1570354/1‑70 (MQ=255)
                                                 cATAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAAc        <  1:58240/71‑1 (MQ=255)
                                                  aTAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACg       <  1:3472775/71‑1 (MQ=255)
                                                  aTAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACg       <  1:3477221/71‑1 (MQ=255)
                                                  aTAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACg       <  1:266633/71‑1 (MQ=255)
                                                  aTAGCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACg       <  1:2544741/71‑1 (MQ=255)
                                                   tAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGAt                           <  1:2860058/50‑1 (MQ=255)
                                                   tAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGAt                           <  1:1175628/50‑1 (MQ=255)
                                                   tAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGAt                           <  1:2102799/50‑1 (MQ=255)
                                                   tAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGAt                           <  1:2436431/50‑1 (MQ=255)
                                                   tAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGAt                           <  1:3567778/50‑1 (MQ=255)
                                                    aGCCCGGAAGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGAt                           >  1:1862285/1‑49 (MQ=255)
                                                     gCCCGGATGGGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCtt    >  1:1401666/1‑71 (MQ=255)
                                                     gCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCg                                >  1:1097726/1‑43 (MQ=255)
                                                     gCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCg                                >  1:305699/1‑43 (MQ=255)
                                                     gCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCg                                >  1:236269/1‑43 (MQ=255)
                                                     gCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCtt    <  1:3083377/71‑1 (MQ=255)
                                                     gCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCtt    <  1:2940334/71‑1 (MQ=255)
                                                     gCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCtt    <  1:1854990/71‑1 (MQ=255)
                                                     gCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCtt    >  1:2361926/1‑71 (MQ=255)
                                                     gCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCtt    >  1:194431/1‑71 (MQ=255)
                                                     gCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCt     <  1:1944942/70‑1 (MQ=255)
                                                     gCCCGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATAGTAGATGTTTTGCATCa                >  1:3520061/1‑59 (MQ=39)
                                                      cccGGATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTc   <  1:1115757/71‑1 (MQ=255)
                                                          gATGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAAc        <  1:3291801/62‑1 (MQ=255)
                                                           aTGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTcg  >  1:3310277/1‑67 (MQ=255)
                                                            tGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGt           <  1:3119100/57‑1 (MQ=255)
                                                            tGCGCATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGt           <  1:326810/57‑1 (MQ=255)
                                                             gcgcATCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGt           <  1:2394053/56‑1 (MQ=255)
                                                               gcgTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTgt                    <  1:1605857/45‑2 (MQ=255)
                                                               gcgTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGc                    <  1:3605173/45‑1 (MQ=255)
                                                               gcgTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGc                    <  1:2048358/45‑1 (MQ=255)
                                                                 |                                                            
AGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACGTATTGCGGGGTAAAATCGTAGATGTTTTGCATCAGGCGTAACGCTTCG  >  Exported/1824223‑1824348

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