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.C79.2% 59.7 / 11.1 24intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (5/0);  new base C (8/11);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 4.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGG  >  NC_002947/1126503‑1126781
                                                                                                                                                    |                                                                                                                                        
gggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCGCATGCGGAGGTATAGCCACCAACGCCTTAGCCGGTaaaacccc                                                                                                                                      >  2:20724/1‑143 (MQ=255)
    cgatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGGTTTAACCCGCATTCGGGGGGTTAGCCCCCATAGCATTAGCCGGGTAACCCGCGcc                                                                                                                                   >  1:106327/1‑147 (MQ=255)
         ccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACACCGAAATCGGAGGGTTATCCCACAACGCATTAGCCGGGAAAACCCCACCGacac                                                                                                                              >  6:63431/1‑147 (MQ=255)
           gatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGAGGTTTAGCCACCCTCGCCATAGCGGGGTAACCCCCTCCGACACaga                                                                                                                           >  3:32565/1‑148 (MQ=255)
                           aaGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGCGGTTTAGCCACCATCGCAATAGCGGGGACACCCCCTGCGACACAGAGcccccccccccaccc                                                                                                           >  7:50523/1‑143 (MQ=255)
                                          aTGGGGACCGCCTGCCCGGGTGTGGGGTGGACATTGAGCCCTGTAGGCGCGGGTTAACCCGCGAATGCGGAGGTATAGCCACCATCGCATGCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg                                                                                            <  3:68795/148‑1 (MQ=37)
                                               gACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGCAAGCGGAGGTTTAGCCCCCATCGCATTAGCGGGTAAACCCCCGCCTACACAGAGCCCGCCCCCCCACCCCGGTGGCGGGGCGGGCTATc                                                                                       >  3:153522/1‑148 (MQ=255)
                                                 cctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAAACCCCAATTCGGAGGTATAAACACCAACGCAATAGCCGGTAAACCCCCACCAAAAAAAAGCCCCCCACCCAACCCACCTGGCGGGGACGGCTAtaac                                                                                     >  6:67157/1‑145 (MQ=255)
                                                       gCCAGGTTGTTGGGTGGACATGCACCCCTGCAGCCGCGGGGTTACCCGCGAATGCGGAGGTATCCCCACCATCGCAGGCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTcc                                                                               <  7:37986/148‑1 (MQ=37)
                                                         cAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGGTTTACCCCCGAAGAGGAGATATAACCACCAAACCAATAGCCGGTAAAAACCCGCCAACACAAACCCCCCCCCCCCACCCACGGGGGCGGGAGGGCTATAAAGCgg                                                                                 >  4:122027/1‑143 (MQ=255)
                                                                 tgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGAGCTATCAGGACGCCTGATAGccct                                                                     <  1:20724/148‑1 (MQ=37)
                                                                               gCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCa                                                        <  7:86247/147‑1 (MQ=255)
                                                                                     tAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGcac                                                 <  8:69885/148‑1 (MQ=255)
                                                                                         cgcgGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCa                                             <  5:63431/148‑1 (MQ=255)
                                                                                           cgGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAgg                                           <  7:85951/148‑1 (MQ=255)
                                                                                                     ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCCCGATCGACGGGCAGGCGCACCCCAGGTTGCGGTCg                                  >  3:21213/1‑148 (MQ=255)
                                                                                                      cgcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCt                                                                                          <  8:44497/91‑1 (MQ=255)
                                                                                                                gAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt                       <  5:67157/148‑1 (MQ=255)
                                                                                                                    tatcGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAATCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGTCAGGCGCACACCAGGTTGCGGGCGCCGAACACTTTGTc                    <  8:24023/147‑1 (MQ=255)
                                                                                                                         ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACACg              >  7:98087/1‑148 (MQ=255)
                                                                                                                           accaTCGGATTCGGGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACTGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGc            <  5:95583/148‑1 (MQ=255)
                                                                                                                           accaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGc            >  8:10242/1‑148 (MQ=255)
                                                                                                                                   aTTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGCCcggc    >  5:116801/1‑148 (MQ=255)
                                                                                                                                    ttCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcg   >  3:51416/1‑148 (MQ=255)
                                                                                                                                       gcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcgg  >  8:117029/1‑146 (MQ=255)
                                                                                                                                                    |                                                                                                                                        
GGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGG  >  NC_002947/1126503‑1126781

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: