Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 6,133,865 C→G 10.5% P412R (CCT→CGT)  cusB → Cation efflux system protein CusB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24406,133,8650CG10.5% 123.6 / 6.0 38P412R (CCT→CGT) cusBCation efflux system protein CusB
Reads supporting (aligned to +/- strand):  ref base C (13/21);  new base G (2/2);  total (15/23)
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.94e-01

TCGCAATCTTGAAAGGCCTGCAGGCCGGGCAGAAAGTCGTGGCATCCGGGCAATTCCTGATCGACTCCGAAGCAAGCCTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCGAAGGCACCCAGCTGACCATCGACCACGGGCCATTCGTCACCTTGGGAATGCCTGGGATGACCATGACCTTTCCCGTCGCTGATCCC  >  KT2440/6133733‑6133995
                                                                                                                                    |                                                                                                                                  
tCGCAATCTTGAACGGCCGGCCGGCCGGGCAGAAAGTCGTGGCATCCGGGCAATTCCTGATCGACTCCGAAGCAAGCCTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTg                                                                                                                                  <  8:103742/135‑1 (MQ=255)
             aGGCCTGCAGGCCGGGCAGAAAGTCGTGGCATCCGGGCAATTCCTGATCGACTCCGAAGCAAGCCTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGc                                                                                                                      <  6:305328/134‑1 (MQ=255)
                   gCAGGCCGGGCAGAAAGTCGTGGCATCCGGGCAATTCCTGATCGACTCCGAAGCAAGCCTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGc                                                                                                               >  4:271980/1‑135 (MQ=255)
                     aGGCCGGGCAGAAAGTCGTGGCATCCGGGCAATTCCTGATCGACTCCGAAGCAAGCCTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCgg                                                                                                             >  6:85118/1‑135 (MQ=255)
                     aGGCCGGGCAGAAAGTCGTGGCATCCGGGCAATTCCTGATCGACTCCGAAGCAAGCCTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCgg                                                                                                             >  5:327450/1‑135 (MQ=255)
                      ggCGGGGCAGAATGTCGTGGCATCCGGGCAATTCCTGATCGACTCCGAAGCAAGCCTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGa                                                                                                            <  1:285941/135‑1 (MQ=255)
                            gCAGAAAGTCGTGGCATCCGGGCAATTCCTGATCGACTCCGAAGCAAGCCTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACGTGAACCCGAACTGCACCAAGCGGATGGcca                                                                                                      >  4:65929/1‑134 (MQ=255)
                            gCAGAAAGTCGTGGCATCCGGGCAATTCCTGATCGACTCCGAAGCAAGCCTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACGTGAACCCGAACTGCACCAAGCGGATGGcca                                                                                                      >  4:61267/1‑134 (MQ=255)
                                aaaGTCGTGGCATCCGGGCAATTCCTGATCGACTCCGAAGCAAGCCTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTc                                                                                                  <  3:271980/135‑1 (MQ=255)
                                     cGTGGCTTCCGGGCAATTCCTGATCGACTCCGAAGCAAGCCTCAAGGGTATGGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACa                                                                                             <  4:112455/135‑1 (MQ=255)
                                        ggCATCCGGGCAATTCCTGATCGACTCCGAAGCAAGCCTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAAt                                                                                          <  3:215512/135‑1 (MQ=255)
                                          cATCCGGGCAATTCCTGATCGACTCCGAAGCAAGCCTCACGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCg                                                                                        <  4:124327/135‑1 (MQ=255)
                                          cATCCGGGCAATTCCTGATCGACTCCGAAGCAAGCCTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCg                                                                                        <  2:46323/135‑1 (MQ=255)
                                                      tCCTGATCGACTCCGAAGCAAGCCTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCGAAGGCACCCAGc                                                                            <  3:181352/135‑1 (MQ=255)
                                                        cTGATCGACTCCGAAGCAAGCCTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCGAAGGCACCCAGCTg                                                                          >  3:235808/1‑135 (MQ=255)
                                                         tGATCGACTCCGAAGCAAGCCTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCGAAGGCACCCAGCTGa                                                                         <  1:127726/135‑1 (MQ=255)
                                                              gACTCCGAAGCAAGCCTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCGAAGGCACCCAGCTGACCATc                                                                    >  2:125634/1‑135 (MQ=255)
                                                                     aagcaagcCTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCGAAGGCACCCAGCTGACCATCGACCACg                                                             <  3:124365/135‑1 (MQ=255)
                                                                            ccTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCGAAGGCACCCAGCTGACCATCGACCACGGGCCAtt                                                      <  8:281923/135‑1 (MQ=255)
                                                                             cTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCGAAGGCACCCAGCTGACCATCGACCACGGGCCATTc                                                     >  1:299838/1‑135 (MQ=255)
                                                                             cTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCGAAGGCACCCAGCTGACCATCGACCACGGGCCATTc                                                     >  8:312005/1‑135 (MQ=255)
                                                                                          aGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGaa                                                                                                                          >  5:338561/1‑53 (MQ=255)
                                                                                          aGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGaa                                                                                                                          <  6:338561/53‑1 (MQ=255)
                                                                                           ggCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCGAAGGCACCCAGCTGACCATCGACCACGGGCCATTCGTCACCTTGGGAAt                                       <  6:327450/135‑1 (MQ=255)
                                                                                                ccacAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCGAAGGCACCCAGCTGACCATCGACCACGGGCCATTCGTCACCTTGGGAATGCCTg                                  <  7:312005/135‑1 (MQ=255)
                                                                                                ccacAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCGAAGGCACCCAGCTGACCATCGACCACGGGCCATTCGTCACCTTGGGAATGCCTg                                  >  2:235637/1‑135 (MQ=255)
                                                                                                 cacaGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCaa                                                                                                                 <  7:84137/55‑1 (MQ=255)
                                                                                                 cacaGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCaa                                                                                                                 >  8:84137/1‑55 (MQ=255)
                                                                                                  acaGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCGAAGGCACCCAGCTGACCATCGACCACGGGCCATTCGTCACCTTGGGAATGCCTggg                                >  4:295294/1‑135 (MQ=255)
                                                                                                  acaGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCGAAGCCACCCAGCTGACCATCGACCACGGGCCATTAGTCACCTTGGGAATGCCTggg                                >  3:16807/1‑135 (MQ=255)
                                                                                                     gCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCGAAGGCACCCAGCTGACCATCGACCACGGGCCATTCGTCACCTTGGGAATGCCTGGGatg                             <  7:7558/135‑1 (MQ=255)
                                                                                                       aTCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCGAAGGCACCCAGCTGACCATCGACCACGGGCCATTCGTCACCTTGGGAATGCCTGGGatgac                           <  3:295294/135‑1 (MQ=255)
                                                                                                           gATTCACCGCCTGCCTCGCAAGCACGTGAACCCGAACTGCACCAAGCGGATGGCCATGTCGTACAAATCGAAGGCACCCAGCTGACCATCGACCACGGGCCATTCGTCACCTTGGGAATGCCTGGGatgaccatg                       <  3:65929/135‑1 (MQ=255)
                                                                                                           gATTCACCGCCTGCCTCGCAAGCACGTGAACCCGAACTGCACCAAGCGGATGGCCATGTCGTACAAATCGAAGGCACCCAGCTGACCATCGACCACGGGCCATTCGTCACCTTGGGAATGCCTGGGatgaccatg                       <  3:61267/135‑1 (MQ=255)
                                                                                                                 ccGCCTGGCTCGCAAGGACCTGAACCCGAACTGCACCAAGCGGATGGGCGTGTCGTACAAATCGCAGGCACCCAGCTGACCATCGACCACGGGCCATTCGTCACCTTGGGAATGCCTGGGATGAGCATGACCttt                 <  4:235808/135‑1 (MQ=255)
                                                                                                                         ctCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCGAAGGCACCCAGCTGACCATCGCCCACGGGCCATTCGTCACCTTGGGAATGCCAGGGATGACCATCACCTTTCCCGTc           >  4:230503/1‑133 (MQ=255)
                                                                                                                             cAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCGAAGGCACCCAGCTGACCATCGACCACGGGCCATTCGTCACCTTGGGAATGCCTGGGATGACCATGACCTTTCCCGTCGCTg       <  5:85118/133‑1 (MQ=255)
                                                                                                                                gCACCTGAACCCGAACTGCACCAAGCGGATGGGCGTGTCGTACATATCGAAGGCAGCCAGCTGACCTTCGACCACGGGCCATTGGTCAGCTTGGGAATGCCTGGGATGACCATGACCTTTCCCGTCGCTGATccc  <  6:22916/135‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                  
TCGCAATCTTGAAAGGCCTGCAGGCCGGGCAGAAAGTCGTGGCATCCGGGCAATTCCTGATCGACTCCGAAGCAAGCCTCAAGGGTATTGAGGCCACCACAGCATCGGATTCACCGCCTGCCTCGCAAGCACCTGAACCCGAACTGCACCAAGCGGATGGCCGTGTCGTACAAATCGAAGGCACCCAGCTGACCATCGACCACGGGCCATTCGTCACCTTGGGAATGCCTGGGATGACCATGACCTTTCCCGTCGCTGATCCC  >  KT2440/6133733‑6133995

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