Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,359,716 T→C 39.7% E47E (GAA→GAG yehC ← predicted periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,359,7160TC39.7% 39.8 / 64.7 93E47E (GAA→GAGyehCpredicted periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base T (31/25);  new base C (19/18);  total (50/43)
Fisher's exact test for biased strand distribution p-value = 8.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

AATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCATGTGCAGGAAGGACCATAGAAAA  >  minE/1359650‑1359784
                                                                  |                                                                    
aaTCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAtttt                                                                     >  1:204522/1‑68 (MQ=255)
aaTCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAtttt                                                                     >  1:2879289/1‑68 (MQ=255)
aaTCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAtttt                                                                     >  1:1403941/1‑68 (MQ=255)
aaTCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTc                                                                    >  1:1169537/1‑69 (MQ=255)
aaTCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTc                                                                    >  1:244046/1‑69 (MQ=255)
aaTCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTc                                                                    >  1:90685/1‑69 (MQ=255)
aaTCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTc                                                                    >  1:698333/1‑69 (MQ=255)
aaTACACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTc                                                                    >  1:416988/1‑69 (MQ=255)
  tCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCt                                                                   <  1:1887216/68‑1 (MQ=255)
  tCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCt                                                                   <  1:1362492/68‑1 (MQ=255)
  tCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCt                                                                   <  1:2601634/68‑1 (MQ=255)
    cACGCCTGCAGCAGCTAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAttt                                                                      >  1:1932592/1‑63 (MQ=255)
    cACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAttt                                                                      >  1:1484451/1‑63 (MQ=255)
    cACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGcc                                                                >  1:3027096/1‑69 (MQ=255)
    cACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGcc                                                                >  1:2094495/1‑69 (MQ=255)
    cACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGcc                                                                >  1:2069513/1‑69 (MQ=255)
    cACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGcc                                                                >  1:2275211/1‑69 (MQ=255)
    cACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGcc                                                                >  1:2905006/1‑69 (MQ=255)
    cACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGcc                                                                >  1:1560037/1‑69 (MQ=255)
    cACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGcc                                                                >  1:2530418/1‑69 (MQ=255)
    cACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGcc                                                                >  1:985391/1‑69 (MQ=255)
                    aagaaCGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTGTGCTCTGCCGGGTAAATAATAcgcg                                                <  1:548162/69‑1 (MQ=255)
                    aagaaCGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTGTGCTCTGCCGGGTAAATAATAcgcg                                                <  1:3288472/69‑1 (MQ=255)
                    aagaaCGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTGTGCTCTGCCGGGTAAATAATAcgcg                                                <  1:578916/69‑1 (MQ=255)
                    aagaaCGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTGTGCTCTGCCGGGTAAATAATAcgcg                                                <  1:3092483/69‑1 (MQ=255)
                    aagaaCGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTGTGCTCTGCCGGGTAAATAATAcgcg                                                <  1:2122658/69‑1 (MQ=255)
                    aagaaCGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTGTGCTCTGCCGGGTAAATAATAcgcg                                                <  1:2655429/69‑1 (MQ=255)
                    aagaaCGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTGTGCTCTGCCGGGTAAATAATAcgcg                                                <  1:2898996/69‑1 (MQ=255)
                    aagaaCGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTGTGCTCTGCCGGGTAAATAATAcgcg                                                <  1:332925/69‑1 (MQ=255)
                    aagaaCGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTGTGCTCTGCCGGGTAAATAATAcgcg                                                <  1:911756/69‑1 (MQ=255)
                    aagaaCGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTGTGCTCTGCCGGGTAAATAATAcgcg                                                <  1:1924127/69‑1 (MQ=255)
                    aagaaCGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTGTGCTCTGCCGGGTAAATAATAcgcg                                                <  1:127003/69‑1 (MQ=255)
                    aagaaCGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTGTGCTCTGCCGGGTAAATAATAcgcg                                                <  1:131029/69‑1 (MQ=255)
                    aagaaCGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTGTGCTCTGCCGGGTAAATAATAcgcg                                                <  1:135001/69‑1 (MQ=255)
                     agaaCGGTTTGCCTGGTTCATCAACTGCACCATCACTTCTTTGTGCTCTGCCGGGTAAATAATAcgcg                                                <  1:1392413/68‑1 (MQ=255)
                     agaaCGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTGTGCTCTGCCGGGTAATTAATAcgcg                                                <  1:1989209/68‑1 (MQ=255)
                     agaaCGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTGTGCTCTGCCGGGTAAATAATAcgcg                                                <  1:775827/68‑1 (MQ=255)
                     agaaCGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTGTGCTCTGCCGGGTAAATAATAcgcg                                                <  1:889278/68‑1 (MQ=255)
                     agaaCGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTGTGCTCTGCCGGGTAAATAATAcgcg                                                <  1:536266/68‑1 (MQ=255)
                           gTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTa                                         <  1:2096407/69‑1 (MQ=255)
                           gTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTa                                         <  1:287325/69‑1 (MQ=255)
                           gTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTa                                         <  1:1099209/69‑1 (MQ=255)
                           gTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTa                                         <  1:1431142/69‑1 (MQ=255)
                           gTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTa                                         <  1:2517420/69‑1 (MQ=255)
                                  tGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAAc                                   <  1:195591/68‑1 (MQ=255)
                                  tGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAAc                                   <  1:1616914/68‑1 (MQ=255)
                                  tGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAAc                                   <  1:887343/68‑1 (MQ=255)
                                  tGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAAc                                   <  1:595272/68‑1 (MQ=255)
                                  tGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAAc                                   <  1:2970110/68‑1 (MQ=255)
                                      tcatcaACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAAc                                   >  1:452147/1‑64 (MQ=255)
                                      tcatcaACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAAc                                   >  1:741699/1‑64 (MQ=255)
                                          caACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCgg                           <  1:1375628/68‑1 (MQ=255)
                                          caACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCgg                           <  1:2416402/68‑1 (MQ=255)
                                           aaCTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCa                         >  1:2902060/1‑69 (MQ=255)
                                           aaCTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCa                         >  1:3127385/1‑69 (MQ=255)
                                           aaCTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCa                         >  1:549915/1‑69 (MQ=255)
                                           aaCTGCACAATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGc                          >  1:2518633/1‑68 (MQ=255)
                                           aaCTGCACATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCTTa                                         <  1:3108028/52‑1 (MQ=255)
                                               gCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCAtg                       >  1:2841305/1‑67 (MQ=255)
                                               gCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCAtg                       >  1:3345353/1‑67 (MQ=255)
                                               gCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCAtg                       >  1:2772161/1‑67 (MQ=255)
                                               gCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCAtg                       >  1:2419513/1‑67 (MQ=255)
                                               gCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCAtg                       >  1:678259/1‑67 (MQ=255)
                                               gCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCAtg                       >  1:3250080/1‑67 (MQ=255)
                                               gCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCAtg                       >  1:1451263/1‑67 (MQ=255)
                                                      cACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCATGTGCaggaag              <  1:3080128/69‑1 (MQ=255)
                                                       aCTTCTTTGTGCTCTGCCGGGTAAATAATACGCGTCCCGTAGATAACaa                                 >  1:2665076/1‑49 (MQ=37)
                                                       aCTTCTTTGTGCTCTGCCGGGTAAATAATACGCGTCCCGTAGATAACAATAcc                             >  1:35341/1‑53 (MQ=255)
                                                       aCTTCTTTGTGCTCTGCCGGGTAAATAATACGCGTCCCGTAGATAACAATAcc                             >  1:436999/1‑53 (MQ=255)
                                                       aCTTCTTTGTGCTCTGCCGGGTAAATAATACGCGTCCCGTAGATAACAATAcc                             >  1:986446/1‑53 (MQ=255)
                                                       aCTTCTTTGTGCTCTGCCGGGTAAATAATACGCGTCCCGTAGATAACAATAcc                             >  1:2709908/1‑53 (MQ=255)
                                                       aCTTCTTTGTGCTCTGCCGGGTAAATAATACGCGTCCCGTAGATAACAATAcc                             >  1:2272617/1‑53 (MQ=255)
                                                       aCTTCTTTGTGCTCTGCCGGGTAAATAATACGCGTCCCGTAGATAACAATAcc                             >  1:1822103/1‑53 (MQ=255)
                                                       aCTTCTTTGTGCTCTGCCGGGTAAATAATACGCGTCCCGTAGATAACAATAcc                             >  1:1785451/1‑53 (MQ=255)
                                                       aCTTCTTTGTGCTCTGCCGGGTAAATAATACGCGTCCCGTAGATAACAATAcc                             >  1:191225/1‑53 (MQ=255)
                                                       aCTTCTTTGTGCTCTGCCGGGTAAATAATACGCGTCCCGTAGATAACAATAcc                             >  1:1948694/1‑53 (MQ=255)
                                                       aCTTCTTTGTGCTCTGCCGGGTAAATAATACGCGTCCCGTAGATAACAATAcc                             >  1:2187782/1‑53 (MQ=255)
                                                       aCTTCTTTGTGCTCTGCCGGGTAAATAATACGCGTCCCGTAGATAACAATAcc                             >  1:2237236/1‑53 (MQ=255)
                                                       aCTTCTTTGTGCTCTGCCGGGTAAATAATACGCGTCCCGTAGATAACAATAcc                             >  1:2242891/1‑53 (MQ=255)
                                                       aCTTCTTTGTGCTCTGCCGGGTAAATAATACGCGTCCCGTAGATAACAATAcc                             >  1:2366736/1‑53 (MQ=255)
                                                       aCTTCTTTGTGCTCTGCCGGGTAAATAATACGCGTCCCGTAGATAACAATAcc                             >  1:2383504/1‑53 (MQ=255)
                                                       aCTTCTTTGTGCTCTGCCGGGTAAATAATACGCGTCCCGTAGATAACAATAcc                             >  1:2778086/1‑53 (MQ=255)
                                                       aCTTCTTTGTGCTCTGCCGGGTAAATAATACGCGTCCCGTAGATAACAATAcc                             >  1:2802822/1‑53 (MQ=255)
                                                       aCTTCTTTGTGCTCTGCCGGGTAAATAATACGCGTCCCGTAGATAACAATAcc                             >  1:2803284/1‑53 (MQ=255)
                                                       aCTTCTTTGTGCTCTGCCGGGTAAATAATACGCGTCCCGTAGATAACAATAcc                             >  1:2865533/1‑53 (MQ=255)
                                                       aCTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTa                                         <  1:91890/41‑1 (MQ=255)
                                                       aCTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCATGTGCaggaagg             >  1:1235336/1‑69 (MQ=255)
                                                       aCTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCATGTGCaggaagg             >  1:285979/1‑69 (MQ=255)
                                                       aCTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCATGTGCaggaagg             >  1:1342583/1‑69 (MQ=255)
                                                       aCTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCATGTGCaggaagg             >  1:2924633/1‑69 (MQ=255)
                                                         ttcttTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCATGTGCAGGAAGGAc           <  1:1107913/69‑1 (MQ=255)
                                                         ttcttTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCATGTGCAGGAAGGAc           <  1:1675215/69‑1 (MQ=255)
                                                         ttcttTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCATGTGCAGGAAGGAc           <  1:3220594/69‑1 (MQ=255)
                                                         ttcttTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCATGTGCAGGAAGGAc           <  1:75621/69‑1 (MQ=255)
                                                         ttcttTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCATGTGCAGGAAGGAc           <  1:1851866/69‑1 (MQ=255)
                                                         ttcttTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCATGTGCAGGAAGGAc           <  1:711363/69‑1 (MQ=255)
                                                         ttcttTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCATGTGCAGGAAGGAc           <  1:3354443/69‑1 (MQ=255)
                                                           cttTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCATGTGCAGGAAGGAc           >  1:659065/1‑67 (MQ=255)
                                                                  ttCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCATGTGCAGGAAGGACCATAGaaaa  <  1:1213525/69‑1 (MQ=255)
                                                                  ttCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCATGTGCAGGAAGGACCATAGaaaa  <  1:2975909/69‑1 (MQ=255)
                                                                  |                                                                    
AATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGTCCCGTAGATAACAATTCCGGCATGTGCAGGAAGGACCATAGAAAA  >  minE/1359650‑1359784

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: