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

GAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCC  >  NC_002947/1126529‑1126772
                                                                                                                      |                                                                                                                                
gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGAGGTATAGGCACCATCGCATTAGCGGGGAAACCCGCTCCTACACAGAGCCCGCCCCCCCAg                                                                                                    >  6:69721/1‑147 (MQ=255)
           gTGACTTGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCggg                                                                                        <  2:116811/148‑1 (MQ=255)
               cATGTGGACCTCCTGCCAGGTTGTGGGGGGTACTTTGAGCCCTGTAGGCGCGGGTTTACCCGCTAATGCTGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAg                                                                                     <  8:103556/147‑1 (MQ=37)
                  gtgGACCTCCTGCCAGGTTGTTGGGGGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCt                                                                                 <  8:82699/148‑1 (MQ=255)
                         tcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCGAATGCGGAGGTATATCCACCATCGCATTAGCGGGTTAACCCCCGCCGACACACACCCCCCCCCCCAACCCCCGTGGCGGGGCGGGCTATCAGGc                                                                          >  3:33913/1‑148 (MQ=255)
                                tggTTGTTGGGTGGACATTGAGCCCTTTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGGGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGa                                                                   <  8:83529/147‑1 (MQ=255)
                                aGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGa                                                                   <  3:69632/148‑1 (MQ=255)
                                      ttgGGGGGACATTGAGCCCTGTCGGCGCGGGTTTCCCCGCGAATGCGTCGGTATAGCCACCATCGCAGTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGct                                                              <  2:113621/147‑1 (MQ=37)
                                          ggggccATTGCGCGCTGTAGGCGCGGGGTTCCCCGCGAAGGCGGAGGTATACCCACCACCGCATTCGCGGTAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTAGCAGGCGTCCTGATAGCTCtcga                                                         <  5:8497/143‑1 (MQ=255)
                                          gTGGACATTGAGCCCTGTAGGCGCGGGTTTTAACCCGAAAGCGGAGGGATAACCACCATACCAATAGCCGGGAAACACCCGCCTAAAAAAAGCCCCCCCCCCAACCCACGGGGCGGGGTGGGCTATAAGGCGGCCAGATaagtctag                                                          >  2:30063/1‑140 (MQ=255)
                                               cttggcgCCCGGTTGGGGCGGTTTTACCCGCGAATGCGCAGCTATAGCCACCAGCCCATTCGCGGGAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGAc                                                    <  2:64160/142‑1 (MQ=255)
                                                    acgccTGTAGGCGCGGGTTTAGCCGCGAATGCGGAGGTATAGCCCCCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACggg                                                 <  7:30017/143‑1 (MQ=38)
                                                       ccTGTAGGCGCGGGTTTACCCCCGAATGCGGAGGGATAGCCACCCTCGCATTAGCGGGTAAACCCCCTCCTACACAGAGCCCCCCCCCCCACCCACGGGGCGGGGAGGGCTATCAGGCGGCCGGCTAGCGGGCGAACGACGGGCCGgc                                            >  6:68641/1‑148 (MQ=255)
                                                                              taaTGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGa                     <  2:113758/147‑1 (MQ=255)
                                                                                    ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGt               >  7:59357/1‑148 (MQ=255)
                                                                                    ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGt               >  3:18908/1‑148 (MQ=255)
                                                                                    ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGt               >  2:39484/1‑148 (MQ=255)
                                                                                     gAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt              <  4:33913/148‑1 (MQ=255)
                                                                                                 ccaCCGCATGCGGGGTAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGcc  <  5:102745/148‑1 (MQ=255)
                                                                                                                      |                                                                                                                                
GAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCC  >  NC_002947/1126529‑1126772

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: