Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 469,429 A→G K187E (AAG→GAG)  mdlA → putative multidrug ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913469,4290AG100.0% 154.4 / NA 43K187E (AAG→GAG) mdlAputative multidrug ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/22);  total (21/22)

ATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGG  >  NC_000913/469300‑469562
                                                                                                                                 |                                                                                                                                     
aTGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGcac                                                                                                                              >  1:1098616/1‑139 (MQ=255)
   ggCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAgg                                                                                                                           >  2:636228/1‑139 (MQ=255)
        cgcTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGgcggcg                                                                                                                      >  1:935337/1‑139 (MQ=255)
        cgcTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGgcggcg                                                                                                                      <  1:963462/139‑1 (MQ=255)
             tgtTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                                 >  1:247087/1‑139 (MQ=255)
             tgtTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                                 >  2:318229/1‑139 (MQ=255)
                    tATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCtt                                                                                                          <  1:103188/139‑1 (MQ=255)
                       gatgTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAAt                                                                                                       >  1:234354/1‑139 (MQ=255)
                       gatgTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAAt                                                                                                       <  2:1098616/139‑1 (MQ=255)
                           tCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAcc                                                                                                   <  2:234354/139‑1 (MQ=255)
                           tCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAcc                                                                                                   <  1:636228/139‑1 (MQ=255)
                            cTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCg                                                                                                  <  2:935337/139‑1 (MQ=255)
                              aggCAAATTAGCTGGCAGTTGACCTTACTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTGTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCa                                                                                                >  2:680431/3‑139 (MQ=255)
                                              aGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                                <  2:1085802/139‑1 (MQ=255)
                                                     cTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATc                                                                         <  1:235689/139‑1 (MQ=255)
                                                               cTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAcc                                                                                                   >  2:592923/1‑103 (MQ=255)
                                                               cTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAcc                                                                                                   <  1:592923/103‑1 (MQ=255)
                                                                   tGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                                >  1:986305/1‑118 (MQ=255)
                                                                   tGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                                <  2:986305/118‑1 (MQ=255)
                                                                          gCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGt                                                    >  1:330591/1‑139 (MQ=255)
                                                                            cAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                                 <  1:224802/76‑1 (MQ=255)
                                                                            cAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                                 >  2:224802/1‑76 (MQ=255)
                                                                                         cATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGcac                                                                                                                              >  1:548486/1‑50 (MQ=255)
                                                                                         cATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGcac                                                                                                                              <  2:548486/50‑1 (MQ=255)
                                                                                                 aGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAc                                                                                                    <  2:935285/68‑1 (MQ=255)
                                                                                                 aGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAc                                                                                                    >  1:935285/1‑68 (MQ=255)
                                                                                                   cgcAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAgg                                                                                          >  1:1130362/1‑76 (MQ=255)
                                                                                                   cgcAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAgg                                                                                          <  2:1130362/76‑1 (MQ=255)
                                                                                                   cgcAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGt                                                                            <  2:276285/90‑1 (MQ=255)
                                                                                                   cgcAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGt                                                                            >  1:276285/1‑90 (MQ=255)
                                                                                                      aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGcc                                                                                    >  1:407057/1‑79 (MQ=255)
                                                                                                      aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGcc                                                                                    <  2:407057/79‑1 (MQ=255)
                                                                                                      aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTg                        >  1:525123/1‑139 (MQ=255)
                                                                                                      aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTg                        >  2:761760/1‑139 (MQ=255)
                                                                                                          gCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCcgc                    <  2:525123/139‑1 (MQ=255)
                                                                                                                cTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCcgcg                   <  2:858535/134‑1 (MQ=255)
                                                                                                                cTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCcgcg                   >  1:858535/1‑134 (MQ=255)
                                                                                                                          aCGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCATGAAAGCTTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATAcc    <  2:718111/139‑1 (MQ=255)
                                                                                                                          aCGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCt                                                                                   <  2:799328/60‑1 (MQ=255)
                                                                                                                          aCGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCt                                                                                   >  1:799328/1‑60 (MQ=255)
                                                                                                                          aCGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCtt                                                        >  1:384015/1‑87 (MQ=255)
                                                                                                                          aCGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCtt                                                        <  2:384015/87‑1 (MQ=255)
                                                                                                                            gCTTTGAGCTGGCACAGGCGGCGTTTTCGAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCgg  <  1:680431/139‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                     
ATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGG  >  NC_000913/469300‑469562

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