Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 299,635 C→T 9.1% G434S (GGT→AGT)  paoC ← PaoABC aldehyde oxidoreductase, Moco‑containing subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913299,6350CT 9.1% 126.6 / 3.7 44G434S (GGT→AGT) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base C (22/18);  new base T (2/2);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

TTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTCAACCTGAG  >  NC_000913/299489‑299781
                                                                                                                                                  |                                                                                                                                                  
ttCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGt                                                                                                                                                  <  2:5074/149‑1 (MQ=255)
 tCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTg                                                                                                                                                 <  2:508346/149‑1 (MQ=255)
    gCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGcc                                                                                                                                              <  2:206264/149‑1 (MQ=255)
       gCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGAc                                                                                                                                           <  2:100446/149‑1 (MQ=255)
                 ccAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTcc                                                                                                                                 >  4:410507/1‑149 (MQ=255)
                 ccAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTcc                                                                                                                                 >  2:201105/1‑149 (MQ=255)
                      aTGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTc                                                                                                                            <  1:82077/149‑1 (MQ=255)
                                tGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGt                                                                                                                  >  4:400883/1‑149 (MQ=255)
                                tGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGt                                                                                                                  >  2:118918/1‑149 (MQ=255)
                                                  tCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTgcgc                                                                                                >  4:248698/1‑149 (MQ=255)
                                                   cATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCgct                                                                                               <  4:177047/149‑1 (MQ=255)
                                                     tGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCgctgc                                                                                             <  1:46970/149‑1 (MQ=255)
                                                       tCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCtt                                                                                           >  3:380390/1‑149 (MQ=255)
                                                       tCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCtt                                                                                           >  1:283045/1‑149 (MQ=255)
                                                       tCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCtt                                                                                           >  1:254507/1‑149 (MQ=255)
                                                       tCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCtt                                                                                           >  3:104652/1‑149 (MQ=255)
                                                                aCGATAACGGTGCTGTTTTTTCCGAGGCGAACCCGAGCACCCGTTGTTTCCAACAGATCATTGCGACAGCCCGCCGCTACGCCAGGGCCGACTAGCCACTCCCCGTCGCGCATCTGTCCTTGTGTGGCTTTTCGCTGCCTCCAGCCaaa                                                                                  <  2:375531/149‑1 (MQ=255)
                                                                 cGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTTCGCTGCTTCCAGCCAGAt                                                                                 <  4:551671/149‑1 (MQ=255)
                                                                      aCGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATc                                                                            <  3:464427/149‑1 (MQ=255)
                                                                               ttttGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGt                                                                   >  2:89374/1‑149 (MQ=255)
                                                                               ttttGTTCGAGGTGAACCCGAGCAACCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGt                                                                   >  2:268621/1‑149 (MQ=255)
                                                                                    ttCGAGGTGAACCCGAGCACCCGATTTTTCCAGGAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGGGTGGCGTTGCGCTGCTTCCAGCCAAAGTTATCCGCTCCGGTGCGCa                                                              >  4:301617/1‑149 (MQ=255)
                                                                                    ttCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCa                                                              >  1:283655/1‑149 (MQ=255)
                                                                                    ttCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCa                                                              >  1:197401/1‑149 (MQ=255)
                                                                                           tGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCTAATTTATCCGCTCCGGTGCGCAAGCACTc                                                       <  2:232394/149‑1 (MQ=255)
                                                                                                ccGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATaa                                                  <  4:104652/149‑1 (MQ=255)
                                                                                                  gAGCACCCGATTTTTCCAGCAGATTATTGCGAACGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGc                                                <  3:61779/149‑1 (MQ=255)
                                                                                                    gCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTg                                              <  4:213438/149‑1 (MQ=255)
                                                                                                        ccGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACgg                                          >  3:440932/1‑149 (MQ=255)
                                                                                                        ccGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACgg                                          >  3:538819/1‑149 (MQ=255)
                                                                                                        ccGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACgg                                          >  1:280619/1‑149 (MQ=255)
                                                                                                          gATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACTGTGGCCGACTAGc                                                                                                                                       >  2:53198/1‑54 (MQ=255)
                                                                                                          gATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACTGTGGCCGACTAGc                                                                                                                                       <  3:17739/54‑1 (MQ=255)
                                                                                                          gATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACTGTGGCCGACTAGc                                                                                                                                       <  1:53198/54‑1 (MQ=255)
                                                                                                          gATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACTGTGGCCGACTAGc                                                                                                                                       >  4:17739/1‑54 (MQ=255)
                                                                                                                  cagcagATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAagc                                >  3:566135/1‑149 (MQ=255)
                                                                                                                  cagcagATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAagc                                >  1:232043/1‑149 (MQ=255)
                                                                                                                   agcagATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAagca                               <  4:528014/149‑1 (MQ=255)
                                                                                                                     cagATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAagcagc                             <  4:98461/149‑1 (MQ=255)
                                                                                                                                     cccGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCggg             <  1:201105/149‑1 (MQ=255)
                                                                                                                                        gccgcAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTc          >  3:215052/1‑149 (MQ=255)
                                                                                                                                           gcAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTCAAc       >  4:301264/1‑149 (MQ=255)
                                                                                                                                           gcAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTCAAc       <  2:522460/149‑1 (MQ=255)
                                                                                                                                                acCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTCAACCTGAg  >  2:182862/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCAACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGATAAGCTGACGGCGAGAGAAGCAGCGCGTCGGGTCGGCGGGGTCAACCTGAG  >  NC_000913/299489‑299781

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