Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 359,429 T→C 10.0% Y220H (TAT→CAT)  sfmD → predicted outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported359,4290TC10.0% 198.3 / 11.5 70Y220H (TAT→CAT) sfmDpredicted outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base T (36/27);  new base C (5/2);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCAT  >  Exported/359362‑359496
                                                                   |                                                                   
tttttttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTAtt                                                                  >  1:3290794/1‑71 (MQ=255)
  tttttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATAGTTCCAGGAACCAttg                                                                 <  1:1843514/70‑2 (MQ=38)
  tttttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATAGTTCCAGGAACCAttg                                                                 <  1:2501335/70‑2 (MQ=38)
    tttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTc                                                               >  1:1482658/1‑70 (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:690946/1‑71 (MQ=255)
                   gggATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCa                                               <  1:1835197/71‑1 (MQ=255)
                   gggATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCa                                               <  1:939532/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:1708034/71‑1 (MQ=255)
                                  ccGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaa                                >  1:1326734/1‑71 (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:1260541/1‑70 (MQ=255)
                                    gTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaat                               >  1:829994/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:2754335/71‑1 (MQ=255)
                                                  acaaTGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTgg                <  1:605177/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:1624151/36‑1 (MQ=255)
                                                         ttCCTGGAACTATTTTCGCGGAAATGGATATCATTc                                            <  1:3371800/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:3311146/1‑55 (MQ=255)
                                                          tCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTg                        >  1:3573846/1‑55 (MQ=255)
                                                          tCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTg                        >  1:2635623/1‑55 (MQ=255)
                                                          tCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTg                        >  1:763290/1‑55 (MQ=255)
                                                          tCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAgcg        >  1:1110694/1‑71 (MQ=255)
                                                          tCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAgcg        >  1:873384/1‑71 (MQ=255)
                                                          tCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAgcg        >  1:2561081/1‑71 (MQ=255)
                                                          tCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAgcg        >  1:647062/1‑71 (MQ=255)
                                                          tCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTAGCACCTGGGTACAgcg        >  1:3144641/1‑71 (MQ=255)
                                                           ccTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGAGGAATAATATTGGCACCTGGGTACAgcgc       >  1:158811/1‑71 (MQ=255)
                                                               gAAATATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaataat                            >  1:252287/1‑46 (MQ=255)
                                                                aaCTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCat  <  1:2529697/71‑1 (MQ=255)
                                                                  cTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATAtt                         >  1:1503105/1‑46 (MQ=255)
                                                                   tATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGc                      <  1:3385251/48‑1 (MQ=255)
                                                                   |                                                                   
TTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCAT  >  Exported/359362‑359496

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: