Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,159,354 C→T R323R (CGC→CGT mdtC → multidrug efflux system, subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,159,3540CT97.9% 158.7 / ‑4.2 48R323R (CGC→CGTmdtCmultidrug efflux system, subunit C
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (26/21);  total (27/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 = 9.15e-01

AATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCT  >  NC_000913/2159259‑2159448
                                                                                               |                                                                                              
aaTATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTccc                                                                                            >  2:1247104/1‑100 (MQ=255)
    tgatCCAGACGGTGGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCccac                                                                                        <  1:130063/98‑1 (MQ=255)
      atCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCccacca                                                                                      >  2:80794/1‑100 (MQ=255)
      atCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCccacca                                                                                      >  2:1048616/1‑100 (MQ=255)
        ccAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCAtt                                                                                    >  2:320198/1‑100 (MQ=255)
              ggTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGcc                                                                              >  2:1405096/1‑100 (MQ=255)
                ttGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTc                                                                            <  2:711099/100‑1 (MQ=255)
                     aGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTgg                                                                       <  1:370922/100‑1 (MQ=255)
                      gCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGga                                                                      >  2:869913/1‑100 (MQ=255)
                         tcccagCAAAATTACCGGATTTGCAGGAAACCTTACCGGCAGCGACTGATCTCCACATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGgaag                                                                    <  2:152286/94‑1 (MQ=255)
                         tCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGgaaga                                                                   >  2:1384341/1‑100 (MQ=255)
                                aaaaTTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAg                                                            <  1:863639/100‑1 (MQ=255)
                                aaaaTTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAg                                                            <  1:727854/100‑1 (MQ=255)
                                 aaaTTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAg                                                            <  2:1146216/99‑1 (MQ=255)
                                  aaTTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCa                                                          >  2:316531/1‑100 (MQ=255)
                                     tACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAAc                                                       >  2:628295/1‑100 (MQ=255)
                                      aCCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACg                                                      <  1:1549113/100‑1 (MQ=255)
                                         ggAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTg                                                   <  2:1411753/100‑1 (MQ=255)
                                          gAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGa                                                  <  1:1416125/100‑1 (MQ=255)
                                           aGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGAt                                                 <  1:1346692/100‑1 (MQ=255)
                                           aGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCACGCCTCGCTGGAAGAAGTCGAGCAAACGCTg                                                   <  2:1066699/98‑1 (MQ=255)
                                            gTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGAtt                                                >  2:304892/1‑100 (MQ=255)
                                               gCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTAtc                                             >  2:1122080/1‑100 (MQ=255)
                                                cAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTAtct                                            >  2:444500/1‑100 (MQ=255)
                                                cAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTAtct                                            >  1:1520065/1‑100 (MQ=255)
                                                cAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTAtct                                            >  2:1168678/1‑100 (MQ=255)
                                                 aGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTAtctc                                           <  2:1020785/100‑1 (MQ=255)
                                                 aGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTAtctc                                           >  1:711135/1‑100 (MQ=255)
                                                           tCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTg                                 >  1:616896/1‑100 (MQ=255)
                                                           tCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTg                                 >  1:402811/1‑100 (MQ=255)
                                                             cggcggCGATTGATCTGCAAATTGCCCAGGATCGTTCCTCCACCATTCGCGCCTCGCTCGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGt                               <  2:859916/100‑1 (MQ=255)
                                                              ggcggcGATTGATCTGCACATTGCCCAGGGTCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTg                              <  2:556817/100‑1 (MQ=255)
                                                                cggcgATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCACTGGTGAt                            >  1:1183534/1‑100 (MQ=255)
                                                                   cgATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCt                         <  1:673232/100‑1 (MQ=255)
                                                                    gATTGATCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCtg                        >  1:165237/1‑100 (MQ=255)
                                                                          tCTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGgtggtg                  <  1:284784/100‑1 (MQ=255)
                                                                           cTGCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGgtgt                 >  1:480690/1‑100 (MQ=255)
                                                                             gCAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGttt               >  2:1184960/1‑100 (MQ=255)
                                                                              caccgtGCCCACGATCGTCCCCCCCCCCGTCGCGTCTCGTTGGAAGAAGTCAAGCCAGCGCTGATTATCTCGGTGGCGCTCGTGATTCTGGTGGTGtttt              <  2:1290397/95‑1 (MQ=255)
                                                                              cAAATTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGtttt              <  1:1430308/100‑1 (MQ=255)
                                                                                 aTTGCCCAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTAt           <  1:1169071/100‑1 (MQ=255)
                                                                                     cccAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGGTGGTGATTCTGGAAGTAATTTTATTCCt       >  2:1169199/1‑100 (MQ=255)
                                                                                      ccAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTg      >  1:872348/1‑100 (MQ=255)
                                                                                       cAGGATCGTTCTCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTgc     >  1:1295860/1‑100 (MQ=255)
                                                                                       cAGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTgc     >  1:1582095/1‑100 (MQ=255)
                                                                                        aGGATCGTTCCCCGACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTgc     <  2:168964/99‑1 (MQ=255)
                                                                                        aGGATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTgcg    >  1:578192/1‑100 (MQ=255)
                                                                                          gATCGTTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCt  >  1:1498705/1‑100 (MQ=255)
                                                                                               |                                                                                              
AATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCT  >  NC_000913/2159259‑2159448

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: