Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,532,229 T→A 9.5% K15M (AAG→ATG) ‡ braE ← high‑affinity branched‑chain amino acid ABC transporter permease BraE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,532,2290TA 9.5% 146.0 / 4.5 52K15M (AAG→ATG) ‡braEhigh‑affinity branched‑chain amino acid ABC transporter permease BraE
Reads supporting (aligned to +/- strand):  ref base T (21/26);  new base A (2/3);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGAAGATAAGGATCATCACCAGCAGCGAGAAGCTGAACACATCCTTGTAGTCGGAGTTGATCA  >  NC_002947/5532083‑5532375
                                                                                                                                                  |                                                                                                                                                   
cacCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTg                                                                                                                                                   >  5:355997/1‑149 (MQ=255)
   cgccgACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑tt                                                                                                                                               >  2:12012/1‑149 (MQ=255)
     ccgACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTc                                                                                                                                             >  8:289244/1‑149 (MQ=255)
       gACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGc                                                                                                                                           >  2:291450/1‑149 (MQ=255)
       gACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCAAGGG‑TTTCga                                                                                                                                           <  6:70798/149‑2 (MQ=255)
       gACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCAAGGG‑TTTCga                                                                                                                                           <  2:346993/149‑2 (MQ=255)
       gACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCAAGGG‑TTTCga                                                                                                                                           <  2:221056/149‑2 (MQ=255)
            gccagGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATa                                                                                                                                      <  7:89726/149‑1 (MQ=255)
             ccagGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTg                                                                                                                                                   <  8:97080/136‑1 (MQ=255)
             ccagGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTg                                                                                                                                                   >  7:97080/1‑136 (MQ=255)
              cagGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACTGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATaga                                                                                                                                    <  1:286590/149‑1 (MQ=255)
                     cgcggcgcggCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGtt                                                                                                                             <  7:289244/149‑1 (MQ=255)
                        ggcgcggcTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTg                                                                                                                          <  7:50864/149‑1 (MQ=255)
                               cTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACAACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGGTTTTCGCTAATAGAGGCAGTTTTGGCAACg                                                                                                                    <  1:3983/149‑1 (MQ=255)
                                  ggCATTGAAGGTGTAGCCGTCGAGCACCACACAAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGa                                                                                                                  <  5:218058/147‑1 (MQ=255)
                                     aTTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGt                                                                                                             >  2:297674/1‑149 (MQ=255)
                                      ttGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGtt                                                                                                            <  2:123001/149‑1 (MQ=255)
                                         aaGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTca                                                                                                          >  8:275880/1‑148 (MQ=255)
                                                           accacaccaACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCggg                                                                                       <  1:12012/149‑1 (MQ=255)
                                                            ccacaccaACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGc                                                                                      <  6:355997/149‑1 (MQ=255)
                                                              acaccaACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGa                                                                                    >  6:309816/1‑149 (MQ=255)
                                                                  caACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGAcc                                                                                  >  2:364172/1‑147 (MQ=255)
                                                                     cGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCcagca                                                                             <  7:275880/149‑1 (MQ=255)
                                                                          ggACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCt                                                                        >  7:234787/1‑149 (MQ=255)
                                                                           gACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCtt                                                                       >  7:287644/1‑149 (MQ=255)
                                                                               agacacGATAAGTGCCAGCAGGCCAGCGACACTGGTCTCCAGCAGGCGGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCggg                                                                                       <  6:10687/127‑1 (MQ=255)
                                                                                aaCACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGa                                                                                                                  <  2:350128/101‑1 (MQ=255)
                                                                                aaCACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGa                                                                                                                  >  1:350128/1‑101 (MQ=255)
                                                                                   acGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCgg                                                               <  6:41184/149‑1 (MQ=255)
                                                                                     gATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCgg                                                               >  1:239506/1‑147 (MQ=255)
                                                                                     gATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGaa                                                             <  1:297674/149‑1 (MQ=255)
                                                                                        aaGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGAAGa                                                           >  2:323915/1‑148 (MQ=255)
                                                                                          gTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTgg                                                                   <  4:128272/138‑1 (MQ=255)
                                                                                          gTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTgg                                                                   >  3:128272/1‑138 (MQ=255)
                                                                                            gCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGa                                                                                    >  1:319499/1‑119 (MQ=255)
                                                                                            gCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGa                                                                                    <  2:319499/119‑1 (MQ=255)
                                                                                               agcagGCCAGCGACAATGGTCTCCAGCATGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGAAGATAAGGatc                                                   <  5:274475/149‑1 (MQ=255)
                                                                                               agcagGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGAAGATAAGGatc                                                   <  8:287644/149‑1 (MQ=255)
                                                                                                    gCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGAAGATAACGATCATCAc                                              >  4:353606/1‑149 (MQ=255)
                                                                                                      cAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGAAGATAAGGATCATCAc                                              >  3:16485/1‑147 (MQ=255)
                                                                                                                 gTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGAAGATAAGGATCATCACCAGCAGCGAGaa                                  <  4:16485/148‑1 (MQ=255)
                                                                                                                    tcCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGAAGATAAGGATCATCACCAGCAGCgag                                    <  3:73367/143‑1 (MQ=255)
                                                                                                                    tcCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGAAGATAAGGATCATCACCAGCAGCgag                                    >  4:73367/1‑143 (MQ=255)
                                                                                                                    tcCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGAAGATAAGGATCATCACCAGCAGCGAGAAGCt                               >  7:48359/1‑148 (MQ=255)
                                                                                                                      cagcagGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGAAGATAAGGATCATCACCAGCAGCGAGAAGCTg                              <  5:309816/147‑1 (MQ=255)
                                                                                                                      cagcagGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGAAGATAAGGATCATCACCAGCAGCGAGAAGCTGaa                            <  1:364172/149‑1 (MQ=255)
                                                                                                                             cTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGAAGATAAGGATCATCACCAGCAGCGAGAAGCTGAACACATcc                     <  7:224880/149‑1 (MQ=255)
                                                                                                                                 gTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTGGTGGGCGGAAGATAAGGATCATCACCAGCAGCGAGAAGCTGAACACAGCCtgg                  >  8:267186/1‑146 (MQ=255)
                                                                                                                                        gTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGAAGATAAGGATCATCACCAGCAGCGAGAAGCTGAACACATCCTTGTAGTCGGa          <  2:239506/149‑1 (MQ=255)
                                                                                                                                         gcGAAACCAAGGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGAAGATAAGGATCATCACCAGCAGCGAGAAGCTGAACACATCCTTGTAGTCGGAg         >  8:39794/2‑149 (MQ=255)
                                                                                                                                         gcGAAACCAAGGG‑ATTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGAAGATAAGGATCATCACCAGCAGCGAGAAGCTGAACACATCCTTGTAGTCGGAg         >  2:19323/2‑149 (MQ=255)
                                                                                                                                                ccTTGG‑TTTCGCTAATATAGGCAGTTTTGGCAACGGACATGTTCACACATTAGCCACGAGCGGGCGACCCAGCAGGCCTTGGGGGGGGATGATAAGGATCATCACCAGCAGCGAGAAGCTGAACACATCCTTGTAGTCGGAGTTGATCa  <  1:291450/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGAAGATAAGGATCATCACCAGCAGCGAGAAGCTGAACACATCCTTGTAGTCGGAGTTGATCA  >  NC_002947/5532083‑5532375

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: ‑