Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,915,796 C→G 9.4% Y255* (TAC→TAG)  PP_3455 → multidrug RND transporter membrane fusion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,915,7960CG 9.4% 158.7 / 7.1 53Y255* (TAC→TAG) PP_3455multidrug RND transporter membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base C (20/28);  new base G (3/2);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

CCCGATCCATGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCC‑CTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCGGGTACGCCTGGAGCAAGCGGTG  >  NC_002947/3915661‑3915937
                                                                                                                                        |                                                                                                                                             
cccGATCCATGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCC‑CTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGc                                                                                                                                   >  5:42756/1‑148 (MQ=255)
 ccGATCCATGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCCGCTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGc                                                                                                                                   >  1:66902/1‑148 (MQ=255)
 ccGATCCATGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCC‑CTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGGGTAGGCGCTGCCgt                                                                                                                                     >  8:304080/1‑144 (MQ=255)
 ccGATCCATGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCC‑CTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGGGTAGGCGCTGCCgt                                                                                                                                     <  7:304080/145‑2 (MQ=255)
    aTCCATGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCC‑CTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAgct                                                                                                                              >  7:176274/1‑149 (MQ=255)
        aTGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCC‑CTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGctgctg                                                                                                                          <  3:54825/149‑1 (MQ=255)
        aTGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCC‑CTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGctgctg                                                                                                                          <  7:154450/149‑1 (MQ=255)
        aTGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCC‑CTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGctgctg                                                                                                                          <  2:66902/149‑1 (MQ=255)
             gATGTCACCCAGTCTACCCGCGAGCTCAACGCC‑CTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGACCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTtctc                                                                                                                     >  1:14167/1‑149 (MQ=255)
             gATGTCACCCAGTCTACCCGCGAGCTCAACGCC‑CTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGACCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTtctc                                                                                                                     <  2:14167/149‑1 (MQ=255)
                       aGTCTACCCGCGAGCTCAACGCC‑CTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGc                                                                                                           <  7:226952/149‑1 (MQ=255)
                               cgcgAGCTCAACGCC‑CTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGAccc                                                                                                   <  4:151985/149‑1 (MQ=255)
                               cgcgAGCTCAACGCC‑CTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGAccc                                                                                                   <  4:147074/149‑1 (MQ=255)
                                 cgAGCTCAACGCC‑CTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCa                                                                                                 <  6:42756/149‑1 (MQ=255)
                                        aaCGCC‑CTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCaccaa                                                                                          >  4:281366/1‑149 (MQ=255)
                                         aCGCC‑CTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCaccaac                                                                                         <  3:80685/149‑1 (MQ=255)
                                          cGCC‑CTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGAc                                                                                                     <  8:67634/136‑1 (MQ=255)
                                          cGCC‑CTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGAc                                                                                                     >  7:67634/1‑136 (MQ=255)
                                                        gccCCGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTgcg                                                                           >  1:95399/1‑149 (MQ=255)
                                                          ccTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGctgctg                                                                                                                          <  6:250534/100‑1 (MQ=255)
                                                          ccTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGctgctg                                                                                                                          >  5:250534/1‑100 (MQ=255)
                                                           cTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAg                                                                        >  6:103073/1‑149 (MQ=255)
                                                           cTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAg                                                                        >  2:14417/1‑149 (MQ=255)
                                                             gcgcgcCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCaccaac                                                                                         <  1:273721/130‑1 (MQ=255)
                                                             gcgcgcCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCaccaac                                                                                         >  2:273721/1‑130 (MQ=255)
                                                                       gAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGa                                                                                     <  8:253115/124‑1 (MQ=255)
                                                                       gAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGa                                                                                     >  7:253115/1‑124 (MQ=255)
                                                                       gAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTccccga                                                            >  8:91474/1‑147 (MQ=255)
                                                                        aaTTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAAc                                                           <  2:89317/149‑1 (MQ=255)
                                                                         aTTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAAcc                                                          <  2:314710/149‑1 (MQ=255)
                                                                             cagcagGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCgg                                                       >  4:90002/1‑148 (MQ=255)
                                                                             cagcagGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCgg                                                       <  3:90002/148‑1 (MQ=255)
                                                                                    tGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCGACCCGGACCTCGAc                                               <  7:91474/149‑1 (MQ=255)
                                                                                         gACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACctgct                                          >  2:210097/1‑149 (MQ=255)
                                                                                           cGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCtgctgc                                        <  5:103073/149‑1 (MQ=255)
                                                                                            ggCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCtgctgc                                        >  2:293730/1‑148 (MQ=255)
                                                                                            ggCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCtgctgc                                        >  8:150905/1‑148 (MQ=255)
                                                                                               cAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACctgct                                          >  3:108275/1‑143 (MQ=255)
                                                                                               cAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACctgct                                          <  4:108275/143‑1 (MQ=255)
                                                                                                     cgcgCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTa                              <  1:293730/149‑1 (MQ=255)
                                                                                                     cgcgCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTa                              <  8:174724/149‑1 (MQ=255)
                                                                                                      gcgcCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTAc                             <  6:121136/149‑1 (MQ=255)
                                                                                                           aCGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCg                        >  6:38093/1‑149 (MQ=255)
                                                                                                            cGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCgg                       <  3:135686/149‑1 (MQ=255)
                                                                                                             gCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATc                                                                                                              <  3:253964/61‑1 (MQ=255)
                                                                                                             gCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATc                                                                                                              >  4:253964/1‑61 (MQ=255)
                                                                                                                    cAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCGGGTACGCCt               <  5:38093/149‑1 (MQ=255)
                                                                                                                     aGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCGGGTACGCCTg              <  1:14417/149‑1 (MQ=255)
                                                                                                                       gATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCGGGTACGCCTGGa            <  1:210097/149‑1 (MQ=255)
                                                                                                                        tgcccGGCAGCGGGTAGGCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACct                                             <  2:223335/111‑1 (MQ=255)
                                                                                                                        tgcccGGCAGCGGGTAGGCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACct                                             >  1:223335/5‑115 (MQ=255)
                                                                                                                           ccGGCAGCGGGTAGGCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCGGGTACGCCTGGAGCaa        >  8:124632/2‑149 (MQ=255)
                                                                                                                                 gcgcCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCGGGTACGCCTGGAGCAAGCGGTg  <  7:36556/149‑1 (MQ=255)
                                                                                                                                        cccGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCGGGTACGCCt               >  5:72285/1‑129 (MQ=255)
                                                                                                                                        cccGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCGGGTACGCCt               <  6:72285/129‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                             
CCCGATCCATGCCGATGTCACCCAGTCTACCCGCGAGCTCAACGCC‑CTGCGCCGCGCCCTGCGCGCCGGCGAATTGCAGCAGGTGGGCGACGGCCAGGCACGCGCCACGCTGATCCAGGATGACGGCAGCGCCTACCCGCTGCCGGGCAAGCTGCTGTTCTCCGACATCAGCGTCGACCCCAGCACCAACCAGATCACCCTGCGCAGCGAGTTCCCCAACCCGGACCTCGACCTGCTGCCAGGCAGCTACGTGCGGGTACGCCTGGAGCAAGCGGTG  >  NC_002947/3915661‑3915937

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