Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 138,172 A→G 100% N665N (AAT→AAC htrE ← predicted outer membrane usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE138,1720AG84.8% 306.6 / 40.7 125N665N (AAT→AAChtrEpredicted outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base A (18/1);  new base G (36/70);  total (54/71)
Fisher's exact test for biased strand distribution p-value = 5.23e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 125 total aligned reads displayed.

CCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCT  >  minE/138113‑138219
                                                           |                                               
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:1656700/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:1344200/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:2047622/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:1429059/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:327932/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:1502314/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:1504186/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:2275440/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:1628798/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:2205021/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:476897/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:1669213/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:1686330/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:1989097/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:1704501/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:1988581/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:1722405/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:466658/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:818785/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:2195804/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:1088324/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:1142772/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:2131139/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:2259045/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:908994/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:188318/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:900052/66‑1 (MQ=255)
ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:876205/66‑1 (MQ=255)
          aTACAAAACGACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAcc                                >  1:794005/1‑67 (MQ=255)
          aTACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAcc                                >  1:577963/1‑67 (MQ=255)
          aTACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAcc                                >  1:1536446/1‑67 (MQ=255)
          aTACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAcc                                >  1:1967407/1‑67 (MQ=255)
          aTACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAcc                                >  1:1280402/1‑67 (MQ=255)
          aTACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAcc                                >  1:1459796/1‑67 (MQ=255)
          aTACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAcc                                >  1:2045496/1‑67 (MQ=255)
          aTACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAcc                                >  1:299270/1‑67 (MQ=255)
          aTACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAcc                                >  1:2150401/1‑67 (MQ=255)
          aTACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAcc                                >  1:1413243/1‑67 (MQ=255)
          aTACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAcc                                >  1:1055559/1‑67 (MQ=255)
          aTACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAcc                                >  1:1040093/1‑67 (MQ=255)
          aTACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAcc                                >  1:1755267/1‑67 (MQ=255)
            aCAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:2278878/54‑1 (MQ=255)
                          ggTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                           <  1:628634/40‑1 (MQ=255)
                            tGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATg               <  1:1202060/66‑1 (MQ=255)
                                    gaAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCg                        <  1:1976157/49‑1 (MQ=255)
                                    gaAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCg                        <  1:2026875/49‑1 (MQ=255)
                                    gaAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCg                        <  1:226832/49‑1 (MQ=255)
                                    gaAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCg                        <  1:1948157/49‑1 (MQ=255)
                                    gaAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCg                        <  1:1028997/49‑1 (MQ=255)
                                    gaAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCg                        <  1:1713034/49‑1 (MQ=255)
                                    gaAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCg                        <  1:169870/49‑1 (MQ=255)
                                    gaAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCg                        <  1:495772/49‑1 (MQ=255)
                                    gaAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCg                        <  1:513056/49‑1 (MQ=255)
                                    gaAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCg                        <  1:76494/49‑1 (MQ=255)
                                    gaAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCg                        <  1:769239/49‑1 (MQ=255)
                                    gaAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCg                        <  1:851138/49‑1 (MQ=255)
                                    gaAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCg                        <  1:1170762/49‑1 (MQ=255)
                                    gaAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCg                        <  1:1100044/49‑1 (MQ=255)
                                    gaAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCg                        <  1:1071460/49‑1 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:282113/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:376497/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:433609/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:952324/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:996211/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:2258413/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:1540243/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:1031215/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:1198293/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:1199453/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:1272139/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:1280932/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:1282025/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:1325150/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:145425/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:1502077/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:2241476/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:1733504/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:1831587/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:2101561/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:2025690/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:1861433/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:1771607/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa     >  1:1774326/1‑67 (MQ=255)
                                     aaaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCACCATGGTGACTCATa     >  1:542634/1‑67 (MQ=255)
                                      aaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGACCCTGCCAGCGtt                      >  1:1598514/1‑49 (MQ=255)
                                      aaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGtt                      >  1:661866/1‑49 (MQ=255)
                                      aaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGtt                      >  1:1515529/1‑49 (MQ=255)
                                      aaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGtt                      >  1:858249/1‑49 (MQ=255)
                                      aaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGtt                      >  1:2311275/1‑49 (MQ=255)
                                      aaCTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGtt                      >  1:1660144/1‑49 (MQ=255)
                                        cTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCt  <  1:613069/67‑1 (MQ=255)
                                        cTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCt  <  1:330799/67‑1 (MQ=255)
                                        cTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCt  <  1:1355777/67‑1 (MQ=255)
                                        cTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCt  <  1:1341749/67‑1 (MQ=255)
                                        cTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCt  <  1:870510/67‑1 (MQ=255)
                                        cTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCt  <  1:194284/67‑1 (MQ=255)
                                        cTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCt  <  1:2093843/67‑1 (MQ=255)
                                        cTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCt  <  1:428349/67‑1 (MQ=255)
                                        cTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCCGCCAGCGTTCCCCATGGTGACTCATAGCt  <  1:505733/67‑1 (MQ=255)
                                                    tATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCt  <  1:1877648/55‑1 (MQ=255)
                                                           |                                               
CCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCT  >  minE/138113‑138219

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