Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,963,420 G→A 7.7% S171L (TCG→TTG)  PP_2593 ← ferric siderophore ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,963,4200GA 7.7% 148.7 / 2.3 52S171L (TCG→TTG) PP_2593ferric siderophore ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (21/27);  new base A (2/2);  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 = 7.71e-01

GCGCACGTTGCGCCGGTAGTCCACGCCCAGTGAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGCCATGACCATGGTCAGCACCAGGCCGACCAGCAGCAGGAAGTACACGTTGTTTTT  >  NC_002947/2963278‑2963565
                                                                                                                                              |                                                                                                                                                 
gcgcACGTTGCGCCGGTAGTCCACGCCCAGTGAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCAGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCAATGccc                                                                                                                                             >  1:272465/1‑149 (MQ=255)
    aacgtGCGCCGGTAGTCCACGCCCAGTGAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGtt                                                                                                                                         >  1:49693/5‑149 (MQ=255)
           gcCGGTAGTCCACGCCCAGTGAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGaa                                                                                                                                  <  4:237172/149‑1 (MQ=255)
           gcCGGTAGTCCACGCCCAGTGAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGaa                                                                                                                                  >  4:120863/1‑149 (MQ=255)
                  gTCCACGCCCAGTGAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACGGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTggg                                                                                                                           >  1:192993/1‑149 (MQ=255)
                  gTCCACGCCCAGTGAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTgg                                                                                                                            <  3:120863/148‑1 (MQ=255)
                   tCCACGCCCAGTGAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGgt                                                                                                                          <  6:126803/149‑1 (MQ=255)
                   tCCACGCCCAGTGAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGgt                                                                                                                          <  1:233758/149‑1 (MQ=255)
                    ccACGCCCAGTGAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACACCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGgtg                                                                                                                         >  6:125196/1‑149 (MQ=255)
                        gCCCAGTGAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGc                                                                                                                     >  1:136609/1‑149 (MQ=255)
                        gCCCAGTGAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACACCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGc                                                                                                                     <  5:125196/149‑1 (MQ=255)
                          ccAGTGAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTa                                                                                                                       <  5:171441/145‑1 (MQ=255)
                          ccAGTGAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTa                                                                                                                       >  6:171441/1‑145 (MQ=255)
                            aGTGAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGc                                                                                                                                >  7:32015/1‑134 (MQ=255)
                            aGTGAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGc                                                                                                                                <  8:32015/134‑1 (MQ=255)
                             gTGAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCcagc                                                                                                                <  2:192993/149‑1 (MQ=255)
                             gTGAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCcagc                                                                                                                <  5:223996/149‑1 (MQ=255)
                               gAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCTGCACGGCCGAATAGAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCcagca                                                                                                               <  4:176868/148‑1 (MQ=255)
                                     gCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCga                                                                                                        <  2:236607/149‑1 (MQ=255)
                                       cTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCgaga                                                                                                      <  6:163416/149‑1 (MQ=255)
                                          gTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATAGATaa                                                                                                     <  7:199787/147‑1 (MQ=255)
                                             gcgTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGc                                                                                                >  4:104170/1‑149 (MQ=255)
                                               gTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGTGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCg                                                                                              >  1:165230/1‑149 (MQ=255)
                                                    aGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCAATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCt                                                                                         <  2:272465/149‑1 (MQ=255)
                                                          aGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACttt                                                                                   <  3:194268/149‑1 (MQ=255)
                                                               gTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTgg                                                                                                                            <  7:265571/103‑1 (MQ=255)
                                                               gTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTgg                                                                                                                            >  8:265571/1‑103 (MQ=255)
                                                               gTCCAGCACCGGTACGGCCCGGGACAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCATGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCt                                                                              >  8:106631/1‑149 (MQ=255)
                                                                   aGCACCGGTACGGCCCGGGCCAGCAGCACACATACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGAc                                                                          >  7:85827/1‑149 (MQ=255)
                                                                          ttACGGCCCGGTCAATCAGCACACACACGCCCGCCACCAGCAGGTCCGAATTCAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTATCCCAGCAGCATCGAGAACTCTCCGGGGCTGACTTTCAGCTGGACGAACTgg                                                                   <  7:216864/148‑1 (MQ=255)
                                                                           tACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTggg                                                                  <  2:136609/149‑1 (MQ=255)
                                                                                    ggCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTgg                                                         <  8:25022/149‑1 (MQ=255)
                                                                                      ccAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGc                                                        >  6:310674/1‑148 (MQ=255)
                                                                                        agcagcACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGCCAt                                                     <  3:104170/149‑1 (MQ=255)
                                                                                              acacacacGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGCCATGACCAt                                               >  3:206904/1‑149 (MQ=255)
                                                                                               cacacacGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGCCATGACCATg                                              <  2:165230/149‑1 (MQ=255)
                                                                                                acacacGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCAATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGACCTGGGTGAACGTGGCCATGACCATgg                                             >  5:31775/1‑149 (MQ=255)
                                                                                                  acacGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGCCATGACCATg                                              >  6:100323/1‑146 (MQ=255)
                                                                                                  acacGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGCCATGACCATg                                              <  5:100323/146‑1 (MQ=255)
                                                                                                      gcccgccACCAGCACGGCCGAATACAGCTGCAACTGCCGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCCTCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGCCATGACCATGGTCAGCa                                       >  1:19732/1‑149 (MQ=255)
                                                                                                      gcccgccACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGcc                                                       >  8:298776/1‑133 (MQ=255)
                                                                                                      gcccgccACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGcc                                                       <  7:298776/133‑1 (MQ=255)
                                                                                                          gccACCAGCACGGCCGAATACAGCAGCAACTGCGGCAATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGCCATGACCATGGTCAGCACCAg                                   <  6:31775/149‑1 (MQ=255)
                                                                                                               caGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGCCATGACCATGGTCAGCACCAGGCCGa                              >  7:168869/1‑149 (MQ=255)
                                                                                                                            tACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTTGCCATGACCATGGTCAGCACCAGGCCGACCAGCAGCAGGaa                 <  8:85827/149‑1 (MQ=255)
                                                                                                                              cagcagCAACTGCGGCGATGCCCGGTTTAAAGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGCCATGACCATGGt                                            <  8:117986/120‑1 (MQ=255)
                                                                                                                              aagcagcAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGCCATGACCATGGt                                            >  7:117986/2‑120 (MQ=255)
                                                                                                                               agcagcAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGCCATGACCATGGTCAGCACCAGGCCGACCAGCAGCAGGAAGTa              >  8:223646/1‑149 (MQ=255)
                                                                                                                                gcagcaACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGCCATGACCATGGTCAGCAc                                      <  6:297545/124‑1 (MQ=255)
                                                                                                                                gcagcaACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGCCATGACCATGGTCAGCAc                                      >  5:297545/1‑124 (MQ=255)
                                                                                                                                gcagcaACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGCCATGACCATGGTCAGCACCAGGCCGACCAGCAGCAGGAAGTac             <  2:49693/149‑1 (MQ=255)
                                                                                                                                   gcaACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGCCATGACCATGGTCAGCACCAGGCCGACCAGCAGCAGGAAGTACACg          <  4:206904/149‑1 (MQ=255)
                                                                                                                                           ggcgATGCCCGGTTGAACGTAGCCTGGGTGTAGGCCAGCAGCATCGAGAAGTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGCCATGACCATGGTCAGCACCAGGCCGACCAGCAGCAGGAAGTACACGTTGttttt  <  7:106631/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GCGCACGTTGCGCCGGTAGTCCACGCCCAGTGAAATGGCCTGGTCGCGTCCGAGTGCCAGCACGTCCAGCACCGGTACGGCCCGGGCCAGCAGCACACACACGCCCGCCACCAGCACGGCCGAATACAGCAGCAACTGCGGCGATGCCCGGTTGAACGAAGCCTGGGTGTAGCCCAGCAGCATCGAGAACTCGCCGGGGCTGACTTTCAGCTGGACGAACTGGGTGAACGTGGCCATGACCATGGTCAGCACCAGGCCGACCAGCAGCAGGAAGTACACGTTGTTTTT  >  NC_002947/2963278‑2963565

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