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.C86.0% 146.4 / 18.3 50intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (6/1);  new base C (20/23);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAG  >  NC_002947/1126500‑1126785
                                                                                                                                                        |                                                                                                                                             
gAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAc                                                                                                                                              <  7:273690/147‑1 (MQ=255)
gAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGGCTATACCTCCGCATTCGCGGGTAAAc                                                                                                                                              <  3:85868/147‑1 (MQ=25)
        gatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGCATGGGGAGGGTTAGGCCCCATCGCATTCGCGGGGTGACCCGCTccgtc                                                                                                                                    >  1:50531/1‑146 (MQ=255)
             cgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                 >  7:304956/1‑147 (MQ=255)
                        aGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGATGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAGTGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCAccc                                                                                                                    <  7:117256/149‑1 (MQ=255)
                          gTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCCCGAAAGCGGAGGTTTAGCCCCCAAAGCATTAGCGGGTTAACCCGCGCCTACACAGAGCCCGccccccca                                                                                                                  >  7:116660/1‑147 (MQ=255)
                             gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCCCGAAATCGGAGGTATATCCAACATCGCATTAGCGGGTAAACCCCCAACTACACACAGCCCCCCACCCaaac                                                                                                                >  1:270395/1‑146 (MQ=255)
                                  cGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCCCCCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAgg                                                                                                            <  1:307343/148‑1 (MQ=38)
                                   gCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGCATGCGGAGGTTTTGCCCCCCTAGCATTAGCGGGTAAAACCGCTCCTACACACAGCCCGCCCCCCCAGCCACGTg                                                                                                          >  5:144920/1‑148 (MQ=255)
                                        gTGACATGTGGACCTCCTGCCTGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCg                                                                                                       <  5:1931/146‑1 (MQ=255)
                                        gTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCgggg                                                                                                    <  2:44341/149‑1 (MQ=255)
                                           aCATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAgg                                                                                                 <  7:270510/149‑1 (MQ=255)
                                            cATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg                                                                                                <  8:304956/149‑1 (MQ=255)
                                             aTGTGCACCTCCTGCCCGGTTGTGGGGTGGACATTGAGCCCTGTCGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg                                                                                                <  5:9976/148‑1 (MQ=38)
                                                     ctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCCCATGAGCAAGTATAGCCAAAAAACCCATAGCCGGTAAACACCCACATACACACAGCCCCCCACCCAACCCAGCTGGCGGGGAGGGCTAtaaag                                                                                        >  6:280974/1‑144 (MQ=255)
                                                        ctGCCCGGTGGTGGGGTGCACTTTGAGCCCTGTCGGCGCGGGTTTACCCGCTAAGGCGGCGGTATACCCACCATCGCATTCGCGGTAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTc                                                                                    <  2:378483/149‑1 (MQ=255)
                                                                  gttgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTATCCCGCGAATGCGGAGGTATAGCCACCATCGCATTAGCGGGTAAACCCGCTGCTACACAGAGCCCGCCCCCCAAGCCAGGTGGCGGGGGGGGCTATCAGGCGTCCTGATAGct                                                                           >  1:64689/1‑149 (MQ=255)
                                                                     gggTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCGAAAGCCGAGGGATAGCCACCAAAGCAATAGCGGGTAAAACCCCGACTACACACAGCCCCCCCCCCCAACCACGTGGCGGGGAGGGCTATAAAGCGGCCGGATAGctct                                                                         >  6:92775/1‑147 (MQ=255)
                                                                           aCATTGCGCCCGGTCGGCGCGGGATTACCCGCGAATGCGGAGGTATCGCCACCATCGAAGTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGGGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGAc                                                                 <  7:321885/149‑1 (MQ=38)
                                                                           aCATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGAc                                                                 <  2:270395/149‑1 (MQ=255)
                                                                              ttGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACggg                                                              <  2:64689/149‑1 (MQ=255)
                                                                                   cccTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgc                                                         <  2:50531/149‑1 (MQ=255)
                                                                                    ccTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgc                                                         <  5:232555/148‑1 (MQ=255)
                                                                                     gtGTAGGCGGGGTTTTACCCGCTAAGGGGTAGTTATACCAACCACCGCATGCGGGGTAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgcgc                                                       <  1:195894/148‑1 (MQ=255)
                                                                                                         cgcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc                                    <  7:344877/149‑1 (MQ=255)
                                                                                                          gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc                                    >  1:111861/1‑148 (MQ=255)
                                                                                                          gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGACCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc                                    >  4:83137/1‑148 (MQ=255)
                                                                                                             aaTGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGa                                  >  5:28640/1‑147 (MQ=255)
                                                                                                               tGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAaca                               <  6:144920/148‑1 (MQ=255)
                                                                                                                gCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAacac                              >  6:138389/1‑148 (MQ=255)
                                                                                                                  ggAGGTATAGCCACGATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACg                             <  5:280974/147‑1 (MQ=255)
                                                                                                                  ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCACGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt                           >  1:385066/1‑149 (MQ=255)
                                                                                                                   gAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGGTGCGGTCGCCGAACACGTTg                          >  4:99481/1‑149 (MQ=255)
                                                                                                                         taGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAc                     >  6:291736/1‑148 (MQ=255)
                                                                                                                             ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCACGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg                 >  1:379862/1‑149 (MQ=255)
                                                                                                                             ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg                 >  4:359321/1‑149 (MQ=255)
                                                                                                                             ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCg                  >  6:233400/1‑148 (MQ=255)
                                                                                                                              caccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGc                <  7:32032/149‑1 (MQ=255)
                                                                                                                              cacattCGCATTCGGGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGc                <  3:7050/144‑1 (MQ=255)
                                                                                                                               accaACGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAACCCAGGTGGCGGGGAGGGCTATCAGGCCTCCTGATAGCACACGAACGACGGGCAGGCGCACACAAGGTTGCGGTCGCCGAACACCTTCTCGACCCGGcc               >  8:171094/1‑149 (MQ=255)
                                                                                                                                  aTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGAc            >  1:164085/1‑149 (MQ=255)
                                                                                                                                   tCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGAcc           >  4:324336/1‑149 (MQ=255)
                                                                                                                                    cGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcg          >  7:272078/1‑149 (MQ=255)
                                                                                                                                    cGCATTCGCGGGTAAACACGCTCCTACACAGACCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGCTAGCGGTCGCCGAACCCGTTGTCGACCCGGCCGAcc           >  6:13734/1‑148 (MQ=255)
                                                                                                                                       aTTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggc        <  5:291736/148‑1 (MQ=255)
                                                                                                                                        ttCGCGGGTATACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcgg      <  2:164085/149‑1 (MQ=255)
                                                                                                                                            cgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCg                  <  2:140629/133‑1 (MQ=255)
                                                                                                                                            cgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCg                  >  1:140629/1‑133 (MQ=255)
                                                                                                                                            cgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAg  <  8:116660/149‑1 (MQ=255)
                                                                                                                                            cgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAg  >  1:294725/1‑149 (MQ=255)
                                                                                                                                            cgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAg  <  2:294725/149‑1 (MQ=255)
                                                                                                                                              gggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAg  >  2:360827/1‑148 (MQ=255)
                                                                                                                                                        |                                                                                                                                             
GAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAG  >  NC_002947/1126500‑1126785

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: