Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,926,722 G→A 9.9% A410V (GCG→GTG)  PP_5195 ← iron ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,926,7220GA 9.9% 103.6 / 2.5 41A410V (GCG→GTG) PP_5195iron ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (16/19);  new base A (2/2);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

ACCTTGAAAAACAATCTTGCGCCGCCAACACCCAGGCTGCGCGAGGCTTCAGGGAGTGAAGGGCGGATGCGCTCCAGGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGAAGGCCAGCATGATCGATACCGCCAACACCGAACCGGGCAGGGCATACCCCAGATTGGC  >  NC_002947/5926584‑5926867
                                                                                                                                          |                                                                                                                                                 
acCTTGAAAAACAATCTTGCGCCGCCAACACCCAGGCTGCGCGAGGCTTCAGGGAGTGAAGGGCGGATGCGCTCCAGGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTa                                                                                                                                        <  1:1117583/150‑1 (MQ=255)
  cTTGAAAAACAATCTTGCGCCGCCAACACCCAGGCTGCGCGAGGCTTCAGGGAGTGAAGGGCGGATGCGCTCCAGGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTAcc                                                                                                                                      >  2:903002/1‑150 (MQ=255)
    tGAAAAACAATCTTGCGCCGCCAACACCCAGGCTGCGCGAGGCTTCAGGGAGTGAAGGGCGGATGCGCTCCAGGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTa                                                                                                                                    <  1:109710/150‑1 (MQ=255)
       aaaaCAATCTTGCGCCGCCAACACCCAGGCTGCGCGAGGCTTCAGGGAGTGAAGGGCGGATGCGCTCCAGGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCCCCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTagca                                                                                                                                 >  1:414025/1‑150 (MQ=255)
        aaaCAATCTTGCGCCGCCAACACCCAGGCTGCGCGAGGCTTCAGGGAGTGAAGGGCGGATGCGCTCCAGGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTagcag                                                                                                                                >  1:81012/1‑150 (MQ=255)
              atttGCGCCGCCAACACCCAGGCTGCGCGAGGCTTCAGGGAGTGAAGGGCGGATGCGCTCCAGGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCgg                                                                                                                          >  1:1062452/3‑150 (MQ=255)
                   cgccgcCAACACCCAGGCTGCGCGAGGCTTCAGGGAGTGAAGGGCGGATGCGCTCCAGGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGc                                                                                                                     <  1:1024705/150‑1 (MQ=255)
                    gccgccAACACCCAGGCTGCGCGAGGCTTCAGGGAGTGAAGGGCGGATGCGCTCCAGGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGcc                                                                                                                    <  2:162830/150‑1 (MQ=255)
                       gccAACACCCAGGCTGCGCGAGGCTTCAGGGAGTGAAGGGCGGATGCGCTCCAGGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGgc                                                                                                                 >  2:780534/1‑150 (MQ=255)
                               ccAGGCTGCGCGAGGCTTCAGGGAGTGAAGGGCGGATGCGCTCCAGGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCa                                                                                                         >  2:98038/1‑150 (MQ=255)
                                      gcgcGAGGCTTCAGGGAGTGAAGGGCGGATGCGCTCCAGGCTGCTCTCCAGCGGCTCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACt                                                                                                  <  2:919778/150‑1 (MQ=255)
                                         cgAGGCTTCAGGGAGTGAAGGGCGGATGCGCTCCAGGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGa                                                                                               <  1:768688/150‑1 (MQ=255)
                                               ttCAGGGAGTGAAGGGCGGATGCGCTCCAGGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGCAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACgg                                                                                         >  2:215557/1‑150 (MQ=255)
                                                  aGGGAGTGAAGGGCGGATGCGCTCCAGGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGAt                                                                                      <  2:839183/150‑1 (MQ=255)
                                                  aGGGAGTGAAGGGCCGATGCGCTCCAAGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACCGGACAACGGGAt                                                                                      <  2:503638/150‑1 (MQ=255)
                                                          aaGGGCGGATGCGCTCCAGGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAggt                                                                              <  1:98038/150‑1 (MQ=255)
                                                           aGGGCGGATGCGCTCCAGGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAggtg                                                                             >  2:501929/1‑150 (MQ=255)
                                                                    gcgcTCCAGGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAg                                                                    <  2:1039868/150‑1 (MQ=255)
                                                                    gcgcTCCAGGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGCTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAg                                                                    <  2:7935/150‑1 (MQ=255)
                                                                             gctgctCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGaa                                                           >  2:545407/1‑150 (MQ=255)
                                                                             gctgctCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACACCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGaa                                                           >  1:728262/1‑150 (MQ=255)
                                                                             gctgctCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACACCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGaa                                                           >  1:729642/1‑150 (MQ=255)
                                                                                    tcCAGCGGCCCGTAGGCCACCGCGATCAAGTTCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGAAGGCCAGc                                                    >  1:203823/1‑150 (MQ=255)
                                                                                        gCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGAAGGCCAGCATGa                                                <  2:81012/150‑1 (MQ=255)
                                                                                              cGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGc                                                                                                       >  1:100430/1‑89 (MQ=255)
                                                                                              cGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGc                                                                                                       <  2:100430/89‑1 (MQ=255)
                                                                                                       ccGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGTACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGAAGGCCAGCATGATCGATACCGCCAaca                                 <  1:501929/150‑1 (MQ=255)
                                                                                                        cgcgATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGAAGGCCAGCATGATCGATACCGCCAacac                                >  1:640246/1‑150 (MQ=255)
                                                                                                              gAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTaa                                                                >  2:165568/1‑112 (MQ=255)
                                                                                                              gAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTaa                                                                <  1:165568/112‑1 (MQ=255)
                                                                                                                   gcACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGAAGGCCAGCATGATCGATACCGCCAACACCGAACCgggaa                     <  2:856905/150‑3 (MQ=255)
                                                                                                                       cAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGAAGGCCAGCATGATCGATACCGCCAACACCGAACCgggcagggc                 <  1:903002/150‑1 (MQ=255)
                                                                                                                              gCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGAAGGCCAGCATGATc                                              <  2:584365/114‑1 (MQ=255)
                                                                                                                              gCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGAAGGCCAGCATGATc                                              >  1:584365/1‑114 (MQ=255)
                                                                                                                              gCCAGCAGCAACACCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGAAGGCCAGCATGATCGATACCGCCAACACCGAACCGGGCAGGGCATAcccc          <  2:728262/150‑1 (MQ=255)
                                                                                                                              gCCAGCAGCAACACCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGAAGGCCAGCATGATCGATACCGCCAACACCGAACCGGGCAGGGCATAcccc          <  2:729642/150‑1 (MQ=255)
                                                                                                                                  gcagcaACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGAAGGCCAGCATGATCGATACCGCCAACACCGAACCGGGCAGGGCATACCCCAGAt      >  1:738946/1‑150 (MQ=255)
                                                                                                                                   cagcaACTCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGAAGGCCAGCATGATCGATACCGCCAACAcc                               <  1:846059/124‑1 (MQ=255)
                                                                                                                                   cagcaACTCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGAAGGCCAGCATGATCGATACCGCCAACAcc                               >  2:846059/1‑124 (MQ=255)
                                                                                                                                   cagcaACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGAAGGCCAGCATGATCGATACCGCCAACACCGAACCGGGCAGGGCATACCCCAGAtt     >  2:123176/1‑150 (MQ=255)
                                                                                                                                      caACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGAAGGCCAGCATGATCGATACCGCCAACACCGAACCGGGCAGGGCATACCCCAGATTggc  <  1:215557/150‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
ACCTTGAAAAACAATCTTGCGCCGCCAACACCCAGGCTGCGCGAGGCTTCAGGGAGTGAAGGGCGGATGCGCTCCAGGCTGCTCTCCAGCGGCCCGTAGGCCACCGCGATGAAGCGCACCAGGTAGGCCAGCAGCAACGCCGAGAGGCTACCTAGCAGCAGCGGCTTGCCGGCGCCGCCCAGCCAACTGGACAACGGGATGACCAGGTGGTTGTCGAGGTAACTGAAGGCCAGCATGATCGATACCGCCAACACCGAACCGGGCAGGGCATACCCCAGATTGGC  >  NC_002947/5926584‑5926867

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