Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,835,133 T→A 7.0% intergenic (+88/+40) norW → / ← hypF NADH:flavorubredoxin oxidoreductase/carbamoyl phosphate phosphatase and [NiFe] hydrogenase maturation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,835,1330TA 7.0% 203.5 / 4.3 57intergenic (+88/+40)norW/hypFNADH:flavorubredoxin oxidoreductase/carbamoyl phosphate phosphatase and [NiFe] hydrogenase maturation protein
Reads supporting (aligned to +/- strand):  ref base T (29/24);  new base A (2/2);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGT  >  NC_000913/2835004‑2835280
                                                                                                                                 |                                                                                                                                                   
aTGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGtt                                                                                                                                                     >  1:258483/1‑130 (MQ=255)
aTGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGtt                                                                                                                                                     <  2:258483/130‑1 (MQ=255)
                       gAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATGCGGCGTAAAcgcctta                                                                                                                                            >  1:170612/1‑113 (MQ=25)
                       gAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATGCGGCGTAAAcgcccta                                                                                                                                            <  2:170612/116‑4 (MQ=25)
                                     tGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATGCGGCGTAAAcgccttat                                                                                                                                           <  2:260509/103‑5 (MQ=21)
                                               cTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                       >  3:1063433/1‑97 (MQ=18)
                                                                       cGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCg                                                                                                       <  2:555392/105‑1 (MQ=255)
                                                                       cGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCg                                                                                                       >  1:555392/1‑105 (MQ=255)
                                                                         aCGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTg                                                                                                  >  2:981427/1‑108 (MQ=255)
                                                                         aCGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTg                                                                                                  <  1:981427/108‑1 (MQ=255)
                                                                               gcgcGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGcccc                                                    >  4:862830/1‑148 (MQ=255)
                                                                                 gcgTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTa                                                                                                           >  1:1027106/1‑91 (MQ=255)
                                                                                 gcgTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTa                                                                                                           <  2:1027106/91‑1 (MQ=255)
                                                                                       tATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTa                                                                                                           >  4:285027/1‑85 (MQ=255)
                                                                                       tATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTa                                                                                                           <  3:285027/85‑1 (MQ=255)
                                                                                        aTCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACaaa                                           <  3:862830/148‑1 (MQ=255)
                                                                                         tCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGgtgt                                                                                                                              >  4:263416/1‑64 (MQ=255)
                                                                                         tCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGgtgt                                                                                                                              >  4:718074/1‑64 (MQ=255)
                                                                                         tCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGgtgt                                                                                                                              <  3:718074/64‑1 (MQ=255)
                                                                                         tCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGgtgt                                                                                                                              <  3:263416/64‑1 (MQ=255)
                                                                                          ccGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGgtg                                                                                                                               <  2:897339/62‑1 (MQ=255)
                                                                                          ccGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGgtg                                                                                                                               >  1:897339/1‑62 (MQ=255)
                                                                                          ccGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGcc                                                      <  4:352646/135‑1 (MQ=255)
                                                                                          ccGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGcc                                                      >  3:352646/1‑135 (MQ=255)
                                                                                                tACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGt                                                                                                               <  3:478677/72‑1 (MQ=255)
                                                                                                tACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGt                                                                                                               >  4:478677/1‑72 (MQ=255)
                                                                                                tACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATc                                   >  3:740232/1‑148 (MQ=255)
                                                                                                 aCGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAac                                                                                                                   <  3:385189/67‑1 (MQ=255)
                                                                                                 aCGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAac                                                                                                                   >  4:385189/1‑67 (MQ=255)
                                                                                                      gACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCg                                                                                       <  1:1113407/90‑1 (MQ=255)
                                                                                                      gACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCg                                                                                       >  2:1113407/1‑90 (MQ=255)
                                                                                                          cATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGt                         >  4:764485/1‑148 (MQ=255)
                                                                                                            tgtgTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGa                                                                                                <  3:188637/75‑1 (MQ=255)
                                                                                                            tgtgTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGa                                                                                                >  4:188637/1‑75 (MQ=255)
                                                                                                            tgtgTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCa                                                                                           <  2:756837/80‑1 (MQ=255)
                                                                                                            tgtgTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCa                                                                                           >  1:756837/1‑80 (MQ=255)
                                                                                                            tgtgTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAAcgcg                                                                           <  2:680015/96‑1 (MQ=255)
                                                                                                            tgtgTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAAcgcg                                                                           >  1:680015/1‑96 (MQ=255)
                                                                                                            tgtgTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCagag                                                <  2:504663/123‑1 (MQ=255)
                                                                                                            tgtgTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCagag                                                >  1:504663/1‑123 (MQ=255)
                                                                                                                tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTc                                                                                            >  4:652799/1‑75 (MQ=255)
                                                                                                                tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTc                                                                                            <  3:652799/75‑1 (MQ=255)
                                                                                                                   gCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGt                                                                                                               >  2:603762/1‑53 (MQ=255)
                                                                                                                   gCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGt                                                                                                               <  1:603762/53‑1 (MQ=255)
                                                                                                                    ccGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACccgcc                            <  4:222004/135‑1 (MQ=255)
                                                                                                                    ccGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACccgcc                            >  3:222004/1‑135 (MQ=255)
                                                                                                                       gATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAgaga                                               >  3:630474/1‑113 (MQ=255)
                                                                                                                       gATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAgaga                                               <  4:630474/113‑1 (MQ=255)
                                                                                                                       gATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGaa            >  3:547471/1‑148 (MQ=255)
                                                                                                                        aTAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGcccc                                                    >  3:740055/1‑107 (MQ=255)
                                                                                                                        aTAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGcccc                                                    <  4:740055/107‑1 (MQ=255)
                                                                                                                        aTAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGaaa           >  1:217260/1‑148 (MQ=255)
                                                                                                                              cGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGTCTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCa      >  4:727093/1‑147 (MQ=255)
                                                                                                                              cGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCaa                                                                               <  3:704785/74‑1 (MQ=255)
                                                                                                                              cGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCaa                                                                               >  4:704785/1‑74 (MQ=255)
                                                                                                                              cGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACccgccg                           >  1:448978/1‑126 (MQ=255)
                                                                                                                              cGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACccgccg                           <  2:448978/126‑1 (MQ=255)
                                                                                                                              cGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCaa     >  2:306489/1‑148 (MQ=255)
                                                                                                                              cGTTTACGCCGCATCCGGCAATGGTGTCAAAATGCAACACGTTTTATCCGTTCTGGTCTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCa      <  3:727093/147‑1 (MQ=255)
                                                                                                                                tttACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCaa                                                                               <  3:3539/72‑1 (MQ=255)
                                                                                                                                tttACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCaa                                                                               >  4:3539/1‑72 (MQ=255)
                                                                                                                                tttACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAAcccc                                                           <  1:186762/92‑1 (MQ=255)
                                                                                                                                tttACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAAcccc                                                           >  2:186762/1‑92 (MQ=255)
                                                                                                                                tttACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAAtg   <  4:547471/148‑1 (MQ=255)
                                                                                                                                tttACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAAtg   >  2:688562/1‑148 (MQ=255)
                                                                                                                                 ttACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAAtgt  <  3:194820/148‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                                   
ATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGT  >  NC_000913/2835004‑2835280

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: