Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 322,879 A→G 8.5% L455P (CTG→CCG) ‡ PP_0267 ← ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947322,8790AG 8.5% 131.8 / 2.9 47L455P (CTG→CCG) ‡PP_0267ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base A (21/22);  new base G (2/2);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

CCGGGTCTCTTCCGGGTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCGAGAACCTGCCACTGATCGTTCAGGCGAATCTGGTCCTGCAGATAAACC  >  NC_002947/322737‑323019
                                                                                                                                              |                                                                                                                                            
ccGGGTCTCTTCCGGGTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTgg                                                                                                                                       <  1:537820/150‑1 (MQ=255)
 cGGGTCTCTTCCGGGTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTggg                                                                                                                                      >  1:796643/1‑150 (MQ=255)
       tctTCCGGGTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGtg                                                                                                                                >  2:431532/1‑150 (MQ=255)
            cGGGTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGtgctgtt                                                                                                                           >  1:38271/1‑150 (MQ=255)
               gTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGtgctgttgct                                                                                                                        >  2:758331/1‑150 (MQ=255)
               gTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGtg                                                                                                                                <  2:334587/142‑1 (MQ=255)
               gTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGtg                                                                                                                                >  1:334587/1‑142 (MQ=255)
                           ggTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGAAGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTc                                                                                                            <  1:130660/150‑1 (MQ=255)
                                    cagttgcccTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGcc                                                                                                   >  2:596028/7‑150 (MQ=255)
                                    cTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGcc                                                                                                   >  2:341467/1‑150 (MQ=255)
                                      gtgCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGt                                                                                                                                 >  2:309236/1‑118 (MQ=255)
                                      gtgCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGt                                                                                                                                 <  1:309236/118‑1 (MQ=255)
                                          cccTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCtt                                                                                             >  2:595470/1‑150 (MQ=255)
                                            cTGGGGGGATGCCGATCCAGCCGCCTCCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCttgt                                                                                           <  2:781720/150‑1 (MQ=255)
                                                 gTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTtggtg                                                                                      <  1:200580/150‑1 (MQ=255)
                                                      gCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTc                                                                                 >  2:742479/1‑150 (MQ=255)
                                                         gATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCAc                                                                              >  1:220589/1‑150 (MQ=255)
                                                          aTCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCTGACGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACt                                                                             >  1:209526/1‑150 (MQ=255)
                                                           tCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACtt                                                                            <  2:488400/150‑1 (MQ=255)
                                                               gccgccACCCACGGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCTTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGaa                                                                        <  2:52046/150‑1 (MQ=255)
                                                               gccgccACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGaa                                                                        <  1:742479/150‑1 (MQ=255)
                                                                ccgccACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAAc                                                                       >  2:476679/1‑150 (MQ=255)
                                                                  gccACCCACCGGTGGGAAAGTCTTGCTAAAGGCGGCGTAAACCTAAGGGTCCCGCCATGGCGTTTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTg                                                                     <  2:209526/150‑1 (MQ=255)
                                                                     acccaccGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTc                                                                  >  2:593973/1‑150 (MQ=255)
                                                                                  aaaGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCTGGGCGTggtc                                                                                                                                     >  1:368356/1‑68 (MQ=255)
                                                                                  aaaGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCTGGGCGTggtc                                                                                                                                     <  2:368356/70‑3 (MQ=255)
                                                                                  aaaGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCg                                                     >  1:63153/1‑150 (MQ=255)
                                                                                      tCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCga                                                 <  2:796643/150‑1 (MQ=255)
                                                                                      tCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCCGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCga                                                 <  1:245496/150‑1 (MQ=255)
                                                                                           cTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCGAGAAcc                                            <  1:596028/150‑1 (MQ=255)
                                                                                              cACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCGAGAACCTGc                                         <  2:38271/150‑1 (MQ=255)
                                                                                                   ggCGTAGATCGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCCGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCGAGAACCTGCCACTg                                    <  1:431532/150‑1 (MQ=255)
                                                                                                    gCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCGAGAACCTGCCACTGa                                   >  1:168560/1‑150 (MQ=255)
                                                                                                       tAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCGAGAACCTGCCACTGATCg                                >  1:372250/1‑150 (MQ=255)
                                                                                                           aCGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCGAGAACCTGCCACTGATCGTTCa                            <  2:168560/150‑1 (MQ=255)
                                                                                                                 ggTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCGAGAACCTGCCACTGATCGTTCAGGCGaa                      <  2:220589/150‑1 (MQ=255)
                                                                                                                  gTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCACGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCGAGAACCTGCCACTGATCGTTCAGGCGAAt                     <  1:787428/150‑1 (MQ=255)
                                                                                                                     gcgcCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTg                                                                     <  1:310895/99‑1 (MQ=255)
                                                                                                                     gcgcCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTg                                                                     >  2:310895/1‑99 (MQ=255)
                                                                                                                     gcgcCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCGAGAACCTGCCACTGATCGTTCAGGCGAATCTg                  <  1:593973/150‑1 (MQ=255)
                                                                                                                            ggCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCGAGAACCTGCCACTGATCGTTCAGGCGAATCTGGTCCTGc           <  1:341467/150‑1 (MQ=255)
                                                                                                                                 gtAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCGAGAACCTGCCACTGATCGTTCAGGCGAATCTGGTCCTGCAGATa      >  2:349707/1‑150 (MQ=255)
                                                                                                                                  cagcccACGCCTGGGCGTGGTCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCGAGAACCTGCCACTGATc                                 >  1:459139/5‑122 (MQ=255)
                                                                                                                                  cagcccACGCCTGGGCGTGGTCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCGAGAACCTGCCACTGATc                                 <  2:459139/118‑1 (MQ=255)
                                                                                                                                     aCCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTAGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCGAGAACCTGCCACTGATCGTTCAGGCGAATCTGGTCCTGCAGATAAAcc  >  2:578068/1‑150 (MQ=255)
                                                                                                                                        aCGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCGAGaa                                              >  2:115199/1‑103 (MQ=255)
                                                                                                                                        aCGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCGAGaa                                              <  1:115199/103‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                            
CCGGGTCTCTTCCGGGTCGGTTTTATTGGTGTTGCCCTGTGCCCCTGGGGTGATGCCGATCAAGCCGCCACCCACCGGTGAGAAAGTCTTGCTCCACGAGGCGTAGAACGAATGGTCGCGCCATGGCGTGTAGACCACGCCCAGGCGTGGGCTGGTGCTGTTGCTGTCCTGCGGTTCGCGCAAGCCATCCTTGTTGGTGGTTTCCACTTCGAACTGGTCAAAGCGCACGCCGGCGAGAACCTGCCACTGATCGTTCAGGCGAATCTGGTCCTGCAGATAAACC  >  NC_002947/322737‑323019

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