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

GGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCA  >  NC_002947/1126504‑1126784
                                                                                                                                                 |                                                                                                                                            
ggCCGATGTCGATCTTTAACTGGTCTAAGAGGGTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGGTAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccgctc                                                                                                                                       <  8:164042/149‑6 (MQ=37)
 gccgatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTAACCCGAATGGGGAGGGGTAGCCACCAACGCCATAGCGGGTTAAccccctc                                                                                                                                       >  1:366969/1‑143 (MQ=255)
 gccgatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGAGGTATAGCCACCATAGCATTAGCGGGTAAACCCGCTcc                                                                                                                                       >  5:104694/1‑149 (MQ=255)
        ccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                 <  8:106198/148‑1 (MQ=255)
           atCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTTGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCGCGCATGCGGGGGTGTAGCCACCAGGGCATTCGCGGGTGAGCCCGGGCCTGCACAGAGc                                                                                                                            >  8:386536/1‑149 (MQ=255)
                gAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACagag                                                                                                                              <  1:254097/143‑1 (MQ=255)
                gAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGccc                                                                                                                          <  2:119446/146‑1 (MQ=255)
                gAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGc                                                                                                                        <  5:126948/148‑1 (MQ=255)
                      gTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGCAAGCGGAGGTATAGCCACCATAGCAATCGCGGGTAAACCCCCGACTACACACAGCCCGCCAcccca                                                                                                                 >  6:318098/1‑147 (MQ=255)
                           aGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGGGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAg                                                                                                            <  1:364386/149‑1 (MQ=255)
                                   aGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCggg                                                                                                    <  7:252079/149‑1 (MQ=255)
                                   aGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCg                                                                                                      <  7:68693/147‑1 (MQ=255)
                                       aCATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGAGGTATAGCCACCATCGCATTAGCGGGTAAACCCGCTCCTACACAGAGCCCGCCCCCCAAGCCAGGGGGCGgggggg                                                                                                >  6:183800/1‑149 (MQ=38)
                                       aCATGTGGACCTCCTGCCAGGTTGTTGGGTGCACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAg                                                                                                 <  3:108074/148‑1 (MQ=255)
                                              gACATCCTGCCAGGTTGTTGGGGGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTAAAGCCACCATCGCTTTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCa                                                                                         <  3:391739/149‑1 (MQ=38)
                                               aCCTCCTGCCAGGTGGTTGGGGGGACATTGAGCCCTGTAGGCGCGGGTGTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAg                                                                                        <  5:318098/149‑1 (MQ=38)
                                                cctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGACCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAg                                                                                        <  1:279254/149‑1 (MQ=255)
                                                  tcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAc                                                                                                                                             <  5:170377/93‑1 (MQ=255)
                                                  tcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAc                                                                                                                                             >  6:170377/1‑93 (MQ=255)
                                                      gCCTGGTTGTTGGGTGCACATTGCCCCCTGTCGGCGGGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATGCGCGGGAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCt                                                                                 <  1:153448/149‑1 (MQ=38)
                                                      gCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAg                                                                                        <  3:331746/142‑1 (MQ=255)
                                                      gCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGAGGTATATCCACCAAAGCATTAGCGGGTAAACCCGCTCCTAAACAGAGCCCCCCCCCCAAGCCAGGTGGCGGGGAGGGCTATCAggagatag                                                                                 >  4:331746/1‑143 (MQ=255)
                                                         agtttgtggGGTTGACATTCAGCACTGTAGGCGCGGGTTTACCCGCTAATGCGGAGGTATAGCCACCATCGCATTCGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGa                                                                               <  3:225216/145‑1 (MQ=255)
                                                         aGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGGTTTACCCCCAAAGCGGAGGGTTTTGCCCCCAACCCATAGACGGGTAAAACCCGACGACACACACCCCCCCCCCCCAACCAGGGGGGGGGGAGGGCGATCAGGCGGCCGGa                                                                               >  4:225216/1‑148 (MQ=255)
                                                          ggTTGTTGGGGGGACATGGCGCCCTGTCGCCGCGGGTTTACCCGCGAATGCGTAGTTATACCCACCATCGCTTTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATa                                                                             <  2:156891/149‑1 (MQ=37)
                                                               tttggTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctct                                                                        <  7:386536/146‑1 (MQ=255)
                                                                tgGGTGGACTTTGAGCCTTGTAGGGGCGGGTTTACCCGCGAAGGCGGAGGTAAAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctct                                                                        <  2:231181/148‑1 (MQ=38)
                                                                 gggTGGCCATTGAGCCCTGTGGGCGCGGGTTTCCCCGCGAATGGGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCCCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCtctc                                                                       <  6:2184/148‑1 (MQ=38)
                                                                   gTGCACATTGCCCCCTGTAGGCGCGGGTTTCCCCGCGAATGCGCAGGTATCCCCACCACCGCATTCGCGGGTAAACCCGCGCCTACACAGAGCCCGCCACCCAAGCCAGGTGGTGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgat                                                                    <  1:61340/149‑1 (MQ=37)
                                                                       aCATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGAc                                                                <  7:66488/149‑1 (MQ=255)
                                                                                       gcgcGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACAcc                                                 <  8:341533/148‑1 (MQ=255)
                                                                                         gcgGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAgg                                              >  5:319854/1‑149 (MQ=255)
                                                                                         gcgGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAgg                                              >  2:280310/1‑149 (MQ=255)
                                                                                                  aCCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGt                                       <  7:267855/148‑1 (MQ=255)
                                                                                                    ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgatc                                                                   <  8:117414/118‑1 (MQ=255)
                                                                                                    ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgatc                                                                   >  4:354743/1‑118 (MQ=255)
                                                                                                    ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgatc                                                                   <  3:354743/118‑1 (MQ=255)
                                                                                                    ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgatc                                                                   >  7:117414/1‑118 (MQ=255)
                                                                                                    ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGAc                                                                <  1:397587/121‑1 (MQ=255)
                                                                                                    ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGAc                                                                >  2:397587/1‑121 (MQ=255)
                                                                                                    ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACAcc                                                 <  5:358715/136‑1 (MQ=255)
                                                                                                    ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACAcc                                                 >  6:358715/1‑136 (MQ=255)
                                                                                                    ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                    <  5:183800/149‑1 (MQ=255)
                                                                                                    ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                    >  2:118982/1‑149 (MQ=255)
                                                                                                                     ataGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTc                                      <  1:162647/131‑1 (MQ=255)
                                                                                                                      taGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAcc                   >  6:411840/1‑149 (MQ=255)
                                                                                                                       aGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAccc                  >  5:165659/1‑149 (MQ=255)
                                                                                                                         ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg                <  1:118982/149‑1 (MQ=255)
                                                                                                                                cGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcg         >  5:260865/1‑149 (MQ=255)
                                                                                                                                 gCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcgg        >  1:20875/1‑149 (MQ=255)
                                                                                                                                   aTTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcg      <  6:319854/149‑1 (MQ=255)
                                                                                                                                       gcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCa  <  7:54861/149‑1 (MQ=255)
                                                                                                                                         gggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAccc                  >  3:95903/1‑131 (MQ=255)
                                                                                                                                         gggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAccc                  <  4:95903/131‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
GGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCA  >  NC_002947/1126504‑1126784

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: