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.4% 157.8 / ‑1.2 54pseudogene (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 (1/2);  new base A (26/25);  total (27/27)
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

CTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCG  >  NC_000913/2197268‑2197456
                                                                                                 |                                                                                           
cTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTAcc                                                                                           <  1:1276280/100‑1 (MQ=255)
cTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTAcc                                                                                           >  2:1187003/1‑100 (MQ=255)
  aCTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGc                                                                                         >  2:1461204/1‑100 (MQ=255)
   cTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCa                                                                                        >  2:501672/1‑100 (MQ=255)
   cTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCa                                                                                        >  2:1200527/1‑100 (MQ=255)
     tGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAAt                                                                                      >  2:1322297/1‑100 (MQ=255)
      ggAGTAATATTCCAGATTTGAATGACGCCTACGTAACCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAAtt                                                                                     <  1:308580/100‑1 (MQ=255)
      ggAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAAtt                                                                                     >  1:63438/1‑100 (MQ=255)
               ttCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATAtt                                                                            >  2:934545/1‑100 (MQ=255)
                tCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTc                                                                           >  2:523689/1‑100 (MQ=255)
                 ccAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTc                                                                           >  1:62338/1‑99 (MQ=255)
                  ccgATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCa                                                                         >  1:1072422/3‑100 (MQ=255)
                  cAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCa                                                                         <  2:1025616/100‑1 (MQ=255)
                   aGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAg                                                                        <  2:658117/100‑1 (MQ=255)
                         gAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggc                                                                  <  1:1095503/100‑1 (MQ=255)
                          aaTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggcg                                                                 >  1:1407946/1‑100 (MQ=255)
                          aaTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggcg                                                                 <  1:469820/100‑1 (MQ=255)
                          aaTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggcg                                                                 >  2:172720/1‑100 (MQ=255)
                           aTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAATAACGTATTACCGCAATTCCACATATTCCAGAATggcg                                                                 >  2:858789/1‑99 (MQ=255)
                           aTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggcgg                                                                <  1:568172/100‑1 (MQ=255)
                           aTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggcgg                                                                <  1:430219/100‑1 (MQ=255)
                                gCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAAcc                                                           <  2:62338/100‑1 (MQ=255)
                                ccgatcctAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAAcc                                                           <  1:25524/93‑1 (MQ=255)
                                 ccTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCa                                                          <  1:1247094/100‑1 (MQ=255)
                                   tACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCAtt                                                        <  1:1199510/100‑1 (MQ=255)
                                   tACGAAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAAGTCCACATATTCCAGAATGGCGGCAACCAtn                                                        <  2:1441562/100‑2 (MQ=255)
                                     cGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTaa                                                      <  1:996105/100‑1 (MQ=255)
                                       tAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCagatcgg                                                    >  2:25524/1‑94 (MQ=255)
                                       tAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAAtt                                                    >  1:983453/1‑100 (MQ=255)
                                          aCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCt                                                 <  2:1402011/100‑1 (MQ=255)
                                          aCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCt                                                 <  1:173841/100‑1 (MQ=255)
                                              ttATCACTGACGTAATAGGCTCTCTGTCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATGTTCCAGAATGGCGGCAACCATTAATTGCttn                                               <  2:333946/98‑2 (MQ=255)
                                              ttATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACTTATTCCAGAATGGCGGCAACCATTAATTGCTTTGt                                             <  2:995295/100‑1 (MQ=255)
                                                aTCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGtta                                           <  2:669603/100‑1 (MQ=255)
                                                 tCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGttat                                          >  2:1087292/1‑100 (MQ=255)
                                                      gACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAccc                                     <  1:897383/100‑1 (MQ=255)
                                                             gagTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGTCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAAc                              >  2:1068438/1‑100 (MQ=255)
                                                               gTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCa                            >  1:1267022/1‑100 (MQ=255)
                                                                 ctctctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGa                          >  1:1336765/1‑100 (MQ=255)
                                                                 ctctctGGAAAAAAAAAACAAAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGa                          <  2:1537452/100‑1 (MQ=255)
                                                                  tctctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGaa                         <  2:531350/100‑1 (MQ=255)
                                                                     ctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAAtt                      <  1:867987/100‑1 (MQ=255)
                                                                       ggCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGc                    <  2:1336758/100‑1 (MQ=255)
                                                                       ggCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGc                    >  1:1194764/1‑100 (MQ=255)
                                                                         cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                  >  2:1578019/1‑100 (MQ=255)
                                                                         cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                  >  2:142441/1‑100 (MQ=255)
                                                                         cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                  >  2:1004171/1‑100 (MQ=255)
                                                                                  gCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTACTTCGTTATTACCCGACAAAACAGAAATTGCACAGGaactcg         >  2:697278/1‑96 (MQ=255)
                                                                                   cacaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGc        >  2:189422/1‑100 (MQ=255)
                                                                                   cacaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGc        <  1:845380/100‑1 (MQ=255)
                                                                                     caAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGccc      >  2:1183505/1‑100 (MQ=255)
                                                                                     caAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGccc      <  2:584578/100‑1 (MQ=255)
                                                                                       aaaaCTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAg    <  1:1087297/100‑1 (MQ=255)
                                                                                         aaCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCg  >  2:131363/1‑100 (MQ=255)
                                                                                                 |                                                                                           
CTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCG  >  NC_000913/2197268‑2197456

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