Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 528,660 C→T 36.1% T98T (ACG→ACA ybhN ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE528,6600CT36.1% 36.6 / 32.0 58T98T (ACG→ACAybhNconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (23/14);  new base T (7/14);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 5.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

CAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGAC  >  minE/528600‑528724
                                                            |                                                                
cAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGt                                                           >  1:2931676/1‑68 (MQ=255)
        gACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGt                                                           >  1:165983/1‑60 (MQ=255)
        gACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGt                                                  >  1:1969600/1‑69 (MQ=255)
        gACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGt                                                  >  1:2917398/1‑69 (MQ=255)
        gACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGt                                                  >  1:3167137/1‑69 (MQ=255)
        gACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGt                                                  >  1:3355994/1‑69 (MQ=255)
        gACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGt                                                  >  1:2363009/1‑69 (MQ=255)
        gACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGt                                                  >  1:732151/1‑69 (MQ=255)
        gACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGt                                                  >  1:1662065/1‑69 (MQ=255)
        gACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGt                                                  >  1:1524660/1‑69 (MQ=255)
        gACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGt                                                  >  1:981758/1‑69 (MQ=255)
             gagTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCa                                              <  1:3390623/68‑1 (MQ=255)
             gagTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCa                                              <  1:906110/68‑1 (MQ=255)
             gagTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCa                                              <  1:1432893/68‑1 (MQ=255)
             gagTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCa                                              <  1:2398113/68‑1 (MQ=255)
             gagTACAAACGATAGCGCATAACAATGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCa                                              <  1:538444/68‑1 (MQ=255)
              agTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCa                                              <  1:2089984/67‑1 (MQ=255)
                tACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCAGATa                                          <  1:1293443/69‑1 (MQ=255)
                tACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCAGATa                                          <  1:1412232/69‑1 (MQ=255)
                tACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCAGATa                                          <  1:397178/69‑1 (MQ=255)
                tACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCAGATa                                          <  1:2754722/69‑1 (MQ=255)
                tACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCAGATa                                          <  1:522225/69‑1 (MQ=255)
                tACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCAGATa                                          <  1:874771/69‑1 (MQ=255)
                tACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCAGATa                                          <  1:1643152/69‑1 (MQ=255)
                tACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCAGATa                                          <  1:1609766/69‑1 (MQ=255)
                 aCAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATa                                          >  1:1733817/1‑68 (MQ=255)
                  cAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATaa                                         >  1:984349/1‑68 (MQ=255)
                    aaCGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAAc                                       <  1:27100/68‑1 (MQ=255)
                    aaCGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAAc                                       <  1:518557/68‑1 (MQ=255)
                    aaCGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAAc                                       <  1:1726542/68‑1 (MQ=255)
                     aCGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGa                                     <  1:950910/69‑1 (MQ=255)
                     aCGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGa                                     <  1:2652091/69‑1 (MQ=255)
                     aCGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGa                                     <  1:1430977/69‑1 (MQ=255)
                      cGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGac                                    <  1:372888/69‑1 (MQ=255)
                      cGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGac                                    <  1:3357152/69‑1 (MQ=255)
                      cGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGac                                    >  1:1202780/1‑69 (MQ=255)
                      cGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGac                                    <  1:1217616/69‑1 (MQ=255)
                      cGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGac                                    <  1:1398014/69‑1 (MQ=255)
                      cGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGac                                    <  1:2453994/69‑1 (MQ=255)
                      cGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGac                                    >  1:1938707/1‑69 (MQ=255)
                      cGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGac                                    <  1:2168850/69‑1 (MQ=255)
                                  ccAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCAcc                        >  1:2831653/1‑69 (MQ=255)
                                  ccAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCAcc                        >  1:2550851/1‑69 (MQ=255)
                                     aTGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCAGATAAACGACTCCAGCATCACCTgg                     >  1:225107/1‑69 (MQ=255)
                                     aTGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTgg                     >  1:3111794/1‑69 (MQ=255)
                                     aTGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTgg                     >  1:1917833/1‑69 (MQ=255)
                                     aTGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTgg                     >  1:1986507/1‑69 (MQ=255)
                                     aTGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTgg                     >  1:2255311/1‑69 (MQ=255)
                                     aTGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTgg                     >  1:2343929/1‑69 (MQ=255)
                                     aTGCCGCCGACCCAGGTACTGAGTGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTg                      >  1:3380469/1‑68 (MQ=255)
                                        ccgccgACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGgcgc                  >  1:593264/1‑69 (MQ=255)
                                        ccgccgACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGgcgc                  >  1:2945436/1‑69 (MQ=255)
                                            cgACCCAGGTACTTAGCGTCAGGTTGAAGGCGTAGCAGATAAAc                                       <  1:641173/44‑1 (MQ=255)
                                                    gTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGtt        <  1:1451688/67‑1 (MQ=255)
                                                         gAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCAt                            >  1:2097351/1‑42 (MQ=255)
                                                         gAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGAc  >  1:504036/1‑68 (MQ=255)
                                                         gAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGAc  >  1:1011873/1‑68 (MQ=255)
                                                         gAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAACATCACCTGGCGCTTCGCCAGTTTGTGAc  >  1:3307228/1‑68 (MQ=255)
                                                            |                                                                
CAACCCCAGACGAGAGTACAAACGATAGCGCATACCAATGCCGCCGACCCAGGTACTGAGCGTCAGGTTGAAGGCGTAGCAGATAAACGACACCAGCATCACCTGGCGCTTCGCCAGTTTGTGAC  >  minE/528600‑528724

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: