Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 2,589,664 A→T intergenic (+300/‑239) acrD → / → yffB aminoglycoside/multidrug efflux system/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,589,6640AT93.3% 249.9 / 7.1 89intergenic (+300/‑239)acrD/yffBaminoglycoside/multidrug efflux system/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (5/1);  new base T (18/65);  total (23/66)
Fisher's exact test for biased strand distribution p-value = 4.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTACCTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGTTAA  >  W3110S.gb/2589602‑2589723
                                                              |                                                           
aCTACCTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAg                                                      >  1:1074056/1‑70 (MQ=255)
aCTACCTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAg                                                      >  1:2889731/1‑70 (MQ=255)
aCTACCTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAg                                                      >  1:3547726/1‑70 (MQ=255)
    ccTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCaa                                                   <  1:2028614/69‑1 (MQ=255)
    ccTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACa                                                 <  1:2096951/71‑1 (MQ=255)
    ccTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACa                                                 <  1:2101592/71‑1 (MQ=255)
    ccTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACa                                                 <  1:2678581/71‑1 (MQ=255)
        cGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACAtt                                               <  1:875089/69‑1 (MQ=255)
        cGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACAtt                                               <  1:1784917/69‑1 (MQ=255)
        cGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACAtt                                               <  1:2013498/69‑1 (MQ=255)
        cGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACAtt                                               <  1:3971296/69‑1 (MQ=255)
        cGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACAtt                                               <  1:2361260/69‑1 (MQ=255)
          ggggACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCga                                           <  1:812764/71‑1 (MQ=255)
           gggACAAATCTTATGTAAATACTATTGTCCTGCGGTGATGATTTGTATGTGTTACACAGCaa                                                   <  1:1832177/62‑1 (MQ=255)
           gggACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCgag                                          <  1:1140426/71‑1 (MQ=255)
           gggACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCgag                                          <  1:2145034/71‑1 (MQ=255)
           gggACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCgag                                          <  1:1174221/71‑1 (MQ=255)
                          tAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCa                            <  1:2641511/70‑1 (MQ=255)
                          tAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCa                            <  1:3730581/70‑1 (MQ=255)
                          tAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCa                            <  1:2142617/70‑1 (MQ=255)
                          tAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCa                            <  1:3899074/70‑1 (MQ=255)
                          tAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCa                            <  1:429560/70‑1 (MQ=255)
                          tAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCa                            <  1:472819/70‑1 (MQ=255)
                          tAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCa                            <  1:557413/70‑1 (MQ=255)
                          tAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCa                            <  1:1464272/70‑1 (MQ=255)
                          tAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCa                            <  1:1676093/70‑1 (MQ=255)
                          tAAATACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCa                            <  1:1539797/70‑1 (MQ=255)
                            aaTACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCt                         <  1:1725117/71‑1 (MQ=255)
                            aaTACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCt                         <  1:2740124/71‑1 (MQ=255)
                            aaTACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCt                         <  1:1436743/71‑1 (MQ=255)
                            aaTACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCt                         <  1:1023495/71‑1 (MQ=255)
                            aaTACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCt                         <  1:3008810/71‑1 (MQ=255)
                            aaTACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCt                         <  1:2453334/71‑1 (MQ=255)
                            aaTACTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCt                         <  1:4057757/71‑1 (MQ=255)
                                cTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAAt                      >  1:2318393/1‑70 (MQ=255)
                                cTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAAt                      >  1:788835/1‑70 (MQ=255)
                                cTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATa                     >  1:4020187/1‑71 (MQ=255)
                                cTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATa                     >  1:3861197/1‑71 (MQ=255)
                                cTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATa                     >  1:3798896/1‑71 (MQ=255)
                                cTATGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATa                     >  1:2124974/1‑71 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATTTGTTACACAGCAACATTTCGAGATAtt                                     <  1:2031483/52‑1 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGTAACATTTCGAGATAtt                                     <  1:1395581/52‑1 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATAtt                                     <  1:1263565/52‑1 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATAtt                                     <  1:516799/52‑1 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATAtt                                     <  1:3988014/52‑1 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATAtt                                     <  1:1161024/52‑1 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATAtt                                     <  1:993416/52‑1 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATAtt                                     <  1:1908240/52‑1 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATAtt                                     <  1:2644037/52‑1 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATAtt                                     <  1:2269302/52‑1 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCa                                   <  1:2885260/54‑1 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCa                                   <  1:2951925/54‑1 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCa                                   <  1:2521457/54‑1 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCa                                   <  1:164602/54‑1 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCa                                   <  1:3308624/54‑1 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCa                                   <  1:1548108/54‑1 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCa                                   <  1:2166457/54‑1 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCa                                   <  1:1341126/54‑1 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCa                                   <  1:148903/54‑1 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATc                          >  1:1144179/1‑63 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATc                          >  1:728614/1‑63 (MQ=255)
                                   tGGTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATc                          >  1:3859182/1‑63 (MQ=255)
                                    ggTCCTGCGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCa                                   <  1:2660148/53‑1 (MQ=255)
                                      tCCTGCGGTGATGATTTGTATGTGTTACACAGCAACa                                                 <  1:3670704/37‑1 (MQ=255)
                                        cTGCGGTGAGGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCa                                   <  1:1353582/49‑1 (MQ=255)
                                           cGGTGATGATTTGTATGTGTTACACAGCAACATTTCGGGATATTCATACGGCATCTAATACTGATTTAAtt          >  1:3975865/1‑71 (MQ=255)
                                           cGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTaa            <  1:1423658/69‑1 (MQ=255)
                                           cGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTaa            <  1:1228443/69‑1 (MQ=255)
                                           cGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTaa            <  1:2152268/69‑1 (MQ=255)
                                           cGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTaa            <  1:3480684/69‑1 (MQ=255)
                                           cGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTaa            <  1:922795/69‑1 (MQ=255)
                                           cGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTaa            <  1:3539850/69‑1 (MQ=255)
                                           cGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAAtt          >  1:3354966/1‑71 (MQ=255)
                                           cGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAAtt          >  1:3320461/1‑71 (MQ=255)
                                           cGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAAt           >  1:2317019/1‑70 (MQ=255)
                                           cGGTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATACTCATACGGCATCTAATACTGATTTAAt           >  1:3197095/1‑70 (MQ=255)
                                            ggTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTaa            <  1:313704/68‑1 (MQ=255)
                                             gTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATAtt                                     <  1:1718435/42‑1 (MQ=255)
                                             gTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATAtt                                     <  1:12338/42‑1 (MQ=255)
                                             gTGATGATTTGTATGTGTTACACAGCAACATTTCGAGATAtt                                     <  1:581127/42‑1 (MQ=255)
                                             gTGATGATTTGTATGTGTTACACAGCAACATATCGAGATAtt                                     <  1:1723391/42‑1 (MQ=38)
                                                  gatTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGtt    >  1:3323035/1‑70 (MQ=255)
                                                  gatTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGtt    >  1:2267172/1‑70 (MQ=255)
                                                  gatTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGtt    >  1:3379086/1‑70 (MQ=255)
                                                  gatTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGtt    >  1:2712701/1‑70 (MQ=255)
                                                  gatTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGtt    >  1:3599937/1‑70 (MQ=255)
                                                   atTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGTTaa  <  1:956846/71‑1 (MQ=255)
                                                   atTTGTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGTTa   >  1:576413/1‑70 (MQ=255)
                                                       gTATGTGTTACACAGCAACATTTCGAGATATTCATACGGCa                            <  1:2360919/41‑1 (MQ=255)
                                                              |                                                           
ACTACCTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGGTGATGATTTGTATGTGATACACAGCAACATTTCGAGATATTCATACGGCATCTAATACTGATTTAATTCTGGTTAA  >  W3110S.gb/2589602‑2589723

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