Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,556,929 T→C V123V (GTA→GTG)  eutC ← ethanolamine ammonia‑lyase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,556,9290TC96.5% 93.8 / NA 29V123V (GTA→GTG) eutCethanolamine ammonia‑lyase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (7/21);  minor base G (1/0);  total (8/21)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCATCAGAAATAACGACCTGTACGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGGTATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGAC  >  NC_000913/2556793‑2557056
                                                                                                                                        |                                                                                                                               
ggCCATCAGAAATAACGACCTGTACGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACt                                                                                                                              >  2:120987/1‑140 (MQ=255)
 gCCATCAGAAATAACTGCCTGTTCGTCCGGATTGGCACCACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACtt                                                                                                                             <  1:321744/140‑1 (MQ=255)
 gCCATCAGAAATAACGACCTGTACGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACtt                                                                                                                             <  1:811475/140‑1 (MQ=255)
 gCCATCAGAAATAACAACCTGTACGGCCGCATTGGCACCACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACtt                                                                                                                             <  1:395718/140‑1 (MQ=255)
              aCGACCTGTACGTCCGGCTTGGTAACACACTCCGCATTCCGCGCTTCATCAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACTTCCAGCAAGCCCTg                                                                                                                <  2:364214/140‑1 (MQ=255)
                gACCTGTACGTCCGGATGGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACTTCCAGCAAGCCCTgcg                                                                                                              <  1:715123/140‑1 (MQ=255)
                        cGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCAcc                                                                                                      <  2:27131/140‑1 (MQ=255)
                           ccGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCAc                                                                                                   <  2:663632/140‑1 (MQ=255)
                            cGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCAct                                                                                                  >  1:575372/1‑140 (MQ=255)
                            cGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGGCGCTAATTGCAGGGCGGACTTCCAGGAAGCTCTGCGCTTTCACCCAct                                                                                                  >  2:256183/1‑140 (MQ=255)
                                ttGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTc                                                                                              <  2:575372/140‑1 (MQ=255)
                                       cacacTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCAct                                                                                       <  2:494315/140‑1 (MQ=255)
                                            tGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCttt                                                                                  >  1:481959/1‑140 (MQ=255)
                                                 tttCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTa                                                                             <  2:604366/140‑1 (MQ=255)
                                                         cTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGGTATCt                                                                     >  1:721138/1‑140 (MQ=255)
                                                          ttCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGGTATCtt                                                                    <  2:850180/140‑1 (MQ=255)
                                                                 gCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGGTATCTTTCGAACg                                                             <  2:681562/140‑1 (MQ=255)
                                                                         acacaGGCGGCGGCCAATATCCCGGCGAGCCAGGTAGTGATTTTTGTCGCTAATTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGGTATCTTTCGAACGAGAGTGAt                                                     <  1:256183/140‑1 (MQ=255)
                                                                                 ggcggcCCATATCCGAGCGAGTAAGTTAGAGATTTTTTTCGCTAATTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGGTATCTTTCGAACGAGAGTGATCGGCCAgg                                             <  1:323094/140‑1 (MQ=255)
                                                                                         atatCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGGTATCTTTCGAACGAGAGTGATCGGCCAGGAAACgcag                                     <  2:126488/140‑1 (MQ=255)
                                                                                         atatCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGGTATCTTTCGAACGAGAGTGATCGGCCAGGAAACgcag                                     <  1:513335/140‑1 (MQ=255)
                                                                                                 gCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGGTATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCt                             >  2:548839/1‑140 (MQ=255)
                                                                                                     gTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGGTATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGt                         >  1:836331/1‑140 (MQ=255)
                                                                                                                        cGCTAATTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGGTATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCa      <  1:120987/140‑1 (MQ=255)
                                                                                                                         gCTAATTTCAGAGCGCCTTTCCCGCAAGCCCTGCGTTTTTACCCACTCATCCGGCACTTCTTTCAGTACGGTATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCAc     <  2:398911/140‑1 (MQ=255)
                                                                                                                          cTAATTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGGTATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACg    <  2:721138/140‑1 (MQ=255)
                                                                                                                           tAATTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGGTATCTTTCGAACGAGAGTGATCGGCCAGGAAACTCAGCAGCACCTACGCACGTGGTCGCGGACCGGCACGa   >  1:398911/1‑140 (MQ=255)
                                                                                                                            aaTTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGGTATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGAc  <  2:47044/140‑1 (MQ=255)
                                                                                                                            aaTTTCAGAGCGCACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGGTATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGAc  <  2:481959/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                               
GGCCATCAGAAATAACGACCTGTACGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGGTATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGAC  >  NC_000913/2556793‑2557056

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