Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,091,823 C→A 38.2% P85P (CCG→CCT sufA ← Fe‑S cluster assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,091,8230CA38.2% 36.4 / 35.3 47P85P (CCG→CCTsufAFe‑S cluster assembly protein
Reads supporting (aligned to +/- strand):  ref base C (14/15);  new base A (3/15);  total (17/30)
Fisher's exact test for biased strand distribution p-value = 3.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCGGGACAAACAGCTTCGCGCCGTCGTGTTCAAACAGCAGATCGTCTTTG  >  minE/1091756‑1091886
                                                                   |                                                               
ggTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAAc                                                                 <  1:2201419/68‑1 (MQ=255)
ggTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAAc                                                                 <  1:3384277/68‑1 (MQ=255)
ggTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAAc                                                                 <  1:3356558/68‑1 (MQ=255)
ggTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAAc                                                                 <  1:3320195/68‑1 (MQ=255)
ggTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAAc                                                                 <  1:300005/68‑1 (MQ=255)
ggTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAAc                                                                 <  1:2799609/68‑1 (MQ=255)
 gTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAAc                                                                 <  1:1554615/67‑1 (MQ=255)
  ttGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGc                                                              <  1:2988388/69‑1 (MQ=255)
  ttGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGc                                                              <  1:2254775/69‑1 (MQ=255)
  ttGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGc                                                              <  1:2349277/69‑1 (MQ=255)
  ttGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGc                                                              <  1:2379195/69‑1 (MQ=255)
  ttGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGc                                                              <  1:2778886/69‑1 (MQ=255)
  ttGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGc                                                              <  1:1941636/69‑1 (MQ=255)
  ttGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGc                                                              <  1:2897300/69‑1 (MQ=255)
  ttGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGc                                                              <  1:1484412/69‑1 (MQ=255)
    gtgAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCAt                                                            >  1:2796387/1‑69 (MQ=255)
    gtgAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCAt                                                            >  1:2656702/1‑69 (MQ=255)
    gtgAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCAt                                                            >  1:629240/1‑69 (MQ=255)
    gtgAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCAt                                                            >  1:852575/1‑69 (MQ=255)
    gtgAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCAt                                                            >  1:940686/1‑69 (MQ=255)
      gAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCg                                                          >  1:377584/1‑69 (MQ=255)
                  tCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCggg                                               >  1:682174/1‑68 (MQ=255)
                        tAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCGGGACAAAc                                         >  1:538617/1‑68 (MQ=255)
                              cTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCGGGACAAACAGCTTcg                                  >  1:634565/1‑69 (MQ=255)
                              cTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCGGGACAAACAGCTTcg                                  >  1:1525523/1‑69 (MQ=255)
                              cTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCGGGACAAACAGCTTc                                   >  1:2107098/1‑68 (MQ=255)
                                tcgcgaacgaaATCGACTTCCGTGCCATCAATAAAAGGCATCGCTTGCAGCGGGACAAACAGCTTcgcg                                <  1:3239881/66‑1 (MQ=255)
                                tcgcgaacgaaATCGACTTCCGTGCCATCAATAAAAGGCATCGCTTGCAGCGGGACAAACAGCTTcgcg                                <  1:1094274/66‑1 (MQ=255)
                                tcgcgaacgaaATCGACTTCCGTGCCATCAATAAAAGGCATCGCTTGCAGCGGGACAAACAGCTTcgcg                                <  1:786062/66‑1 (MQ=255)
                                tcgcgaacgaaATCGACTTCCGTGCCATCAATAAAAGGCATCGCTTGCAGCGGGACAAACAGCTTcgcg                                <  1:760754/66‑1 (MQ=255)
                                tcgcgaacgaaATCGACTTCCGTGCCATCAATAAAAGGCATCGCTTGCAGCGGGACAAACAGCTTcgcg                                <  1:1191591/66‑1 (MQ=255)
                                tcgcgaacgaaATCGACTTCCGTGCCATCAATAAAAGGCATCGCTTGCAGCGGGACAAACAGCTTcgcg                                <  1:3393885/66‑1 (MQ=255)
                                tcgcgaacgaaATCGACTTCCGTGCCATCAATAAAAGGCATCGCTTGCAGCGGGACAAACAGCTTcgcg                                <  1:3351432/66‑1 (MQ=255)
                                tcgcgaacgaaATCGACTTCCGTGCCATCAATAAAAGGCATCGCTTGCAGCGGGACAAACAGCTTcgcg                                <  1:2174457/66‑1 (MQ=255)
                                tcgcgaacgaaATCGACTTCCGTGCCATCAATAAAAGGCATCGCTTGCAGCGGGACAAACAGCTTcgcg                                <  1:3041736/66‑1 (MQ=255)
                                tcgcgaacgaaATCGACTTCCGTGCCATCAATAAAAGGCATCGCTTGCAGCGGGACAAACAGCTTcgcg                                <  1:19102/66‑1 (MQ=255)
                                tcgcgaacgaaATCGACTTCCGTGCCATCAATAAAAGGCATCGCTTGCAGCGGGACAAACAGCTTcgcg                                <  1:1549622/66‑1 (MQ=255)
                                tcgcgaacgaaATCGACTTCCGTGCCATCAATAAAAGGCATCGCTTGCAGCGGGACAAACAGCTTcgcg                                <  1:2737329/66‑1 (MQ=255)
                                tcgcgaacgaaATCGACTTCCGTGCCATCAATAAAAGGCATCGCTTGCAGCGGGACAAACAGCTTcgcg                                <  1:2632102/66‑1 (MQ=255)
                                tcgcgaacgaaATCGACTTCCGTGCCATCAATAAAAGGCATCGCTTGCAGCGGGACAAACAGCTTcgcg                                <  1:2126262/66‑1 (MQ=255)
                                       cgaaATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCGGGACAAACAGCTTCGCGCcgtcgt                         <  1:1778219/69‑1 (MQ=255)
                                            tCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCGGGACAAACAGCTTCGCGCCGTCGTGTTCa                    <  1:1987038/69‑1 (MQ=255)
                                            tCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCGGGACAAACAGCTTCGCGCCGTCGTGTTCa                    <  1:1173458/69‑1 (MQ=255)
                                                cTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCGGGACAAACAGCTTCGCGCCGTCGTGTTCAAAca                <  1:2878389/69‑1 (MQ=255)
                                                cTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCGGGACAAACAGCTTCGCGCCGTCGTGTTCAAAca                <  1:2652110/69‑1 (MQ=255)
                                                       gCCATCAATAAACGGCATCGCTTGCAGCGGGACAAACAGCTTCGCGCCGTCGTGTTCAAACAGCAGATc         <  1:1597249/69‑1 (MQ=255)
                                                         cATCAATAAAAGGCATCGCTTGCAGCGGGACAAACAGCTTcgcg                                <  1:3386500/44‑1 (MQ=255)
                                                             aaTAAACGGCATCGCTTGCAGCGGGACAAACAGCTTCGCGCCGTCGTGTTCAAACAGCAGATCGTCttt   >  1:142557/1‑69 (MQ=255)
                                                             aaTAAACGGCATCGCTTGCAGCGGGACAAACAGCTTCGCGCCGTCGTGTTCAAACAGCAGATCGTCttt   >  1:515879/1‑69 (MQ=255)
                                                             aaTAAACGGCATCGCTTGCAGCGGGACAAACAGCTTCGCGCCGTCGTGTTCAAACAGCAGATCGTCttt   >  1:2320048/1‑69 (MQ=255)
                                                              aTAAAAGGCATCGCTTGCAGCGGGACAAACAGCTTCGCGCCTTCGTGTTCAAACAGCAGATCGTCTTTg  >  1:2316535/1‑69 (MQ=255)
                                                              aTAAAAGGCATCGCTTGCAGCGGGACAAACAGCTTCGCGCCTTCGTGTTCAAACAGCAGATCGTCTTTg  >  1:2024479/1‑69 (MQ=255)
                                                              aTAAAAGGCATCGCTTGCAGCGGGACAAACAGCTTCGCGCCTTCGTGTTCAAACAGCAGATCGTCTTTg  >  1:1308429/1‑69 (MQ=255)
                                                                 aaCGGCATCGCTTGCAGCGGGACAAACAGCTTCGCGCCGTCGTGTTCAAACAGCAGATc         <  1:1008788/59‑1 (MQ=255)
                                                                   |                                                               
GGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGTGCCATCAATAAACGGCATCGCTTGCAGCGGGACAAACAGCTTCGCGCCGTCGTGTTCAAACAGCAGATCGTCTTTG  >  minE/1091756‑1091886

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