Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,587,781 T→C 5.8% F511L (TTT→CTT)  acrD → aminoglycoside/multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,587,7810TC 5.8% 195.2 / 3.7 69F511L (TTT→CTT) acrDaminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base T (38/27);  new base C (2/2);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAA  >  W3110S.gb/2587712‑2587847
                                                                     |                                                                  
cTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtt                                                                   <  1:3151330/71‑1 (MQ=255)
   aCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCttttt                                                                >  1:708777/1‑71 (MQ=255)
   aCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCttttt                                                                >  1:113821/1‑71 (MQ=255)
   aCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtttt                                                                 <  1:1733117/70‑1 (MQ=255)
   aCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtttt                                                                 <  1:2045502/70‑1 (MQ=255)
   aCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtttt                                                                 <  1:3111402/70‑1 (MQ=255)
    cTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtttttt                                                               <  1:1905397/71‑1 (MQ=255)
     tCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTg                                                              >  1:3951154/1‑71 (MQ=255)
       cGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGcc                                                            >  1:4828931/1‑71 (MQ=255)
       cGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGcc                                                            <  1:5017457/71‑1 (MQ=255)
       cGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGc                                                             <  1:1798630/70‑1 (MQ=255)
        ggCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTATGCCt                                                           >  1:3924473/1‑71 (MQ=255)
         gCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTg                                                          >  1:3394237/1‑71 (MQ=255)
         gCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTg                                                          >  1:2269036/1‑71 (MQ=255)
          ctctGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCt                                                                    >  1:603672/1‑60 (MQ=255)
               tgtgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTa                                                     <  1:295056/70‑1 (MQ=255)
                 tgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTTGTTTAAc                                                   <  1:811102/70‑1 (MQ=255)
                 tgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAc                                                   >  1:1071659/1‑70 (MQ=255)
                  gCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAc                                                   >  1:4633337/1‑69 (MQ=255)
                   ccACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAg                                                <  1:4568354/71‑1 (MQ=255)
                     acacTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGAt                                              <  1:4305998/71‑1 (MQ=255)
                     acacTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGAt                                              >  1:4659448/1‑71 (MQ=255)
                     acacTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGAt                                              <  1:3558190/71‑1 (MQ=255)
                      cacTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATg                                             >  1:1033274/1‑71 (MQ=255)
                      cacTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATg                                             >  1:898521/1‑71 (MQ=255)
                      cacTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATg                                             <  1:755347/71‑1 (MQ=255)
                       actactTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGGTTTTTTGCCTGGTTTAACCAGATGt                                            >  1:2782449/1‑71 (MQ=255)
                       actactTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGt                                            >  1:4214645/1‑71 (MQ=255)
                       actactTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGt                                            >  1:4208733/1‑71 (MQ=255)
                        ctactTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtttttt                                                               <  1:1548640/51‑1 (MQ=255)
                        ctactTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCttttt                                                                <  1:2080586/50‑1 (MQ=255)
                         tactTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGtt                                           >  1:3141591/1‑70 (MQ=255)
                          actTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtt                                                                   <  1:2893112/45‑1 (MQ=255)
                           ctTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGttt                                          <  1:3516230/69‑1 (MQ=255)
                           ctTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGttt                                          <  1:3618210/69‑1 (MQ=255)
                           ctTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGttt                                          <  1:61790/69‑1 (MQ=255)
                              aaGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTg                                                          >  1:435640/1‑50 (MQ=255)
                              aaGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTg                                                          >  1:1815891/1‑50 (MQ=255)
                                gCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTa                                                     <  1:34189/53‑1 (MQ=255)
                                gCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGc                                    >  1:2792804/1‑70 (MQ=255)
                                  cACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAAc                                 >  1:3654260/1‑71 (MQ=255)
                                    cTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACGAGATGTTTAACCGCAACGc                               >  1:3757809/1‑71 (MQ=255)
                                     tGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCa                                                 >  1:715317/1‑52 (MQ=255)
                                      gAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGcc                              >  1:4225466/1‑70 (MQ=255)
                                       aaaaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCg                             <  1:4534051/70‑1 (MQ=255)
                                       aaaaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCg                             >  1:1760199/1‑70 (MQ=255)
                                                gAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACg                    >  1:1680047/1‑70 (MQ=255)
                                                gAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAAACGCAACGCCGAACGCTACg                    >  1:4418823/1‑70 (MQ=255)
                                                gAGCATCATGGGCAAAAAGGCTTTTTTGACTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACg                    >  1:4025989/1‑70 (MQ=255)
                                                 aGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGaa                  <  1:5009732/71‑1 (MQ=255)
                                                 aGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGa                   >  1:3752489/1‑70 (MQ=255)
                                                     tcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGttt                                          <  1:4908598/43‑1 (MQ=255)
                                                     tcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGttt                                          <  1:1694109/43‑1 (MQ=255)
                                                     tcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAg              >  1:691496/1‑71 (MQ=255)
                                                     tcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAg              <  1:2079975/71‑1 (MQ=255)
                                                      catGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATg                                             <  1:5051138/39‑1 (MQ=255)
                                                      catGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCg                             <  1:905833/55‑1 (MQ=255)
                                                       atGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgg             <  1:2062610/70‑1 (MQ=255)
                                                       atGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgg             <  1:598481/70‑1 (MQ=255)
                                                       atGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAg              >  1:2434722/1‑69 (MQ=255)
                                                         aggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGaaa                 <  1:3891644/63‑1 (MQ=38)
                                                         aggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGaaa                 <  1:4942688/63‑1 (MQ=38)
                                                             aaaaaaGCCTTTTTGCCTGGTTTAACCAGATGTTTAAcc                                      >  1:1258710/1‑39 (MQ=37)
                                                             aaaaaaGCCTTTTTGCCTGGTTTAACCAGATGTTTAAcc                                      >  1:3680009/1‑39 (MQ=37)
                                                             aaaaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGaa                           <  1:2145003/50‑1 (MQ=255)
                                                              aaaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACg                                <  1:3313626/44‑1 (MQ=255)
                                                              aaaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTg         >  1:1218310/1‑67 (MQ=255)
                                                              aaaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCg      >  1:1999193/1‑70 (MQ=255)
                                                               aaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACg                    >  1:3013108/1‑55 (MQ=255)
                                                               aaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTg         >  1:1708826/1‑66 (MQ=255)
                                                               aaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaa    <  1:2911888/71‑1 (MQ=255)
                                                               aaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaa    <  1:661183/71‑1 (MQ=255)
                                                               aaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGa     >  1:428385/1‑70 (MQ=255)
                                                                aaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCa                                   >  1:1792780/1‑39 (MQ=255)
                                                                aaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGcc                              >  1:3722287/1‑44 (MQ=255)
                                                                aaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaa    >  1:4091433/1‑70 (MQ=255)
                                                                 aGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGaaaaa               >  1:2253051/1‑58 (MQ=255)
                                                                 aGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaaaa  >  1:408980/1‑71 (MQ=255)
                                                                   gCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGTACGCTACGAAAAAGGGGTGGCg      <  1:4360091/65‑1 (MQ=255)
                                                                   gCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGaaa                 <  1:537365/54‑1 (MQ=255)
                                                                     |                                                                  
CTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAA  >  W3110S.gb/2587712‑2587847

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