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.C84.0% 68.1 / 8.4 25intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (4/0);  new base C (5/16);  total (9/16)
Fisher's exact test for biased strand distribution p-value = 9.96e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGG  >  NC_002947/1126500‑1126781
                                                                                                                                                     |                                                                                                                                         
gAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCGCCAATGGGGACGTTTAGCCCCCATCGCCTTAGCGGGTaaaac                                                                                                                                          >  4:110642/1‑146 (MQ=255)
         atgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTAACGCGAATTCGGTGGGATATGCACCATCGCATTAGCCGGTAAAACCGCTCCTac                                                                                                                                 >  7:131238/1‑148 (MQ=255)
            ccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCCGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGTTTTACCCGCTAATGCGGAGGTATAGCCACCATCGCATTCGCGGTTAAACCCGCTCCTAcaca                                                                                                                              <  6:96690/148‑1 (MQ=255)
                     aaCCGGGCGCAGACGCTCTGGGACATGTGGACGTCCTGCCAGGTTGTGGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATCGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACaga                                                                                                                            <  8:110386/141‑1 (MQ=25)
                     aaCAGGGCGAAGACGCTCAGGGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACaga                                                                                                                            <  7:74399/141‑1 (MQ=255)
                                   gCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTAACCCGCAAGGGGAGGTATATGCCCCAAAGCAATAGCCGGTAAAACCCCGGCTACACAACGCCCCCCCCcccacccacggg                                                                                                      >  4:127285/1‑144 (MQ=255)
                                     tCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGGTTTACCCCCAATGCGGGGGGTTAGCCCACATCGCCCTAGCGGGGTAACCCCCGGCGACACCCAGCCCCCCCCCCCACCCACGTGGc                                                                                                    >  5:89783/1‑148 (MQ=255)
                                     tCAGTGACATGTGGACCTCCTGCCAGGTTGTGGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGc                                                                                                    <  5:39766/148‑1 (MQ=255)
                                                      tcctGCGAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTCCCCGCGAATGCGGCGGTATAGCCACCATCGCATTCGGGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGc                                                                                   <  2:102024/148‑1 (MQ=38)
                                                          gCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCGCGAATTAGGTGGTATAGCCCCCATAGCATTAGCCGGTAAAACCGCGACAAAACAAAGCCCCCACCCCCAAGCACGTGGCGGGGTGGGGGATCAGGCGt                                                                                 >  4:82830/1‑146 (MQ=255)
                                                            cagggggttgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTg                                                                             <  3:127285/142‑1 (MQ=255)
                                                               gttgttGGGTGGACATTGAGCCCTGTCGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATa                                                                          <  3:110642/148‑1 (MQ=255)
                                                               gttgtgGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAAGGCGGAGGTATAGCCACCATCGCATTCGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGCGGAGGGCTATCAGGCGTCCTGATa                                                                          <  6:139316/148‑1 (MQ=39)
                                                                    tgGGGGCACATTGAGCCCGGTAGGCGCGGGTTTACCCGCGAAGGCGGAGGTATCGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGGGGCGGGGAGGGCTATCAGGCGTCCTGATAGctct                                                                     <  7:55380/148‑1 (MQ=37)
                                                                            catgccgCCCTGACGGCGGGGTTTTACCCGCTAAGGCGGAGGAATCCCCACCATCCCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGAc                                                             <  5:8321/142‑1 (MQ=255)
                                                                            cATTGCGCCCTGTCGGCGCGGGTTTACCCGCGCAGGCGGAGGTATCCCCACCACCGCAGTCGCGGTAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCGCTCGATCGAc                                                             <  8:62825/148‑1 (MQ=255)
                                                                             aTTGAGCCCTGTAGGCGAGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGTGTAAACCCGCTCCTAAACAGAGCCCGCCACCAAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACg                                                            <  8:29988/148‑1 (MQ=255)
                                                                                                    ttccccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                      <  5:122028/145‑1 (MQ=255)
                                                                                                    ttACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                      <  4:157611/148‑1 (MQ=255)
                                                                                                        ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGCTCGACGGGCAGGCGCACACCAGGTTGCGGTCg                                  >  1:94329/1‑148 (MQ=255)
                                                                                                        ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGCCGGGCAGGCGCACACCAGGTTGCGGTCg                                  >  6:78442/1‑148 (MQ=255)
                                                                                                                  ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGt                        >  7:154926/1‑148 (MQ=255)
                                                                                                                      gTATAGCCACCATCGCATTCGCGGTTAAACCCGCTCCTACACAGAGCCGGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACTCCAGGTTGCGGTCGGCGAACACGTTGt                     <  8:131238/147‑1 (MQ=255)
                                                                                                                          cgCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAcc                <  2:115337/147‑1 (MQ=255)
                                                                                                                                caTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCg          >  1:6188/1‑148 (MQ=255)
                                                                                                                                          gcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcgg  >  2:178817/1‑146 (MQ=255)
                                                                                                                                                     |                                                                                                                                         
GAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGG  >  NC_002947/1126500‑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: