Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 465,464 T→C S380P (TCG→CCG)  mdlA → SmdA family multidrug ABC transporter permease/ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273465,4640TC90.5% 51.0 / ‑0.1 21S380P (TCG→CCG) mdlASmdA family multidrug ABC transporter permease/ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (9/10);  total (10/11)
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.61e-01

GCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTG  >  NZ_CP009273/465377‑465549
                                                                                       |                                                                                     
gcgcTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCg                                                                                     >  2:127044/1‑90 (MQ=255)
 cgcTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGc                                                                                    >  2:446462/1‑90 (MQ=255)
 cgcTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGc                                                                                    <  2:275414/90‑1 (MQ=255)
     ggAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCAt                                                                                >  1:57851/1‑90 (MQ=255)
     ggAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCAt                                                                                >  1:383758/1‑90 (MQ=255)
                              cccGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAg                                                       <  1:142299/90‑1 (MQ=255)
                                            gggTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCtt                                         >  2:566013/1‑90 (MQ=255)
                                            gggTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCtt                                         >  1:61424/1‑90 (MQ=255)
                                             ggTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCttt                                        <  2:57851/90‑1 (MQ=255)
                                                aTCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCAt                                     >  2:690539/1‑90 (MQ=255)
                                                aTCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCAt                                     <  2:383758/90‑1 (MQ=255)
                                                 tCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATg                                    >  2:118574/1‑90 (MQ=255)
                                                    gCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGata                                 >  2:440243/1‑90 (MQ=255)
                                                     cGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGatat                                <  1:566013/90‑1 (MQ=255)
                                                     cGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGatat                                <  1:148415/90‑1 (MQ=255)
                                                      gggCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCtt                                         <  2:446210/80‑1 (MQ=255)
                                                      gggCCGACTGGTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCtt                                         >  1:446210/1‑80 (MQ=255)
                                                               ggTTCCGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACg                      <  2:114084/90‑1 (MQ=255)
                                                                    cGGCAAAAGTACCCTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGtt                 <  1:753512/90‑1 (MQ=255)
                                                                                 cTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGc    <  2:370624/90‑1 (MQ=255)
                                                                                 cTGTTGCCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGc    <  1:479459/90‑1 (MQ=255)
                                                                                   gttgcCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGTTATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTg  >  2:211823/1‑90 (MQ=255)
                                                                                       |                                                                                     
GCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTG  >  NZ_CP009273/465377‑465549

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