Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,749,309 C→G 24.4% V301V (GTC→GTG macB → Macrolide export ATP‑binding/permease proteinMacB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,749,3090CG24.4% 58.7 / 2.0 39V301V (GTC→GTGmacBMacrolide export ATP‑binding/permease proteinMacB
Reads supporting (aligned to +/- strand):  ref base C (9/19);  new base G (8/2);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 2.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGAAAGGTGAACTGCTCGAATCGCTGCAGGCGGCATGGCGAGTGATGTGGATCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTCGGTGGTGGTGATGTTGGCTGTCGGTGAAGGCAGCAAACGCCAGGTAATGGCACAGATGGCCGCCTTCGGCTCCAACATCCTATACCTCAATGGCAGCCCCCCGACCCTGCGCGAACCGGCTGGGCGCATCACCCTGGACGATGTT  >  KT2440/4749180‑4749432
                                                                                                                                 |                                                                                                                           
cTGGAAAGGTGAACTGCTCGAATCGCTGCAGGCGGCATGGCGAGTGATGTGGATCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTGGGTGGGGGTGTTTTTTGGTGTGGGt                                                                                                                          >  2:138985/1‑133 (MQ=255)
cTGGAAAGGTGAACTGCTCGAATCGCTGCAGGCGGCATGGCGAGTGATGTGGATCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTCGGTGGTGGTGATGTTGGGTGTCGGTg                                                                                                                         >  4:10580/1‑134 (MQ=255)
cTGGAAAGGTGAACTGCTCGAATCGCTGCAGGCGGCATGGCGAGTGATGTGGATCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTCGGTGGTGGTGATGTTGGCTGTCGGTg                                                                                                                         >  4:161524/1‑134 (MQ=255)
cTGGAAAGGTGAACTGCTCGAATCGCTGCAGGCGGCATGGCGAGTGATGTGGATCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTCGGTGGTGGTGATGTTGGCTGTCGGTg                                                                                                                         >  4:80632/1‑134 (MQ=255)
cTGGAAAGGTGAACTGCTCGAATCGCTGCAGGCGGCATGGCGAGTGATGTGGATCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTCGGTGGTGGTGAGGTTGGCTGTCGGt                                                                                                                          >  8:283199/1‑133 (MQ=255)
cTGGAAAGGTCAACTGCTCGAATCGCTGCAGGCGGCATGTCGAGTGATGTGGATCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTCGGTGGTGGGGATGTTGGGTGTCGGCGa                                                                                                                        >  2:241817/1‑135 (MQ=255)
   taaaGGTGAACTGCTCTAATCGCTGCAGGGGGAATGGCGAGTGAGGTGGATCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTCGGTGGTGGTGATGTTGGCTGTCGGTGAAGgc                                                                                                                    <  4:356496/135‑1 (MQ=255)
          gAACTGCTCGAATCGCTGCAGGCGGCATGGCGAGTGATGTGGATCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTCGGTGGTGGTGATGTTGGCTGTCGGTGAAGGCAGCaaa                                                                                                              <  5:151429/135‑1 (MQ=255)
           aaCTGCTCGAATCGCTGCAGGCGGCATGGAGAGTGATGTGGATCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGGGGGGGGTGGGGATGTTGGCAGGGGGGGACGGCAGCAAAc                                                                                                             >  2:173951/1‑135 (MQ=255)
             cTGCTCGAATCGCTGCAGGCGGCATGGCGAGTGATGTGGATCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTCGGGGGTGGTGATGTTGGCTGTCGGTGAAGGCAGCAAACg                                                                                                            <  8:150790/134‑1 (MQ=255)
               gCTCGAACCGCTGCAGGCGGCAGGGCGAGGGAGGTGGAGCACCCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCAGCGGAGTCGCGTCGGTGGTGGGGATGTTGGCTGTCGGTGAAGGCAGCAAACGCCAg                                                                                                        <  3:76637/136‑1 (MQ=255)
                   gAATCGCTGCAGGCGGCATGGCGAGTGATGTGGATCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTCGGTGGTGGTGATGTTGGGTGTCGGTGCAGGCAGCAAACGCCAGGTa                                                                                                     >  6:19643/1‑135 (MQ=255)
                     aTCGCTGCAGGCGGCATGGCGAGTGATGTGGATCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTCGGTGGTGGTGATGTTGGCTGTCGGTGAAGGCAGCAAACGCCAGGTAATg                                                                                                  <  3:80632/136‑1 (MQ=255)
                     aTCGCTGCAGGCGGCATGGCGAGTGATGTGGATCAACCGCTTCCGCACCGCACTCACCCTGCTCGGAATCATCATCGGAGTCGCGGCGGTGGGGGAGATTTTGGCTGTCGGGGAATGGAGGAAACGCCAGGGAAt                                                                                                   >  8:157275/1‑135 (MQ=255)
                                   cATGGCGAGTGAGGTGGAGCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTCGGTGGTGGTGATGTTGGCTGTCGGTGAAGGCAGCAAACGCCAGGTAATGGCACAGATGgccg                                                                                     <  1:201244/135‑1 (MQ=255)
                                             gaggtgGATCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTCGGTGGTGGTGATGTTGGCGGTCGGTGAAGGCAGCAAACGCCAGGTAATGGCACAGATGGCCGCCTTCGGCTc                                                                           <  5:19643/132‑1 (MQ=255)
                                             gATGTGGATCACCCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTCGGTGGTGGTGATGTTGGCTGTCGGTGAAGGCAGCAAACGCCAGGTAATGGCACAGATGGCCGCCTTCGGCTc                                                                           <  1:138985/135‑1 (MQ=255)
                                             gATGTGGATCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCAGCGGAGTCGCGTCGGTGGTGGTGATGTGGGCGGGCGGTGAAGGCAGCAAACGCCAGGTAATGGCACAGATGGCCGCCTTCGGCTc                                                                           <  2:119583/135‑1 (MQ=255)
                                             gATGTGGATCAACCGCTTCCGCACGGCACTGACCCGGCTCGGAATCATCATCGGAGTCGCGGCGGGGGTGGTGATGTTGGCTGTCGGTGAAGGCAGCAAACGCCAGGTAATGGCACAGATGGCCGCCTTCGGCt                                                                            <  8:246382/134‑1 (MQ=255)
                                              aggtgGATCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTCGGTGGTGGTGATGTTGGCTGTCGGTGAAGGCAGCAAACGCCAGGTAATGGCACAGATGGCCGCCTTCGGCTcc                                                                          <  4:106107/133‑1 (MQ=255)
                                              aTGTGGATCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGGCGGGGGGGGGGATGTTGTGTGTGGGTGGAGGCCCCCAACCCCCGGTCATGGCACAGCTGGCCCCCGATGGCTcc                                                                          >  2:163746/1‑135 (MQ=255)
                                               ggtgGATCCACCGCTTCCGCACGGCACTGACCCTGCTCGGATGCATCCGCGGAGTCGCGTCGGGGGTGGGGATGTTGGCTGTCGGTGAAGGCAGCAAACGCCAGGTAATGGCACAGATGGCCGCCTTCGGCTCCa                                                                         <  1:129870/134‑1 (MQ=255)
                                                   gATCAACGGCTTGCGCACGGAACTGACCCTGCTCGTAATCATCATCGGAGGCGCGTCGGTGGGGGTGATGTTGGCTGTGGGTTAAGGCAGCAAACGCCAGGTTATTGCACAGATGGCCGCCTTCGGCTCCAACATc                                                                    <  3:103100/136‑1 (MQ=255)
                                                    aTCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGGCGGTGGTGGTGATGTTGGGTGTGGGGGAAGGCAGCCAACGCCAGGTAATGGCACAGATGGCCGCCGTCGGCTCCCACATcc                                                                   >  3:230944/1‑136 (MQ=255)
                                                           gCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGGCGGTGGTGGTGATGTTGGCTGTGGGGGCAGGCAGCAAACGCCAGGGTATGGCACAGATGGCCGCCTTCGGCTCCAACATCCTATACCt                                                            >  3:347033/1‑136 (MQ=255)
                                                            cTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTGGGGGGGGGGGTTGTTGGGGGTGGGGGAAGGGAGCCAACGCCCGGGAATGGGCCAGATTGCCCCCGTCGGCGCCACCATCCTATAcc                                                             >  6:248272/1‑134 (MQ=255)
                                                               ccGCACGGCACTCACCCTGCTCGGAATCATCATCGGAGTCGCGTGGGGGGTGGTGATGTTGGCTGGCGGTGAAGGCAGCAAACGCCAGGTAATGGCACAGATGGCCGCCTTCGGCTCCAACATCCTATACCTCAAt                                                        <  3:10580/136‑1 (MQ=255)
                                                                cGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGGGTCGGTGGTGGTGATGTTGGCTGTCGGTGAAGGCAGCAAACGCCAGGTAATGGCACAGATGGCCGCCTTCGGCTCCAACATCCTATACCTCAAt                                                        <  5:48338/135‑1 (MQ=255)
                                                                  cACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTCGGTGGGGGTGATGTTGGCTGTCGGTGAAGGCAGCAAACGCCAGGTAATGGCACAGATGGCCGCCTTCGGCTCCAACATCCTATACCTCAATgg                                                      <  7:236117/135‑1 (MQ=255)
                                                                     ggCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTCGGTGGTGGTTATGTTGGCTGTCGGTGCAGGCAGCAAAAGCCAGGTAATTGCACAGATGGCCGCCTTCGGCTCCACAAACCTAGACCTAAATCGCac                                                   >  7:257821/1‑134 (MQ=255)
                                                                               cTGCTCGGAATCATCATCGGAGTCGCGGCGGTTGTGGTGATTTTGGGTGTGGGTGAAGGCAGCAAACACCCCGTAATGGCACAGATGGCCCCCGTCGGCGCCAACATACGAAACCGCCAATGCCGCCCCCCCAcc                                         >  7:151181/1‑135 (MQ=255)
                                                                               cTGCTCGGAATCATCATCGGAGTCGCGGCGGGGGGGGTGATGTTGGCGGTGGGGTGAGGCAGCAAACGCCCGGTTATGGGACAGATTGCCCCCTTACGCGCCAACAACCAAAAACTAAATGGCCGCCCCCCCAccc                                        >  4:256179/1‑136 (MQ=255)
                                                                                   tCGGAATCATCATCGGAGTCGCGTCGGGGGTGGTGAGGTTGGCTGTCGGTGAAGGCAGCAAACGCCAGGTAATGGCACAGATGGC‑‑CCTTCGGCTCCAACATCCTATACCTCAATGGCAGCCCCCCGACCCTgc                                     <  2:237527/133‑1 (MQ=255)
                                                                                    cGGAATCATCATCGGAGTCGCGTCGGTGGGGGTGATGTTGGCTGTCGGTGAAGGCAGCAAACGCCAGGTAATGGCACAGATGGCCGCCTTCGGCTCCAACATCCTATACCTCAATGGCAGCCCCCCGCCCCTgcg                                    >  8:63656/1‑135 (MQ=255)
                                                                                         tcatcatcGGAGTCGCGTCGGTGGTGGTGATGTTGGCGGGAGGGGAAGGCAGCAAACGCCAGGTAAATGCACAGATTGCCGCCTTCGGCTCCACCATACAATAACTCAATGGCCGACCCCCCACCCAgcccccac                               >  7:245585/1‑129 (MQ=255)
                                                                                            tcatcGGAGTCGCGTCGGTGGTGGTGATGTTGGCTGTCGGTGAAGGCAGCAAACGCCAGGTAATGGCACAGATGGCCGCCTTCGGCTCCAACATCCTATACCTCAATGGCAGCCCCCCGACCCTGCGCGAACCGGc                           <  4:230944/136‑1 (MQ=255)
                                                                                                   tggcgcgTCGGGGGTGGGGATGTTGGCTGTCGGGGACGGCAGCAAACGCCAGGTAATCGCACAGATGGCCGCCTTCGGCTCCAACATCCTATACCTCAATGGCAGCCCCCCGACCCTGCGCGAACCGGCTGGgcg                     <  1:275450/132‑1 (MQ=255)
                                                                                                            gggtgtggtgAGGTGGGCGGGGGGTGAAGGCAGAAAACGCCAGGTAATGGCACAGATGGCCGCCTTCGGCTCCAACATCCTATACCTCAATGGCAGCCCCCCGACCCTGCGCGAACCGGCTGGGCGCATCACCCt            <  1:163746/131‑1 (MQ=255)
                                                                                                                      aTGTTGGCTGTCGGTGAAGGCAGCAAACGCCAGGTAATGGCACAGATGGCCGCCTTCGGCTCCAACATCCTATACCTCAATGGCAGCCCCCCGACCCTGCGCGAACCGGCTGGGCGCATCACCCTGGACGATGtt  <  5:248272/135‑1 (MQ=255)
                                                                                                                                 |                                                                                                                           
CTGGAAAGGTGAACTGCTCGAATCGCTGCAGGCGGCATGGCGAGTGATGTGGATCAACCGCTTCCGCACGGCACTGACCCTGCTCGGAATCATCATCGGAGTCGCGTCGGTGGTGGTGATGTTGGCTGTCGGTGAAGGCAGCAAACGCCAGGTAATGGCACAGATGGCCGCCTTCGGCTCCAACATCCTATACCTCAATGGCAGCCCCCCGACCCTGCGCGAACCGGCTGGGCGCATCACCCTGGACGATGTT  >  KT2440/4749180‑4749432

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