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.C82.0% 105.7 / 20.7 39intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (7/0);  new base C (13/19);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 8.32e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGC  >  NC_002947/1126499‑1126782
                                                                                                                                                      |                                                                                                                                          
ggAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCGCGAAAGCGGAGGTTTTGGCACCCTCGCATTAGCGGGTaaaac                                                                                                                                           >  2:287794/1‑147 (MQ=255)
     ggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGATTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccgctc                                                                                                                                    <  1:27311/149‑6 (MQ=255)
             ccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGAGGGATAGCCCCCATCGCATTAGCGGGGTAACCCGCGCCGaca                                                                                                                                >  2:54245/1‑146 (MQ=255)
               gatCTTGAACAGGGCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACaga                                                                                                                            <  3:221597/148‑1 (MQ=255)
                     gAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCGAAAGGGGGGGTATAGGCACCCTAGCCTTAGCGGGGTAACCCCCGCCGACCCAGAGCCCGcc                                                                                                                     >  5:14218/1‑149 (MQ=255)
                     gAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACagag                                                                                                                           <  7:364447/143‑1 (MQ=255)
                     gAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCCAATGCGGAGGTATAGCCACCAACGCAATAGCCGGTAAACCCCCGACGACACAGACCCCGcc                                                                                                                     >  5:557/1‑149 (MQ=255)
                           gTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGAGGTATAGCCACCATCGCATTAGCGGGTAAACACGCTCCGACACAGAGCCCGCCACCCaa                                                                                                               >  7:256968/1‑149 (MQ=255)
                            tCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCCCGAATGCGGGGGTATAACCACCCACGCCATAGCGGGGAAAACCCCGCCGACAAAAAACCCCCCAccccac                                                                                                              >  6:119312/1‑146 (MQ=255)
                                        agggACATGTGGACCTCCTGCCAGGTTGTTGGGGGGACATTCAGCCCTGTAGGCGCGGGTTAACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCg                                                                                                    <  3:315543/144‑1 (MQ=38)
                                          tGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCCAATGCGGAGGGATAGCCACCAAAGCAATAGCGGGGAAAACCCCGCCAACAAAAAGCCCCCCCCCCCACCCAcggggcggggg                                                                                                >  4:116102/1‑148 (MQ=255)
                                           ctcATGTGGACCTCCTGCCAGGTTGTTGGGGGCACATTGAGCCCTGAAGGCGCGCGTTTACCCGCGAATGCGGAGGTATACCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAg                                                                                               <  7:46307/147‑1 (MQ=38)
                                                    aCCTCCTCCCATGTGTGTGGGTGGACATTGAGCCCTGAAGGCCCTGGTTTACCCGCGAATGCGGTGGTATAGCCTCCACCGCATTCTCGGGTAAACCCGCGCCAACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAg                                                                                      <  6:356675/149‑1 (MQ=255)
                                                       tcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGc                                                                                            <  7:307886/140‑1 (MQ=255)
                                                        cctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGAGGTATAGCCACCATCGCATTAGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGTGGGGAGGGCTATCAGGCGt                                                                                  >  1:275072/1‑149 (MQ=38)
                                                         ctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGAGGTATAGCCACCATAGCATTAGCGGGTAAAACCGCTCCTACACAGAGCCCGCCACCCAAACCAGGTGGCGGGGAGGGCTATAAGGCGTc                                                                                 >  7:83745/1‑149 (MQ=255)
                                                           gCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGCATGCGGGAGGATAGACCCCAACGCATTGCACGGGTAAACACCGAAAACAAAGACCCCCCCCCCCAACCCACCAGGACGGGACGGCTATCATGccgact                                                                               >  6:304986/1‑144 (MQ=255)
                                                                gttgttGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGACGCCTGATAg                                                                          <  3:101588/149‑1 (MQ=38)
                                                                      gggTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCCGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCt                                                                                           <  7:165037/126‑1 (MQ=255)
                                                                           gACATTGAGCCCTGTAGGCGCGGTTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGa                                                                             <  1:263555/135‑1 (MQ=255)
                                                                                   gCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAgg                                                       <  5:119312/149‑1 (MQ=255)
                                                                                     ccTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgcg                                                     <  1:54245/149‑1 (MQ=255)
                                                                                                   tttACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                       <  1:287794/149‑1 (MQ=255)
                                                                                                       cccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgatcg                                                                >  7:294483/1‑120 (MQ=255)
                                                                                                       cccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACgg                                                            >  3:248598/1‑124 (MQ=255)
                                                                                                       cccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACgg                                                            <  1:34040/124‑1 (MQ=255)
                                                                                                       cccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACgg                                                            >  2:34040/1‑124 (MQ=255)
                                                                                                       cccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACgg                                                            <  4:248598/124‑1 (MQ=255)
                                                                                                        ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGAACACCAGGTTGCGGTCGc                                  >  7:319117/1‑149 (MQ=255)
                                                                                                                gCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgatc                                                                 <  3:110995/110‑1 (MQ=255)
                                                                                                                gCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgatc                                                                 >  4:110995/1‑110 (MQ=255)
                                                                                                                 cGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACg                          <  3:116102/148‑1 (MQ=255)
                                                                                                                  ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAacac                           <  6:202049/146‑1 (MQ=255)
                                                                                                                  ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAacac                           >  5:202049/1‑146 (MQ=255)
                                                                                                                  ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt                        >  1:52175/1‑149 (MQ=255)
                                                                                                                            ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg              >  7:168561/1‑149 (MQ=255)
                                                                                                                               ccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGcc            >  4:166763/1‑148 (MQ=255)
                                                                                                                                   cGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcg       >  6:142020/1‑149 (MQ=255)
                                                                                                                                      aTTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcg    <  5:304986/149‑1 (MQ=255)
                                                                                                                                          cgcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc  <  7:285816/148‑1 (MQ=255)
                                                                                                                                              gggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgcgc                                                    >  6:108329/1‑95 (MQ=255)
                                                                                                                                                      |                                                                                                                                          
GGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGC  >  NC_002947/1126499‑1126782

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