Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,011,322 T→G 27.5% 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,3220TG27.5% 62.5 / 9.3 45V406G (GTT→GGT) U712_04945putative multidrug resistance ABC transporter ATP‑binding/permease protein YheH
Reads supporting (aligned to +/- strand):  ref base T (5/27);  new base G (10/3);  total (15/30)
Fisher's exact test for biased strand distribution p-value = 1.75e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGCGGCGGGGATTGTTTCTATCGGTGTGCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAAGAGGTATTAAAGCATAT  >  CP006881/1011189‑1011451
                                                                                                                                     |                                                                                                                                 
aaTGCGGCGGGGATTGTTTCTATCGGTGTGCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGGTGTAGCTTGCAAGGGTCTCCGCCGGGCGTGgttttt                                                                                                                              >  2:350553/1‑138 (MQ=255)
  tGCGGCGGGGATTGTTTCTATCGGTGTGCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGaa                                                                                                                            <  6:231622/139‑1 (MQ=255)
  tGCGGCGGGGATTGTTTCTATCGGTGTGCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGaa                                                                                                                            <  5:216352/139‑1 (MQ=255)
  tGCGGCGGGGATTGTTTCTATCGGTGTGCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGaa                                                                                                                            <  5:121996/139‑1 (MQ=255)
           gATTGTTTCTATCGGTGTGCTGTACGCGTTTGTCGATTATTTGATCCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTtaa                                                                                                                   <  1:218531/139‑1 (MQ=255)
               gTTTCTATCGGTGTGCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTgaagaa                                                                                                                <  7:5451/138‑1 (MQ=255)
                tttCTATCGGTGTGCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGGCTCCGCCGGGCGTGTTTTTGAACTGCTTGaagaaaa                                                                                                              <  6:79809/139‑1 (MQ=255)
                     aTCGGTGTGCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGGTTTTTGTCTGCGTGAAGATAAGAAt                                                                                                          >  5:52531/1‑138 (MQ=255)
                     aTCGGTGTGCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGGTTTTTAACGTAGTGaagaaaagaaa                                                                                                          >  6:37603/1‑137 (MQ=255)
                      tCGGTGTGCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGGGTTTGTACTGCTTGAATAAAAAAAAAc                                                                                                        >  2:93896/1‑139 (MQ=255)
                          tgtgCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCATTCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTTAACTGCTTTAAGAAAAGAATAAgg                                                                                                      >  3:193604/1‑137 (MQ=255)
                             gCTGTACGCGTTTGTCGATTATTTGAACCGGGTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCGTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGgaagaa                                                                                                 <  5:159781/139‑1 (MQ=255)
                              cTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGgaagaa                                                                                                 <  5:159783/138‑1 (MQ=255)
                                   cgcgTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGGTTTTGGACTGCGTGAAGAAAAGAATACGGAAGAAGCGGGt                                                                                           >  2:261090/1‑139 (MQ=255)
                                   cgcgTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAAGCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGTACTGCGTGGAGCAAAGAATAAGGTAGCAGCGGGt                                                                                           >  6:369530/1‑139 (MQ=255)
                                    gcgTTTGTCGATTATTTGAAACGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGGTTTTTTAATGCTTGAAGAAAATAATACGGAAGAAGCgggcg                                                                                          >  6:356023/1‑137 (MQ=255)
                                      gTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTgca                                                                                        >  3:389810/1‑137 (MQ=255)
                                       tttGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAAc                                                                                       <  1:6029/139‑1 (MQ=255)
                                         tGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGGTTTTTTAAAGCTTGAATAAAAAAAAAAGCCACAAACGGGTGATccc                                                                                     >  6:362146/1‑138 (MQ=255)
                                         tGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGGTTTTGAACAGGATGAAGAAAATAATAAGGAAGAAGCCGGGGAACCg                                                                                     >  3:411518/1‑139 (MQ=255)
                                                  tttGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGGTGTTTTCGTGGTGGAAGAAAAGAATACGGAAGAAGCGGGGGCACCCGCGCACGAg                                                                            >  2:368111/1‑139 (MQ=255)
                                                        gaaagtATTTCAGCCGATCACCGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAgcgcgcg                                                                      <  5:300488/134‑1 (MQ=255)
                                                                        aTCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGGTTTTTTAATTCTATAAGAAAATAAAAAGGAACTCGCGGGTGGACCCGCGAAAGACGGCGCGCGGGGGCGGGtgtgtg                                                      >  6:300488/1‑134 (MQ=255)
                                                                                   cGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTTTACTGCTTGAAGAAAAGAATAAGGAAGAAGCGGGTGCACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTGGGGAtgtg                                           >  7:148499/1‑139 (MQ=255)
                                                                                   cGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGGTTTTTAACTGGTTGAAGAAAATAATAAGGAAGAAGCGGGGGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGGTGGTTTTAGTGTtgtg                                           >  4:322223/1‑139 (MQ=255)
                                                                                                  aaaGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCa                             <  3:349204/138‑1 (MQ=255)
                                                                                                  aaaGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGGTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCaa                            <  1:330753/139‑1 (MQ=255)
                                                                                                  aaaGCTGGAGCTTGCAAGGGGCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCa                             <  6:52531/138‑1 (MQ=255)
                                                                                                   aggCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCaag                           <  3:104352/137‑1 (MQ=255)
                                                                                                   aaGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTTAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCaag                           <  4:193604/139‑1 (MQ=255)
                                                                                                    aGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGGTTTTGAACTGCTTGAAGAAAAGAATACGGCAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCATATCaaga                          <  6:266806/139‑1 (MQ=255)
                                                                                                    aGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGGTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCaaga                          <  8:148499/139‑1 (MQ=255)
                                                                                                            cTTGCAAGGGTCTGCGCCGGGCGTGTTTTTGAACTGCATGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAag                  <  5:362146/139‑1 (MQ=255)
                                                                                                            cTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAATAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAag                  <  5:381780/139‑1 (MQ=255)
                                                                                                            cTTGCAAGGGGCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAag                  <  3:322223/139‑1 (MQ=255)
                                                                                                               ggcaGGGGCTCCGCCGGGCGTGTTTTTGACCTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAagag                <  2:381460/135‑1 (MQ=255)
                                                                                                                    gggctcCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAAGAGGTATTa          <  1:350553/136‑1 (MQ=255)
                                                                                                                     gTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAAGAgg               <  2:235480/133‑1 (MQ=255)
                                                                                                                     gTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAAGAgg               >  1:235480/1‑133 (MQ=255)
                                                                                                                          cgccgGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGGGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAAGAGGTATTAAAGCat    <  4:287440/139‑1 (MQ=255)
                                                                                                                          cgccgGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAAGAGGTATTAAAGCat    <  5:37603/139‑1 (MQ=255)
                                                                                                                          cgccgGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAAGAGGTATTAAAGCat    <  3:75733/139‑1 (MQ=255)
                                                                                                                          cgccgGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAAGAGGTATTAAAGCat    <  1:261090/139‑1 (MQ=255)
                                                                                                                          cgccgGGCGTGTGTTTGAACTGCGTGAAGAAAAGAATCCGGAAGAAGCGGGCGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAAGAGGTATTAAAGCat    <  5:356023/139‑1 (MQ=255)
                                                                                                                          cgccgGGCGAGTTTTTTCAGTGCTTGAAGACAAGAATACGGAAGAAGCGGGGGAACCGGCGTAGGAGCGCGCGCTGGGTCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAAGAGGTATTAAAGCat    <  6:62383/139‑1 (MQ=255)
                                                                                                                            ccgGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAAGAGGTATTAAAGCatat  <  1:93896/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                 
AATGCGGCGGGGATTGTTTCTATCGGTGTGCTGTACGCGTTTGTCGATTATTTGAACCGGTTATTTCAGCCGATCACAGGCATCGTCAATCAGTTTTCAAAGCTGGAGCTTGCAAGGGTCTCCGCCGGGCGTGTTTTTGAACTGCTTGAAGAAAAGAATACGGAAGAAGCGGGTGAACCGGCGAAGGAGCGCGCGCTGGGGCGGGTTGAATTTCGTGATGTGTCATTTGCGTATCAAGAAGGTGAAGAGGTATTAAAGCATAT  >  CP006881/1011189‑1011451

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