Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,011,322 T→G 25.0% V406G (GTT→GGT)  U712_04945 → putative multidrug resistance ABC transporter ATP‑binding/permease protein YheH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,011,3220TG25.0% 55.6 / 8.6 43V406G (GTT→GGT) U712_04945putative multidrug resistance ABC transporter ATP‑binding/permease protein YheH
Reads supporting (aligned to +/- strand):  ref base T (11/21);  new base G (8/3);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 3.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGGGGATTGTTTCTATCGGTGTGCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAAGAGGTATTAAAGCATATTTC  >  CP006881/1011194‑1011454
                                                                                                                                |                                                                                                                                    
ggcggGGATTGTTTCTATCGGTGTGCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGATCTg                                                                                                                            >  6:137546/1‑139 (MQ=255)
ggcggGGATTGTTTCTATCGGTGTGCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTg                                                                                                                            >  2:82320/1‑139 (MQ=255)
ggcggGGATTGTTTCTATCGGTGTGCTGTACGCGTTTGTCGATTATTTGAACAGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGGTTTTGAACTg                                                                                                                            >  2:87750/1‑139 (MQ=255)
          gTTTCTATCGGTGTGCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTgaagaa                                                                                                                   <  5:291229/138‑1 (MQ=255)
           tttCTATCGGTGTGCTGTACGCGTTTCTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTTAACTGCTTGaagaaa                                                                                                                  >  3:343924/1‑138 (MQ=255)
                aTCGGTGTGCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGaa                                                                                                                               <  7:275451/120‑1 (MQ=255)
                aTCGGTGTGCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGGTTTTGaa                                                                                                                               >  8:275451/1‑120 (MQ=255)
                     tgtgCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACgg                                                                                                         >  1:67439/1‑137 (MQ=255)
                     tgtgCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGgaa                                                                                                       <  6:144929/139‑1 (MQ=255)
                       tgCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGGTTTTTTAGTGGGTGAATAAAAAAATAAGgcag                                                                                                      >  3:220029/1‑138 (MQ=255)
                         cTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGgaagaa                                                                                                    <  7:16571/138‑1 (MQ=255)
                         cTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGaagaag                                                                                                   >  6:126286/1‑139 (MQ=255)
                                  tttGTCGATTATTAGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCACAGCGGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAAc                                                                                          <  1:188348/139‑1 (MQ=255)
                                    tGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCg                                                                                        >  8:246659/1‑139 (MQ=255)
                                    tGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGGTTTTTAAATGCGTGTATAAAATAATAACGAATAAGCGGGGGCAccc                                                                                        >  2:231282/1‑138 (MQ=255)
                                      tCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAAATGCTTGAAGAAAAGAATAAGGAAGAAGCGGGTGAACCGGc                                                                                      >  3:188287/1‑139 (MQ=255)
                                                gAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAgcg                                                                             <  1:82320/138‑1 (MQ=255)
                                                gAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGGTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAgcg                                                                             <  1:87750/138‑1 (MQ=255)
                                                         aTTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTggg                                                                    <  2:67439/138‑1 (MQ=255)
                                                           ttCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGATCTGCGTGTAGAAAAGAAAAAGGAAGAAGGCGGGGCACCGGCGTAGGCGAGCGCGCTGGGGCg                                                                 >  6:331534/1‑139 (MQ=255)
                                                                ccGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTTGACAGCATGAAGAAAATCAAAACGCAGAAGCGGGGGCAACGAGGAAGGCGCGCGCGCGGGGGGGGGGTg                                                            >  6:30797/1‑139 (MQ=255)
                                                                 cGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGa                                                           <  5:126286/139‑1 (MQ=255)
                                                                       caGGCATCGTCAATCAGTTTTCAAAGCTGGAGCGTGCAAGGGTCGCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCg                                                     <  2:360653/139‑1 (MQ=255)
                                                                            aTCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGAtg                                                <  4:343924/139‑1 (MQ=255)
                                                                             tCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGAtgt                                               <  3:59661/139‑1 (MQ=255)
                                                                                    taaGTCTTCAAAGATGCAGGATGCCAGGGTCGCCGCCGGGCGTGGTTATGAACTGCTTGAAGAAAAGAATACGTCAGACGCGGGTGACCCAGCGAAGGAGCGCGCGCTGGGGCGGGGTGAATTTCGTGATGTGTCAttt                                        <  4:74143/137‑1 (MQ=255)
                                                                                    tCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCAttt                                        <  4:338584/139‑1 (MQ=255)
                                                                                    tCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGGTTTTTAAATGGTTGAAGAAAAGAAAAACGAAGAAGCGGGTGAACCCGCGAAGGAGCGCGCGGTGGGGCGGGTTGATTTTCGTGATGTGTCAttt                                        >  8:368260/1‑139 (MQ=255)
                                                                                        ttttCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCGTGAAAAAAAGAATACGGCAGAAGCGGGGGAACCGGCGAAGGAGCGCGCGCGGGGGGGGGGTGAATTTAGGGATGTGTCAttgtcgt                                    >  6:143078/1‑134 (MQ=255)
                                                                                             aaaGCTGGAGCTTGCACGGGGCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCa                                <  1:191439/138‑1 (MQ=255)
                                                                                             aaaGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCa                                <  2:65137/138‑1 (MQ=255)
                                                                                               aGCTGGAGCTTTCAAGGGTCTGCGCCGGTCGTGTTTTTGAAATGCTTGAAGAAAAGAATACGGTAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCaaga                             <  5:25525/139‑1 (MQ=255)
                                                                                                    gAGCTTGCAAGGGTCTCCGCCGGGCGTGTGTTTTAAATGGTTGAAAAAAATAAATAGGAACAAACCGGTGTACCGGCCAAGGACAGCGCGCGTGGGGGGGTTGAAATTTATGAAGTGTCATTTAGGGATCaataacggg                        >  4:140173/1‑135 (MQ=255)
                                                                                                       cTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGACGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAag                     <  4:220029/139‑1 (MQ=255)
                                                                                                         tGCAAGGGTCTCCGCCGGGCGTGGTTTTTTACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGGGAACCGGCGAAGGAGCGCGCGCGGGGGGGGGGTGATTTTAGTGTTGTGTGCTTTGCGTATCAAGAAGGTGTagag                   >  8:217189/1‑139 (MQ=255)
                                                                                                          gCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAaga                    <  7:246659/137‑1 (MQ=255)
                                                                                                           cAAGGGTCTCCGCCGGGCGTGGTTTTTAACTGCGTGAAGAAAAGAATACGGAAGAAGCGGGGGGACCGGCGCAGGGGCGCGCGCGGGGGCGGGGGGAATTTCGTGATGTGTGATTTGCGTTTCAAGAAGGTGAAGAGGt                 >  3:278449/1‑139 (MQ=255)
                                                                                                               ggTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAAGAGGTATTa             <  3:91495/139‑1 (MQ=255)
                                                                                                                gTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAAAACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAAGAGGTATTaa            <  5:116207/139‑1 (MQ=255)
                                                                                                                     cgccgGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGGGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAAGAGGTATTAAAGCat       <  1:283932/139‑1 (MQ=255)
                                                                                                                     cgccgGGCGTGGTTTTGAACTGCTTGAAGAAAAGAATAAGGAAGAAGCGGGGGAACCGGCGAAGGATCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAAGAGGTATTAAAGCat       <  5:30797/139‑1 (MQ=255)
                                                                                                                      gccgGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAAGAGGTATTAAAGCata      <  1:231282/139‑1 (MQ=255)
                                                                                                                          ggCGTGGTTTTGAACTGCGTGAAGAAAAGAAGACGGAAGCAGCGGGTGAGCCGGCGAAGGAGCGCGCGCTGGGGCGGGTGGATTTTCGTGAGGTGTCATTTGCGTATCAAGAAGGTGAAGAGGTTGTTAAGCATATTTc  >  8:109269/1‑139 (MQ=255)
                                                                                                                                |                                                                                                                                    
GGCGGGGATTGTTTCTATCGGTGTGCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAAGAGGTATTAAAGCATATTTC  >  CP006881/1011194‑1011454

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: ‑