Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,415,880 A→T 7.0% P331P (CCA→CCT acrF → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,415,8800AT 7.0% 228.0 / 8.8 85P331P (CCA→CCTacrFmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (47/32);  new base T (2/4);  total (49/36)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

AAACTGGCGGAATTACAGCCATTCTTCCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTATGCTGGTGTTCCTGGTG  >  W3110S.gb/3415812‑3415949
                                                                    |                                                                     
aaaCTGGCGGAATTACAGCCATTCTTCCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCAtt                                                                     >  1:3328075/1‑71 (MQ=255)
 aaCTGGCGGAATTACAGCCATTCTTCCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTc                                                                    >  1:1558527/1‑71 (MQ=255)
  aCTGGCGGAATTACAGCCATTCTTCCCGCAGGGAATTAAGGTTCTCTACCCTTATGACACCACGCCATTCg                                                                   >  1:32242/1‑71 (MQ=255)
  aCTGGCGGAATTACAGCCATTCTTCCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCg                                                                   >  1:4002204/1‑71 (MQ=255)
  aCTGGCGGAATTACAGCCATTCTTCCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCg                                                                   >  1:3450605/1‑71 (MQ=255)
   cTGGCGGAATTACAGCCATTCTTCCCGCAGGGAATGAAGGTTCTTTACCCTTATGACACCACGCCTTTCGt                                                                  <  1:1291274/71‑1 (MQ=255)
   cTGGCGGAATTACAGCCATTCTTCCCGCAGGGAATGAAGGTTCTTTACCCTTATGACACCACGCCTTTCGt                                                                  <  1:3363154/71‑1 (MQ=255)
   cTGGCGGAATTACAGCCATTCTTCCCGCAGGGAATGAAGGTTCTTTACCCTTATGACACCACGCCTTTCGt                                                                  <  1:3184732/71‑1 (MQ=255)
       cGGAATTACAGCCATTCTTCCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAg                                                              <  1:3818479/71‑1 (MQ=255)
        ggAATTACAGCCATTCTTCCCGCAGGGAATGAAGGTTCTTTACCCTTATGACACCACGCCTTTCGt                                                                  <  1:252198/66‑1 (MQ=255)
         gAATTACAGCCATTCTTCCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCt                                                            >  1:1031002/1‑71 (MQ=255)
         gAATTACAGCCATTCTTCCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCt                                                            <  1:3922977/71‑1 (MQ=255)
         gAATTACACCCATTCTTCCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCt                                                            >  1:222664/1‑71 (MQ=255)
           aTTACAGCCATTCTTCCCGCAGGGAATGAAGGTTCTCTACCCTTATGACAGCACGCCATTCGTCCAGCttt                                                          >  1:3565234/1‑71 (MQ=255)
           aTTACAGCCATTCTTCCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCttt                                                          >  1:27047/1‑71 (MQ=255)
           aTTACAGCCATTCTTCCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCttt                                                          >  1:46904/1‑71 (MQ=255)
            ttACAGCCATTCTTCCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTc                                                         >  1:3537065/1‑71 (MQ=255)
            ttACAGCCATTCTTCCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTc                                                         >  1:221881/1‑71 (MQ=255)
             tACAGCCATTCTTCCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTc                                                         >  1:3896638/1‑70 (MQ=255)
              aCAGCCATTCTTCCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTa                                                       >  1:3376114/1‑71 (MQ=255)
               cAGCCATTCTTCCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTa                                                       <  1:2523770/70‑1 (MQ=255)
                aGCCATTCTTCCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTAtt                                                     >  1:2737829/1‑71 (MQ=255)
                    aTTCTTCCCGCAGGGAATGAAGGTTCTTTACCCTTATGACACCACGCCTTTCGTCCAGc                                                             >  1:4066446/1‑59 (MQ=255)
                       cttcCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGc                                                             >  1:2246471/1‑56 (MQ=255)
                       cttcCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAg                                              >  1:79575/1‑71 (MQ=255)
                       cttcCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAg                                              >  1:2936774/1‑71 (MQ=255)
                       cttcCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAg                                              >  1:1707132/1‑71 (MQ=255)
                        ttcCCGGAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACg                                                 >  1:1370101/1‑67 (MQ=255)
                          cccGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGt                                             <  1:33680/69‑1 (MQ=255)
                          cccGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTg                                            <  1:1495425/70‑1 (MQ=255)
                            cGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGt                                          >  1:3541505/1‑70 (MQ=255)
                              cAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTaaa                                       <  1:3402330/71‑1 (MQ=255)
                               aGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTaaa                                       <  1:1238784/70‑1 (MQ=255)
                                   aaTGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTaaaaa                                     >  1:2132862/1‑68 (MQ=255)
                                   aaTGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTaaaaa                                     >  1:3515428/1‑68 (MQ=255)
                                     tGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTc                                                         <  1:2139963/46‑1 (MQ=255)
                                     tGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTc                                                         <  1:3368633/46‑1 (MQ=255)
                                     tGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGt                                             <  1:190565/58‑1 (MQ=255)
                                       aaGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCttt                                                          >  1:3910781/1‑43 (MQ=255)
                                       aaGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGt                                             >  1:4280633/1‑56 (MQ=255)
                                       aaGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGt                                             >  1:4218249/1‑56 (MQ=255)
                                       aaGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGtt                              >  1:4284497/1‑71 (MQ=255)
                                       aaGGTTCTCTACCCTTATGACACCACGCCATACGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGtt                              >  1:146667/1‑71 (MQ=255)
                                         ggTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCg                            <  1:1941964/71‑1 (MQ=255)
                                         ggTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCg                            <  1:495293/71‑1 (MQ=255)
                                         ggTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCg                            <  1:256324/71‑1 (MQ=255)
                                          gTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACg                                                 <  1:3552915/49‑1 (MQ=255)
                                          gTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACg                                                 <  1:3397562/49‑1 (MQ=255)
                                            tctcTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTg                                            >  1:189749/1‑52 (MQ=255)
                                             ctttACCCTTATGACACCACGCCTTTCGTCCAGCTTTCTATCCAGGAAGTGGTGAAAACGCTGTTCGAAGc                        >  1:1005542/4‑71 (MQ=21)
                                             ctctACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGt                                             <  1:3360699/50‑1 (MQ=255)
                                              tctACCCTTATGACACCACGCCATTCGTCCAGCTTTCTa                                                       >  1:2486125/1‑39 (MQ=255)
                                              tctACCCTTATGACACCACGCCATTCGTCCAGCTTTCTa                                                       >  1:3647204/1‑39 (MQ=255)
                                              tctACCCTTATGACACACGCCATTCGTCCAGCTTTCTa                                                       >  1:429580/1‑38 (MQ=38)
                                               ctACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGtt                              >  1:1352145/1‑63 (MQ=255)
                                               ctACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGtt                              >  1:2436003/1‑63 (MQ=255)
                                               ctACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCa                      >  1:2669585/1‑71 (MQ=255)
                                                tACCCTTATGACACCACGCCATTCGTCCAGCTTTc                                                         <  1:923944/35‑1 (MQ=255)
                                                  cccTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCa                      <  1:3951127/68‑1 (MQ=255)
                                                  cccTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTa                   >  1:3803666/1‑71 (MQ=255)
                                                   ccTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCAt                     <  1:1997156/68‑1 (MQ=255)
                                                    cTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTg                                <  1:4161939/56‑1 (MQ=255)
                                                     ttATGACACCACGCCATTCGTCCAGCTTTCTATTCACGa                                                <  1:4343462/39‑1 (MQ=255)
                                                     ttATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTg                                >  1:2684529/1‑55 (MQ=255)
                                                      tATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTa                   >  1:4243694/1‑67 (MQ=255)
                                                      tATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTATGCt               >  1:3258488/1‑71 (MQ=255)
                                                      tATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCACTa                   >  1:3012146/1‑67 (MQ=255)
                                                      tATGACACAACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTa                   >  1:281696/1‑67 (MQ=255)
                                                       aTGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTa                   >  1:3921516/1‑66 (MQ=255)
                                                       aTGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTATGc                <  1:2322665/69‑1 (MQ=255)
                                                       aTGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTATGc                <  1:2297027/69‑1 (MQ=255)
                                                       aTGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTATGc                <  1:2133089/69‑1 (MQ=255)
                                                        tGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTATGCTgg             <  1:825074/71‑1 (MQ=255)
                                                         gACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTATGCTgg             >  1:4318840/1‑70 (MQ=255)
                                                         gACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTATGCTGgt            >  1:420004/1‑71 (MQ=255)
                                                          acacCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTATGCTGgtg           <  1:2345050/71‑1 (MQ=255)
                                                           caccacGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTATGCTGgtgt          <  1:1492256/71‑1 (MQ=255)
                                                              cacGCCATTCGTCCAGCTTTCTATTAACGAAGTGGTAAAAACGCTGTTCGAAGCCATTATGCTGGTGTTcc       >  1:3616489/1‑71 (MQ=255)
                                                                cGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTa                   >  1:1531778/1‑57 (MQ=255)
                                                                cGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTATGCTGGTGTTCCTg     >  1:4410774/1‑71 (MQ=255)
                                                                  ccATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTATGCTGGTGTTCCTGGt   <  1:2056033/71‑1 (MQ=255)
                                                                  ccATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTATGCTGGTGTTCCTGGt   <  1:1774983/71‑1 (MQ=255)
                                                                  ccATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTATGCTGGTGTTCCTGGt   <  1:3250232/71‑1 (MQ=255)
                                                                   cATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTATGCTGGTGTTCCTGGTg  <  1:840143/71‑1 (MQ=255)
                                                                   cATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTATGCTCGTGTTCCTGGTg  <  1:813956/71‑1 (MQ=255)
                                                                    aTTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCAt                     <  1:4071159/51‑1 (MQ=255)
                                                                    |                                                                     
AAACTGGCGGAATTACAGCCATTCTTCCCGCAGGGAATGAAGGTTCTCTACCCTTATGACACCACGCCATTCGTCCAGCTTTCTATTCACGAAGTGGTAAAAACGCTGTTCGAAGCCATTATGCTGGTGTTCCTGGTG  >  W3110S.gb/3415812‑3415949

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: