Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,370,069 T→A 7.8% E54D (GAA→GAT)  ramA ← (R)‑stereoselective amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,370,0690TA 7.8% 156.0 / 4.2 51E54D (GAA→GAT) ramA(R)‑stereoselective amidase
Reads supporting (aligned to +/- strand):  ref base T (28/19);  new base A (2/2);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

TCCCTTGCGGGGTGACCAGCAAGGAGGTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCGGTATCGCCATCGCGACCGGAGAGTTG  >  NC_002947/4369927‑4370212
                                                                                                                                              |                                                                                                                                               
tCCCTTGCGGGGTGACCAGCAAGGAGGTGTTGTAAACGGTGCCGGAGTCGTTCTTAGCCAGACCAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGa                                                                                                                                          <  2:1063272/150‑1 (MQ=255)
  ccTTGCGGGGTGACCAGCAAGGAGGTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATAACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGc                                                                                                                                        >  2:394979/1‑150 (MQ=255)
   cTTGCGGGGTGACCAGCAAGGAGGTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCt                                                                                                                                       >  1:564162/1‑150 (MQ=255)
      gCGGGGTGACCAGCAAGGAGGTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGTCTTTTTCTGCTAGCTGtt                                                                                                                                    <  2:160578/150‑1 (MQ=255)
       gggggggACCAGCAAGGAGGTGTTGTAAACGGTGCCGGAGTGGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTg                                                                                                                                   <  2:188810/149‑1 (MQ=255)
        ggggTGACCAGCAAGGAGGTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTgg                                                                                                                                  <  1:288747/150‑1 (MQ=255)
                cAGCAAGGAGGTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTggg                                                                                                                                 >  1:13612/1‑143 (MQ=255)
                cAGCAAGGAGGTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTggg                                                                                                                                 <  2:13612/143‑1 (MQ=255)
                 aGCAAGGAGGTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCa                                                                                                                         >  1:554136/1‑150 (MQ=255)
                                  aaCGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGt                                                                                                        <  2:330569/150‑1 (MQ=255)
                                         ccGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTggg                                                                                                 >  1:304216/1‑150 (MQ=255)
                                               tCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATc                                                                                           >  1:700170/1‑150 (MQ=255)
                                                 gTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATcag                                                                                         >  2:748625/1‑150 (MQ=255)
                                                  ttCTCAGCCAGACCAATCACCGCTGCGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTCAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCCGCTGGGTCTCTGGGAAAATcagc                                                                                        <  2:554136/150‑1 (MQ=255)
                                                     tcAGCCAGACCAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGa                                                                                     >  2:736115/1‑150 (MQ=255)
                                                        gCCAGACCAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTc                                                                                                                                               <  1:361363/89‑1 (MQ=255)
                                                        gCCAGACCAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTc                                                                                                                                               >  2:361363/1‑89 (MQ=255)
                                                             aCCAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCCGACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATc                                                                             <  2:579305/150‑1 (MQ=255)
                                                               cAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCa                                                                           >  1:463658/1‑150 (MQ=255)
                                                               cAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCa                                                                           <  1:1022091/150‑1 (MQ=255)
                                                                aaTCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCGGGGTCTCTGGGAAAATCAGCAGAGCGGTATCCaa                                                                          >  2:2473/1‑149 (MQ=255)
                                                                 aTCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGCCGAAATCGCCGAGc                                                                                                                                        >  2:500079/1‑87 (MQ=255)
                                                                 aTCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGCCGAAATCGCCGAGc                                                                                                                                        <  1:500079/87‑1 (MQ=255)
                                                                   cACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGgc                                                                       >  2:501825/1‑150 (MQ=255)
                                                                    aCCGCTGAGACGTTTCGTTCGTGTCAGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGgcg                                                                      >  1:704695/1‑150 (MQ=255)
                                                                      cGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGAGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCa                                                                    >  1:782631/1‑150 (MQ=255)
                                                                       gCTGAGACGTTTCGGTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAg                                                                   <  2:20563/150‑1 (MQ=255)
                                                                          gagaCGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTgg                                                                                                                           <  2:383320/91‑1 (MQ=255)
                                                                          gagaCGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTgg                                                                                                                           >  1:383320/1‑91 (MQ=255)
                                                                                tttCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTa                                                                <  1:686050/144‑1 (MQ=255)
                                                                                tttCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTa                                                                >  2:686050/1‑144 (MQ=255)
                                                                                tttCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATg                                                          <  2:463658/150‑1 (MQ=255)
                                                                                tttCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATg                                                          >  2:925696/1‑150 (MQ=255)
                                                                                    gttcgtGTACGGCGCGTTGCACGGCCTGGACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAGAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCt                                                      >  1:18134/1‑150 (MQ=255)
                                                                                        gtgtACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGAAACAGCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCaa                                                  >  1:1039919/1‑150 (MQ=255)
                                                                                                   gTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCa                                       <  2:564162/150‑1 (MQ=255)
                                                                                                   gTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCa                                       >  1:902604/1‑150 (MQ=255)
                                                                                                         cGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTtat                                 >  1:841289/1‑150 (MQ=255)
                                                                                                            ccTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAg                              >  1:295556/1‑150 (MQ=255)
                                                                                                              tGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCg                            <  1:501825/150‑1 (MQ=255)
                                                                                                                   acTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCGGTATc                       <  2:782631/150‑1 (MQ=255)
                                                                                                                        gCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCGGTATCGCCAt                  >  2:166302/1‑150 (MQ=255)
                                                                                                                          ccGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCGGTATCGCCATcg                >  1:387734/1‑150 (MQ=255)
                                                                                                                            gTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCGGTATCGCCATcgcg              >  1:553391/1‑150 (MQ=255)
                                                                                                                                    gCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCGGTATCGCCATCGCGACCGgaga      >  2:566896/1‑150 (MQ=255)
                                                                                                                                      tCGGCGATATCGGCGAGCTGTTGGGCGCTGGCAAATCCCATc                                                                                                                <  1:703297/42‑1 (MQ=255)
                                                                                                                                      tCGGCGATATCGGCGAGCTGTTGGGCGCTGGCAAATCCCATc                                                                                                                >  2:703297/1‑42 (MQ=255)
                                                                                                                                       cGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCACTGTACGTTCCAGGTTATAAGCgg                           <  2:556027/126‑1 (MQ=255)
                                                                                                                                       cGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCgg                           >  1:556027/1‑126 (MQ=255)
                                                                                                                                        ggCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCGGTATCGCCATCGCGACCGGAGAGTTg  >  2:224888/1‑150 (MQ=255)
                                                                                                                                        ggCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCGGTATCGCCATCGCGACCGGAGAGTTg  <  1:166302/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
TCCCTTGCGGGGTGACCAGCAAGGAGGTGTTGTAAACGGTGCCGGAGTCGTTCTCAGCCAGACCAATCACCGCTGAGACGTTTCGTTCGTGTACGGCGCGTTGCACGGCCTGCACACTCGGCCCGTTAACGGGCTCGGCGATTTCGGCGAGCTGTTGGGCGCTGGCAAATCCCATCAGCTGGGTCTCTGGGAAAATCAGCAGATCGGTATCCACGGCGCAGGTAGCGATGGCCTGCAATGTACGTTCCAGGTTATAAGCGGTATCGCCATCGCGACCGGAGAGTTG  >  NC_002947/4369927‑4370212

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