Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,126,645:1 +C 100% intergenic (‑99/+50) 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.9% 94.9 / 15.9 35intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (6/0);  new base C (11/18);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 7.62e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAG  >  NC_002947/1126498‑1126785
                                                                                                                                                         |                                                                                                                                             
tGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCGCGAATGCGGAGGTATAGCCCACATGGCATTCGCGGGGaaaa                                                                                                                                               >  5:282168/1‑148 (MQ=255)
  gAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAccc                                                                                                                                           <  5:341016/149‑2 (MQ=255)
   aGGGGCCGATGCCGATCTTGAACCGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccg                                                                                                                                          <  2:352809/149‑3 (MQ=255)
      ggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccgctc                                                                                                                                       <  2:77122/149‑6 (MQ=255)
               cgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAc                                                                                                                                              >  6:65157/1‑134 (MQ=255)
               cgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAc                                                                                                                                              <  5:65157/134‑1 (MQ=255)
                    ttGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCCAAATCGGATGTTTTTCCACCATAGCATTACCCGGGAAAACCCCTCCGAAAAAAAGCCCg                                                                                                                          >  7:194014/1‑149 (MQ=255)
                       aaCAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGCATGCGGAGGTATAGCCACCATCGCATTAGCGGGTAAACCCGCTCCTACACCGAGCCCGccc                                                                                                                       >  6:269869/1‑148 (MQ=38)
                            gTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGGTTTACCCCGAATGGGGGGGGTTAGGCCCCCTCGCCTTAGCGGGGGAACCCCCGCCTGCACACAGccccccccccca                                                                                                                  >  2:31213/1‑147 (MQ=255)
                             tCGAAGACGAGCAGTTACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTAACCCGCTAATGCGGAGGTATAGCCACCATCGCATTCGGGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCaa                                                                                                                  <  8:194014/148‑1 (MQ=255)
                                   aCGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGAGGTATAGCCACCATAGCATTAGCGGGTAAACCCCCTACGACACAAAGCCCCCCACCCAACCCACGt                                                                                                           >  5:73916/1‑149 (MQ=255)
                                    cGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTg                                                                                                          <  5:269869/149‑1 (MQ=255)
                                       tCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTCCCCCGAATTCGGAGGTATAGCCACCATCGCATTAGCGGGTAAACCCGCTACTACACACAGCCCGCCACCCAAGCCAGGGGGCg                                                                                                       >  4:112426/1‑149 (MQ=255)
                                             aCATGTGGACCTCCTGCCAGGTTGTTGGGGGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACACAGCCCGCCACCCAAGCCAGGTGGCGGGGAg                                                                                                  <  1:31213/148‑1 (MQ=255)
                                              cATGTGCACCTCCTGCCCGGTTCTGGGCTGCACAATGCGCCCGGTAGCCGCGGGTTAACCCGCGAATGCGCAGGTAAACCCACCACCGCATTCGGGGTAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg                                                                                                <  3:6481/149‑1 (MQ=255)
                                                tgtgGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGAGGTATAGCCACCATCGCATTAGCGGGGAAACACGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGAGc                                                                                               >  7:341669/1‑148 (MQ=255)
                                                  tgGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTa                                                                                             <  6:172528/148‑1 (MQ=255)
                                                   ggACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCCAATTCGGTGGTATAGCCACCATAGCATTAGCGGGTTAACCCGCTCCTACACAGAGCCCCCCACCCAACCCAGGTGGCGGGGAGGGCTAt                                                                                            >  3:62145/1‑148 (MQ=255)
                                                    gACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCa                                                                                          <  5:141315/149‑1 (MQ=255)
                                                       ctcctcCCTGGTGGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTCCCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGc                                                                                       <  5:206198/149‑1 (MQ=38)
                                                       ctcctGCCAGGTTGTGGGGGGGACATGGCGCCCTGTAGGCGCGGGTTTACCCGCGAAGGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGc                                                                                       <  5:211342/149‑1 (MQ=37)
                                                           tGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCCCGCATGCGGAGGTATAGCCCCCAAAGCCTTAGCGGGGTAAACCGCGCCGACCCAGAGCCCCCCCCCCAACCCACGGGGCGGGGAGGGCAAAAAGGCGtgc                                                                                   >  8:109484/1‑147 (MQ=255)
                                                             cccggtggtTGGGTGCACATCGAGCCCTGTAGGCGGGGTTTTCCCCGCGAATGCGGCGGTATAGCCACCATCGCATTCGCGGGAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTcc                                                                                   <  4:177028/144‑1 (MQ=255)
                                                             ccAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCt                                                                                  <  1:92890/148‑1 (MQ=255)
                                                                  ttgttgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGc                                                                            <  2:205090/149‑1 (MQ=255)
                                                                    gttgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCGAAAGCGGGGGGATAGCCACCATCGCAATAGCGGGGAAAACCCCGCCGACACACAGCCCCCCCCCCCAGCCCGGGGGGGGGGGGGGCGATGACGCGGCCGGATAg                                                                             >  8:67913/1‑146 (MQ=255)
                                                                          tGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAg                                                                                         <  1:8498/128‑1 (MQ=255)
                                                                            taCATTGAGCCCTGTAGGGGGGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgatcga                                                                  <  2:91709/148‑1 (MQ=255)
                                                                                            ggCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACAcc                                                  <  6:73916/149‑1 (MQ=255)
                                                                                                      tttACCCGCTAATGCGGAGGTAAAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                          <  6:282168/149‑1 (MQ=255)
                                                                                                           ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                     >  7:59639/1‑149 (MQ=255)
                                                                                                           ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                     >  3:90196/1‑149 (MQ=255)
                                                                                                             gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc                                    >  3:142495/1‑148 (MQ=255)
                                                                                                                     ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt                           <  7:67913/149‑1 (MQ=255)
                                                                                                                               ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGACCCCGTTGTCGACCCgg                 >  8:333351/1‑149 (MQ=255)
                                                                                                                               ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg                 >  1:20450/1‑149 (MQ=255)
                                                                                                                                  ccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGcc               >  3:72238/1‑148 (MQ=255)
                                                                                                                                        cATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggc        >  3:11689/1‑149 (MQ=255)
                                                                                                                                         aTTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcg       >  2:38589/1‑149 (MQ=255)
                                                                                                                                             cgcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc     <  4:72238/148‑1 (MQ=255)
                                                                                                                                               cgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAg  <  4:62145/149‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                             
TGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAG  >  NC_002947/1126498‑1126785

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: