Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,532,228 C→A 9.8% K15N (AAG→AAT) ‡ 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,2280CA 9.8% 148.8 / 7.3 51K15N (AAG→AAT) ‡braEhigh‑affinity branched‑chain amino acid ABC transporter permease BraE
Reads supporting (aligned to +/- strand):  ref base C (21/25);  new base A (2/3);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.92e-01

TCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGAAGATAAGGATCATCACCAGCAGCGAGAAGCTGAACACATCCTTGTAGTCGGAGTTGATCA  >  NC_002947/5532080‑5532375
                                                                                                                                                    |                                                                                                                                                    
tCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAAccc                                                                                                                                                      <  6:341067/149‑1 (MQ=255)
   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                                                                                                                                           <  2:346993/149‑2 (MQ=255)
          gACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCAAGGG‑TTTCga                                                                                                                                           <  6:70798/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)
                                                                                                                                                    |                                                                                                                                                    
TCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGGCCAGCGACAATGGTCTCCAGCAGGCTGCGTTTAAGGTCGAAACCCTTGG‑TTTCGCTAATAGAGGCAGTTTTGGCAACGGACATGTTCACACCTTAGCCACGAGCGGGCGACCCAGCAGGCCTTGTGGGCGGAAGATAAGGATCATCACCAGCAGCGAGAAGCTGAACACATCCTTGTAGTCGGAGTTGATCA  >  NC_002947/5532080‑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: ‑