Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,197,365 G→A pseudogene (892/2847 nt) yehH → DUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,197,3650GA94.5% 231.6 / 3.0 73pseudogene (892/2847 nt)yehHDUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (37/32);  total (39/34)
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.86e-01

CTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTC  >  NC_000913/2197268‑2197458
                                                                                                 |                                                                                             
cTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTAcc                                                                                             <  2:1731969/100‑1 (MQ=255)
 tACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCg                                                                                            <  1:1415206/100‑1 (MQ=255)
  aCTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGc                                                                                           >  1:1888542/1‑100 (MQ=255)
      ggAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAAtt                                                                                       >  1:749669/1‑100 (MQ=255)
      ggAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAAtt                                                                                       <  1:1916111/100‑1 (MQ=255)
      ggAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAAtt                                                                                       <  1:1913421/100‑1 (MQ=255)
       gAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGATGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTc                                                                                      >  2:1248623/1‑100 (MQ=255)
        aGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTcc                                                                                     >  2:1258747/1‑100 (MQ=255)
         gtaaTATTTTAGATTTGAAAGAGGCCTACGCAAAAATTTATCACTGACGGAATAGACTCTCTGACAAAAATATGCACAAAAACTTATTACAGCAATTCca                                                                                    >  1:1352941/3‑100 (MQ=255)
         gCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCca                                                                                    >  2:112799/1‑100 (MQ=255)
         gCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCca                                                                                    <  1:1494040/100‑1 (MQ=255)
          cAATATTGCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCACAAATGTGCACAAAAACTTATTACCGCAATTCcac                                                                                   >  1:1023033/1‑100 (MQ=255)
           aaTATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCcaca                                                                                  >  1:1631165/1‑100 (MQ=255)
           aaTATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCcaca                                                                                  <  2:1053487/100‑1 (MQ=255)
           aaTATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCcaca                                                                                  <  2:1162508/100‑1 (MQ=255)
              atTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACatat                                                                               >  2:604747/1‑100 (MQ=255)
               ttCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATAtt                                                                              >  1:1748966/1‑100 (MQ=255)
               ttCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATAtt                                                                              >  1:657565/1‑100 (MQ=255)
                  cAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCa                                                                           >  1:1444699/1‑100 (MQ=255)
                  cAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCa                                                                           >  1:1247468/1‑100 (MQ=255)
                     aTTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGaa                                                                        >  2:879233/1‑100 (MQ=255)
                      tttGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAAt                                                                       <  2:379429/100‑1 (MQ=255)
                       ttGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATg                                                                      >  1:951687/1‑100 (MQ=255)
                         gAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggc                                                                    <  2:7712/100‑1 (MQ=255)
                          aaTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggcg                                                                   <  2:113862/100‑1 (MQ=255)
                          aaTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggcg                                                                   >  1:1440055/1‑100 (MQ=255)
                           aTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggcgg                                                                  <  2:1383201/100‑1 (MQ=255)
                              aCGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCaa                                                               <  1:1503290/100‑1 (MQ=255)
                                 ccTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCa                                                            <  1:992953/100‑1 (MQ=255)
                                 ccTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCa                                                            >  1:858294/1‑100 (MQ=255)
                                   tACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCAtt                                                          >  1:157714/1‑100 (MQ=255)
                                     cGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTaa                                                        >  1:1033689/1‑100 (MQ=255)
                                      gTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAAt                                                       >  2:774567/1‑100 (MQ=255)
                                        aaaCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTg                                                     >  2:41556/1‑100 (MQ=255)
                                           ttccgatctcTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCtt                                                  <  2:91578/91‑1 (MQ=255)
                                             tttATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTg                                                <  2:506077/100‑1 (MQ=255)
                                              ttATCACTGATGGAAGAGAATATCTAGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATGTTTCAGAATGGCGGCAACCATTAATTGCTTTGt                                               <  1:1278302/100‑1 (MQ=255)
                                               tATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGtt                                              >  2:1644040/1‑100 (MQ=255)
                                               tATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGtt                                              <  2:902890/100‑1 (MQ=255)
                                                aTCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGtta                                             <  1:218630/100‑1 (MQ=255)
                                                 tCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGttat                                            <  2:988836/100‑1 (MQ=255)
                                                   aCTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGttatta                                          >  2:249409/1‑100 (MQ=255)
                                                    cTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCAAGAATGGAGGCAACCATTAATTGCttagatcggaa                                         >  1:91578/1‑91 (MQ=255)
                                                       aCGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCg                                      >  2:1249034/1‑100 (MQ=255)
                                                          gAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGaaa                                   >  2:858685/1‑100 (MQ=255)
                                                              agTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAAcc                               <  2:1401602/100‑1 (MQ=255)
                                                              agTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAAcc                               <  2:422988/100‑1 (MQ=255)
                                                                 ctctctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGa                            <  2:1219746/100‑1 (MQ=255)
                                                                   ctctGGCAAAAATGTGCACACAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCCTTAATTGCTTTGTTATTACCCGAAAAACAAGaaa                          >  1:873174/1‑100 (MQ=255)
                                                                   ctctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGaaa                          >  1:769199/1‑100 (MQ=255)
                                                                   ctctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCAAATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGaaa                          >  1:134047/1‑100 (MQ=255)
                                                                    tctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAAt                         <  2:564688/100‑1 (MQ=255)
                                                                        gCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGca                     <  1:279189/100‑1 (MQ=255)
                                                                        gCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGca                     >  1:332879/1‑100 (MQ=255)
                                                                         cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                    >  1:1489472/1‑100 (MQ=255)
                                                                         cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                    >  1:593918/1‑100 (MQ=255)
                                                                         cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                    >  1:1603996/1‑100 (MQ=255)
                                                                         cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                    <  2:1611673/100‑1 (MQ=255)
                                                                         cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATAATCCAGAATGGCGGCAACCATTAATTGCTTTGAAAAAAAAAAAAAAACCAGAAATTGca                     <  1:431096/99‑1 (MQ=255)
                                                                                gtgCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAAt             <  2:764638/100‑1 (MQ=255)
                                                                                 tgCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAAtt            >  1:119768/1‑100 (MQ=255)
                                                                                 tgCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAAtt            >  1:730732/1‑100 (MQ=255)
                                                                                  gCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTg           >  2:1923992/1‑100 (MQ=255)
                                                                                    acaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGcc         <  1:518073/100‑1 (MQ=255)
                                                                                     caAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGccc        <  2:1631175/100‑1 (MQ=255)
                                                                                     caAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGccc        >  1:393284/1‑100 (MQ=255)
                                                                                     caAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGccc        >  1:378263/1‑100 (MQ=255)
                                                                                     caAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGccc        <  2:1263673/100‑1 (MQ=255)
                                                                                     caAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGACAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACAGGAAATTGccc        >  1:1279806/1‑100 (MQ=255)
                                                                                       aaaaCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAg      <  2:1497208/100‑1 (MQ=255)
                                                                                         aaCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCg    <  1:41556/100‑1 (MQ=255)
                                                                                          aCTAATTGCCGAACTTCCAAATATTCCAGCCTGGCGGCAAACATTCATTGCTTTGTTATTACCCGCAAACCCAGAACTTGCACCTGAAATTGcaaag      <  1:771089/97‑5 (MQ=255)
                                                                                           cTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTc  <  2:1230360/100‑1 (MQ=255)
                                                                                                 |                                                                                             
CTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTC  >  NC_000913/2197268‑2197458

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