Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,003,645 C→A 9.4% G98C (GGC→TGC) ‡ amaB ← L‑piperidine‑6‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,003,6450CA 9.4% 157.4 / 6.9 53G98C (GGC→TGC) ‡amaBL‑piperidine‑6‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (24/24);  new base A (3/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.63e-01

ACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGG  >  NC_002947/6003503‑6003787
                                                                                                                                              |                                                                                                                                               
aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCt                                                                                                                                            <  1:197650/148‑1 (MQ=255)
aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTg                                                                                                                                           >  5:73228/1‑149 (MQ=255)
 cGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTg‑g                                                                                                                                         <  2:175149/149‑1 (MQ=255)
  gCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTg‑gg                                                                                                                                        >  7:104851/1‑149 (MQ=255)
  gCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTg‑g                                                                                                                                         <  4:148585/148‑1 (MQ=255)
  gCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTg‑g                                                                                                                                         >  3:148585/1‑148 (MQ=255)
   ccGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGt                                                                                                                                       <  7:195815/149‑1 (MQ=255)
      aCCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGAt                                                                                                                                    <  2:224439/149‑1 (MQ=255)
          cacCCAGCGGGTGCCAGGTTTCACGCATGTGGTTGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCACGCCTTCCTG‑GGTGATCTTg                                                                                                                                <  4:232979/149‑1 (MQ=255)
            cccAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGcc                                                                                                                              <  6:112508/149‑1 (MQ=255)
             ccAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGt                                                                                                                                       <  3:35358/139‑1 (MQ=255)
             ccAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGt                                                                                                                                       >  4:35358/1‑139 (MQ=255)
             ccAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGcc                                                                                                                              >  6:124576/1‑148 (MQ=255)
             ccAGCGGGTGCCAGGTTTCACGCATGAGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTCGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCg                                                                                                                             >  5:24482/1‑149 (MQ=255)
                       ccAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTtcgatc                                                                                                                    <  6:24482/148‑1 (MQ=255)
                        cAGGTTTCACGCATGTGGTGGGCCGGGCGCTCGGAAGCTATGGTCAGGCCTTACAGCTGGCGCGACAGGCCAACTGCGATGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGTCCTTCCTG‑GGTGATCTTGCCGGCTtcgatcga                                                                                                                  <  8:176678/149‑1 (MQ=255)
                            tttCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGCGGTGATCTTGCCGGCTTCGATCGAGAc                                                                                                               <  1:12211/149‑1 (MQ=255)
                              tCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGGCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCa                                                                                                             <  6:266719/148‑1 (MQ=255)
                                 cGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTc                                                                                                         <  2:249510/149‑1 (MQ=255)
                                            gCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCgg                                                                                              >  1:157633/1‑149 (MQ=255)
                                             cccGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGc                                                                                             <  8:104851/149‑1 (MQ=255)
                                                gggCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGGCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCtt                                                                                          <  7:174044/149‑1 (MQ=255)
                                                    gcTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGAAGGCCTG‑GGTGaa                                                                                                                                    <  3:3197/103‑2 (MQ=255)
                                                    gcTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGAAGGCCTG‑GGTGaa                                                                                                                                    >  4:3197/1‑102 (MQ=255)
                                                               aTGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGa                                                                                                    >  3:74091/1‑124 (MQ=255)
                                                               aTGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGa                                                                                                    <  4:74091/124‑1 (MQ=255)
                                                                  gTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTc                                                                        >  7:54548/1‑149 (MQ=255)
                                                                     aGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGc                                                                      <  8:54548/148‑1 (MQ=255)
                                                                       gCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTc                                                                                    <  8:89158/132‑1 (MQ=255)
                                                                       gCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTc                                                                                    >  7:89158/1‑132 (MQ=255)
                                                                               gCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGgcg                                                            >  1:69082/1‑148 (MQ=255)
                                                                                 tGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACgcagc                                                                             <  6:139185/129‑1 (MQ=255)
                                                                                 tGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACgcagc                                                                             >  5:139185/1‑129 (MQ=255)
                                                                                   gcgcGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCg                                                                                                     >  7:241908/1‑103 (MQ=255)
                                                                                     gcgACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACgcag                                                                              <  4:189123/124‑1 (MQ=255)
                                                                                     gcgACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACgcag                                                                              >  3:189123/1‑124 (MQ=255)
                                                                                        aCAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCg                                                  <  2:69082/149‑1 (MQ=255)
                                                                                              cAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACgg                                            >  6:185099/1‑149 (MQ=255)
                                                                                                        gTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCg                                  <  1:263517/149‑1 (MQ=255)
                                                                                                              gATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCacac                            <  5:185099/149‑1 (MQ=255)
                                                                                                                    gATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTgcgc                       >  1:74738/1‑148 (MQ=255)
                                                                                                                    gATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTgcgc                       >  7:24570/1‑148 (MQ=255)
                                                                                                                     aTCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTgcgc                       >  3:29099/1‑147 (MQ=255)
                                                                                                                         tttCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCt                 >  4:51805/1‑149 (MQ=255)
                                                                                                                          ttCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCa                     >  4:244274/1‑144 (MQ=255)
                                                                                                                          ttCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCa                     <  3:244274/144‑1 (MQ=255)
                                                                                                                              tGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCg              >  2:28135/1‑147 (MQ=255)
                                                                                                                                  gatCACCCAGGAAGGCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGgcg        >  7:239062/3‑149 (MQ=255)
                                                                                                                                   atCACCCAGGAAGGCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTtgt                                                                                        <  3:112198/67‑1 (MQ=37)
                                                                                                                                   atCACCCAGGAAGGCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTtgt                                                                                        >  4:112198/2‑68 (MQ=37)
                                                                                                                                    tCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCg                                                  >  2:189645/1‑105 (MQ=255)
                                                                                                                                       cccAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCAACAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTg   >  8:225186/1‑149 (MQ=255)
                                                                                                                                        ccAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTgg  >  5:151146/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                               
ACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTG‑GGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGG  >  NC_002947/6003503‑6003787

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑