Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 300,585 G→A A117V (GCC→GTC)  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,5850GA90.9% 21.5 / ‑3.4 12A117V (GCC→GTC) paoCaldehyde dehydrogenase: molybdenum cofactor‑binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (3/7);  total (3/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

AGACGGCCATCGAGGCATGCGGCTCCATCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGTCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCAT  >  NC_000913/300340‑300786
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                         
aGACGGCCATCGAGGCATGCGGCTCCATCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGTCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGGGTGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGACCACt                                                                                                                                                                                                      <  1:253425/251‑1 (MQ=255)
   ctgCCATCGAGGCATGCGGCTCCATCGACATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCCCCGCAGCGGCGGTGAAACCCCCGTCAAAGTCACCGACGTTTTTGTCGGGCGTGTCTTACGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGTTTAAGCTCCTTTATTACGGCGATAGCGTGCCTCCGCCAGCGAGTCCGCCGCTCGCGCCTGTTCGAAGGTCTCGACCACTAcc                                                                                                                                                                                                   >  1:11429/3‑251 (MQ=255)
                                                                        tCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGTCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGACCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGACGGTGTTTTTGTCg                                                                                                                                 <  1:343721/248‑1 (MQ=255)
                                                                                  ggAGGTGAAAGCCCCGTCAAAGTCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGTCGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGACCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCt                                                                                                                                                         <  1:440551/214‑1 (MQ=255)
                                                                                        gAAAGCCCCGTCAAAGTCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCgn                                                                                                                                            <  1:107704/221‑2 (MQ=255)
                                                                                                                                      ggcggcTGATTGACGGCCTGTTTTTCGTCNGCCAGGGAGNNAGNTCCTGNANNNCNNNGANAGTGTGCCTGCACCAGCGANGCNNCCGCTCNCGCCTGTTCGAAGGTCTCGACCACTa                                                                                                                                                                                                     <  2:99856/118‑1 (MQ=255)
                                                                                                                                      ggcggcTGATTGAAGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGACCACTa                                                                                                                                                                                                     >  1:99856/1‑118 (MQ=255)
                                                                                                                                                     gCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGANNGTCTCGNnc                                                                                                                                                                                                         >  2:5631/1‑99 (MQ=255)
                                                                                                                                                                               aGCTCCTTTATTACGGCGATAGGGTGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGACCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCGAACAGCCTGACGGTGTTTTTGTGGCCTTTGCCGAGTGCCCCGTCGGTACTGGCGGTAATGACAGCCAGGACGGCCGGCGGTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATgg                       <  1:400573/251‑1 (MQ=255)
                                                                                                                                                                                   ccTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGACCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGACGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCg                     <  1:5631/249‑1 (MQ=255)
                                                                                                                                                                                                    gtgtGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGACCACTACCAGCGCACTGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGACGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAAAGACAGCCCGTACACCCCGCGCTTTTTTCGCGGCGGCCGTATCCACGGCGGTGTGGCGTCCTTTGGCAATGGCGGAACCGACGGTATTGccct  >  1:222140/1‑249 (MQ=255)
                                                                                                                                                                                                                                   gcCTGTTCGAAGGTCTCGACCACTACCAGCGCAATGGCCTGATGAGAGTGCTc                                                                                                                                                                         <  1:176442/53‑1 (MQ=255)
                                                                                                                                                                                                                                      tGNTCGNANNTNNCGCNCACTACCAGCGCAATGGCCTGNNGATAGTGCTNAATAGTGGGGCCGCCTAACAGCCTGACGGTGTTTTTGTCGCCTTTGCCGAGTGNCCCGGCGTTNNTGNCGGTANTNNNANNNAGNACGCCCGGCGCTTTTTGCGNGGNNTCCGTANCAAGGGCGGTGAGGCGTCCTTTGGCa                           <  2:96526/192‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                         
AGACGGCCATCGAGGCATGCGGCTCCATCGCCATATGGCTCTGGTCCGGGGTCGTGTAGGTAGCATCAATCTTCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGTCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCAT  >  NC_000913/300340‑300786

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