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

CCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAC  >  NC_002947/1126506‑1126788
                                                                                                                                                |                                                                                                                                                 
ccgatgccgatCTTGACCCGGGCGAATACGCGCAGTGACATGTGGACCTCCTGCCTGGTTGTTGGGTGTACTTTGAGCCCTGTTGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATGCGCGGGTAAACCCGCTCCt                                                                                                                                          <  8:312354/149‑1 (MQ=37)
             tGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCCCCAAAGCGGAGGTATAGCACCCAACGCAATCGCGGGTAAAACCGCGACGACACAGAGCCCGc                                                                                                                            >  2:82334/1‑149 (MQ=255)
             tGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCCCCAATGCGGAGGTATACGCACCATCGCATTAGCGGGTAAAAACGCTACAACACAGAGGCCGc                                                                                                                            >  2:87563/1‑149 (MQ=255)
               aaCAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCa                                                                                                                          <  1:249124/149‑1 (MQ=255)
                      cGAACACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGGGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGc                                                                                                                   <  7:132750/149‑1 (MQ=255)
                                 aGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCggg                                                                                                        <  1:275695/149‑1 (MQ=255)
                                  gTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCgggg                                                                                                       >  5:363252/1‑149 (MQ=255)
                                      cATGTGGCCCCCAGGCCCGGTTGTTGGGGGTCCTTTGCGCCCTGTCGGCGCGGGTTTCCCGGCTAATGCGGAGGTAACCCCACCATCGCATCCGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg                                                                                                   <  8:190771/149‑1 (MQ=255)
                                      cATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg                                                                                                   <  5:314992/149‑1 (MQ=255)
                                       aTGTGGACCCCCTGCCCGGTTGTTGGGTGGACATTGCGCCCGGTAGGCGCGGGTTTACCCGCTAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGc                                                                                                  <  8:218798/149‑1 (MQ=38)
                                         gtgGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCCCGAATGCGGAGGTATAGCCCACATAGCATTAGCGGGGTAAACCGCTGCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTa                                                                                                >  5:286697/1‑149 (MQ=255)
                                               ctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTCGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCGAGGTGGCGGGGAGGGCTATCAGGc                                                                                          <  1:87563/149‑1 (MQ=255)
                                               ctcctGCCAGGTTGTTGGGGGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGc                                                                                          <  1:82334/149‑1 (MQ=255)
                                                 cctGCCAGGTGGTGGGGTGGACCTTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATACCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGt                                                                                        <  5:210975/149‑1 (MQ=38)
                                                    tgcctgtttttGGGGGGAACATTCAGCCCTTTCGGCGCGGTTTTACCCGCCAATGCGTAGTTATACCCACAACCGAATTCGGGTTAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTcc                                                                                      <  2:166139/144‑1 (MQ=255)
                                                       aGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGGGGATGTATTTCCACCATCGCATTAGCGGGTTAAACCGCACCTACACAGAGACCGCCACCCAACCCAGGTGGCGGGGAGGGCTATAACGCGtacaga                                                                                   >  5:108652/1‑143 (MQ=37)
                                                           tgttgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATTCGGAGGTATAGCCACCATCGCATTAGCGGGTTAACACGCTCCTAAACAGAGCCCGCCACCCAAGCCAGGGGGCGGGGAGGGCTATCAGGCGTCCTGATAg                                                                                >  7:252097/1‑147 (MQ=255)
                                                                             gCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGAATAGCCCCCATCGCATTCGCGGGTAAACCCGCTCCTACAGAGAGCCCGCCACCCAAGCCTGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAg                                                              <  6:89512/148‑1 (MQ=38)
                                                                              cccTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgc                                                            <  6:363252/149‑1 (MQ=255)
                                                                                          gggTTTACCCGCGAATGGGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGtt                                                <  5:1221/149‑1 (MQ=255)
                                                                                          gggTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCa                                                    <  7:22118/145‑1 (MQ=255)
                                                                                          gggTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCa                                                    >  8:22118/1‑145 (MQ=255)
                                                                                                 cccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTc                                          <  3:26397/148‑1 (MQ=255)
                                                                                                 cccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCg                                         >  2:206534/1‑149 (MQ=255)
                                                                                                   cgcgAATGCGGAGGTATAGCCGCCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc                                       >  2:373042/1‑149 (MQ=255)
                                                                                                   cgcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc                                       <  6:286697/149‑1 (MQ=255)
                                                                                                                gTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTc                           >  3:316434/1‑148 (MQ=255)
                                                                                                                  ataGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGa                                                                                   >  1:183352/1‑90 (MQ=255)
                                                                                                                  ataGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGa                                                                                   <  2:183352/90‑1 (MQ=255)
                                                                                                                  ataGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGAt                                                                                  <  5:363924/91‑1 (MQ=255)
                                                                                                                  ataGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGAt                                                                                  >  6:363924/1‑91 (MQ=255)
                                                                                                                     gccgccATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAccc                      <  1:373042/148‑1 (MQ=255)
                                                                                                                     gCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCg                     <  6:108652/149‑1 (MQ=255)
                                                                                                                      ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg                    >  1:340610/1‑149 (MQ=255)
                                                                                                                      ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg                    >  5:326079/1‑149 (MQ=255)
                                                                                                                       caccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGc                   >  3:21029/1‑149 (MQ=255)
                                                                                                                               cATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAgg                                                             <  5:289172/99‑1 (MQ=255)
                                                                                                                               cATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAgg                                                             >  6:289172/1‑99 (MQ=255)
                                                                                                                               cATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggc           <  3:64311/149‑1 (MQ=255)
                                                                                                                                 ttCGCGGGTAAACCCGCTCCTACACAGAGCCTCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcg          >  1:39601/1‑149 (MQ=255)
                                                                                                                                 ttCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATTAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcgg         >  4:130664/1‑149 (MQ=255)
                                                                                                                                  tCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt                              <  6:326190/127‑1 (MQ=255)
                                                                                                                                  tCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt                              >  5:326190/1‑127 (MQ=255)
                                                                                                                                   cgcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATTAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc        <  3:130664/148‑1 (MQ=255)
                                                                                                                                   cgcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc        <  4:21029/148‑1 (MQ=255)
                                                                                                                                    gcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                             >  1:316254/1‑110 (MQ=255)
                                                                                                                                    gcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                             <  2:316254/110‑1 (MQ=255)
                                                                                                                                    gcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTc                           >  3:206709/1‑128 (MQ=255)
                                                                                                                                    gcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTc                           <  4:206709/128‑1 (MQ=255)
                                                                                                                                         gTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAc  <  4:303505/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAC  >  NC_002947/1126506‑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: