Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 359,417 G→A 13.5% G216S (GGT→AGT)  sfmD → predicted outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE359,4170GA13.5% 161.5 / 14.9 59G216S (GGT→AGT) sfmDpredicted outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base G (26/25);  new base A (6/2);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

GTAACGATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGG  >  minE/359349‑359482
                                                                    |                                                                 
gTAACGATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATgg                                                                  <  1:65021/70‑1 (MQ=255)
gTAACGATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGt                                                                 <  1:3421123/71‑1 (MQ=255)
   aCGATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTc                                                               <  1:321655/70‑1 (MQ=255)
    cGATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCt                                                             >  1:939703/1‑71 (MQ=255)
    cGATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCt                                                             <  1:47867/71‑1 (MQ=255)
     gATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTg                                                            >  1:3454605/1‑71 (MQ=255)
      aTAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGATCCTg                                                            >  1:79117/1‑70 (MQ=255)
             tttttttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTAtt                                                    >  1:3290794/1‑71 (MQ=255)
               tttttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATAGTTCCAGGAACCAttg                                                   <  1:2501335/70‑2 (MQ=38)
               tttttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATAGTTCCAGGAACCAttg                                                   <  1:1843514/70‑2 (MQ=38)
                 tttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTgg                                                           >  1:2324704/1‑60 (MQ=255)
                 tttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTc                                                 >  1:1482658/1‑70 (MQ=255)
                 tttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATAGTTCCAgg                                                           >  1:282727/1‑60 (MQ=255)
                 tttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATAGTTCCAGGAACCAttg                                                   >  1:2527780/1‑67 (MQ=38)
                 tttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATAGTTCCAGGAACCAttg                                                   >  1:265835/1‑67 (MQ=38)
                  ttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTcg                                                <  1:1751062/70‑1 (MQ=255)
                  ttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTcg                                                <  1:3017912/70‑1 (MQ=255)
                        aGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGaaa                                          <  1:3528974/70‑1 (MQ=255)
                        aGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAAt                                         <  1:1328461/71‑1 (MQ=255)
                         gCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATg                                        <  1:1385374/71‑1 (MQ=255)
                         gCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATg                                        >  1:2060088/1‑71 (MQ=255)
                              gCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATg                                        <  1:3262091/66‑1 (MQ=255)
                               cGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGatat                                   >  1:3182397/1‑70 (MQ=255)
                               cGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGatat                                   >  1:2430152/1‑70 (MQ=255)
                                gggATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCa                                 <  1:939532/71‑1 (MQ=255)
                                gggATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCa                                 >  1:690946/1‑71 (MQ=255)
                                gggATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCa                                 >  1:1416540/1‑71 (MQ=255)
                                gggATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCa                                 <  1:1835197/71‑1 (MQ=255)
                                   aTTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTc                              <  1:2561643/71‑1 (MQ=255)
                                      aaTATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAg                            <  1:1291712/70‑1 (MQ=255)
                                      aaTATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAg                            >  1:1800028/1‑70 (MQ=255)
                                         atTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGaa                          >  1:3283148/1‑69 (MQ=255)
                                          ttGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTc                                                 >  1:3544505/1‑45 (MQ=255)
                                          ttGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAAc                         <  1:1600020/69‑1 (MQ=255)
                                          ttGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAg                       >  1:3273700/1‑71 (MQ=255)
                                              cccGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGa                   <  1:3145295/71‑1 (MQ=255)
                                               ccGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaa                  >  1:1326734/1‑71 (MQ=255)
                                               ccGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaa                  <  1:1708034/71‑1 (MQ=255)
                                                cGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaat                 <  1:784800/71‑1 (MQ=255)
                                                cGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaat                 <  1:772477/71‑1 (MQ=255)
                                                 gTGGCGTTTACGCAACAATGGTTCTTGGAACTAtttt                                                  >  1:1514068/1‑37 (MQ=255)
                                                 gTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaat                 >  1:1404443/1‑70 (MQ=255)
                                                 gTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaat                 >  1:829994/1‑70 (MQ=255)
                                                 gTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaat                 >  1:1260541/1‑70 (MQ=255)
                                                 gTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaat                 >  1:2534419/1‑70 (MQ=255)
                                                 gTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaa                  >  1:287779/1‑69 (MQ=255)
                                                  tGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTcgcg                                              <  1:415405/40‑1 (MQ=255)
                                                  tGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaataa               <  1:3075361/71‑1 (MQ=255)
                                                    gCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATTGATATCATTCAGAACAGTTGAataata             <  1:3328808/71‑1 (MQ=255)
                                                    gCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaataat              <  1:349478/70‑1 (MQ=255)
                                                    gCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaataat              <  1:729915/70‑1 (MQ=255)
                                                     catttCCGCGAAAATAGTTCCAGGAACCATTGTTGCGGAAATGGATATCATTCAGAACAGTGGAATAatat            >  1:2465091/3‑71 (MQ=9)
                                                         tACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATAtt           >  1:1559657/1‑68 (MQ=255)
                                                          aCGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCa       >  1:1140651/1‑71 (MQ=255)
                                                              aacaaTGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTg   >  1:2494508/1‑71 (MQ=255)
                                                               acaaTGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTgg  <  1:605177/71‑1 (MQ=255)
                                                               acaaTGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTgg  <  1:2754335/71‑1 (MQ=255)
                                                                aaaTAGTTCCAGGAACCATTGTTGCGTAAATGGATATCATTCAGAACAGTGGaataat              >  1:3106890/2‑58 (MQ=17)
                                                                aaaTAGTTCCAGGAACCATTGTTGCGTAAATGGATATCATTCAGAACAGTGGaataat              >  1:1187370/2‑58 (MQ=17)
                                                                 aaTGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGa                   >  1:28939/1‑52 (MQ=255)
                                                                    |                                                                 
GTAACGATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGG  >  minE/359349‑359482

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: