Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,131,842 G→A 8.8% R38R (CGC→CGT)  PP_3635 ← sulfonate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,131,8420GA 8.8% 124.9 / 2.6 46R38R (CGC→CGT) PP_3635sulfonate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (17/24);  new base A (2/2);  total (20/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 = 7.58e-01

CTACCAGCAAGCCCAACGGCACGCCCATCAGCAGCGCCAGCAGCAGCCCAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCGACCAGCCCAAGGCCGGCATGTACGACATGGGTGCGCATCAGACCGCAAACTCCTCGCTACGCTCGAAGTTGTC  >  NC_002947/4131697‑4131990
                                                                                                                                                  |                                                                                                                                                    
cTACCAGCAAGCCCAACGGCACGCCCATCAGCAGCGCCAGCAGCAGCCCAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAgcgc                                                                                                                                                   <  4:57252/149‑1 (MQ=255)
cTACCAGCAAGCCCAACGGCACGCCCATCAGCAGCGCCAGCAGCAGCCCAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAgcgc                                                                                                                                                   <  8:346183/149‑1 (MQ=255)
cTACCAGCAAGCCCAACGGCACGCCCATCAGCAGCGCCAGCAGCAGCCCAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAgcgc                                                                                                                                                   <  5:249140/149‑1 (MQ=255)
 tACCAGCAAGCCCAACGGCACGCCCATCAGCAGCGCCAGCAGCAGCCCAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAgcgcg                                                                                                                                                  >  6:370626/1‑149 (MQ=255)
 tACCAGCAAGCCCAACGGCACGCCCATCAGCAGCGCCAGCAGCAGCCCAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCACCTCGACCAGGCTGGCCAGGGTTGCCGCCAGCGAAAAgcgcg                                                                                                                                                  >  4:260610/1‑149 (MQ=255)
   ccAGCAAGCCCAACGGCACGCCCATCAGCAGCGCCAGCAGCAGCCCAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAACGCGCGcc                                                                                                                                                >  5:198331/1‑149 (MQ=255)
         aGCCCAACGGCACGCCCATCAGCAGCGCCAGCAGCAGCCCAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAgg                                                                                                                                          >  7:310048/1‑149 (MQ=255)
          gCCCAACGGCACGCCCATCAGCAGCGCCAGCAGCAGCCCAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GACCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGCTTGCCGCCGGCGAAAAGCGCGCCGACAGGc                                                                                                                                         >  6:194743/1‑149 (MQ=255)
               aCGGCACGCCCATCAGCAGCGCCAGCAGCAGCCCAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCg                                                                                                                                    <  2:62978/149‑1 (MQ=255)
                  gCACGCCCATCAGCAGCGCCAGCAGCAGCCCAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGcc                                                                                                                                 <  5:321477/149‑1 (MQ=255)
                          aTCAGCAGCGCCAGCAGCAGCCCAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGaa                                                                                                                         <  5:370626/149‑1 (MQ=255)
                           tCAGCAGCGCCAGCAGCAGCCCAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCt                                                                                                                                <  5:368529/141‑1 (MQ=255)
                           tCAGCAGCGCCAGCAGCAGCCCAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCt                                                                                                                                >  6:368529/1‑141 (MQ=255)
                             atcagcGCCAGCAGCAGCCCAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACa                                                                                                                       <  2:321714/146‑1 (MQ=255)
                              gcagcGCCAGCAGCAGCCCAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGt                                                                                                                     >  2:230730/1‑149 (MQ=255)
                                  cgcCAGCAGCAGCCCAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCAACAGGCCGTCGGCCTGGCCGAACAGTGCTa                                                                                                                 <  8:271506/149‑1 (MQ=255)
                                           cagcCCAATCAGGATGCGCTTGAGGCTGACCCAGAGATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGGTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGccc                                                                                                        <  7:159351/149‑1 (MQ=255)
                                              cccAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCccac                                                                                                      >  2:293138/1‑148 (MQ=255)
                                                cAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCAt                                                                                                   <  5:170130/149‑1 (MQ=255)
                                                     aGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACa                                                                                               <  1:230730/148‑1 (MQ=255)
                                                     aGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACa                                                                                               <  6:198331/148‑1 (MQ=255)
                                                     aGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACa                                                                                               <  5:247498/148‑1 (MQ=255)
                                                     aGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACa                                                                                               <  8:310048/148‑1 (MQ=255)
                                                      ggATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGc                                                                                                                                         >  2:260523/1‑105 (MQ=255)
                                                      ggATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGc                                                                                                                                         <  1:260523/105‑1 (MQ=255)
                                                      ggATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGc                                                                                                                   >  4:255810/1‑127 (MQ=255)
                                                      ggATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGc                                                                                                                   <  3:255810/127‑1 (MQ=255)
                                                            gcTTGAGGCTGACCCAGATATGCCCATACACCTC‑TCCCAGACCCAGCAGCTCGACCACGCTGGCCAGGGTTGCCGCAGGCAAAAAGACAGCCGACACGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCCTAACACccacacg                                                                                       >  2:265889/1‑147 (MQ=255)
                                                                   gCTGACCCAGATATGCCCATAGACCTCGGCCATCCCCCAGCACCTCGACCAGGCTGGCCAGGTTTCCCGCCGGCGAAAACCCCCCCGACAGCCCTTCGGCCTGCCCGAACAGTCATAAACCGGCCCACCATAACCGCAACAAGGCCAg                                                                                  >  1:198929/1‑148 (MQ=255)
                                                                         ccAGATATGCCCATAGACCTC‑GCCCTGCCCCACCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCGa                                                                          >  3:159727/1‑149 (MQ=255)
                                                                          cAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCGAc                                                                         <  7:209826/149‑1 (MQ=255)
                                                                                  cccATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGcc                                                                                    >  5:201288/1‑130 (MQ=255)
                                                                                  cccATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGcc                                                                                    <  6:201288/130‑1 (MQ=255)
                                                                                     attgACCTA‑GCCCTGCCCCAGCAGGTCGACCAGGCTGGCGAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGATCAGTGCTACAGCTGCACACCATATCAGCAACAAGGCCAGCAACGGTACCTGCCCAAGGc                                                              <  4:159727/146‑1 (MQ=255)
                                                                                           cTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCGACCAGCCCAAGGCCGGCAt                                                        >  4:25358/1‑149 (MQ=255)
                                                                                               gCCCTGGCCCAGCAGCTCGACCAGGCTGGCCATGGTTGCCGCCTGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGTACAGTGCTACACCGGCCCACTATAACAGCAACATGGCCAGCAACCCGACCAGCCCAATGCCGGCATGTa                                                     <  1:265889/149‑1 (MQ=255)
                                                                                                 ccTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCGACCAGCCCAAGGCCGGCATGTAc                                                    >  7:150394/1‑148 (MQ=255)
                                                                                                       ccAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCGGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACAATAACAGCAAGAAGGCCAGCAACCCGACCAGCGCAAGGCCGGCATGTACGACATgg                                             <  3:260610/149‑1 (MQ=255)
                                                                                                           cagcTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCGACCAGCCCAAGGCCGGCATGTACGACATGGGTgc                                         <  3:25358/149‑1 (MQ=255)
                                                                                                             gcTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCGACCAGCCCAAGGCCGGCATGTACGACATGGGTgcg                                        >  2:52817/1‑148 (MQ=255)
                                                                                                              cTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCGACCAGCCCAAGGCCGGCATGTACGACATGGGTGCGCa                                      >  7:308592/1‑149 (MQ=255)
                                                                                                                   ccAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGTCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCGACCAGCCCAAGGCCGGCATGTACGACATGGGTGCGCATCAGa                                 <  5:194743/149‑1 (MQ=255)
                                                                                                                    cAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCGACCAGCCCAAGGCCGGCATGTACGACATGGGTGCGCATCAGAc                                <  1:52817/149‑1 (MQ=255)
                                                                                                                    cAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCGACCAGCCCAAGGCCGGCATGTACGACATGGGTGCGCATCAGAc                                >  8:311346/1‑149 (MQ=255)
                                                                                                                           gCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCGACCAGCCCAAGGCCGGCATGTACGACATGGGTGCGCATCAGACCGCAAAc                         >  7:119931/1‑149 (MQ=255)
                                                                                                                                     ccgccgGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCGACCAGCCCAAGGCCGGCATGTACGACATGGGTGCGCATCAGACCGCAAACTCGTCGCTAc               <  2:198929/149‑1 (MQ=255)
                                                                                                                                            cGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCTACCAGCCCAAGGCCGGCATGTACGACATGGGTGCGCATCAGACCGCAAACTCCTCGCTACGCTCGaa        <  1:293138/149‑1 (MQ=255)
                                                                                                                                             gAAAAACGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCGACCAGCCCAAGGCCGGCATGTACGACATGGGTGCGCAt                                     <  2:91855/119‑1 (MQ=255)
                                                                                                                                             gAAAAACGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCGACCAGCCCAAGGCCGGCATGTACGACATGGGTGCGCAt                                     <  2:96572/119‑1 (MQ=255)
                                                                                                                                             gAAAAACGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCGACCAGCCCAAGGCCGGCATGTACGACATGGGTGCGCAt                                     >  1:96572/1‑119 (MQ=255)
                                                                                                                                             gAAAAACGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCGACCAGCCCAAGGCCGGCATGTACGACATGGGTGCGCAt                                     >  1:91855/1‑119 (MQ=255)
                                                                                                                                                  gcgcgcCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCGACCAGCCCAAGGCCGGCATGTACGACATGGGTGCGCATCAGACCGCAAACTCCTCGCTACGCTCGAAGTTGTc  <  8:308592/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTACCAGCAAGCCCAACGGCACGCCCATCAGCAGCGCCAGCAGCAGCCCAATCAGGATGCGCTTGAGGCTGACCCAGATATGCCCATAGACCTC‑GCCCTGCCCCAGCAGCTCGACCAGGCTGGCCAGGGTTGCCGCCGGCGAAAAGCGCGCCGACAGGCCGTCGGCCTGGCCGAACAGTGCTACACCGGCCCACCATAACAGCAACAAGGCCAGCAACCCGACCAGCCCAAGGCCGGCATGTACGACATGGGTGCGCATCAGACCGCAAACTCCTCGCTACGCTCGAAGTTGTC  >  NC_002947/4131697‑4131990

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