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.C97.8% 146.8 / ‑3.2 45intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (20/24);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

GGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGG  >  NC_002947/1126502‑1126778
                                                                                                                                                 |                                                                                                                                      
ggggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccgc                                                                                                                                   <  5:40779/149‑4 (MQ=255)
          ccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTaca                                                                                                                            >  6:255780/1‑146 (MQ=255)
           cgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCGAATGCGGAGGTATAGCCACCATAGCCATAGCCGGTAAACCCGCTCCAAcaca                                                                                                                          >  6:96408/1‑147 (MQ=255)
                  gAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGcc                                                                                                                     <  7:63049/145‑1 (MQ=255)
                        ggcGAAGACGCGCAGTGAAATGTGGACCTCGTGCCTGGTTGTGGGGGTGACTTTGAGCCCTGTCGGCGCGGGTTTACGCGCGAATGCGGAGTTATAGCCACAATCGAATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCaa                                                                                                           <  6:28882/147‑1 (MQ=25)
                         tCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCGCGAATGCGGAGGTATAGCCCCCATCGCCTTAGCGGGGAAACCCGCGCCGACCCAGCGCCCGCCCCCCCAg                                                                                                          >  6:253892/1‑149 (MQ=255)
                           gAAGACGCTCATGGACATGTGGACCTCCTGCCAGGTTGTGGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGcc                                                                                                        <  8:41483/149‑1 (MQ=38)
                            aaGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGGGGAGGTATTGCCACCATCGCATTAGCGGGTTAACCCGCTCCTAAACAGAGCCCCCCCCCCAAGc                                                                                                         >  5:146526/1‑147 (MQ=255)
                                       tGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGa                                                                                            <  5:255780/149‑1 (MQ=255)
                                       tGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGa                                                                                            <  3:89347/149‑1 (MQ=255)
                                                 aCCTCCTGCCAGGTTGTTGGGGGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAg                                                                                  <  3:273983/149‑1 (MQ=255)
                                                   ctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGc                                                                                >  5:238985/1‑149 (MQ=255)
                                                         ccAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTcc                                                                            <  2:83047/147‑1 (MQ=255)
                                                           aGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGATTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCt                                                                                                                                <  1:26428/94‑1 (MQ=255)
                                                             gttgtgGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGACGTATAGCCACCATCCCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAg                                                                      <  1:283828/149‑1 (MQ=38)
                                                                 ttgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctct                                                                  >  1:57275/1‑149 (MQ=255)
                                                                         aCATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGAc                                                          <  5:253892/149‑1 (MQ=255)
                                                                         aCATTGAGCCCTGTAGGCGCGGGTTAACCCGCGAATGCGCAGGAAAAGCCACCATCGCATGCGCGGGAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGAc                                                          <  5:204046/149‑1 (MQ=255)
                                                                              gcgCCCTGTCGGGGCGGGTTTACCCGCGAATGCGGCGGTATAGCCACCATCCCATTCGCTGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGc                                                      <  4:228900/146‑1 (MQ=37)
                                                                                   cTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgcgc                                                <  6:238985/149‑1 (MQ=255)
                                                                                          cgcgGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACCGGCAGGCGCACACCAg                                         >  4:40096/1‑149 (MQ=255)
                                                                                                ttaaCCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                   <  5:174000/146‑1 (MQ=255)
                                                                                                     ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                              <  3:40096/149‑1 (MQ=255)
                                                                                                     ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                              <  5:96408/149‑1 (MQ=255)
                                                                                                      cgcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTATCACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                              >  1:245340/1‑149 (MQ=255)
                                                                                                      cgcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTATCACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                              >  1:249519/1‑149 (MQ=255)
                                                                                                       gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCg                            <  2:57275/149‑1 (MQ=255)
                                                                                                       gcgAATGCGGAGGTATAGCCACCATCGAATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc                             >  5:79684/1‑148 (MQ=255)
                                                                                                         gAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGaa                          >  3:227951/1‑149 (MQ=255)
                                                                                                         gAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACACAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGaa                          >  8:253121/1‑149 (MQ=255)
                                                                                                                  ggTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTc                 <  6:146526/149‑1 (MQ=255)
                                                                                                                   gTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCg                            <  4:192432/137‑1 (MQ=255)
                                                                                                                   gTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCg                            >  3:192432/1‑137 (MQ=255)
                                                                                                                      taGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                              >  4:3960/1‑132 (MQ=255)
                                                                                                                      taGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                              <  3:3960/132‑1 (MQ=255)
                                                                                                                      taGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAcc             >  7:189442/1‑149 (MQ=255)
                                                                                                                      taGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAcc             >  7:194275/1‑149 (MQ=255)
                                                                                                                      taGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAcc             >  1:197504/1‑149 (MQ=255)
                                                                                                                         ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg          >  8:66279/1‑149 (MQ=255)
                                                                                                                             caTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGa      >  3:306298/1‑149 (MQ=255)
                                                                                                                              aTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTc                 <  6:281246/137‑1 (MQ=255)
                                                                                                                              aTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTc                 >  5:281246/1‑137 (MQ=255)
                                                                                                                               tCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGa               >  6:50899/1‑138 (MQ=255)
                                                                                                                               tCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGa               <  5:50899/138‑1 (MQ=255)
                                                                                                                                 gCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcgg  <  8:189442/149‑1 (MQ=255)
                                                                                                                                 gCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcgg  <  8:194275/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                      
GGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGG  >  NC_002947/1126502‑1126778

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: