Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 299,256 G→A A560V (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_000913299,2560GA93.8% 35.6 / ‑3.5 16A560V (GCC→GTC) paoCaldehyde dehydrogenase: molybdenum cofactor‑binding subunit
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (7/8);  total (8/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-01

CCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTC  >  NC_000913/299014‑299491
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           
cctccATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGACCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGACTTCATGt                                                                                                                                                                                                                                     >  1:191559/1‑251 (MQ=255)
                     tgcccatAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAATGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGACTTCATGTAGCGTgggg                                                                                                                                                                                                                            <  1:415483/239‑3 (MQ=255)
                       cccatAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGACCTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGTCTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGACTTCa                                                                                                                                                                                                                                        >  1:448926/1‑225 (MQ=255)
                                   gcgcCAATGCCCTGGTTGCGCGCAGTTTTCGGAGTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGGCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTTCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTTCCGCCTGCGGTGACTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGc                                                                                                                                                                                                   <  1:317587/250‑1 (MQ=255)
                                                             tttCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGACTTCATg                                                                                                                                                                                                                                      <  1:147179/189‑1 (MQ=255)
                                                                                         cacacaGCGTGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGACTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTc                                                                                                                                            >  1:95479/1‑251 (MQ=255)
                                                                                                cGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAACTGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGACTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAGACTGCGACTGCTCAGGATCAAACCCGAATGCCGAGGCAATCATTTCGCGAAGc                                                                                                                                     >  1:338602/1‑251 (MQ=255)
                                                                                                     nnaCNNNNNACCNGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGACTTc                                                                                                                                                                                                                                         <  2:268536/144‑1 (MQ=255)
                                                                                                                                gcgcTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGACTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTGGCCGTc                                                                                                                                                                                               <  1:120342/161‑1 (MQ=255)
                                                                                                                                       gCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCg                                                                                                                                                                >  1:215767/1‑185 (MQ=255)
                                                                                                                                                            gCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCGCAGTGTTCCGAAGTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGACTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAAGCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTcn                                                                                             <  1:51560/231‑2 (MQ=255)
                                                                                                                                                                          gNCTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGACTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAggc                                                                                                                          >  1:348367/1‑188 (MQ=255)
                                                                                                                                                                              gctgGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCCGTCTGCCGCCTGCGGTGACTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAAc                                                                                                                                                                   >  1:176698/1‑143 (MQ=255)
                                                                                                                                                                                                                                 ntCTGCCGCCTGCGGTGACTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTACTCTTGACGTCCGCAAACTGAGACTGCTAAGGATCAAAACCGACCGCCGAGGCAATCATTTCGCGAACCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGGATTCGCGCCCCATTGCCCACCAGAACCCGCAGAAACCGGGAAACTGGAATAGCCGAGGTGAACCGCCACCTGCTCCAGCGGTACGCCAAGCAtt    >  1:53734/2‑251 (MQ=255)
                                                                                                                                                                                                                                 gTCTGCCGCCTGCGGTGACTTCATGTAGCNTCGNNCTNNGGNTACCg                                                                                                                                                                                                                <  2:256002/47‑1 (MQ=25)
                                                                                                                                                                                                                                     gccgccTGCGGTGACTTCATGTAGCGTGGCGCTTCGGTTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCAtttn  <  1:55366/249‑2 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           
CCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACCAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTC  >  NC_000913/299014‑299491

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