Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 756,321 G→C 12.2% E150Q (GAA→CAA)  lolD → outer membrane‑specific lipoprotein transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE756,3210GC12.2% 138.0 / 8.7 49E150Q (GAA→CAA) lolDouter membrane‑specific lipoprotein transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (30/13);  new base C (4/2);  total (34/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GAGATGTTAAAAGCGGTGGGGCTGGATCATCGTGCGAATCACCGCCCATCTGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGTGAATAACCCGCGCCTGGTACTGGCGGATGA  >  minE/756255‑756385
                                                                  |                                                                
gagaTGTTAAAAGCGGTGGGGCTGGATCATCTTGCGAATCACCGCCCATCTGAACTTTCTGGCGGCGAACg                                                              >  1:4551754/1‑71 (MQ=255)
gagaTGTTAAAAGCGGTGGGGCTGGATCATCGTGCGAATCACCGCCCATCTGAACTTTCTGGCGGCGAACg                                                              >  1:891801/1‑71 (MQ=255)
gagaTGTTAAAAGCGGTGGGGCTGGATCATCGTGCGAATCACCGCCCATCTGAACTTTCTGGCGGCGAACg                                                              >  1:181060/1‑71 (MQ=255)
gagaTGTTAAAAGCGGTGGGGCTGGATCATCGTGCGAATCACCGCCCATCTGAACTTTCTGGCGGCGAACg                                                              >  1:2222086/1‑71 (MQ=255)
gagaTGTTAAAAGCGGTGGGGCTGGATCATCGTGCGAATCACCGCCCATCTGAACTTTCTGGCGGCGAACg                                                              >  1:2795022/1‑71 (MQ=255)
    tGTTAAAAGCGGTGGGGCTGGATCATCGTGCGAATCACCGCCCATCTGAACTTTCTGGCGGCGAACGCCAg                                                          <  1:3312325/71‑1 (MQ=255)
        aaaaGCGGTGGGGCTGGATCATCGTGCGAATCACTGCCCATCTGAACTTTCTGGCGGCGaa                                                                >  1:4113748/1‑61 (MQ=255)
        aaaaGCGGTGGGGCTGGATCATCGTGCGAATCACCGCCCATCTGAACTTTCTGGCGGCGaa                                                                >  1:3010365/1‑61 (MQ=255)
        aaaaGCGGTGGGGCTGGATCATCGTGCGAATCACCGCCCATCTGAACTTTCTGGCGGCGaa                                                                >  1:4936149/1‑61 (MQ=255)
        aaaaGCGGTGGGGCTGGATCATCGTGCGAATCACCGCCCATCTGAACTTTCTGGCGGCGaa                                                                >  1:624727/1‑61 (MQ=255)
        aaaaGCGGTGGGGCTGGATCATCGTGCGAATCACCGCCCATCTGAACTTTCTGGCGGCGAACGCCAGCgtg                                                      >  1:2557694/1‑71 (MQ=255)
         aaaGCGGTGGGGCTGGATCATCGTGCGAATCACCGCCCATCTGAACTTTCTGGCGGCGAACGCCAGCgtgt                                                     >  1:1230763/1‑71 (MQ=255)
            gCGGTGGGGCTGGATCATCGTGCGAATCACCGCCCATCTGAACTTTCTGGCGGCGAACGCCAGCGTGTGGc                                                  <  1:3660592/71‑1 (MQ=255)
                    gCTGGATCATCGTGCGAATCACCGCCCATCTGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGa                                                <  1:1465715/65‑1 (MQ=255)
                    gCTGGATCATCGTGCGAATCACCGCCCATCTGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGccc                                          >  1:4117877/1‑71 (MQ=255)
                     cTGGATCATCGTGCGAATCACCGCCCATCTGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCg                                         >  1:532499/1‑71 (MQ=255)
                          tcatcGTGCGAATCACCGCCCATCTGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTgcgc                                    >  1:4637102/1‑71 (MQ=255)
                             acGTGCACATCACTACCCTTCCGAACTTTCTGGTGGTCAACAGCAGCGTGTGGCGATTGCCCGCGCGCTgg                                 <  1:1602388/70‑1 (MQ=9)
                              cGTGCGAATCACCGCCCATCTGAACTTTCTGGCGGCGAACGCCAGCgtg                                                      >  1:5556916/1‑49 (MQ=255)
                              cGTGCACATCACTACCCTTCCGAACTTTCTGGTGGTCAACAGCAGCGTGTGGCGATTGCCCGCGCGCTggc                                >  1:4738873/1‑70 (MQ=9)
                                   gAATCACCGCCCATCTGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGa                                                >  1:3205053/1‑50 (MQ=255)
                                    caTCACTACCCTTCCGAACTTTCTGGTGGTCAACAGCAGCGTGTGGCGATTGccc                                          >  1:3741188/2‑55 (MQ=9)
                                    caTCACTACCCTTCCGAACTTTCTGGTGGTCAACAGCAGCGTGTGGCGATTGccc                                          >  1:2360772/2‑55 (MQ=9)
                                    aaTCACCGCCCATCTGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGTGAATaa                          <  1:295227/71‑1 (MQ=255)
                                     aTCACTACCCTTCCGAACTTTCTGGTGGTCAACAGCAGCGTGTGGCGATTGCCCgcg                                       >  1:5039457/1‑55 (MQ=9)
                                     aTCACCGCCCATCTGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGTGAATAAc                         <  1:2119147/71‑1 (MQ=255)
                                      tCACCGCCCATCTGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGTGAATAAcc                        >  1:110642/1‑71 (MQ=255)
                                        aCCGCCCATCTGAACTTTCTGGCGGCGAACGCCAGCgt                                                       >  1:799255/1‑38 (MQ=255)
                                        aCCGCCCATCTGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGTGAATAACCcg                      >  1:2491998/1‑71 (MQ=255)
                                         ccGCCCATCTGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGTGAATAACCcgc                     <  1:2408147/71‑1 (MQ=255)
                                          cGCCCATCTGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGTGAATAACCcgcg                    >  1:5257174/1‑71 (MQ=255)
                                          cGCCCATCTGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGTGAATAACCcgc                     >  1:297353/1‑70 (MQ=255)
                                            cccATCTGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCTATTGCCCGTCCGCTGGTGAATAACCCGCGcc                  <  1:747219/71‑1 (MQ=255)
                                            cccATCTGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGTGAATAACCCGCGcc                  <  1:3835882/71‑1 (MQ=255)
                                                  tGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGccc                                          >  1:3982004/1‑41 (MQ=255)
                                                  tGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGt                                >  1:4986793/1‑51 (MQ=255)
                                                  tGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGt                                >  1:449435/1‑51 (MQ=255)
                                                  tGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGt                                >  1:3204423/1‑51 (MQ=255)
                                                   gAACTTTCTGGTGGTCAACAGCAGCGTGTGGCGATTGCCCgcg                                       <  1:2434335/43‑3 (MQ=9)
                                                   gAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGTGAATAACCCGCGCCTgg               <  1:2890172/67‑1 (MQ=255)
                                                    aaCTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGTGAATAACCCGCGCCTGGt              <  1:5209298/67‑1 (MQ=255)
                                                    aaCTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGTGAATAACCCGCGCCTGGt              <  1:4163726/67‑1 (MQ=255)
                                                    aaCTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGTGAATAACCCGCGCCTGGt              <  1:5636335/67‑1 (MQ=255)
                                                    aaCTTTCTGGCGGCGAACGCCAGCCTGTGGCGATTGCCCGTGCGCTGGTGAATAACCCGCGCCTGGTACTg          >  1:1534050/1‑71 (MQ=255)
                                                     aCTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGTGAATAACCCGCGCCTGGTACTgg         >  1:4669592/1‑71 (MQ=255)
                                                      cTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGTGAATAACCCGCGCCTGGTAc            >  1:1533487/1‑67 (MQ=255)
                                                      cTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGACCGTGCGCTGGTGAATAACCCGCGCCTGGTAc            >  1:430245/1‑67 (MQ=255)
                                                         tCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGTGAATAACCCGCGCCTGGTACTGGCgg      <  1:5689915/70‑1 (MQ=255)
                                                         tCTGGCGGCGAACGCCAGCGTGTGGCGATTGACCGTGCGCTGGTGAATAACCCGCGCCTGGTACTGGCGGa     >  1:920414/1‑71 (MQ=255)
                                                             gcggcgAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGTGAATAACCCGCGCCTGGTACTGGCGGATGa  <  1:4567061/70‑1 (MQ=255)
                                                                  |                                                                
GAGATGTTAAAAGCGGTGGGGCTGGATCATCGTGCGAATCACCGCCCATCTGAACTTTCTGGCGGCGAACGCCAGCGTGTGGCGATTGCCCGTGCGCTGGTGAATAACCCGCGCCTGGTACTGGCGGATGA  >  minE/756255‑756385

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: