Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,161,380 T→C Y999H (TAT→CAT)  mdtC → multidrug efflux system, subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,161,3800TC100.0% 106.7 / NA 29Y999H (TAT→CAT) mdtCmultidrug efflux system, subunit C
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/17);  total (12/17)

CCCGATTATGATGACTACCCTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAATTGTGGATTGTGGCTTTCGGCTTCTTTA  >  NC_000913/2161241‑2161524
                                                                                                                                           |                                                                                                                                                
cccGATTATGATGACTACCCTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCg                                                                                                                                       >  1:394515/1‑151 (MQ=255)
  cGATTATGATGACTACCCTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCggt                                                                                                                                     <  2:233804/151‑1 (MQ=255)
        tgatgaCTACCCTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGGGgtgt                                                                                                                                >  1:143075/1‑150 (MQ=255)
             aCTACCCTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATctct                                                                                                                          >  1:513355/1‑151 (MQ=255)
                aCCCTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCtttt                                                                                                                       <  1:354666/151‑1 (MQ=255)
                            cTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGc                                                                                                           >  2:459769/1‑151 (MQ=255)
                                    tGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGt                                                                                                   >  1:114676/1‑151 (MQ=255)
                                          gCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGttttt                                                                                               <  1:203645/149‑1 (MQ=255)
                                          gCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTc                                                                                              <  2:120795/150‑1 (MQ=255)
                                           ccGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTcgc                                                                                            <  2:221861/151‑1 (MQ=255)
                                            cGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTcgc                                                                                            <  1:87075/150‑1 (MQ=255)
                                                ggTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTaa                                                                                        <  2:297994/150‑1 (MQ=255)
                                                   aTTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCt                                                                                    >  1:610510/1‑151 (MQ=255)
                                                   aTTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCt                                                                                    >  1:480442/1‑151 (MQ=255)
                                                      gTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTaa                                                                                  <  2:18044/150‑1 (MQ=255)
                                                      gTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTaa                                                                                  <  1:233912/150‑1 (MQ=255)
                                                                gACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTaa                                                                       <  2:256945/151‑1 (MQ=255)
                                                                gACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTa                                                                        >  2:485321/1‑150 (MQ=255)
                                                                        ggAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTaa                                                               <  1:50399/151‑1 (MQ=255)
                                                                        ggAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTa                                                                >  2:77548/1‑150 (MQ=255)
                                                                                  cAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGAc                                                          >  2:102930/1‑146 (MQ=255)
                                                                                     ccccTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTc                                                  <  2:383670/151‑1 (MQ=255)
                                                                                      cccTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTcc                                                 >  1:364390/1‑151 (MQ=255)
                                                                                                cACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCAc                                        <  1:74637/150‑1 (MQ=255)
                                                                                                              gACTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAattgtg                         <  1:202603/151‑1 (MQ=255)
                                                                                                               aCTGGTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAattgtgg                        >  2:79993/1‑151 (MQ=255)
                                                                                                                   gTAATGAGCCAGCTCCTTACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAattgtgga                       <  1:548835/148‑1 (MQ=255)
                                                                                                                                    tACGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAATTGTGGATTGTGGCTTTCGGcttctt    <  2:21968/150‑1 (MQ=255)
                                                                                                                                     aCGCTGCATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAATTGTGGATTGTGGCTTTCGGCTTCTTTa  <  1:443575/151‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                
CCCGATTATGATGACTACCCTGGCGGCGCTGTTTGGTGCGCTGCCGCTGGTATTGTCGGGCGGCGACGGCTCGGAGCTGCGGCAACCCCTGGGGATCACCATTGTCGGCGGACTGGTAATGAGCCAGCTCCTTACGCTGTATACCACGCCGGTGGTGTATCTCTTTTTCGACCGTCTGCGGCTGCGTTTTTCGCGTAAACCTAAACAAACGGTAACCGAGTAAATGACAGATCTTCCCGACAGCACCCGTTGGCAATTGTGGATTGTGGCTTTCGGCTTCTTTA  >  NC_000913/2161241‑2161524

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