Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,157,507 C→G 8.0% L749V (CTG→GTG)  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,5070CG 8.0% 167.9 / 6.3 50L749V (CTG→GTG) mdtBmultidrug efflux system, subunit B
Reads supporting (aligned to +/- strand):  ref base C (27/19);  new base G (2/2);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

CAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAA  >  NC_000913/2157360‑2157654
                                                                                                                                                   |                                                                                                                                                   
cAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCt                                                                                                                                                    <  3:153222/149‑1 (MQ=255)
       cacaGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGATGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCAGTCTGGGGATCAGAATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTc                                                                                                                                              >  1:288126/1‑148 (MQ=255)
                          ccgCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGcc                                                                                                                           >  6:19903/3‑148 (MQ=255)
                          cAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATCCTCAGGcc                                                                                                                           >  6:210024/1‑148 (MQ=255)
                                   gCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATc                                                                                                                 >  4:330816/1‑149 (MQ=255)
                                   gCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATc                                                                                                                 >  4:242204/1‑149 (MQ=255)
                                    cAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCTGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCAAACCAGTATcg                                                                                                                >  4:264349/1‑149 (MQ=255)
                                     agtaCAAAGGGCTTGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGTGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGGCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATcgc                                                                                                               <  2:288126/146‑1 (MQ=255)
                                          aaaGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCgtggt                                                                                                          >  1:184210/1‑149 (MQ=255)
                                          aaaGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCgtggt                                                                                                          >  6:365370/1‑149 (MQ=255)
                                          aaaGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCAAACCAGTATCGCgtggt                                                                                                          <  3:264349/149‑1 (MQ=255)
                                                   gtggCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGca                                                                                                 >  7:110976/1‑149 (MQ=255)
                                                      gCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCacaa                                                                                              <  5:365370/149‑1 (MQ=255)
                                                            gTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTAt                                                                                                                                        <  2:52887/101‑1 (MQ=255)
                                                            gTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTAt                                                                                                                                        >  1:52887/1‑101 (MQ=255)
                                                            gTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCg                                                                                         >  3:242807/1‑148 (MQ=255)
                                                               aaTGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGaaaa                                                                                     >  8:306815/1‑149 (MQ=255)
                                                               aaTGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGaaaa                                                                                     <  5:19903/149‑1 (MQ=255)
                                                                aTGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAAt                                                                                    <  8:238270/149‑1 (MQ=255)
                                                                 tGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGa                                                                                                    <  6:287928/132‑1 (MQ=255)
                                                                 tGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGa                                                                                                    >  5:287928/1‑132 (MQ=255)
                                                                       tCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAg                                                                             <  2:184210/149‑1 (MQ=255)
                                                                                   cAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGgcgc                                                                 >  3:276071/1‑149 (MQ=255)
                                                                                       cGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGa                                                             <  7:306815/149‑1 (MQ=255)
                                                                                                aTCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGGTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTc                                                     <  8:26081/148‑1 (MQ=255)
                                                                                                aTCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGGTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTc                                                     >  7:26081/1‑148 (MQ=255)
                                                                                                aTCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGGTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTc                                                     <  8:20412/148‑1 (MQ=255)
                                                                                                aTCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGGTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTc                                                     >  7:20412/1‑148 (MQ=255)
                                                                                                            gATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACcagcag                                        <  3:330816/149‑1 (MQ=255)
                                                                                                                tCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCACGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGAc                                    >  8:220282/1‑149 (MQ=255)
                                                                                                                  gATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGAcgg                                  <  8:258473/149‑1 (MQ=255)
                                                                                                                        gCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCgt                            >  5:143652/1‑149 (MQ=255)
                                                                                                                         cccTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAACACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCgtg                           >  8:372077/1‑149 (MQ=255)
                                                                                                                          ccTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCgtgg                          >  3:173653/1‑149 (MQ=255)
                                                                                                                                cAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCa                                                                              <  2:96675/91‑1 (MQ=255)
                                                                                                                                cAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCa                                                                              >  1:96675/1‑91 (MQ=255)
                                                                                                                                 aaCGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGc                    >  5:252657/1‑148 (MQ=255)
                                                                                                                                 aaCGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGc                    >  1:11631/1‑148 (MQ=255)
                                                                                                                                  aCGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCt                   >  1:12438/1‑148 (MQ=255)
                                                                                                                                  aCGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCt                   >  2:6375/1‑148 (MQ=255)
                                                                                                                                   cgcgTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCGGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTaa                 <  5:210024/149‑1 (MQ=255)
                                                                                                                                    gcgTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAg                >  5:64648/1‑149 (MQ=255)
                                                                                                                                    gcgTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAg                >  3:108937/1‑149 (MQ=255)
                                                                                                                                      gTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCgtg                           <  6:195222/136‑1 (MQ=255)
                                                                                                                                      gTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCgtg                           >  5:195222/1‑136 (MQ=255)
                                                                                                                                      gTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCt              <  1:6375/149‑1 (MQ=255)
                                                                                                                                           gTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAAtt         <  4:242807/149‑1 (MQ=255)
                                                                                                                                           gTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAAtt         <  3:242204/149‑1 (MQ=255)
                                                                                                                                            tCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTg        >  8:127579/1‑149 (MQ=255)
                                                                                                                                                  gCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCaaaa  <  6:247994/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAA  >  NC_000913/2157360‑2157654

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