Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,157,974 C→T 7.7% G904G (GGC→GGT)  mdtB → multidrug efflux system, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,157,9740CT 7.7% 160.6 / 3.2 53G904G (GGC→GGT) mdtBmultidrug efflux system, subunit B
Reads supporting (aligned to +/- strand):  ref base C (29/19);  new base T (2/2);  total (31/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.23e-01

CACCCTCGCCTTCCAGTCGGCGCTGGGCAGCACTGTCTGGCTGATTGTCGCGGCGGTGGTGGCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTGAAGAAGAACGCCATCATGATGATCGACTTCGCGCTGGCTGCTGAGCGCGAGCAAGGCATGT  >  NC_000913/2157827‑2158116
                                                                                                                                                   |                                                                                                                                              
caCCCTCGCCTTCCAGTCGGCGCTGGGCAGCACTGTCTGGCTGATTGTCGCGGCGGTGGTGGCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGATGgcg                                                                                                                                               >  3:223918/1‑149 (MQ=255)
   ccTCGCCTTCCAGTCGGCGCTGGGCAGCACTGTCTGGCTGATTGTCGCGGCGGTGGTGGCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCAc                                                                                                                                            >  2:46494/1‑149 (MQ=255)
   ccTCGCCTTCCAGTCGGCGCGGGGCAGCACTGTCTGGCTGATTGTCGCGGCGGTGGGGGCGATGTATATCGTGCTCGGCAGTCTGTACGAGAGCTTTATTCACCCGAGCAGCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCAc                                                                                                                                            >  8:156333/1‑149 (MQ=255)
       gCCTTCCAGTCGGCGCTGGGCAGCACTGTCTGGCTGATTGTCGCGGCGGTGGTGGCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCActgct                                                                                                                                        >  2:161397/1‑149 (MQ=255)
       gCCTTCCAGTCGGCGCTGGGCAGCACTGTCTGGCTGATTGTCGCGGCGGTGGTGGCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCActgct                                                                                                                                        >  1:338692/1‑149 (MQ=255)
       gCCTTCCAGTCGGCGCTGGGCAGCACTGTCTGGCTGATTGTCGCGGCGGTGGTGGCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCActg                                                                                                                                          <  6:104479/147‑1 (MQ=255)
         cTTCCAGTCGGCGCTGGGCAGCACTGTCTGGCTGATTGTCGCGGCGGTGGTGGCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTgg                                                                                                                                      >  4:283442/1‑149 (MQ=255)
          ttCCAGTCGGCGCTGGGCAGCACTGTCTGGCTGATTGTCGCGGCGGTGGAGGCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGc                                                                                                                                     >  5:193822/1‑149 (MQ=255)
               gTCGGCGCTGGGCAGCACTGTCTGGCTGATTGTCGCGGCGGTGGTGGCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGTGCACTGCTGGCGTTg                                                                                                                                 >  8:240470/1‑148 (MQ=255)
               gTCGGCGCTGGGCAGCACTGTCTGGCTGATTGTCGCGGCGGTGGTGGCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGTGCACTGCTGGCGTTg                                                                                                                                 <  7:234609/148‑1 (MQ=255)
               gTCGGCGCTGGGCAGCACTGTCTGGCTGATTGTCGCGGCGGTGGTGGCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGTGCACTGCTGGCGTTg                                                                                                                                 <  7:240470/148‑1 (MQ=255)
               gTCGGCGCTGGGCAGCACTGTCTGGCTGATTGTCGCGGCGGTGGTGGCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGTGCACTGCTGGCGTTg                                                                                                                                 >  8:234609/1‑148 (MQ=255)
                         ggCAGCACTGTCTGGCTGATTGTCGCGGCGGTGGTGGCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTg                                                                                                                       <  4:344598/148‑1 (MQ=255)
                             gcaCTGTCTGGCTGATTGTCGCGGCGGTGGTGGCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGTGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGc                                                                                                                  <  3:236738/149‑1 (MQ=255)
                                cTGTCTGGCTGATTGTCGCGGCGGTGGTGGCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCTCCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGaa                                                                                                               <  4:41731/149‑1 (MQ=255)
                                cTGTCTGGCTGATTGTCGCGGCGGTGGTGGCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGaa                                                                                                               <  6:193822/149‑1 (MQ=255)
                                   tCTGGCTGATTGTCGCGGCGGTGGTGGCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTg                                                                                                            <  3:283442/149‑1 (MQ=255)
                                       gCTGATTGTCGCGGCGGTGGTGGCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGAtg                                                                                                        >  3:148055/1‑149 (MQ=255)
                                         tGATTGTCGCGGCGGTGGTGGCGATGTATATCGTGCTGGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGAtgt                                                                                                       >  7:162256/1‑148 (MQ=255)
                                                      ggtggtGGCGATGTATATCGTGCTGGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCg                                                                                         <  8:162256/149‑1 (MQ=255)
                                                          gtggCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGAGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCat                                                                                     >  4:344605/1‑149 (MQ=255)
                                                                aTGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTAtttt                                                                               <  1:337583/149‑1 (MQ=255)
                                                                 tGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTg                                                                              <  4:223918/149‑1 (MQ=255)
                                                                   tataTCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCt                                                                            >  7:328533/1‑149 (MQ=255)
                                                                   tataTCGTGCCCGTCATTCTGGACGTGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCAGCGCAGGGGTTGGCGCACTGCGGGCGGTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGc                                                                             <  8:104834/148‑1 (MQ=255)
                                                                        cGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCg                                                                       >  2:79651/1‑149 (MQ=255)
                                                                               ggCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATc                                                                  <  2:338692/147‑1 (MQ=255)
                                                                               ggCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGt                                                                >  1:331292/1‑149 (MQ=255)
                                                                               ggCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGAGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCg                                                                 <  1:9038/148‑1 (MQ=255)
                                                                                 cATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTga                                                              <  2:400931/149‑1 (MQ=255)
                                                                                        tACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTgaagaagaa                                                       <  1:159095/149‑1 (MQ=255)
                                                                                           gagaGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGAATCGTGAAGAAGAACGc                                                    >  1:272287/1‑149 (MQ=255)
                                                                                            agagCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTGAAGAAGAACGcc                                                   <  1:79651/149‑1 (MQ=255)
                                                                                                    aTTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTGAAGAAGCACGCCATCatgat                                           >  1:135707/1‑149 (MQ=255)
                                                                                                        aCCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTGAAGAAGAACGCCATCATGATGATc                                       >  2:166332/1‑149 (MQ=255)
                                                                                                        aCCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGAATATCGGCATTATTTTGCTGATCGGTATCGTGAAGAAGAACGCCATCATGATGATc                                       >  4:249749/1‑149 (MQ=255)
                                                                                                          ccGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTGAAGAAGAACGCCATCATGATGATCGa                                     >  4:22450/1‑149 (MQ=255)
                                                                                                          ccGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTGAAGAAGAACGCCATCATGATGATCGa                                     >  8:285714/1‑149 (MQ=255)
                                                                                                          ccGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTGAAGAAGAACGCCATCATGATGATCGa                                     >  1:333457/1‑149 (MQ=255)
                                                                                                            gATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATTGTGAAGAAGAACGCCATCATGATGATCGACt                                   >  2:78193/1‑149 (MQ=255)
                                                                                                                 ccATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTGAAGAAGAACGCCATCATGATGATCGACTTcgcg                              >  5:14578/1‑149 (MQ=255)
                                                                                                                  cATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTgaa                                                             >  1:21423/1‑117 (MQ=255)
                                                                                                                  cATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTgaa                                                             <  2:21423/117‑1 (MQ=255)
                                                                                                                  cATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTGAAGAAGAACGCCATCATGATGATCGACTTCgcgc                             <  4:148055/149‑1 (MQ=255)
                                                                                                                         tcGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTGAAGAAGAACGCCATCATGATGATCGACTTCGCGCTGgctgc                      >  8:358292/1‑149 (MQ=255)
                                                                                                                         tcGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTGAAGAAGAACGCCATCATGATGATCGACTTCGCGCTGgctgc                      >  4:243881/1‑149 (MQ=255)
                                                                                                                           gACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTGAAGAAGAACGCCATCATGATGATCGACTTCGCGCTGgctgct                     >  8:280533/1‑148 (MQ=255)
                                                                                                                                tACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTgaaga                                                           <  3:47339/105‑1 (MQ=255)
                                                                                                                                tACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTgaaga                                                           >  4:47339/1‑105 (MQ=255)
                                                                                                                                   ccACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTGAAGAAGAACGCCATCATGATGATCGACTTCGCGCTGGCTGCTGAGCGCGAg            >  8:108368/1‑149 (MQ=255)
                                                                                                                                   ccACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTGAAGAAGAACGCCATCATGATGATCGACTTCGCGCTGGCTGCTGAGCGCGAg            >  1:234605/1‑149 (MQ=255)
                                                                                                                                       cGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTGAAGAAGAACGCCATCATGATGATCGACTTCGCGCTGGCTGCTGAGCGCGAGCAAg        >  7:338136/1‑149 (MQ=255)
                                                                                                                                         cAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGAATCGTGAAGAAGAACGCCATCATGATGATCGACTTCGCGCTGGCGGCTGAGCGCGAGCAAGGc      <  2:272287/149‑1 (MQ=255)
                                                                                                                                             ggTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTGAAGAAGAACGCCATCATGATGATCGACTTCGCGCTGGCTGCTGAGCGCGAGCAAGGCATGt  <  1:161397/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CACCCTCGCCTTCCAGTCGGCGCTGGGCAGCACTGTCTGGCTGATTGTCGCGGCGGTGGTGGCGATGTATATCGTGCTCGGCATTCTGTACGAGAGCTTTATTCACCCGATCACCATTCTCTCGACGCTACCCACCGCAGGGGTTGGCGCACTGCTGGCGTTGCTGATTGCTGGTAGCGAACTGGATGTGATTGCGATTATCGGCATTATTTTGCTGATCGGTATCGTGAAGAAGAACGCCATCATGATGATCGACTTCGCGCTGGCTGCTGAGCGCGAGCAAGGCATGT  >  NC_000913/2157827‑2158116

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