Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,126,647 +C intergenic (‑101/+48) tdcG‑II ← / ← gcvP‑I L‑serine dehydratase/glycine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,126,6451.C80.0% 48.4 / 9.5 20intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (4/0);  new base C (5/11);  total (9/11)
Fisher's exact test for biased strand distribution p-value = 2.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGA  >  NC_002947/1126501‑1126774
                                                                                                                                                     |                                                                                                                                 
aGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGGGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAc                                                                                                                                  <  3:134281/146‑1 (MQ=255)
  gggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccgc                                                                                                                              <  8:68534/148‑4 (MQ=255)
    gccgatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTAACCCGAATGGGGTGGTATAGCCCCCAAAGCATTCCCGGGGTAACCCGcac                                                                                                                             >  7:76586/1‑146 (MQ=255)
            cgatCTTGAACCGGGCGAAGACGCTCAGTGACCTGGGCACCTCCTGCCCGGTTGTGGGGTGCAAATTGAGCCCTGTCGGCGGGGGTTTACCCGCGAATGGGGAGGTATAGCCACCATCGCATGCGCGGGTAAACCCGCTCCTAcaca                                                                                                                      <  6:94982/147‑1 (MQ=25)
                    aaCAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCGAATTGGGAGGTTTTGCCACCATCGCAATAGCCGGTTAAACCCCGCCCACACACACcccccc                                                                                                             >  1:112175/1‑148 (MQ=255)
                      cAGGTCGAAGACGCTCAGTTACATGTGGACCTCCGGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCAc                                                                                                           <  7:29279/148‑1 (MQ=255)
                                  gCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGCCATTGAGCCCTGTAGGCGCGGGGTTTCACGCGGATTCGGGGGGTTTGGCACCCGCGCATTAGCGGGTTAACCCGGTGCTGGACCTCGCCCGCCACCCAAAGCAGGTg                                                                                               >  8:125642/1‑148 (MQ=255)
                                           cATGGGGACCTCCTGCCAGGTGGGTGGGTGGACATTGCGCCCTGTAGGCGCGGGTTTACCCGCGAAGGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAg                                                                                       <  7:12744/147‑1 (MQ=37)
                                                        ggCCAGGTTGTTGGGTGGACAGTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATACGCGGGTAAACCCGCGCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTc                                                                         <  5:131879/147‑1 (MQ=38)
                                                             ggtggggggggggACATTGAGCCCTGTCGGCGCGGGTTTACCCGCGAAGGCGGAGGAATAGCCACCAGCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGAt                                                                    <  5:103381/141‑1 (MQ=38)
                                                             ggTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCCCGAAATCGGAGGTATAGCCCCCATAGCATTAGCCGGTAAAACCGCTCCTACACCGAGCCCCCCACCCAACCCAGGTGGCGGGGAGGGCTATAAGGCggacgg                                                                      >  4:131483/1‑141 (MQ=255)
                                                               ttgttgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAg                                                                  <  4:133547/148‑1 (MQ=255)
                                                               ttgtggGGGGGACATTGAGCCCTGTCGGCGCGGGTTTACCCGCGAATGCGGAGGTATACCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGACGCCTGATAg                                                                  <  2:123137/148‑1 (MQ=37)
                                                                   ggggTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctct                                                              <  8:118352/147‑1 (MQ=255)
                                                                     ggTGGACATTGAGCCCTGTAGGCGCGGGTTTTAACCCGAATTGGGAGGTATATCCCCCAAAGCATTAGCGGGTAAACCCCCTCCGACACACAGACCCCCCCCCCAGCCACGGGGCGGGGAGGGCTATAAGGCGGCCTGATAGcgatcg                                                            >  2:19947/1‑143 (MQ=255)
                                                                             ttGAGCCCTGTAGGCGCGGGTTTCCCCGCGAATGCGGAGGTATCGCCACCAGCGCATGCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACgg                                                    <  1:42633/148‑1 (MQ=38)
                                                                                   ccTGTAGGCGGGGGTTTACCCGCGAATGCGTAGGTATAGCCACCAGCGCCTGCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGGCCTGATAGCTCTCGATCGACGGGCAGgc                                              <  1:114808/148‑1 (MQ=38)
                                                                                                      ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCg                           >  3:43085/1‑148 (MQ=255)
                                                                                                                 gAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt                >  3:84786/1‑148 (MQ=255)
                                                                                                                   gttacccccaccaCCCCATCCGGGGTAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGGGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGt              <  2:181284/141‑1 (MQ=255)
                                                                                                                            accaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGc     >  1:19271/1‑148 (MQ=255)
                                                                                                                               aTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGa  >  8:42774/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                 
AGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGA  >  NC_002947/1126501‑1126774

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