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

TGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCA  >  NC_002947/1126498‑1126784
                                                                                                                                                        |                                                                                                                                           
tGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATTCGGAGGTATATCCACCATCGCCATAGCGGGTaaa                                                                                                                                              >  7:112284/1‑148 (MQ=255)
tGGAGGGGCCGATGCCGATCTTGAACAGGGCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTaaa                                                                                                                                              <  7:105859/148‑1 (MQ=255)
      ggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccgct                                                                                                                                       <  2:165073/148‑5 (MQ=255)
       gccgatgccgatCTTGAACCGGTCGAAGACGTTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccgctc                                                                                                                                      <  7:1427/148‑6 (MQ=255)
                       aaCAGGTCGAAGACGCTCAGGGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGc                                                                                                                        <  2:31871/147‑1 (MQ=255)
                             ccgacgaCGCGCAGTGACATGGGGACCTCCGGCCAGGTGGGGGGGTGGACATGCAGCCCTGGAGGCGCGGGTTTACCCGCGAAGGCGGAGGTAAAGCCACCACCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCaa                                                                                                                 <  3:16826/147‑1 (MQ=21)
                               gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCGAAAGGGGGGGTTTAACCCCCAAAGCAGTAGCGGGTAAAACCCCGCCGACACAGAGCCCCCCCCCCCAg                                                                                                                >  4:74653/1‑147 (MQ=255)
                                  gACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAg                                                                                                            <  3:132305/148‑1 (MQ=255)
                                   aCGCGCAGTGACATGTGGACCTCCTGCCAGGTTGTGGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAgg                                                                                                           <  7:160887/148‑1 (MQ=255)
                                                tggggACCTCCTGCCAGGTTGTGGGGGGCACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGc                                                                                              <  3:62640/145‑1 (MQ=37)
                                                   ggACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTAt                                                                                           <  6:69813/148‑1 (MQ=255)
                                                              cAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTcc                                                                                  <  3:146609/146‑1 (MQ=255)
                                                                         gTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCGAATGGGGGGGTATAGCCACCAACGCATTAGCGGGTAAACCCCCGCCGACACACAGCCCCCCCCCCCAGCCAGGGGGCGGGGGGGGCGATAAGGCgggcgga                                                                               >  2:121347/1‑132 (MQ=255)
                                                                          tGGACATTGAGCCCTGTAGGCGCGGGTTTTAACCCCCATGGCGGGGGTTAGCCACCCACGCAATAGCGGGTGAACCCCCGCCGACACACACCCCCCCCACCCACCCCGGTGGGGGGGCGGGGTATCACGCGGCCTGGTAGccggcgca                                                                    >  8:129290/1‑141 (MQ=255)
                                                                          gggACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGAt                                                                              <  1:121347/137‑1 (MQ=255)
                                                                                           aGGCGGGGGTTTACCCGCGAATGGGGCGGGAAACCCCCCCACGCAGTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGca                                                     <  3:128993/146‑1 (MQ=255)
                                                                                           aGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGGCCTGATAGCTCTCGATCGACGGGCAGGCGcaca                                                   <  8:112284/148‑1 (MQ=255)
                                                                                               gcgGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAACCCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCCCACCAg                                               >  1:40850/1‑148 (MQ=255)
                                                                                                   gTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTTAACCCGCTGCTACCCAGAGCCCGCCACCCACGCCAGGTGGCGGGGAGGGCTCTCAGGCGTCCTGATCGCTCTCGATCGCCGGGCAGGCGCACCCCAGGTTg                                           >  5:703/1‑148 (MQ=255)
                                                                                                          ccGCGAATGCGGAGGTATAGCCCCCATCGCATTCGCGGGTAAACCCGCTCCTCCACAGAGCCCGCCCCCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGCCGGGCAGGCGCACACCAGGTTGCGGTCg                                     >  4:41298/1‑148 (MQ=255)
                                                                                                          ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCg                                     >  6:103092/1‑148 (MQ=255)
                                                                                                           cgcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAACCCCGCTCCTACACAGAGCCCGCCCCCCACGCCAGGGGGCGGGGAGGGCTATCAGGCGTCCTGCTAGCTCTCGAGCGACGGGCAGGCGCCCCCCAGGTTGCGGgcgc                                    >  5:127406/1‑148 (MQ=255)
                                                                                                                          atccccaccaTCCCATGCGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGa                     <  6:20402/144‑1 (MQ=255)
                                                                                                                              ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGCCGGGCAGGCGCACACCACGTTGCGGTCGCCGAACACGTAGTCGACCCg                 >  2:91364/1‑148 (MQ=255)
                                                                                                                                       cATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACtgg        >  3:46154/1‑148 (MQ=255)
                                                                                                                                         ttCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGGCGCCGAACACGTTGTCGACCCGGCCGACcggcg      <  8:85112/148‑1 (MQ=255)
                                                                                                                                            gcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc    <  7:148138/147‑1 (MQ=255)
                                                                                                                                              gggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGCCCGTCGGCCa  >  8:58148/1‑147 (MQ=255)
                                                                                                                                                        |                                                                                                                                           
TGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCA  >  NC_002947/1126498‑1126784

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: