Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,043,569 C→A 10.5% D76E (GAC→GAA)  pstS → phosphate ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,043,5690CA10.5% 135.1 / 4.9 45D76E (GAC→GAA) pstSphosphate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (18/22);  new base A (2/3);  total (20/25)
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.94e-01

CTGCAGAAAACGTTCGGTTGACCGGTTCGGGCGCCAGCTTCCCTGCACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACATTGCCAAGGTTGGCGAAGGCGTGCAATTGCTGCCAATGACCGCCGGTGAAATCGTCCTGGCTTACAACCTGCCGGGTAACCCTA  >  NC_002947/3043427‑3043711
                                                                                                                                              |                                                                                                                                              
ctgcAGAAAACGTTCGGTTGACCGGTTCGGGCGCCAGCTTCCCTGCACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCt                                                                                                                                           >  6:250964/1‑148 (MQ=255)
ctgcAGAAAACGTTCGGTTGACCGGTTCGGGCGCCAGCTTCCCTGCACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCt                                                                                                                                           <  8:181191/148‑1 (MQ=255)
ctgcAGAAAACGTTCGGTTGACCGGTTCGGGCGCCAGCTTCCCTGCACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCt                                                                                                                                           >  7:181191/1‑148 (MQ=255)
ctgcAGAAAACGTTCGGTTGACCGGTTCGGGCGCCAGCTTCCCTGCACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTg                                                                                                                                          <  1:308474/149‑1 (MQ=255)
     gAAAACGTTCGGTTGACCGGTTCGGGCGCCAGCTTCCCTGCACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACaa                                                                                                                                     <  1:265226/149‑1 (MQ=255)
     gAAAACGTTCGGTTGACCGGTTCGGGCGCCAGCTTCCCTGCACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGAAGTCCTGAAcgc                                                                                                                                     <  1:383333/149‑3 (MQ=255)
     gAAAACGTTCGGTTGACCGGTTCGGGCGCCAGCTTCCCTGCACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGAAGTCCTGAAcgc                                                                                                                                     >  2:383333/1‑147 (MQ=255)
     gAAAACGTTCGGTTGACCGGTTCGGGCGCCAGCTTCCCTGCACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGAAGTCCTGAAcgc                                                                                                                                     <  7:93932/149‑3 (MQ=255)
       aaaCGTTCGGTTGACCGGTTCGGGCGCCAGCTTCCCTGCACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGAAGTCCTGAAcgccc                                                                                                                                   <  4:323737/149‑5 (MQ=255)
            ttCGGTTGACCGGTTCGGGCGCCAGCTTCCCTGCACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTc                                                                                                                              <  5:250964/149‑1 (MQ=255)
               ggTTGACCGGTTCGGGCGCCAGCTTCCCTGCACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGa                                                                                                                            <  4:99534/148‑1 (MQ=255)
               ggTTGACCGGTTCGGGCGCCAGCTTCCCTGCACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGAc                                                                                                                           <  6:320030/149‑1 (MQ=255)
                                  cAGCTTCCCTGCACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCgg                                                                                                        >  8:318305/1‑149 (MQ=255)
                                   aGCTTCCCTGCACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGc                                                                                                       >  5:342003/1‑149 (MQ=255)
                                         ccTGCACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGa                                                                                                  <  6:244918/148‑1 (MQ=255)
                                          cTGCACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGa                                                                                                                                         >  8:71150/1‑108 (MQ=255)
                                          cTGCACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGa                                                                                                                                         <  7:71150/108‑1 (MQ=255)
                                             cACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGa                                                                                              <  1:137029/148‑1 (MQ=255)
                                                 gATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGaa                                                                                             >  5:208396/1‑145 (MQ=255)
                                                 gATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGaa                                                                                             <  6:208396/145‑1 (MQ=255)
                                                       cctgacctgGTTCAAGGACTTCAGCAAGAACACGGCTGGCTTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACATTg                                                                                   <  8:161526/149‑1 (MQ=255)
                                                       cctgacctgGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACAtt                                                                                    <  8:228496/148‑1 (MQ=255)
                                                        ctgacctgGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGa                                                                                                  <  6:152621/133‑1 (MQ=255)
                                                        ctgacctgGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGa                                                                                                  >  5:152621/1‑133 (MQ=255)
                                                             ctgGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACATTGCCAAgg                                                                             >  5:369695/1‑149 (MQ=255)
                                                              tgGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGc                                                                                                       >  3:88474/1‑122 (MQ=255)
                                                              tgGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGc                                                                                                       <  4:88474/122‑1 (MQ=255)
                                                                      ggACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACATTGCCAAGGTTGGCGAAg                                                                    >  2:93289/1‑149 (MQ=255)
                                                                              gCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGTGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACATTGCCAAGGTTGGCGAAGGCGTGCaa                                                            >  5:216643/1‑149 (MQ=255)
                                                                                 aGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACATTGCCAAGGTTGGCGAAGGCGTGCAATtg                                                         <  6:127483/149‑1 (MQ=255)
                                                                                     accGGCTGGCGTCTCCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACATTGCCAAGGTTGGCGAAGGCGTGCAATtgctgc                                                     <  8:331446/147‑1 (MQ=255)
                                                                                        ggCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACATTGCCAAGGTTGGCGAAGGCGTGCAATTGCTGCCaa                                                  <  1:257384/149‑1 (MQ=255)
                                                                                               gTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACATTGCCAAGGTTGGCGAAGGCGTGCAATTGCTGCCAATGAccgc                                           >  2:401993/1‑149 (MQ=255)
                                                                                                    cGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACATTGCCAAGGTTGGCGCAGGCGGGCAATTGCTGCCAATGACCGCCGGt                                       >  3:148479/1‑148 (MQ=255)
                                                                                                     gTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCTTCGACTTCGCCGCGAGCGGCTCGGCCAGGAGCGAAGCAGACATTGACAAGGTTGGCAGAGGCGTGCAATGGCTGCCAATGCCAGCAGGTg                                      >  8:358778/1‑148 (MQ=255)
                                                                                                      tCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACATTGCCAAGGTTGGCGAAGGCGTGCAATTGCTGCCAATGACCGCCGGTGa                                     >  8:203772/1‑148 (MQ=255)
                                                                                                             ccAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACAtt                                                                                    >  6:354921/1‑94 (MQ=255)
                                                                                                             ccAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACAtt                                                                                    <  5:354921/94‑1 (MQ=255)
                                                                                                             ccAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACATTGCCAAGGTTGGCGAAGGCGTGCAATTGCTGCCAATGACCGCCGGTGAAATCGTc                              <  7:318305/148‑1 (MQ=255)
                                                                                                              cAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACATTGCCAAGGTTGGCGAAGGCGTGCAATTGCTGCCAATGACCGCCGGTGAAATCGTCCt                            >  7:18716/1‑149 (MQ=255)
                                                                                                                 tCCAAGGGCAGCGGCGCGGGCGTGCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACATTGCCAAGGTTGGCGAAGGCGTGCAATTGCTGCCAATGACCGCCGGTGAAATCGTCCTGGc                         <  6:342003/149‑1 (MQ=255)
                                                                                                                   cAAGGGCAGCGGCGCGGGCGTTAAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACATTGCCAAGGTTGGCGAAGGCGTGCAATTGCTGCCAATGACCGCCGGTGAAATCGTCCTGGCt                        >  7:59681/1‑148 (MQ=255)
                                                                                                                         cAGCGTCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACATTGCCAAGGTTGGCGAAGGCGTGCAATTGCTGCCAATGACCGCCGGTGAAATCGTCCTGGCTTACAAcc                 <  6:216643/149‑1 (MQ=255)
                                                                                                                                       ttCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACATTGCCAAGGTTGGCGAAGGCGTGCAATTGCTGCCAATGACCGCCGGTGAAATCGTCCTGGCTTACAACCTGCCGGGTAACCCt   <  7:203772/149‑1 (MQ=255)
                                                                                                                                         cAGGAAGTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACATTGCCAAGGTTGGCGAAGGCGTGCAATTGCTGCCAATGACCGCCGGTGAAATCGTCCTGGCTTACAACCTGCCGGGTAACCCTa  >  4:303972/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                              
CTGCAGAAAACGTTCGGTTGACCGGTTCGGGCGCCAGCTTCCCTGCACCGATCTACCTGACCTGGTTCAAGGACTTCAGCAAGAACACGGCTGGCGTCACCGTCGATTACCAGTCCAAGGGCAGCGGCGCGGGCGTTCAGGACTTCCTGAACAAGACCGTCGACTTCGCCGCCAGCGACTCGGCCATGAGCGAAGCAGACATTGCCAAGGTTGGCGAAGGCGTGCAATTGCTGCCAATGACCGCCGGTGAAATCGTCCTGGCTTACAACCTGCCGGGTAACCCTA  >  NC_002947/3043427‑3043711

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