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

CGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTA  >  NC_002947/1126513‑1126787
                                                                                                                                        |                                                                                                                                              
cgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                   >  5:232915/1‑147 (MQ=255)
cgatCTTGAACAGGCCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                   <  7:350828/147‑1 (MQ=255)
cgatCTTGAACAGGCCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGCAAGCGGAGGTATAGCCACCATCGCATTAGCGGGTTAACCCGCTCCTAcaca                                                                                                                                   >  8:350828/1‑147 (MQ=255)
                gAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCGAATTCGGGGGTTTATCCACCATAGCAATAGCGGGTTAACCCGCTCCGAAACAGAGCCCGCCCCCCaaac                                                                                                                  >  7:714/1‑146 (MQ=255)
                      gCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTAAACCCAATAAGGAGGTTTATCCACCCTACCAATAGCCGGGTAAACCGCAACAAAACAGACCCCCCCCCCCAACCCACGTg                                                                                                            >  4:243485/1‑148 (MQ=255)
                      gCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCCCCATCGCCTTATCGGGTAAACCCGCGCCTAAACAGAACCCGCCCCCCAAACCCGGTgg                                                                                                           >  1:107851/1‑149 (MQ=37)
                      gCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGACCCCGCCACCCAAGCCAGGt                                                                                                             <  3:243485/147‑1 (MQ=255)
                            tGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGa                                                                                                     <  6:232915/149‑1 (MQ=255)
                               cATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg                                                                                                  <  1:20118/149‑1 (MQ=255)
                               cATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAggg                                                                                                  <  4:50026/149‑1 (MQ=255)
                                  gtgGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCCCGAATGCGGAGGTATATCCACCATCGCATTAGCGGGTTAACCCCCTCCTACACAGAGCCCCCCCCCCAAACCACGTGGCGGGGAGGGcaa                                                                                               >  3:85285/1‑147 (MQ=255)
                                   tgcccctcctGCCAGGTTGTGGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTAt                                                                                              <  3:279830/145‑1 (MQ=255)
                                     taCCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATc                                                                                             <  7:248530/147‑1 (MQ=255)
                                     gACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTAACGCGAATGCGGAGGTATAGCCCCCATCGCCTTAGCGGGTAAACCCCCGCCGAAACAGAGCCCCCCACCCAAGCCACGGGGCGGGGAGGGCGAt                                                                                              >  6:432259/1‑147 (MQ=255)
                                           ctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTc                                                                                      <  7:210729/149‑1 (MQ=255)
                                              ccAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTcc                                                                                     <  5:432259/147‑1 (MQ=255)
                                               cAGGTTGTGGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGa                                                                                  <  5:152200/149‑1 (MQ=255)
                                                        gggTGGACATTGATCCCTGTAGGGGCGGGTTTACCGGCGAATGCGGAGGTATAGCAACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCtctc                                                                          <  4:85285/148‑1 (MQ=38)
                                                          gTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCGAATGCGGAGGGATAGCCACCATAGCAATAGCCGGGAAAACCGCTCCGACACAGAGCCCCCCACCCCACCCAGGGGGCGGGGGGGGCGATCAGGCGTGCAGATAGcag                                                                            >  6:85292/1‑142 (MQ=255)
                                                                     gCCCTGTAGGCGCGGGTTTACCCGCTAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAgg                                                            <  3:447640/149‑1 (MQ=255)
                                                                                  gggTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGtt                                               >  2:188970/1‑149 (MQ=255)
                                                                                     tttACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGTGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                            >  5:188043/1‑149 (MQ=255)
                                                                                        tcccGCGTATGCGGTTGTAAAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCgg                                           <  7:180727/146‑1 (MQ=39)
                                                                                         cccGCGAATGCGGCGTTAACCCCACCAGCCCATGCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTc                                         <  7:379694/148‑1 (MQ=255)
                                                                                          ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCCCCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGGCCTGATAGCTCTCGAGCCACGGCAAGGCGCCCACCAGGTTCCggac                                         >  8:180727/1‑145 (MQ=38)
                                                                                          ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                       >  6:194818/1‑149 (MQ=255)
                                                                                          ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                       >  6:447195/1‑149 (MQ=255)
                                                                                              gAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGaa                                   <  5:150890/149‑1 (MQ=255)
                                                                                                    ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt                             >  7:401038/1‑149 (MQ=255)
                                                                                                     gAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTg                            >  5:233907/1‑149 (MQ=255)
                                                                                                     gAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTg                            <  2:84666/149‑1 (MQ=255)
                                                                                                      aGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                            >  6:339183/1‑132 (MQ=255)
                                                                                                      aGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                            <  5:339183/132‑1 (MQ=255)
                                                                                                      aGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGt                           <  2:107851/149‑1 (MQ=255)
                                                                                                           taGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGca                                                        <  3:226014/115‑1 (MQ=255)
                                                                                                           taGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGca                                                        >  4:226014/1‑115 (MQ=255)
                                                                                                             gCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCg                    >  6:196322/1‑149 (MQ=255)
                                                                                                              ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg                   <  5:85292/149‑1 (MQ=255)
                                                                                                                accaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACg                               >  4:403035/1‑135 (MQ=255)
                                                                                                                accaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGcc                 >  6:85143/1‑149 (MQ=255)
                                                                                                                 cgaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGTGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCg                <  5:379706/147‑1 (MQ=255)
                                                                                                                  caTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGa               >  7:292644/1‑149 (MQ=255)
                                                                                                                   aTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGc                                                                                                          <  7:399379/57‑1 (MQ=25)
                                                                                                                   aTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGc                                                                                                          >  8:399379/1‑57 (MQ=25)
                                                                                                                   aTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGAc              >  4:32665/1‑149 (MQ=255)
                                                                                                                     cGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGGCCTGATAGCACTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcg            >  7:297959/1‑149 (MQ=255)
                                                                                                                       cATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTg                            <  7:215406/131‑1 (MQ=255)
                                                                                                                       cATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTg                            >  8:215406/1‑131 (MQ=255)
                                                                                                                       cATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggc          <  8:714/149‑1 (MQ=255)
                                                                                                                          tCGCGGGTAAGCCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc       >  3:186728/1‑149 (MQ=255)
                                                                                                                           cgcgTGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc       <  2:301313/148‑1 (MQ=255)
                                                                                                                                ggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTa  <  5:43235/149‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                              
CGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTA  >  NC_002947/1126513‑1126787

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: