Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,673,345 (T)7→8 coding (157/366 nt) mdtJ ← multidrug efflux system transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,673,3381.T96.7% 188.3 / ‑0.9 60coding (164/366 nt)mdtJmultidrug efflux system transporter
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base T (30/28);  total (31/29)
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.29e-01

GGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCC  >  NC_000913/1673219‑1673472
                                                                                                                        |                                                                                                                                       
gggTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAag                                                                                                                      >  1:1220192/1‑139 (MQ=255)
   tggtTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGAtttttttt                                                                                                                                 >  1:766115/1‑124 (MQ=255)
   tggtTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGAtttttttt                                                                                                                                 <  2:766115/125‑2 (MQ=255)
    ggtTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAagagag                                                                                                                  <  2:1220192/139‑1 (MQ=255)
      ttAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGa                                                                                                                         >  1:453104/1‑130 (MQ=255)
      ttAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGa                                                                                                                         <  2:453104/130‑1 (MQ=255)
            cGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGaa                                                                                                                >  2:833164/1‑133 (MQ=255)
            cGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGaa                                                                                                                <  1:833164/133‑1 (MQ=255)
             ggCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGaaa                                                                                                                       >  2:373912/1‑125 (MQ=255)
             ggCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGaaa                                                                                                                       <  1:373912/125‑1 (MQ=255)
                                aaaCTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCAcc                                                                                      <  2:230778/139‑1 (MQ=255)
                                         tcgAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAag                                                                                                                      <  2:449629/98‑1 (MQ=255)
                                         tcgAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAag                                                                                                                      >  1:449629/1‑98 (MQ=255)
                                                    aaCGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAgccgcc                                                                  <  1:1477414/139‑1 (MQ=255)
                                                      cGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCAt                                                                                 >  1:494453/1‑122 (MQ=255)
                                                      cGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCAt                                                                                 <  2:494453/122‑1 (MQ=255)
                                                                ggTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCAcc                                                                                      <  1:1284792/107‑1 (MQ=255)
                                                                ggTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCAcc                                                                                      >  2:1284792/1‑107 (MQ=255)
                                                                     taaataaaATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTAAACGGCGAAAGAGAGAAatata                                                                                                           >  2:453393/1‑81 (MQ=255)
                                                                     taaataaaATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTAAACGGCGAAAGAGAGAAatata                                                                                                           <  1:453393/81‑1 (MQ=255)
                                                                            aataCCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCAc                                                                                       >  2:1283223/1‑94 (MQ=255)
                                                                            aataCCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCAc                                                                                       <  1:1283223/94‑1 (MQ=255)
                                                                              taCCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCAtt                                        >  2:1407015/1‑139 (MQ=255)
                                                                                   aTACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCAt                                         >  2:1511046/1‑133 (MQ=255)
                                                                                   aTACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCAt                                         <  1:1511046/133‑1 (MQ=255)
                                                                                   aTACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCAtt                                   >  1:1133261/1‑139 (MQ=255)
                                                                                   aTACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCAtt                                   >  1:1242587/1‑139 (MQ=255)
                                                                                         tCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGaa                                                                                                                <  2:667259/56‑1 (MQ=39)
                                                                                         tCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGaa                                                                                                                >  1:667259/1‑56 (MQ=39)
                                                                                          cccACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACt                                                      <  1:281183/113‑1 (MQ=255)
                                                                                          cccACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACt                                                      >  2:281183/1‑113 (MQ=255)
                                                                                          cccACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCg                            >  2:1097806/1‑139 (MQ=255)
                                                                                               aGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGa                                 >  2:1128692/1‑129 (MQ=255)
                                                                                               aGCGCATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGa                                 <  1:1128692/129‑1 (MQ=255)
                                                                                                gcgcATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCAcc                                                                                      <  1:131912/75‑1 (MQ=255)
                                                                                                gcgcATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCAcc                                                                                      >  2:131912/1‑75 (MQ=255)
                                                                                                 cgcATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAagaga                                                                                                                   >  2:736566/1‑45 (MQ=38)
                                                                                                 cgcATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAagaga                                                                                                                   <  1:736566/45‑1 (MQ=38)
                                                                                                 cgcATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGaa                                                                                                                >  2:577190/1‑48 (MQ=38)
                                                                                                 cgcATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGaa                                                                                                                <  1:577190/48‑1 (MQ=38)
                                                                                                 cgcATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAgc                                                                      <  2:54197/90‑1 (MQ=255)
                                                                                                 cgcATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAgc                                                                      >  1:54197/1‑90 (MQ=255)
                                                                                                 cgcATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCAt                                         >  1:177105/1‑119 (MQ=255)
                                                                                                 cgcATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCAt                                         <  2:177105/119‑1 (MQ=255)
                                                                                                 cgcATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGt                     >  1:775441/1‑139 (MQ=255)
                                                                                                 cgcATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTa                    >  2:830841/1‑139 (MQ=255)
                                                                                                   cATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAatata                                                                                                           >  1:650573/1‑51 (MQ=38)
                                                                                                   cATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAatata                                                                                                           <  2:650573/51‑1 (MQ=38)
                                                                                                   cATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGaaa                                                                                              <  2:427959/64‑1 (MQ=255)
                                                                                                   cATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGaaa                                                                                              >  1:427959/1‑64 (MQ=255)
                                                                                                   cATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAatc                                                                                            <  2:250218/66‑1 (MQ=255)
                                                                                                   cATAAGCTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAatc                                                                                            >  1:250218/1‑66 (MQ=255)
                                                                                                         cTACGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTg             <  2:1133261/139‑1 (MQ=255)
                                                                                                           aCGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTa           >  2:382757/1‑139 (MQ=255)
                                                                                                            cGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAg                                                                       <  1:574825/78‑1 (MQ=255)
                                                                                                            cGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAg                                                                       >  2:574825/1‑78 (MQ=255)
                                                                                                            cGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTg             >  2:1013182/1‑136 (MQ=255)
                                                                                                            cGCCTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTg             <  1:1013182/136‑1 (MQ=255)
                                                                                                              ccTAAGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCa        <  1:1097806/139‑1 (MQ=255)
                                                                                                                  aGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCAt                                         >  2:179595/1‑102 (MQ=255)
                                                                                                                  aGGCGATTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCAt                                         <  1:179595/102‑1 (MQ=255)
                                                                                                                   ggCGATTTTTTTTAACGGCGAAAGAGAGAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGcc  <  1:830841/139‑1 (MQ=255)
                                                                                                                        |                                                                                                                                       
GGGTGGTTAACCCGGCAATTTTCATCAGCGATAAACTTTCGTCGAATAACAAAACGCTAAACAAGGTAATAAATAAAATACCGATACCTTCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCC  >  NC_000913/1673219‑1673472

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