Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,854,647 C→T 38.4% L361L (TTG→TTA yjeS ← predicted Fe‑S electron transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,854,6470CT38.4% 17.8 / 22.2 39L361L (TTG→TTAyjeSpredicted Fe‑S electron transport protein
Reads supporting (aligned to +/- strand):  ref base C (6/18);  new base T (10/5);  total (16/23)
Fisher's exact test for biased strand distribution p-value = 1.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

CTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTC  >  minE/2854579‑2854715
                                                                    |                                                                    
cTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGc                                                                      <  1:2661660/69‑1 (MQ=255)
cTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGc                                                                      <  1:2862048/69‑1 (MQ=255)
cTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGc                                                                      <  1:442473/69‑1 (MQ=255)
cTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGc                                                                      <  1:3097952/69‑1 (MQ=255)
     gATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGTAATTg                                                                 <  1:67357/69‑1 (MQ=255)
     gATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGTAATTg                                                                 <  1:549688/69‑1 (MQ=255)
     gATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGTAATTg                                                                 <  1:1523988/69‑1 (MQ=255)
     gATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGTAATTg                                                                 <  1:1801820/69‑1 (MQ=255)
     gATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTg                                                                 <  1:2852244/69‑1 (MQ=255)
     gATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTg                                                                 <  1:138280/69‑1 (MQ=255)
      aTGAATCATGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGTAATTg                                                                 <  1:2575860/68‑1 (MQ=255)
        gAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCAc                                                              >  1:2108823/1‑69 (MQ=255)
              ggCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGAc                                                        >  1:468019/1‑69 (MQ=255)
              ggCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGAc                                                        >  1:763814/1‑69 (MQ=255)
                           aaCCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCAttt                                           >  1:3325733/1‑69 (MQ=255)
                                      aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                 <  1:1158844/68‑1 (MQ=255)
                                      aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                 <  1:2552945/68‑1 (MQ=255)
                                      aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                 <  1:794428/68‑1 (MQ=255)
                                      aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                 <  1:123378/68‑1 (MQ=255)
                                      aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                 <  1:1603856/68‑1 (MQ=255)
                                      aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                 <  1:2116875/68‑1 (MQ=255)
                                      aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                 <  1:1865837/68‑1 (MQ=255)
                                      aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                 <  1:2418921/68‑1 (MQ=255)
                                      aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                 <  1:1891284/68‑1 (MQ=255)
                                       aTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                 <  1:835764/67‑1 (MQ=255)
                                       aTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                 <  1:1127586/67‑1 (MQ=255)
                                       aTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                 <  1:1586240/67‑1 (MQ=255)
                                        tCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTAttt                              >  1:777116/1‑69 (MQ=255)
                                          aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGAcc                                                       >  1:382908/1‑41 (MQ=39)
                                          aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                            >  1:296791/1‑69 (MQ=255)
                                          aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                            >  1:131832/1‑69 (MQ=255)
                                          aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                            >  1:2766498/1‑69 (MQ=255)
                                          aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                            >  1:2653249/1‑69 (MQ=255)
                                          aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                            >  1:2648339/1‑69 (MQ=255)
                                          aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                            >  1:2642364/1‑69 (MQ=255)
                                          aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                            >  1:2258613/1‑69 (MQ=255)
                                          aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                            >  1:2568991/1‑69 (MQ=255)
                                          aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                            >  1:2472206/1‑69 (MQ=255)
                                          aCCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGAc                                                        <  1:485177/41‑1 (MQ=255)
                                          aCCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGAc                                                        <  1:2585789/41‑1 (MQ=255)
                                               aaCCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                 <  1:2068060/59‑1 (MQ=255)
                                                  cAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAAt                    <  1:1848239/69‑1 (MQ=255)
                                                       gCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGcc                >  1:239483/1‑68 (MQ=255)
                                                                    cAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTc  <  1:1523506/69‑1 (MQ=255)
                                                                    |                                                                    
CTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTC  >  minE/2854579‑2854715

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: