Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 300,444 T→C D164G (GAC→GGC)  paoC ← aldehyde dehydrogenase: molybdenum cofactor‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913300,4440TC97.7% 141.8 / NA 44D164G (GAC→GGC) paoCaldehyde dehydrogenase: molybdenum cofactor‑binding subunit
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (23/20);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GAGGTCCAGAGGGTAAGCTTATTTCCATCCCAGACGGCCATCGAGGCATGCGGCTCCATCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGTCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGT  >  NC_000913/300309‑300580
                                                                                                                                       |                                                                                                                                        
gAGGTCCAGAGGGTAAGCTTATTTCCATCCCAGACGGCCATCGAGGCATGCGGCTCCATCGCCACATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCAcc                                                                                                                                      <  1:162263/140‑1 (MQ=255)
   gTCCAGAGGGTAAGCTTATTTCCATCCCAGACGGCCATCGAGGCATGCGGCTCCATCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGAc                                                                                                                                   >  2:752570/1‑140 (MQ=255)
      cAGAGGGTAAGCTTATTTCCATCCCAGACGGCCATCGAGGCATGCGGCTCCATCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCAGCGCCGCGGAGGTGAAAGCCCCTTCAAAGCCACCGACGtt                                                                                                                                >  2:273213/1‑140 (MQ=255)
                gCTTATTTCCATCCCAGACGGCCATCGAGGCATGCGGCTCCATCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGc                                                                                                                      <  1:123750/140‑1 (MQ=255)
                gCTTATTTCCATCCCAGACGGCAATCGAGACATGCGGCTCCATCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGc                                                                                                                      <  1:372236/140‑1 (MQ=255)
                     tttCCATCCCAGACGGCCATCGAGGCATGCGGCTCCATCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTc                                                                                                                 >  2:734204/1‑140 (MQ=255)
                     tttCCATCCCAGACGGCCATCGAGGCATGCGGCTCCATCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTc                                                                                                                 >  2:435665/1‑140 (MQ=255)
                     tttCCATCCCAGACGGCCATCGAGGCATGCGGCTCCATCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTc                                                                                                                 >  1:726731/1‑140 (MQ=255)
                                gACGGCCATCGAGGCATGCGGCTCCATCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCcggcggtt                                                                                                      <  1:752570/140‑3 (MQ=255)
                                       aTCGAGGCATGCGGCTCCATCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGAc                                                                                               >  2:744094/1‑140 (MQ=255)
                                         cGAGGCATGCGGCTCCATCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACgg                                                                                             >  1:492973/1‑140 (MQ=255)
                                         cGAGGCATGCGGCTCCATCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACgg                                                                                             >  1:520779/1‑140 (MQ=255)
                                                     cTCCATCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGt                                                                                 <  1:435665/140‑1 (MQ=255)
                                                       ccATCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTcc                                                                               >  1:476152/1‑140 (MQ=255)
                                                        cATCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTATCTTCCGGCGGCTGGTTGACGACCTGTTTTTCGTCCg                                                                              >  1:345511/1‑140 (MQ=255)
                                                        cATCGCCATATGGCTCTGGTCCGGGGCCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTGTGTCGGGCGTGTATTCGGGCGGCTGTTTGACGGCCTGTTTTTCGTccc                                                                              >  2:298722/1‑139 (MQ=255)
                                                         aTCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGc                                                                             >  1:790689/1‑140 (MQ=255)
                                                           cGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCa                                                                           <  2:726731/140‑1 (MQ=255)
                                                           cGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCa                                                                           <  2:37403/140‑1 (MQ=255)
                                                             ccATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCTGGCGTGTCCTCCTGCGGCTGATTGAGGGCCTGTTTTTCGTCCGCCAgg                                                                         >  2:247588/1‑140 (MQ=255)
                                                                   ggCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTaa                                                                   <  2:492973/140‑1 (MQ=255)
                                                                   ggCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTaa                                                                   <  2:520779/140‑1 (MQ=255)
                                                                    gCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAg                                                                  <  2:162033/140‑1 (MQ=255)
                                                                    gCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAg                                                                  <  1:734204/140‑1 (MQ=255)
                                                                     ctctGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGc                                                                 >  1:93637/1‑140 (MQ=255)
                                                                     ctctGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGc                                                                 >  1:574854/1‑140 (MQ=255)
                                                                     ctctGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGAGTAAGCt                                                                >  2:704685/1‑140 (MQ=255)
                                                                       ctGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTc                                                               <  2:591128/140‑1 (MQ=255)
                                                                           tCCGGGGTCGTGTAGATAGCATTAAGCTTTACCGCAGGGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCttt                                                           <  1:344999/140‑1 (MQ=255)
                                                                             cGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTttat                                                         <  2:93637/140‑1 (MQ=255)
                                                                             cGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTttat                                                         <  2:574854/140‑1 (MQ=255)
                                                                                         ggtagCATCAATCTTTACCGCACCGGAGGGGAAAGCCACGTCAAAGTCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAgt                                             <  1:275375/140‑1 (MQ=255)
                                                                                            agCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGtgtg                                          >  2:463341/1‑140 (MQ=255)
                                                                                              cATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGTCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGcc                                        >  1:360610/1‑140 (MQ=255)
                                                                                                    tCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGAGAGTAAGATCCTTTATTACGGGGATAGTGTGCCTGCAcc                                  >  2:211715/1‑140 (MQ=255)
                                                                                                              aGCGGAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGgccg                        >  1:650577/1‑140 (MQ=255)
                                                                                                                  gAGGTGAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCcgccgc                    >  2:433202/1‑140 (MQ=255)
                                                                                                                       gAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTcgcg               <  1:433202/140‑1 (MQ=255)
                                                                                                                       gAAAGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTCgcgc              <  1:704685/140‑1 (MQ=255)
                                                                                                                        aaaGCCCCGTCAAAGCCGCCGAAGTTTTTGTCGGGCGTGGCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTCgcgc              <  2:345511/140‑1 (MQ=255)
                                                                                                                        aaaGCCCCGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTCgcgc              >  2:692857/1‑140 (MQ=255)
                                                                                                                             cccGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTcg                 >  2:33494/1‑132 (MQ=255)
                                                                                                                             cccGTCAAAGCCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTcg                 <  1:33494/132‑1 (MQ=255)
                                                                                                                                    aaGCCACCGACGTTTTTGTCGGCCGCGTCGTCCGGCGCCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGt  <  1:211715/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                        
GAGGTCCAGAGGGTAAGCTTATTTCCATCCCAGACGGCCATCGAGGCATGCGGCTCCATCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGTCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGT  >  NC_000913/300309‑300580

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