Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,907,962:1 +G 7.3% coding (61/1095 nt) spuE ← spermidine/putrescine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,907,9621.G 7.3% 215.8 / 11.1 55coding (61/1095 nt)spuEspermidine/putrescine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (24/27);  new base G (2/2);  total (26/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 = 7.94e-01

AGCAGCTTGCCTTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTCG  >  NC_002947/5907817‑5908107
                                                                                                                                                  |                                                                                                                                                 
agcagcTTGCCTTCCAGGGTTTTGTTGGAGTCGAACACGTTTTGCATAGGCTTGATGCCCGTGGCTTTCTCGAAGTTCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCAt                                                                                                                                               <  2:16598/150‑1 (MQ=255)
    gcTTGCCTTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCa                                                                                                                                           <  1:1848515/150‑1 (MQ=255)
     cTTGCCTTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAg                                                                                                                                          >  1:270641/1‑150 (MQ=255)
      ttGCCTTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGc                                                                                                                                         >  1:493276/1‑150 (MQ=255)
            tCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGGCATGCAGCGTccaa                                                                                                                                   <  1:1053046/150‑4 (MQ=255)
                   tttCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCgcgc                                                                                                                            >  2:122999/1‑150 (MQ=255)
                    ttCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCgcc                                                                                                                           <  2:1540543/150‑1 (MQ=255)
                      cGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATTCCAGCGTCAGGCCCGCgccag                                                                                                                         <  1:1189492/150‑1 (MQ=255)
                       gTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAATTGCACCGTCGGCGCCGCTTGGACGCTGGCATGCAGCGTccaa                                                                                                                                   <  1:343555/139‑4 (MQ=255)
                       gTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAATTGCACCGTCGGCGCCGCTTGGACGCTGGCATGCAGCGTccaa                                                                                                                                   >  2:343555/1‑136 (MQ=255)
                       gTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCgccagc                                                                                                                        >  2:1932721/1‑150 (MQ=255)
                             gTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATc                                                                                                                    >  2:842293/1‑148 (MQ=255)
                             gTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATc                                                                                                                    <  1:842293/148‑1 (MQ=255)
                                gAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGc                                                                                                               >  1:1252782/1‑150 (MQ=255)
                                gAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGc                                                                                                               <  1:122999/150‑1 (MQ=255)
                                gAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGc                                                                                                               <  1:1212688/150‑1 (MQ=255)
                                 aaCACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGcc                                                                                                              <  2:1218777/150‑1 (MQ=255)
                                  acacATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCt                                                                                                             >  1:1756575/1‑150 (MQ=255)
                                      aTCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATGGAAGTGCACCGTCGGCGCCGCTTGGACGCTGc                                                                                                                                                 >  2:757210/1‑110 (MQ=255)
                                      aTCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGc                                                                                                                                                 <  1:757210/110‑1 (MQ=255)
                                      aTCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCg                                                                                                         >  1:1722713/1‑150 (MQ=255)
                                        cTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGa                                                                                                       <  2:1252782/150‑1 (MQ=255)
                                          tGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAAt                                                                                                     <  2:1263750/150‑1 (MQ=255)
                                           gCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGAGGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATa                                                                                                    <  1:1189485/150‑1 (MQ=255)
                                               aGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTg                                                                                                                                                  <  2:970961/100‑1 (MQ=255)
                                               aGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTg                                                                                                                                                  >  1:970961/1‑100 (MQ=255)
                                                  cTTGATGCCCGTGGCGTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTACGGAATACAGAAAt                                                                                             >  1:1050604/1‑150 (MQ=255)
                                                      aTGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCtt                                                                                                            >  2:886450/1‑131 (MQ=255)
                                                      aTGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCtt                                                                                                            <  1:886450/131‑1 (MQ=255)
                                                          ccGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGACGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCCGGCCCGCGCCAGCCAACAAGGCCTTGCGGAATACAGAAATAGACAAGt                                                                                     >  1:1270660/1‑150 (MQ=255)
                                                          ccGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCgcca                                                                                                                          >  1:101558/1‑113 (MQ=255)
                                                          ccGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCgcca                                                                                                                          <  2:101558/113‑1 (MQ=255)
                                                            gTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTgg                                                                                   <  1:1093324/150‑1 (MQ=255)
                                                                 tttCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGc                                                                                                                                  >  1:1519183/1‑98 (MQ=255)
                                                                 tttCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGc                                                                                                                                  <  2:1519183/98‑1 (MQ=255)
                                                                     tcGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTca                                                                          <  2:493276/150‑1 (MQ=255)
                                                                       gAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTgggg                                                                                 >  1:1028145/1‑141 (MQ=255)
                                                                       gAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTgggg                                                                                 <  2:1028145/141‑1 (MQ=255)
                                                                        aaGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACaa                                                                       <  1:1932721/150‑1 (MQ=255)
                                                                             cgcgAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGc                                                                                                          <  2:382477/110‑1 (MQ=255)
                                                                             cgcgAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGc                                                                                                          >  1:382477/1‑110 (MQ=255)
                                                                               cgAGAGTTGTGCGGGCGATTTAGTCGGTCCAGTACTAGAAGTGCTCCGGCGGCGCCGCTTGGACGCTGCATGGTAGCGTCAGGCCCGTGCTAGCTTTCACTGGCTTGGGGAGTACAGAAATAGACAGGTGGGGTTCCTCACAATCAGttc                                                                <  2:270641/150‑3 (MQ=255)
                                                                                              cGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCa                                                                    >  2:999005/1‑131 (MQ=255)
                                                                                              cGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCa                                                                    <  1:999005/131‑1 (MQ=255)
                                                                                              cGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAAt                                                 >  2:397757/1‑150 (MQ=255)
                                                                                                    aGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGcc                                           >  2:171524/1‑150 (MQ=255)
                                                                                                        ggACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGcac                                       >  1:864336/1‑150 (MQ=255)
                                                                                                        ggACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAAACAGTGCCTGGGTAGCGACAATGCTGCCGcac                                       >  1:712695/1‑150 (MQ=255)
                                                                                                           ccAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCaa                                                                                                                  <  2:15415/72‑1 (MQ=255)
                                                                                                           ccAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCaa                                                                                                                  >  1:15415/1‑72 (MQ=255)
                                                                                                                           aCCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCaa                    <  2:864336/150‑1 (MQ=255)
                                                                                                                                    gccgcTTAGACGCTGCATCCCAGCGACAGGCCGGCGCAAGCCATTAAGGCCTTGCGGAATACAAAAATAGACAAAGGGAGGTTCACGCAATCAGGACATGGGTAGCGACAATGCCGCCGCACACGCAAATACGGCGCGCAACTTAACTTc           <  2:1641621/150‑1 (MQ=255)
                                                                                                                                      cgcTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGc         <  1:171524/150‑1 (MQ=255)
                                                                                                                                        cTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAg       >  1:451394/1‑150 (MQ=255)
                                                                                                                                          ctgACGCTGGCATGCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAgat     >  2:937797/3‑148 (MQ=255)
                                                                                                                                           ggACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCtt    <  2:451394/150‑1 (MQ=255)
                                                                                                                                             aCGCTGCATGCCCGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGGCTACAGACATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAACACCGGCGCGCACCTTAACTTCGCAGCTTCg  >  2:515709/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGCAGCTTGCCTTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTCG  >  NC_002947/5907817‑5908107

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