Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 359,423 T→A 11.1% W218R (TGG→AGG)  sfmD → predicted outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported359,4230TA11.1% 198.6 / 13.2 71W218R (TGG→AGG) sfmDpredicted outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base T (37/26);  new base A (6/2);  total (43/28)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGC  >  Exported/359353‑359491
                                                                      |                                                                    
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:1843514/70‑2 (MQ=38)
           tttttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATAGTTCCAGGAACCAttg                                                            <  1:2501335/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:265835/1‑67 (MQ=38)
             tttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATAGTTCCAGGAACCAttg                                                            >  1:2527780/1‑67 (MQ=38)
              ttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTcg                                                         <  1:3017912/70‑1 (MQ=255)
              ttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTcg                                                         <  1:1751062/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:690946/1‑71 (MQ=255)
                            gggATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCa                                          <  1:939532/71‑1 (MQ=255)
                            gggATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCa                                          <  1:1835197/71‑1 (MQ=255)
                            gggATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCa                                          >  1:1416540/1‑71 (MQ=255)
                               aTTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTc                                       <  1:2561643/71‑1 (MQ=255)
                                  aaTATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAg                                     >  1:1800028/1‑70 (MQ=255)
                                  aaTATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAg                                     <  1:1291712/70‑1 (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:772477/71‑1 (MQ=255)
                                            cGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaat                          <  1:784800/71‑1 (MQ=255)
                                             gTGGCGTTTACGCAACAATGGTTCTTGGAACTAtttt                                                           >  1:1514068/1‑37 (MQ=255)
                                             gTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaat                          >  1:2534419/1‑70 (MQ=255)
                                             gTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaat                          >  1:1260541/1‑70 (MQ=255)
                                             gTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaat                          >  1:829994/1‑70 (MQ=255)
                                             gTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaat                          >  1:1404443/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:729915/70‑1 (MQ=255)
                                                gCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaataat                       <  1:349478/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)
                                                                  ttCCTGGAACTATTTTCGCGGAAATGGATATCATTc                                       <  1:3371800/36‑1 (MQ=255)
                                                                  ttCCTGGAACTATTTTCGCGGAAATGGATATCATTc                                       <  1:1624151/36‑1 (MQ=255)
                                                                  ttCCTGGAACTATTTTCGCGGAAATGGATATCATTc                                       <  1:1317486/36‑1 (MQ=255)
                                                                  ttCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGATCAGTGGAATAATATTGGCACCTGGGTACa      >  1:3335786/1‑69 (MQ=255)
                                                                  ttCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAgc    >  1:3382831/1‑71 (MQ=255)
                                                                   tCCTGGAACTATTTTCGCGGAAATGGATATCATTc                                       <  1:754257/35‑1 (MQ=255)
                                                                   tCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTg                   >  1:2635623/1‑55 (MQ=255)
                                                                   tCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTg                   >  1:3573846/1‑55 (MQ=255)
                                                                   tCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTg                   >  1:3311146/1‑55 (MQ=255)
                                                                   tCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTg                   >  1:763290/1‑55 (MQ=255)
                                                                   tCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAgcg   >  1:873384/1‑71 (MQ=255)
                                                                   tCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAgcg   >  1:1110694/1‑71 (MQ=255)
                                                                   tCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAgcg   >  1:647062/1‑71 (MQ=255)
                                                                   tCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAgcg   >  1:2561081/1‑71 (MQ=255)
                                                                   tCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTAGCACCTGGGTACAgcg   >  1:3144641/1‑71 (MQ=255)
                                                                    ccTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGAGGAATAATATTGGCACCTGGGTACAgcgc  >  1:158811/1‑71 (MQ=255)
                                                                      |                                                                    
CGATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGC  >  Exported/359353‑359491

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: