Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,345,321 C→T 36.8% G177S (GGT→AGT)  wza ← lipoprotein required for capsular polysaccharide translocation through the outer membrane

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,345,3210CT36.8% 19.9 / 20.0 38G177S (GGT→AGT) wzalipoprotein required for capsular polysaccharide translocation through the outer membrane
Reads supporting (aligned to +/- strand):  ref base C (7/17);  new base T (6/8);  total (13/25)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

GCGTCGAGAATAGTCAGTGGCACGTTGGTGATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGATATAGGCCTTTTGTGAGCGGAAGGCGGCGATATTAACGTCCACCTGCGGGTCAGCGATG  >  minE/1345260‑1345385
                                                             |                                                                
gCGTCGAGAATGGTCATAGGCACGTTGGTAATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGata                                                           <  1:1877423/69‑1 (MQ=255)
gCGTCGAGAATGGTCAGAGGCACGTTGGTAATCGCCTGTTGACCGGATTTATTCACCTGACCTGAGatt                                                           <  1:6779/69‑2 (MQ=255)
gCGTCGAGAATGGTCAGAGGCACGTTGGTAATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGata                                                           <  1:2395049/69‑1 (MQ=255)
gCGTCGAGAATGGTCAGAGGCACGTTGGTAATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGata                                                           <  1:651214/69‑1 (MQ=255)
gCGTCGAGAATGGTCAGAGGCACGTTGGTAATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGata                                                           <  1:573838/69‑1 (MQ=255)
gCGTCGAGAATGGTCAGAGGCACGTTGGTAATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGata                                                           <  1:506962/69‑1 (MQ=255)
gCGTCGAGAATGGTCAGAGGCACGTTGGTAATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGata                                                           <  1:44546/69‑1 (MQ=255)
gCGTCGAGAATGGTCAGAGGCACGTTGGTAATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGata                                                           <  1:3072901/69‑1 (MQ=255)
gCGTCGAGAATGGTCAGAGGCACGTTGGTAATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGata                                                           <  1:1977678/69‑1 (MQ=255)
gCGTCGAGAATGGTCAGAGGCACGTTGGTAATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGata                                                           <  1:1798335/69‑1 (MQ=255)
gCGTCGAGAATGGTCAGAGGCACGTTGGTAATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGata                                                           <  1:1353096/69‑1 (MQ=255)
gCGTCGAGAATGGTCAGAGGCACGTTGGTAATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGata                                                           <  1:1354011/69‑1 (MQ=255)
          tAGTCAGTGGCACGTTGGTGATCGCCTGTTGACCGGATTTATTCACCTGACTGGAGATATAGGCCtttt                                                 >  1:1520224/1‑69 (MQ=255)
          tAGTCAGTGGCACGTTGGTGATCGCCTGTTGACCGGATTTATTCACCTGACTGGAGATATAGGCCtttt                                                 >  1:1636109/1‑69 (MQ=255)
          tAGTCAGTGGCACGTTGGTGATCGCCTGTTGACCGGATTTATTCACCTGACTGGAGATATAGGCCtttt                                                 >  1:499713/1‑69 (MQ=255)
          tAGTCAGTGGCACGTTGGTGATCGCCTGTTGACCGGATTTATTCACCTGACTGGAGATATAGGCCtttt                                                 >  1:3361740/1‑69 (MQ=255)
          tAGTCAGTGGCACGTTGGTGATCGCCTGTTGACCGGATTTATTCACCTGACTGGAGATATAGGCCtttt                                                 >  1:3221963/1‑69 (MQ=255)
          tAGTCAGTGGCACGTTGGTGATCGCCTGTTGACCGGATTTATTCACCTGACTGGAGATATAGGCCtttt                                                 >  1:2905486/1‑69 (MQ=255)
          tAGTCAGTGGCACGTTGGTGATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGATATAGGCCtttt                                                 >  1:1116152/1‑69 (MQ=255)
          tAGTCAGTGGCACGTTGGTGATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGATATAGGCCtttt                                                 >  1:303649/1‑69 (MQ=255)
            gTCAGTGGCACGTTGGTGATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGATATAGGCCtttt                                                 >  1:1074528/1‑67 (MQ=255)
            gTCAGTGGCACGTTGGTGATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGATATAGGCCtttt                                                 >  1:2909805/1‑67 (MQ=255)
            gTCAGTGGCACGTTGGTGATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGATATAGGCCtttt                                                 >  1:1431375/1‑67 (MQ=255)
                          ggTGATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGATATAGGCCTTTTGTGAGCGGAAggcggc                                 <  1:33222/69‑1 (MQ=255)
                          ggTGATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGATATAGGCCTTTTGTGAGCGGAAggcggc                                 <  1:1997990/69‑1 (MQ=255)
                              aTCGCCTGTTGACCGGATTTATTCACCTGACCGGAGATATAGGCCTTTTGTGAGCGGAAGGCGGCGat                              >  1:2995150/1‑68 (MQ=255)
                                cGCCTGTTGACCGGATTTATTCACCTTACTGGAGATATAGGCCTTTTGTGAGCGGAAGGCGGCGATAtt                           <  1:1414535/69‑1 (MQ=255)
                                cGCCTGTTGACCGGATTTATTCACCTGACTGGAGATATAGGGCTTTTGTGAGCGGAAGGCGGCGATAtt                           <  1:1205296/69‑1 (MQ=255)
                                cGCCTGTTGACCGGATTTATTCACCTGACTGGAGATATAGGCCTTTTGTGAGCGGAAGGCGGCGATAtt                           <  1:2345525/69‑1 (MQ=255)
                                cGCCTGTTGACCGGATTTATTCACCTGACTGGAGATATAGGCCTTTTGTGAGCGGAAGGCGGCGATAtt                           <  1:379086/69‑1 (MQ=255)
                                cGCCTGTTGACCGGATTTATTCACCTGACTGGAGATATAGGCCTTTTGTGAGCGGAAGGCGGCGATAtt                           <  1:1985007/69‑1 (MQ=255)
                                cGCCTGTTGACCGGATTTATTCACCTGACTGGAGATATAGGCCTTTTGTGAGCGGAAGGCGGCGATAtt                           <  1:1150999/69‑1 (MQ=255)
                                cGCCTGTTGACCGGATTTATTCACCTGACTGGAGATATAGGCCTTTTGTGAGCGGAAGGCGGCGATAtt                           <  1:1436200/69‑1 (MQ=255)
                                cGCCTGTTGACCGGATTTATTCACCTGACCGGAGATATAGGCCTTTTGTGAGCGGAAGGCGGCGATAtt                           <  1:74194/69‑1 (MQ=255)
                                cGCCTGTTGACCGGATTTATTCACCTGACTGAGATATAGGCCTTTTGTGAGCGGAAGGCGGCGATAtt                           <  1:209541/68‑1 (MQ=255)
                                     gTTGACCGGATTTATTCACCTGACCGGAGATATAGGCCTTTTGTGAGCGGAAggcggc                                 <  1:766859/58‑1 (MQ=255)
                                       tGACCGGATTTATTCACCTGACTGGAGATATAGGCCTTTTGTGAGCGGAAGGCGGCGATAtt                           <  1:3119726/62‑1 (MQ=255)
                                                         tGACCGGAGATATAGGCCTTTTGTGAGCGGAAGGCGGCGATATTAACGTCCACCTGCGGGTCAGCGATg  >  1:2897049/1‑69 (MQ=255)
                                                             |                                                                
GCGTCGAGAATAGTCAGTGGCACGTTGGTGATCGCCTGTTGACCGGATTTATTCACCTGACCGGAGATATAGGCCTTTTGTGAGCGGAAGGCGGCGATATTAACGTCCACCTGCGGGTCAGCGATG  >  minE/1345260‑1345385

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: