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.C94.7% 118.7 / 0.0 38intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (13/23);  total (15/23)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTA  >  NC_002947/1126502‑1126787
                                                                                                                                                   |                                                                                                                                                 
ggggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccgc                                                                                                                                             <  1:197788/149‑4 (MQ=255)
  ggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccgctc                                                                                                                                           <  6:79983/149‑6 (MQ=255)
  ggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccgctc                                                                                                                                           <  4:163867/149‑6 (MQ=255)
  ggCCGATGCCGATCTTGAACAGGGCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccgctc                                                                                                                                           <  7:255130/149‑6 (MQ=255)
     cgatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTa                                                                                                                                         <  7:275538/149‑1 (MQ=255)
     cgatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTGGGGGGCACATTGAGCCCTGAAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTa                                                                                                                                         <  5:120640/149‑1 (MQ=38)
        tgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTacac                                                                                                                                      <  3:136460/149‑1 (MQ=255)
              tCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCGCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCCCGAATGCGGAGGTATAGCCACCAACGCATTAGCGGGTAAACCCGGACCTAAACAGAGc                                                                                                                                 >  1:88887/1‑149 (MQ=255)
                  gAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGc                                                                                                                                 <  8:203467/144‑1 (MQ=255)
                  gAACAGGTCGAAGACGCTCAGTGACATGTCGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCg                                                                                                                              <  7:121287/147‑1 (MQ=255)
                        gTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcc                                                                                                                                                <  2:332886/122‑1 (MQ=255)
                        gTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCCAATGCGGAGGTATAGCCACCAAAGCATTAGCCGGTAAACACCCACCTAAAACGAGCCCCCCCCCCaa                                                                                                                     >  8:386262/1‑149 (MQ=255)
                        gTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATAGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcc                                                                                                                                                >  1:332886/1‑122 (MQ=255)
                                 gCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTgg                                                                                                            <  1:304910/149‑1 (MQ=255)
                                  cTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTGGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGc                                                                                                           <  2:88887/149‑1 (MQ=255)
                                        gACATGTGCACCTCCTGCCTGGTTGTTGGGGGGACATTGAGCCCTGTCGGCGCGGTTTTACCCGCGAAGGCGGAGGTATCGCCACCATCCCATCCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGGGGCGGGGa                                                                                                      <  3:190712/148‑1 (MQ=37)
                                         aCATGTGGACCTCGTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAgg                                                                                                    <  7:108750/149‑1 (MQ=255)
                                                 aCCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGc                                                                                                                   <  6:135051/126‑1 (MQ=255)
                                                 aCCTCCTGCCAGGTTGTTGGGGGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAg                                                                                            <  1:99795/149‑1 (MQ=255)
                                                    tcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTacac                                                                                                                                      <  5:395085/105‑1 (MQ=255)
                                                      ctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCTAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTc                                                                                       <  1:356284/149‑1 (MQ=255)
                                                       tGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCCCGAATGCGGAGGTATAGCCACCATCGCAATAGCGGGTAAACCCGCTCCTAACAAGAGCCCCCCACCCCACCCACGTGGCCGGGAGGGCTATCAcgcggc                                                                                       >  2:60470/1‑146 (MQ=25)
                                                           cAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGGTTTACCCCGCAAGGGGGGGTATATGCACCATAGCATTAGCCGGGTAAACCGCGCCGAAAAACACCCCCCCCCCCCCGCCCCGGGGCGGGGGGGGCTATGAGGCGTccgg                                                                                    >  6:241035/1‑146 (MQ=255)
                                                            aGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTAAAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGAt                                                                                  <  6:82409/149‑1 (MQ=255)
                                                            aGGTGGTTGTGTGGACATTGAGCCATGTAGGCGCGGGTTAACCCGCGAATGCGGAGTTATAGCCACAATCGCATTAGCGTGTAAACCCGCTCATACACAGAGCCCGCAACCCAAGCCAGGTGTCGGTGAGGGCTATCAGGCGTCCTGAt                                                                                  <  6:82410/149‑1 (MQ=255)
                                                                      gTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGTATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCTCCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgat                                                                        >  2:271245/1‑149 (MQ=255)
                                                                      gTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGTATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCTCCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgat                                                                        >  2:277362/1‑149 (MQ=255)
                                                                          aCATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAACCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGAc                                                                   <  7:181792/149‑1 (MQ=255)
                                                                                    cTGTAGGCGCGGGTTTAAACCCGCATGAGGAGGTATAGCCCCCAACGCATTAGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCCAGCCAGGGGGCGGGGGGGGCGATAAGGCGGCCGGATAGCGCTCGATCGACCGGCACgcgc                                                         >  2:26944/1‑149 (MQ=255)
                                                                                         ggCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACgg                                                                 <  5:444504/136‑1 (MQ=255)
                                                                                                            aaTGCGGAGGTATAGCCACCATCGTATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCTCCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACCGGGCAGGCGCACACCAGGTTGCGGTCGCCGaa                                   <  1:271245/149‑1 (MQ=39)
                                                                                                            aaTGCGGAGGTATAGCCACCATCGTATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCTCCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACCGGGCAGGCGCACACCAGGTTGCGGTCGCCGaa                                   <  1:277362/149‑1 (MQ=39)
                                                                                                            aaTGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAac                                  <  7:386262/149‑1 (MQ=255)
                                                                                                               gCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACg                               >  2:389566/1‑149 (MQ=255)
                                                                                                                 ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACTGCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGt                              >  1:342734/1‑149 (MQ=39)
                                                                                                                      tataGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGGCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGa                        >  8:307049/1‑149 (MQ=255)
                                                                                                                       ataGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAc                       >  3:320894/1‑149 (MQ=255)
                                                                                                                            caccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGc                  <  7:191922/149‑1 (MQ=255)
                                                                                                                                 tCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                            <  6:35235/118‑1 (MQ=255)
                                                                                                                                 tCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                            >  5:35235/1‑118 (MQ=255)
                                                                                                                                        cgcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc       <  1:300172/148‑1 (MQ=255)
                                                                                                                                            gggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAg    >  6:335050/1‑148 (MQ=255)
                                                                                                                                            gggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAg    >  8:310352/1‑148 (MQ=255)
                                                                                                                                             ggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTa  >  6:443700/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTA  >  NC_002947/1126502‑1126787

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: