Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 799,720 (A)6→5 coding (856/954 nt) ybhD ← putative DNA‑binding transcriptional regulator YbhD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913799,7150A.93.3% 158.4 / 7.4 45coding (861/954 nt)ybhDputative DNA‑binding transcriptional regulator YbhD
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base . (19/23);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CGCTGATTTACCCGGCGCAGTCTCTCCTGCGCCGGTGTATTAACCTATCTCCTGTAACGCGTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCACCAATGGGATACCAATAATCGAAGAATAGGGCGCATGAGGCATTGCCA  >  NC_000913/799583‑799842
                                                                                                                                    |                                                                                                                               
cGCTGATTTACCCGGCGCAGTCTCTCCTGCGCCGGTGTATTAACCTATCTCCTGTAACGCGTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCa                                                                                                                                 >  1:449082/1‑133 (MQ=255)
  cTGATTTACCCGGCGTAGTCTCGCCTGCGCTGGTGTATTAACCTATCTCCTGTAACGCGTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCaaa                                                                                                                               <  1:253347/133‑1 (MQ=255)
       ttACCCGGCGCAGTCTCTCCTGCGCCGGTGTATTAACCTATCTCCTGTAACGCGTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCaaaaagcg                                                                                                                          <  1:104599/133‑4 (MQ=255)
          cccGGCGCAGTCTCTCCTGCGCCGGTGTATTAACCTATCTCCTGTAACGCGTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTGATCGGCGAGGCAAGAAGCGTTc                                                                                                                      >  1:185420/1‑133 (MQ=255)
             ggCGCAGTCTCTCCTGCGCCGGTGTATTAACCTATCTCCTGTAACGCGTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATAAGCCAGGCAAAAAACGTTctgc                                                                                                                   >  1:258328/1‑133 (MQ=255)
               cgcAGTCTCTCCTGCGCCGGTGTATTAACCTATCTCCTGTAACGCGTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTctgctg                                                                                                                 <  1:156369/133‑1 (MQ=255)
                       ctcCTGCGCCGGTGTCTTAACCTACCTCTTGTCACGCGTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGaa                                                                                                         <  2:185420/133‑1 (MQ=255)
                         ccTGCGCCGGTGTATTAACCTATCTCCTGTAACGCGTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAg                                                                                                       >  1:397004/1‑133 (MQ=255)
                              gcCGGTGTATTAACCTATCTCCTGTAACGCGTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACgg                                                                                                  >  1:77739/1‑133 (MQ=255)
                                cGGTGTATTAACCTATCTCCTGTAACGCGTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCa                                                                                                >  1:603926/1‑133 (MQ=255)
                                          aCCTATCTCCTGTAACGCGTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGc                                                                                      <  2:570476/133‑1 (MQ=255)
                                               tctcCTGTAACGCGTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGc                                                                                      >  2:574652/1‑128 (MQ=255)
                                               tctcCTGTAACGCGTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGc                                                                                      <  1:574652/128‑1 (MQ=255)
                                               tctcCTGTAACGCGTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAatt                                                                                 >  1:571123/1‑133 (MQ=255)
                                               tctcCTGTAACGCGTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAatt                                                                                 >  1:101760/1‑133 (MQ=255)
                                                    tGTAACGCGTGTCTCTGGCGTTCGATGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTcc                                                                            >  2:433342/1‑133 (MQ=255)
                                                    tGTAACGCGTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTc                                                                             >  1:568642/1‑133 (MQ=255)
                                                    tGTAACGCGCGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTcc                                                                            >  2:463663/1‑133 (MQ=255)
                                                           cgTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGtt                                                                     <  1:579185/133‑1 (MQ=255)
                                                           cgTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGATCTGCTGCTCGAGAAAGTCCGGCATCGTTACGGCGAATTAGTCCCAATGtt                                                                     >  2:252159/1‑133 (MQ=255)
                                                             tgtCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCg                                                                   >  2:429083/1‑133 (MQ=255)
                                                             tgtCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCg                                                                   <  2:603926/133‑1 (MQ=255)
                                                             tgtCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCAACTTAACGGCGAATTATTCCCAATGTTCg                                                                   >  2:276121/1‑133 (MQ=255)
                                                              gtCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGa                                                                  <  2:101760/133‑1 (MQ=255)
                                                                 tctgctgTTCGACGATATTGGTCCATAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGt                                                               <  2:397004/127‑1 (MQ=255)
                                                                     gCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAAATAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATaa                                                           <  2:584975/133‑1 (MQ=255)
                                                                       gTTCGACGATATTGGTCCACAATTTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACg                                                         <  1:252159/133‑1 (MQ=255)
                                                                        ttCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGAGAATTATACCCAATGTTCGACGTATAACgg                                                        >  1:180067/1‑133 (MQ=255)
                                                                          cgacgaTATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGt                                                      >  2:375816/1‑133 (MQ=255)
                                                                           gacgaTATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGt                                                      <  2:568642/133‑1 (MQ=255)
                                                                                  ttGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCACCaa                                              <  1:541577/133‑1 (MQ=255)
                                                                                   tGGTCCACAACTTGTCTTTATCGTCAGTCCATCAATTAATCTGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCACCAAt                                             <  2:180067/133‑1 (MQ=255)
                                                                                   tGGTCCACAAATTGTTTTTATCGTCAGTTCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCACCAAt                                             <  2:271882/133‑1 (MQ=255)
                                                                                      tCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCACCAATggg                                          <  1:433342/133‑1 (MQ=255)
                                                                                       ccACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCACCAATGGGa                                         <  2:119340/133‑1 (MQ=255)
                                                                                              ttGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCACCAATGGGATACCaat                                  <  2:77739/133‑1 (MQ=255)
                                                                                              ttGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCACCAATGGGATACCaat                                  <  1:375816/133‑1 (MQ=255)
                                                                                                gTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCACCAATGGGATACCaataa                                <  2:571123/133‑1 (MQ=255)
                                                                                                gTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCACCAATGGGATACCaataa                                >  2:580450/1‑133 (MQ=255)
                                                                                                            aGTCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCACCAATGGGATACCAATAATCGAAGAATAgg                    <  1:340188/133‑1 (MQ=255)
                                                                                                              tCCATAAATTAATCAGCAGGGCCAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCACCAATGGGATACCAATAATCGAAGAATAGGgc                  <  1:272751/133‑1 (MQ=255)
                                                                                                              tCCATAAATTAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCACCAATGGGATACCAATAATCGAAGAATAGGgc                  <  2:515550/133‑1 (MQ=255)
                                                                                                                       tAATCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCACCAATGGGATACCAATAATCGAAGAATAGGGCGCATGAGGc         <  2:449082/133‑1 (MQ=255)
                                                                                                                        aaTCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCACCAATGGGATACCAATAATCGAAGAATAGGGCGCATGAGGCa        >  1:393516/1‑133 (MQ=255)
                                                                                                                        aaTCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCAACAATGGGATACCAATACTCGAAGAAGAGGGCGCATGAGGCa        >  2:272751/1‑133 (MQ=255)
                                                                                                                          tCAGCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCACCAATGGGATACCAATAATCGAAGAATAGGGCGCATGAGGCAtt      >  1:544316/1‑133 (MQ=255)
                                                                                                                             gCAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCACCAATGGGATACCAATAATCGAAGAATAGGGCGCATGAGGCATTgcc   <  1:580450/133‑1 (MQ=255)
                                                                                                                              cAAGGCAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCACCAATGGGATACCAATAATCGAAGAATAGGGCGCATGAGGCATTgcca  >  2:524619/1‑133 (MQ=255)
                                                                                                                                    |                                                                                                                               
CGCTGATTTACCCGGCGCAGTCTCTCCTGCGCCGGTGTATTAACCTATCTCCTGTAACGCGTGTCTCTGGCGTTCGACGATATTGGTCCACAAATTGTCTTTATCGTCAGTCCATAAATTAATCAGCAAGGCAAAAAAGCGTTCTGCTGCTGGAGAAAGTACGGCATCTTTACGGCGAATTATTCCCAATGTTCGACGTATAACGGGTTCCACCAATGGGATACCAATAATCGAAGAATAGGGCGCATGAGGCATTGCCA  >  NC_000913/799583‑799842

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