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.C92.9% 84.3 / 0.5 28intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (13/13);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAC  >  NC_002947/1126502‑1126788
                                                                                                                                                  |                                                                                                                                                
ggggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccgc                                                                                                                                             <  7:4197/149‑4 (MQ=255)
                      aGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTaaa                                                                                                                                                   >  2:84130/1‑122 (MQ=255)
                      aGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCCCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTaaa                                                                                                                                                   <  1:84130/122‑1 (MQ=255)
                         tCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGCTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGACGTGATGGCACGCACCCAATAGCCGGTAAAACCGCTCCAACACAAAGCCCCCCCCCCaa                                                                                                                      >  8:195124/1‑148 (MQ=255)
                          cGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGc                                                                                                                    <  5:270772/149‑1 (MQ=255)
                           gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGAGGTATAGCCACCAGAGCAATAGCGGGTAAACCCGCGCCAACACAGAGCCCGcccccccaccc                                                                                                                   >  6:270772/1‑144 (MQ=255)
                                               ggACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCGAATGCGGAGGTATAGCCACCCTCGCATTAGCGGGTAAACCCGCGACTACACAGAGCCCCCCCCCCAACCCAGGGGGCGGGGGGGGCTAt                                                                                               >  6:308547/1‑148 (MQ=255)
                                                gACCTCCTGCCAGGATGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTAACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATc                                                                                              <  5:92212/148‑1 (MQ=255)
                                                  cctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACagag                                                                                                                                  <  7:84434/111‑1 (MQ=255)
                                                         ccAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCGCGAAAGCGGAGGTATAGCCACCAACGCAATAGCGGGGAAAACCCCGCCGAAACAGAGCCCCCCACCCAAGCCAGGTGGCGGGGCGGGCTATCAAGCGTCCt                                                                                     >  2:92585/1‑148 (MQ=255)
                                                         ccAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGTGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGTCGGGGAg                                                                                                     <  1:8997/132‑1 (MQ=255)
                                                                 ttgGGGGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctct                                                                            <  4:214911/149‑1 (MQ=255)
                                                                           aTTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACgg                                                                  <  4:85829/149‑1 (MQ=255)
                                                                                        ggCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGCGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCacac                                                      <  1:92585/149‑1 (MQ=255)
                                                                                              gttttACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTg                                               <  7:236713/147‑1 (MQ=255)
                                                                                              ggTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTg                                               >  2:300509/1‑149 (MQ=255)
                                                                                                    cccGCGAATGCGGAGGTATAGCCACCAGCGCTTTCGGGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGACGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCg                                         <  1:142335/149‑1 (MQ=255)
                                                                                                      cgcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGCTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc                                       >  8:23387/1‑149 (MQ=255)
                                                                                                               ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt                              >  3:233633/1‑149 (MQ=255)
                                                                                                                   gTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg                                                                                                   >  7:29157/1‑76 (MQ=255)
                                                                                                                   gTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg                                                                                                   <  8:29157/76‑1 (MQ=255)
                                                                                                                     ataGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAc                        >  8:171255/1‑149 (MQ=255)
                                                                                                                         ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg                    >  7:155828/1‑149 (MQ=255)
                                                                                                                         ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg                    >  2:274323/1‑149 (MQ=255)
                                                                                                                          caccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGc                   >  7:46289/1‑149 (MQ=255)
                                                                                                                           accaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGcc                  <  1:274323/149‑1 (MQ=255)
                                                                                                                             caTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGa                <  1:300509/149‑1 (MQ=255)
                                                                                                                              aTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATGGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGAc               >  3:110067/1‑149 (MQ=255)
                                                                                                                              aTCCCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGAc               <  5:36442/149‑1 (MQ=255)
                                                                                                                                        cgGGTAAACCCGATCATACACAGAGCCCGCCACCCAACCCAGGTGGCGGCGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGCCCGGCGGCCAg     >  2:312922/1‑149 (MQ=255)
                                                                                                                                            gTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAc  >  7:243525/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAC  >  NC_002947/1126502‑1126788

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: