Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,195,634 C→T 44.7% S54S (AGC→AGT yhcB → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,195,6340CT44.7% 13.9 / 51.3 67S54S (AGC→AGTyhcBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/11);  new base T (18/12);  total (44/23)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCGTTGCAGTACGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTACTGGATACCATGGCGCACGACTATCGCCAGCTG  >  minE/2195566‑2195691
                                                                    |                                                         
cAGGCGTTGCAGTACGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGc                                                           <  1:1000500/69‑1 (MQ=255)
cAGGCGTTGCAGTACGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGc                                                           <  1:2226125/69‑1 (MQ=255)
cAGGCGTTGCAGTACGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGc                                                           <  1:653791/69‑1 (MQ=255)
cAGGCGTTGCAGTACGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGc                                                           <  1:1413794/69‑1 (MQ=255)
cAGGCGTTGCAGTACGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGc                                                           <  1:1486295/69‑1 (MQ=255)
     gTTGCAGTACGAACTGGAAAAGAATAAAGCTTAACTGGACGAGTATCGCGAAGAGCTGGTTAGc                                                           <  1:545364/64‑1 (MQ=255)
     gTTGCAGTACGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGc                                                           <  1:1421500/64‑1 (MQ=255)
          aGTACGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  <  1:2407541/68‑1 (MQ=255)
          aGTACGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  <  1:1700322/68‑1 (MQ=255)
          aGTACGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGAGAAGAGCTGGTTAGTCACTTTGcc                                                  <  1:1756453/68‑1 (MQ=255)
          aGTACGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGAGAAGAGCTGGTTAGTCACTTTGcc                                                  <  1:2348877/68‑1 (MQ=255)
          aGTACGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGAGAAGAGCTGGTTAGTCACTTTGcc                                                  <  1:631188/68‑1 (MQ=255)
          aGTACGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGAGAAGAGCTGGTTAGTCACTTTGcc                                                  <  1:1595003/68‑1 (MQ=255)
            tACGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGc                                               >  1:1634258/1‑69 (MQ=255)
              ggAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  >  1:2948903/2‑64 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  >  1:451518/1‑64 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  >  1:2564900/1‑64 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  >  1:2459618/1‑64 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  >  1:910741/1‑64 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  >  1:2333778/1‑64 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  >  1:231581/1‑64 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  >  1:3194591/1‑64 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  >  1:2119223/1‑64 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  >  1:207255/1‑64 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  >  1:1447922/1‑64 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  >  1:1176385/1‑64 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  >  1:57659/1‑64 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  >  1:1246993/1‑64 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  >  1:1253718/1‑64 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  >  1:885111/1‑64 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  >  1:732906/1‑64 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGcc                                                  >  1:1815278/1‑64 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCa                                              >  1:871168/1‑68 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAg                                             >  1:170591/1‑69 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAg                                             >  1:1488576/1‑69 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAg                                             >  1:1589514/1‑69 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAg                                             >  1:162584/1‑69 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAg                                             >  1:2888550/1‑69 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAg                                             >  1:1169130/1‑69 (MQ=255)
              cGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAg                                             >  1:3251321/1‑69 (MQ=255)
                          gAATAAAGCTGAACTGGACGAGTATCGCGAATAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTAc                                  >  1:1051717/1‑68 (MQ=255)
                          gAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGCCCGCAGCGCGGAATTACt                                 <  1:1042365/69‑1 (MQ=255)
                          gAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGCCCGCAGCGCGGAATTACt                                 <  1:3395987/69‑1 (MQ=255)
                          gAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGCCCGCAGCGCGGAATTACt                                 <  1:2708417/69‑1 (MQ=255)
                          gAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGCCCGCAGCGCGGAATTACt                                 <  1:2633977/69‑1 (MQ=255)
                          gAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGCCCGCAGCGCGGAATTACt                                 <  1:630983/69‑1 (MQ=255)
                          gAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGTCACTTTGCCCGCAGCGCGGAATTACt                                 <  1:1960211/69‑1 (MQ=255)
                          gAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCGGCAGCGCGGAATTAc                                  >  1:897261/1‑68 (MQ=255)
                          gAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTAc                                  >  1:2551571/1‑68 (MQ=255)
                          gAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTAc                                  >  1:2505155/1‑68 (MQ=255)
                          gAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTAc                                  >  1:2865195/1‑68 (MQ=255)
                          gAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTAc                                  >  1:1737279/1‑68 (MQ=255)
                          gAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTAc                                  >  1:3125913/1‑68 (MQ=255)
                          gAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTAc                                  >  1:224605/1‑68 (MQ=255)
                          gAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTAc                                  >  1:1785965/1‑68 (MQ=255)
                          gAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTAc                                  >  1:3276813/1‑68 (MQ=255)
                           aaTAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTACTg                                <  1:2643466/69‑1 (MQ=255)
                           aaTAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTACTg                                <  1:1121183/69‑1 (MQ=255)
                           aaTAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTACTg                                <  1:1289471/69‑1 (MQ=255)
                           aaTAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTACTg                                <  1:178584/69‑1 (MQ=255)
                           aaTAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTACTg                                <  1:2981650/69‑1 (MQ=255)
                               aaGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAAt                                     <  1:53847/60‑1 (MQ=255)
                                aGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGc                                                           <  1:2313176/37‑1 (MQ=255)
                                 gCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTACTGGATAc                           <  1:1842855/68‑1 (MQ=255)
                                     aaCTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTACTg                                <  1:684967/59‑1 (MQ=255)
                                            cGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTACTGGATACCATGGCGCAc                 <  1:17095/67‑1 (MQ=255)
                                            cGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTACTGGATACCATGGCGCAc                 <  1:1352341/67‑1 (MQ=255)
                                                  tCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTACTGGATACCATGGCGCACGACTa            >  1:498572/1‑66 (MQ=255)
                                                  tCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTACTGGATACCATGGCGCACGACTa            >  1:1841636/1‑66 (MQ=255)
                                                  tCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTACTGGATACCATGGCGCACGACTa            >  1:3007550/1‑66 (MQ=255)
                                                         gagCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTACTGGATACCATGGCGCACGACTATCGCCAGCt   >  1:2887355/1‑68 (MQ=255)
                                                         gagCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTACTGGATACCATGGCGCACGACTATCGCCAGCTg  >  1:187313/1‑69 (MQ=255)
                                                         gagCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTACTGGATACCATGGCGCACGACTATCGCCAGCTg  >  1:607079/1‑69 (MQ=255)
                                                         gagCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTACTGGATACCATGGCGCACGACTATCGCCAGCTg  >  1:2285651/1‑69 (MQ=255)
                                                          agCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTACTg                                <  1:1971536/38‑1 (MQ=255)
                                                                    |                                                         
CAGGCGTTGCAGTACGAACTGGAAAAGAATAAAGCTGAACTGGACGAGTATCGCGAAGAGCTGGTTAGCCACTTTGCCCGCAGCGCGGAATTACTGGATACCATGGCGCACGACTATCGCCAGCTG  >  minE/2195566‑2195691

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