Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,585,535 G→A 20.0% P285S (CCG→TCG)  ttgB ← efflux pump membrane protein TtgB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,585,5350GA20.0% 38.6 / 4.9 20P285S (CCG→TCG) ttgBefflux pump membrane protein TtgB
Reads supporting (aligned to +/- strand):  ref base G (7/9);  new base A (2/2);  total (9/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 = 6.70e-01

CCCGGCGGGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCGTCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCG  >  NC_002947/1585429‑1585659
                                                                                                          |                                                                                                                            
cccGGCGGGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGt                                                                                                 <  2:159943/136‑1 (MQ=255)
     ggggAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGttt                                                                                            <  2:401938/135‑1 (MQ=255)
     cggGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGttt                                                                                            >  1:467979/1‑136 (MQ=255)
          agaaCGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCAcc                                                                                       <  1:416019/136‑1 (MQ=255)
                                    tcGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACg                                                                                                     <  2:465437/96‑1 (MQ=255)
                                    tcGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACg                                                                                                     >  1:465437/1‑96 (MQ=255)
                                        gcAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGTGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAg                                                         <  1:199403/136‑1 (MQ=255)
                                        gcAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAg                                                         <  2:262503/136‑1 (MQ=255)
                                               cTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCAcc                                                  >  1:6164/1‑136 (MQ=255)
                                                 tGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGACTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCtg                                                >  2:143117/1‑136 (MQ=255)
                                                 tGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGACTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCtg                                                >  2:146278/1‑136 (MQ=255)
                                                   gCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCtgtg                                              <  1:287471/136‑1 (MQ=255)
                                                         tCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCGt                                        >  1:290206/1‑136 (MQ=255)
                                                         tCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGCCGCTGACGGCGTAGATTTCACCCCCCAGGCCTACTTGGGCGACAACGCCCAGGCGCACCTGTGAACCGt                                        >  1:415954/1‑136 (MQ=255)
                                                                      cgcCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAg                                                         <  1:160791/106‑1 (MQ=255)
                                                                      cgcCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAg                                                         >  2:160791/1‑106 (MQ=255)
                                                                      cgcCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGACTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCGTCTTTGTTGACCtt                           <  1:146278/136‑1 (MQ=255)
                                                                      cgcCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGACTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCGTCTTTGTTGACCtt                           <  1:143117/136‑1 (MQ=255)
                                                                                  gTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCGTCTTTGTTGACCTTGAGCAGGATAct               >  2:576452/1‑136 (MQ=255)
                                                                                               gCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCGTCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTcg  <  1:211571/136‑1 (MQ=255)
                                                                                                          |                                                                                                                            
CCCGGCGGGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCGTCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCG  >  NC_002947/1585429‑1585659

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