Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,152,494 G→T A87S (GCG→TCG)  metC → cystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,152,4940GT93.8% 42.6 / ‑5.6 16A87S (GCG→TCG) metCcystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (14/1);  minor base C (0/1);  total (14/2)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

AGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAAT‑GGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGG‑CGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACT  >  NC_000913/3152249‑3152679
                                                                                                                                                                                                                                                       |                                                                                                                                                                                         
aGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAAT‑GGAGAGTTGTTCTATGGACGGCGCGGAACGTTAAACCATTTCTCCTTACATCAAGCGATGTGTGAACTGGAAGGTGG‑CGCAGGCTGCGTGCTATTTCCCTGCGGgtcggc                                                                                                                                                                                       >  2:207449/1‑250 (MQ=255)
            aaCTGGTGAATGCAGGACGCAGCAAAAAATACAATCTCGGCGAGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAACACCACGCGACACGTAATCGCGCCAAT‑GGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGG‑CGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAAt                                                                                                                                                                          >  1:207943/1‑251 (MQ=255)
                                  caAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAAT‑GGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGG‑CGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAAc                                                                                                                                                    >  1:55382/1‑251 (MQ=255)
                                        aTACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAAT‑GGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGG‑CGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCg                                                                                                                                              >  1:63158/1‑251 (MQ=255)
                                         gacacTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCACAAAACACGCGACACGTAATCGCGCCAAT‑GGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGG‑AGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGGTGCTACTTCCATTCTTGCTTTTATCGAACAGGGCg                                                                                                                                              >  2:46002/2‑250 (MQ=255)
                                                                       cAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAAT‑GGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGG‑CGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATg                                                                                                                >  1:169241/1‑250 (MQ=255)
                                                                               ttcttcGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAAT‑GGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGG‑CGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTc                                                                                                       >  2:203933/1‑251 (MQ=255)
                                                                               ttcgtggCTGGTCTTTGACAGTGCAGAAGCTAAAAAACCCGCGACACGTAATCGCGCCAAT‑GGAGAGCTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACACGCGATGAGGTCACTGGAAGGGGG‑CGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGATCAACACCGCCTATGAACAGAGTc                                                                                                       >  2:158413/7‑251 (MQ=255)
                                                                                     gCTGGTCTTTGCCAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAAT‑GGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGG‑CGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGAtt                                                                                                 >  1:61362/1‑251 (MQ=255)
                                                                                                         aacccAAAAAACCCGCCACACTTAGTCGCGC‑AATATGAGTGTTGTTCTATGGGCGGCGGGGAGCGTTACCCCATGTCTCCTTCCAACAACCGATGTGT‑AAATTGTAGGTGGGCGAGGGCTGCGTTTTATTTCCCTGCGGGCCGGCTGCGGTTGCTAATTCCATTCTTCCTTTTATCGAACAGGGCATTCATTTTTTGATCACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGc                                                                             <  2:219853/248‑1 (MQ=255)
                                                                                                                       gcgACACGTAATCGCGCCAAT‑GGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGG‑CGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCAATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGAAAAATACTCAGCAAACTGGGCGTaa                                                                >  2:126689/1‑250 (MQ=255)
                                                                                                                                                    tgttCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGG‑CGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGt                                   >  2:84712/1‑251 (MQ=255)
                                                                                                                                                    tgttCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGG‑CGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTAATCAAACAGGGCGATCATGTGTTGATGACAAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTgg                                    <  1:123022/250‑1 (MQ=255)
                                                                                                                                                                        aCGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGG‑CGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCt               >  2:178306/1‑251 (MQ=255)
                                                                                                                                                                          gTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGG‑CGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGc             >  2:107673/1‑251 (MQ=255)
                                                                                                                                                                                    ttCTCATTACAACAAGCGATGTGTGAACTGGAAGGTGG‑CGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAacac   >  1:161401/1‑251 (MQ=255)
                                                                                                                                                                                     tctcCTTACAACAAGCGATGTGTGAACTGGAAGGTGG‑CGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCACACACt  >  2:21913/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                         
AGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAAT‑GGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGG‑CGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACT  >  NC_000913/3152249‑3152679

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