Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 300,645 G→A A97V (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,6450GA83.3% 21.1 / ‑3.1 12A97V (GCC→GTC) paoCaldehyde dehydrogenase: molybdenum cofactor‑binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (2/8);  total (2/10)
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

TCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGTCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTCGATGCGGTCATGGGGACGACCG  >  NC_000913/300412‑300881
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              
tCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGTCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGACCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGACGGTGTTTTTGTCg                                                                                                                                                                                                                                  <  1:343721/248‑1 (MQ=255)
                gAAAGCCCCGTCAAAGTCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCgn                                                                                                                                                                                                                                             <  1:107704/221‑2 (MQ=255)
                                                                                                       aGCTCCTTTATTACGGCGATAGGGTGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGACCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCGAACAGCCTGACGGTGTTTTTGTGGCCTTTGCCGAGTGCCCCGTCGGTACTGGCGGTAATGACAGCCAGGACGGCCGGCGGTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATgg                                                                                                                        <  1:400573/251‑1 (MQ=255)
                                                                                                           ccTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGACCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGACGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCg                                                                                                                     <  1:5631/249‑1 (MQ=255)
                                                                                                                            gtgtGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGACCACTACCAGCGCACTGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGACGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAAAGACAGCCCGTACACCCCGCGCTTTTTTCGCGGCGGCCGTATCCACGGCGGTGTGGCGTCCTTTGGCAATGGCGGAACCGACGGTATTGccct                                                                                                  >  1:222140/1‑249 (MQ=255)
                                                                                                                                                              tGNTCGNANNTNNCGCNCACTACCAGCGCAATGGCCTGNNGATAGTGCTNAATAGTGGGGCCGCCTAACAGCCTGACGGTGTTTTTGTCGCCTTTGCCGAGTGNCCCGGCGTTNNTGNCGGTANTNNNANNNAGNACGCCCGGCGCTTTTTGCGNGGNNTCCGTANCAAGGGCGGTGAGGCGTCCTTTGGCa                                                                                                                            <  2:96526/192‑1 (MQ=255)
                                                                                                                                                                      ggtcgcgaccACTACCAGCTCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGACGGTGTTTTTGTc                                                                                                                                                                                                                                   <  1:201014/73‑1 (MQ=255)
                                                                                                                                                                                aCTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGACGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCgtgccgtgcc                                             <  1:379286/251‑1 (MQ=255)
                                                                                                                                                                                                      gatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGt                       <  1:278104/250‑2 (MQ=255)
                                                                                                                                                                                                               cAATAGTGGGGCCGCCTAACAGCCTGACGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAgg                                                                                                                                        <  1:302761/131‑1 (MQ=255)
                                                                                                                                                                                                                           cgccTAACAGCCTGACGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGGCAGCCATGACGGCCGGCGCTTTTTGTGCGGCGTCCGCATCCAGGGCGGTGAGGCGGCCTTTTGGACTGGACGAAACGCACCTATAGACCTAGGCGGCGGTGGGGGGTTCTTCAATACCTTCGTAAGCGTAGCGGGCCGTTCCGGGAGTTTTTATCGGTCCGTCgcgtggggcatgggggcgaccc  >  1:223785/1‑230 (MQ=255)
                                                                                                                                                                                                                                    gCCTGACGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCGTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATgn                                                                                                                        <  1:47449/126‑2 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              
TCACCGCAGCGGAGGTGAAAGCCCCGTCAAAGTCACCGACGTTTTTGTCGGGCGTGTCTTCCGGCGGCTGATTGACGGCCTGTTTTTCGTCCGCCAGGGAGTAAGCTCCTTTATTACGGCGATAGTGTGCCTGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTCGATGCGGTCATGGGGACGACCG  >  NC_000913/300412‑300881

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: