Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 558,960 A→G 12.2% Q69Q (CAA→CAG cmr → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE558,9600AG12.2% 99.7 / 5.8 41Q69Q (CAA→CAGcmrmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (17/19);  new base G (2/3);  total (19/22)
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.41e-01

ATATCAGGCGGGCATTGATTGGGTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAG  >  minE/558891‑559018
                                                                     |                                                          
atatCAGGCGGGCATTGATTGGGTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAAt                                                           <  1:3480977/71‑1 (MQ=255)
 tatCAGGCGGGCATTGATTGGGTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATg                                                          <  1:770611/71‑1 (MQ=255)
 tatCAGGCGGGCATTGATTGGGTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATg                                                          <  1:5421174/71‑1 (MQ=255)
   tCAGGCGGGCATTGATTGGGTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAAt                                                           >  1:3938526/1‑68 (MQ=255)
                gattgGGTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCt                                           >  1:3735237/1‑71 (MQ=255)
                gattgGGTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCt                                           >  1:3099087/1‑71 (MQ=255)
                gattgGGTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCt                                           >  1:1613794/1‑71 (MQ=255)
                 attgGGTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCt                                           <  1:5664165/70‑1 (MQ=255)
                  ttgGGTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGt                                         >  1:3592948/1‑71 (MQ=255)
                    gggTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGc                                            >  1:2431449/1‑66 (MQ=255)
                     ggTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCg                                       >  1:4834167/1‑70 (MQ=255)
                      gTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGa                                     <  1:5145216/71‑1 (MQ=255)
                        tCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTgg                                                  >  1:4855454/1‑56 (MQ=255)
                        tCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTgg                                                  >  1:5493626/1‑56 (MQ=255)
                          cTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGTATCGt                                 <  1:2509527/71‑1 (MQ=255)
                          cTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCg                                  <  1:3370055/70‑1 (MQ=255)
                          cTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCg                                  <  1:4622537/70‑1 (MQ=255)
                          cTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGt                                 <  1:1168012/71‑1 (MQ=255)
                                cGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTa                                >  1:2171554/1‑66 (MQ=255)
                                cGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTg                             <  1:2587475/69‑1 (MQ=255)
                                cGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTg                             <  1:4481084/69‑1 (MQ=255)
                                cGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGt                           >  1:648688/1‑71 (MQ=255)
                                cGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGt                           >  1:983826/1‑71 (MQ=255)
                                     aCCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGcc                       <  1:4582602/70‑1 (MQ=255)
                                       cgcgTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTc                    >  1:2157784/1‑71 (MQ=255)
                                         tctATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAATTGGCCGCAGGCCg                  >  1:1792202/3‑71 (MQ=255)
                                         cgTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTcc                   <  1:2524799/70‑1 (MQ=255)
                                          gTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCgg                 <  1:3187092/71‑1 (MQ=255)
                                          gTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCgg                 <  1:5559143/71‑1 (MQ=255)
                                             tCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAATTGGCCg                         <  1:4474671/60‑1 (MQ=255)
                                             tCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAATTGGCCGCAGGCCgg                 <  1:3590787/68‑1 (MQ=255)
                                             tCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAATTGGCCGCAGGCCgg                 <  1:2676914/68‑1 (MQ=255)
                                                     gcgggATGTTTTTACAATGGCTGCTGGGGCCGCTGTCgg                                      <  1:5157700/39‑1 (MQ=255)
                                                         gATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAg  <  1:4317998/71‑1 (MQ=255)
                                                            gTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCg                         >  1:3151194/1‑45 (MQ=255)
                                                            gTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCg                         >  1:367359/1‑45 (MQ=255)
                                                             tttttACAATGGCTGCTGGGGCCGCTGTCGGATCg                                  <  1:2850000/35‑1 (MQ=255)
                                                             tttttACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTcc                   <  1:1351681/50‑1 (MQ=255)
                                                               tttACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTg               >  1:2286765/1‑52 (MQ=255)
                                                               tttACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGa              >  1:776816/1‑53 (MQ=255)
                                                                  aCAGTGGCTGCTGGGGCCGCTTTCCGACAGAATTGGCCGCAGGCCg                  >  1:3574839/1‑46 (MQ=255)
                                                                     |                                                          
ATATCAGGCGGGCATTGATTGGGTTCCTACTTCGATGACCGCGTATCTGGCGGGCGGGATGTTTTTACAATGGCTGCTGGGGCCGCTGTCGGATCGTATTGGTCGCCGTCCGGTGATGCTGGCGGGAG  >  minE/558891‑559018

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: