Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 922,060 C→T 18.3% A431V (GCC→GTC)  macB → fused macrolide transporter subunits and ATP‑binding component and membrane component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb922,0600CT18.3% 140.9 / 22.3 60A431V (GCC→GTC) macBfused macrolide transporter subunits and ATP‑binding component and membrane component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base C (25/24);  new base T (4/7);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTA  >  W3110S.gb/921991‑922130
                                                                     |                                                                      
cTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGATGTGATTCTGGTCGGCAATATGCCCGc                                                                        >  1:2417245/1‑70 (MQ=255)
cTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGc                                                                        >  1:3428657/1‑70 (MQ=255)
cTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGc                                                                        >  1:364449/1‑70 (MQ=255)
cTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGc                                                                        >  1:1349983/1‑70 (MQ=255)
cTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGc                                                                        >  1:2428860/1‑70 (MQ=255)
cTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGc                                                                        >  1:796078/1‑70 (MQ=255)
   gccgccAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCa                                                                      >  1:638012/1‑69 (MQ=255)
   gccgccAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCa                                                                      >  1:3199301/1‑69 (MQ=255)
         aGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGTCAGAGTCAt                                                               >  1:460540/1‑70 (MQ=255)
         aGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGTCAGAGTCAt                                                               >  1:1820359/1‑70 (MQ=255)
         aGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCAtt                                                              <  1:567964/71‑1 (MQ=255)
           cTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTgg                                                            <  1:382028/71‑1 (MQ=255)
                     aTAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGTCAGAGTCATTGGTGTGGCgg                                                    >  1:4021684/1‑69 (MQ=255)
                        aaGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGTCAGAGTCATTGGTGTGGCGgaaga                                                >  1:1670817/1‑70 (MQ=255)
                          gCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGaaa                                              >  1:3334677/1‑70 (MQ=255)
                                    ttGGGGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATg                                   <  1:3531227/71‑1 (MQ=255)
                                    ttGGCGAGGTGATTCTGGTCGGCAATATGCCCTCCAGAGTCATTGGTGTGGCGGATGAAAAACAGTCGATg                                   <  1:1023582/71‑1 (MQ=255)
                                      ggCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGtt                                 >  1:2125638/1‑71 (MQ=255)
                                      ggCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGtt                                 >  1:3898008/1‑71 (MQ=255)
                                       gCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTgg                                                            <  1:1927004/43‑1 (MQ=255)
                                         gAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTgg                              <  1:3158789/71‑1 (MQ=255)
                                         gAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTgg                              <  1:2218184/71‑1 (MQ=255)
                                           ggTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTgg                              >  1:3940366/1‑69 (MQ=255)
                                           ggTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTgg                              >  1:2193511/1‑69 (MQ=255)
                                               aTTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATg                                   <  1:1593969/60‑1 (MQ=255)
                                                ttCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGaaa                                              <  1:2318055/48‑1 (MQ=255)
                                                ttCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGaaa                                              <  1:3780950/48‑1 (MQ=255)
                                                ttCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGaaa                                              <  1:3567146/48‑1 (MQ=255)
                                                 tCTGGTCGTCAATATGCCCGTCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTa                            <  1:1007373/65‑1 (MQ=255)
                                                 tCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGt                       >  1:967738/1‑70 (MQ=255)
                                                  cTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGaaaaa                                            >  1:3539547/1‑48 (MQ=255)
                                                  cTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGaaaaa                                            >  1:319464/1‑48 (MQ=255)
                                                  cTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTTGATGTTTGGTAGCAGTa                      >  1:4021634/1‑70 (MQ=255)
                                                  cTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTa                      >  1:285368/1‑70 (MQ=255)
                                                  cTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTa                      >  1:2944807/1‑70 (MQ=255)
                                                  cTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTa                      >  1:1961425/1‑70 (MQ=255)
                                                  cTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTa                      <  1:1580136/70‑1 (MQ=255)
                                                   tGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTgg                              <  1:4028401/61‑1 (MQ=255)
                                                   tGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTgg                              <  1:2909205/61‑1 (MQ=255)
                                                     gTCGGCAATATGCGCGTCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTATAGt                  <  1:402873/71‑1 (MQ=255)
                                                     gTCGGCAATATGCCCGTCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGt                  <  1:2394691/71‑1 (MQ=255)
                                                     gTCGGCAATATGCCCGTCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGt                  <  1:892995/71‑1 (MQ=255)
                                                          cAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAg                                         <  1:745992/43‑1 (MQ=255)
                                                          cAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAg                                         <  1:214002/43‑1 (MQ=255)
                                                            atatGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGAt                                    <  1:1039796/46‑1 (MQ=255)
                                                            atatGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTa                      <  1:3961003/60‑1 (MQ=255)
                                                            atatGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTa                      <  1:3465279/60‑1 (MQ=255)
                                                            atatGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTa                      <  1:1901120/60‑1 (MQ=255)
                                                            atatGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTa                      <  1:3416616/60‑1 (MQ=255)
                                                             tatGCCCGTCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCgtg          <  1:2612312/71‑1 (MQ=255)
                                                             tatGCCCGTCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCgtg          <  1:3353903/71‑1 (MQ=255)
                                                             tatGCCCGCCAGAGTCATTGGTGTGGCGGAAGaaaaa                                            >  1:2965529/1‑37 (MQ=255)
                                                             tatGCCCGCCAGAGTCATTGGTGTGGCGGAAGaaaaa                                            >  1:2154788/1‑37 (MQ=255)
                                                             tatGCCCGCCAGAGTCATTGGTGTGGCGGAAGaaaaa                                            <  1:1454858/37‑1 (MQ=255)
                                                             tatGCCCGCCAGAGTCATTGGTGTGGCGGAAGaaaaa                                            >  1:1282479/1‑37 (MQ=255)
                                                             tatGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACa                                          <  1:2736339/39‑1 (MQ=255)
                                                             tatGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTa                      <  1:943677/59‑1 (MQ=255)
                                                              atGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGt                                  >  1:3519700/1‑46 (MQ=255)
                                                              atGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGt                                  >  1:2046895/1‑46 (MQ=255)
                                                              atGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTATAGTGCTGCgtgt         >  1:1015349/1‑71 (MQ=255)
                                                              atGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCgtgt         >  1:711258/1‑71 (MQ=255)
                                                              atGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCg            >  1:706242/1‑68 (MQ=255)
                                                               tggccGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTc        >  1:2747991/4‑71 (MQ=255)
                                                               tGCCCGTCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTc        <  1:2199108/71‑1 (MQ=255)
                                                               tGCCCGCCAGAGTCATTTGTGTGGCGGAAGAAAAACAGTCGa                                     >  1:871357/1‑42 (MQ=255)
                                                               tGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGa                                     >  1:2748895/1‑42 (MQ=255)
                                                                     ccAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTa  <  1:3045971/71‑1 (MQ=255)
                                                                     |                                                                      
CTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTA  >  W3110S.gb/921991‑922130

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