Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,407,046 T→G 9.4% N124T (AAC→ACC)  PP_4751 ← amino acid ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,407,0460TG 9.4% 154.0 / 8.2 53N124T (AAC→ACC) PP_4751amino acid ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (26/21);  new base G (2/3);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

ATACGTTGCTGCTGGCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGTCGGCGG  >  NC_002947/5406902‑5407194
                                                                                                                                                |                                                                                                                                                    
atACGTTGCTGCTGGCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTgcg                                                                                                                                                  >  7:63708/1‑149 (MQ=255)
        ctgctgGCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCAcc                                                                                                                                           <  5:4435/148‑1 (MQ=255)
         tgctgGCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGt                                                                                                                                         >  5:297771/1‑149 (MQ=255)
           ctgGCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGGTGCTGCGCAACGTCa                                                                                                                                       <  1:125182/149‑1 (MQ=255)
              gCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGt                                                                                                                                    <  8:63708/149‑1 (MQ=255)
              gCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGGTGCTGCGCAACGTCAgc                                                                                                                                     <  2:158206/148‑2 (MQ=255)
                   caGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCgg                                                                                                                                >  6:81344/1‑148 (MQ=255)
                      agcagcGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAAc                                                                                                                            <  1:108575/149‑1 (MQ=255)
                        cagcGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAg                                                                                                                          <  5:81344/149‑1 (MQ=255)
                        cagcGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAg                                                                                                                          <  5:126132/149‑1 (MQ=255)
                         agcGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAgg                                                                                                                         <  8:36624/149‑1 (MQ=255)
                           cgcgcCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGtt                                                                                                                       <  6:297771/149‑1 (MQ=255)
                                         cTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAac                                                                                                          <  1:225120/148‑1 (MQ=255)
                                                   tCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCaa                                                                                                >  3:184593/1‑148 (MQ=255)
                                                      gccacccCAACGCGATCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTACTGCCCAAGCTGCGGATTTGCGCCAGGCTGACGATGCGCAGCACCGTCAGGAGCGGGAACAGCTTGAACGACTAGAACACCAAACAAAAgcg                                                                                            >  8:334753/1‑149 (MQ=255)
                                                       ccagccCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATgcgc                                                                                           >  5:214328/1‑149 (MQ=255)
                                                       ccagccCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGGTGCTGCGCAACGTCAgcc                                                                                                                                    <  3:236503/108‑3 (MQ=255)
                                                       ccagccCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGGTGCTGCGCAACGTCAgcc                                                                                                                                    >  4:236503/1‑106 (MQ=255)
                                                           cccAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGgc                                                                                       >  6:230063/1‑149 (MQ=255)
                                                              aaCGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACAGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAAAGCGCCGgcgc                                                                                     >  4:165238/1‑148 (MQ=255)
                                                               aCGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCaa                                                                                   <  6:214328/149‑1 (MQ=255)
                                                                            gTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCgcagc                                                                                                                                              >  6:190079/1‑77 (MQ=255)
                                                                            gTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCgcagc                                                                                                                                              <  5:190079/77‑1 (MQ=255)
                                                                            gTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCt                                                                      >  8:89111/1‑149 (MQ=255)
                                                                             ttGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTg                                                                     >  4:215783/1‑149 (MQ=255)
                                                                                 aaGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACg                                                                 <  7:343990/149‑1 (MQ=255)
                                                                                    gAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAAc                                                                                                                   <  8:81463/96‑1 (MQ=255)
                                                                                    gAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAAc                                                                                                                   >  7:81463/1‑96 (MQ=255)
                                                                                      gCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATa                                                                                                    >  1:238515/1‑109 (MQ=255)
                                                                                      gCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATa                                                                                                    <  2:238515/109‑1 (MQ=255)
                                                                                       cGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCaa                                                           >  6:275042/1‑149 (MQ=255)
                                                                                        gggCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAAc                                                          >  4:251155/1‑149 (MQ=255)
                                                                                           cGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGcc                                                       <  5:275042/149‑1 (MQ=255)
                                                                                             tCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCtt                                                     <  7:334753/149‑1 (MQ=255)
                                                                                               ggCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTg                                                    <  4:184593/148‑1 (MQ=255)
                                                                                                                 gCCCAAGGTGCGGATTTGCGCCAGGCTGACCTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGc                                 >  1:30407/1‑149 (MQ=255)
                                                                                                                  cccAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGc                                 >  3:288789/1‑148 (MQ=255)
                                                                                                                    cAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGccat                              >  3:338773/1‑149 (MQ=255)
                                                                                                                        gTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGccatccat                          >  8:18307/1‑149 (MQ=255)
                                                                                                                         tGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATg                         >  3:91860/1‑149 (MQ=255)
                                                                                                                         tGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATg                         <  5:230063/149‑1 (MQ=255)
                                                                                                                           cGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGc                        <  3:215783/148‑1 (MQ=255)
                                                                                                                           cGGATTTGCGACTGGCTGACGTTGGGCAGCACCGTCAGGTGCGTGAACTGGTTGAACGACTGGAACTCAATACAAATGCGCCTGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGc                        <  3:165238/148‑1 (MQ=255)
                                                                                                                                   cgcCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGcc               >  5:179764/1‑149 (MQ=255)
                                                                                                                                   cgcCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGcc               >  8:97517/1‑149 (MQ=255)
                                                                                                                                   cccgcaccTGACGGTGCTGCGCAACGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGcc               >  4:167430/8‑149 (MQ=255)
                                                                                                                                    gcCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGccc              >  2:60120/1‑149 (MQ=255)
                                                                                                                                    gcCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGccc              <  4:338773/149‑1 (MQ=255)
                                                                                                                                     ccAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGcccc             >  3:271205/1‑149 (MQ=255)
                                                                                                                                      cAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCt            >  6:133932/1‑149 (MQ=255)
                                                                                                                                          cTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGt        <  1:60120/149‑1 (MQ=255)
                                                                                                                                            gACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGTcg      >  6:163381/1‑149 (MQ=255)
                                                                                                                                               gTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGTcggcg   >  7:10192/1‑149 (MQ=255)
                                                                                                                                                ttGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGTcggcgg  <  7:97517/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ATACGTTGCTGCTGGCCACCAGAGCAGCGCGCCGGATAATGCTGGGCCTTGTCAGCCAGCCCAACGCGTTCGAGCAGTTGCAAGGAGCGGGCGTCGGCTTCTTTGGCGCTGCGGCCCAAGGTGCGGATTTGCGCCAGGCTGACGTTGCGCAGCACCGTCAGGTGCGGGAACAGGTTGAACGACTGGAACACCATACCAATGCGCCGGCGCAACTTCAGCAGTTCTGTACGGGGCAACGCCTTGCCAGCTTCGATGTGCACGCCATCCATGCTGACCCGCCCCTTGGTCGGCGG  >  NC_002947/5406902‑5407194

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