Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,690,443 T→C 37.4% P393P (CCA→CCG hscA ← DnaK‑like molecular chaperone specific for IscU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,690,4430TC37.4% 37.4 / 45.0 72P393P (CCA→CCGhscADnaK‑like molecular chaperone specific for IscU
Reads supporting (aligned to +/- strand):  ref base T (18/27);  new base C (19/8);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

ACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGC  >  minE/1690389‑1690511
                                                      |                                                                    
aCGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCaa                                                         >  1:726919/1‑68 (MQ=255)
   cgGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGc                                                       >  1:2751126/1‑67 (MQ=255)
      aaTCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAg                                                  <  1:932204/69‑1 (MQ=255)
      aaTCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAg                                                  <  1:735318/69‑1 (MQ=255)
          aCTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCAtt                                              >  1:2089211/1‑69 (MQ=255)
          aCTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCAtt                                              >  1:2985316/1‑69 (MQ=255)
          aCTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCAtt                                              >  1:3022060/1‑69 (MQ=255)
          aCTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCAtt                                              >  1:2680402/1‑69 (MQ=255)
          aCTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCAtt                                              >  1:243675/1‑69 (MQ=255)
          aCTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCAtt                                              >  1:231213/1‑69 (MQ=255)
          aCTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCAtt                                              >  1:2107631/1‑69 (MQ=255)
          aCTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCAtt                                              >  1:625111/1‑69 (MQ=255)
          aCTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCAtt                                              >  1:2019645/1‑69 (MQ=255)
          aCTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCAtt                                              >  1:1850426/1‑69 (MQ=255)
          aCTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCAtt                                              >  1:659763/1‑69 (MQ=255)
          aCTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCAtt                                              >  1:708358/1‑69 (MQ=255)
          aCTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCAt                                               >  1:2876435/1‑68 (MQ=255)
                 ccaccaGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTg                                        <  1:1690079/68‑1 (MQ=255)
                 ccaccaGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTg                                        <  1:2935794/68‑1 (MQ=255)
                 ccaccaGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTg                                        <  1:1502747/68‑1 (MQ=255)
                 ccaccaGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTg                                        <  1:2818700/68‑1 (MQ=255)
                 ccaccaGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTg                                        <  1:372941/68‑1 (MQ=255)
                 ccaccaGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTg                                        <  1:3138308/68‑1 (MQ=255)
                 ccaccaGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTg                                        <  1:1331762/68‑1 (MQ=255)
                              cATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTAcc                          <  1:1756529/69‑1 (MQ=255)
                              cATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTAcc                          <  1:2865068/69‑1 (MQ=255)
                              cATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTAcc                          <  1:123153/69‑1 (MQ=255)
                              cATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTAcc                          <  1:2312864/69‑1 (MQ=255)
                               aTCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTAcc                          <  1:2281996/68‑1 (MQ=255)
                               aTCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTAcc                          <  1:852163/68‑1 (MQ=255)
                                  gTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTAccc                         >  1:1007978/1‑66 (MQ=255)
                                  gTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTAccc                         >  1:3312450/1‑66 (MQ=255)
                                  gTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTAccc                         >  1:604020/1‑66 (MQ=255)
                                       gAGGCCCAGCGACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGaa                  <  1:819569/68‑1 (MQ=255)
                                       gAGGCCCAGCGACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGaa                  <  1:2514860/68‑1 (MQ=255)
                                       gAGGCCCAGCGACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGaa                  <  1:2723011/68‑1 (MQ=255)
                                       gAGGCCCAGCGACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGaa                  <  1:760250/68‑1 (MQ=255)
                                       gAGGCCCAGCGACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGaa                  <  1:1201769/68‑1 (MQ=255)
                                       gAGGCCCAGCGACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGaa                  <  1:645925/68‑1 (MQ=255)
                                       gAGGCCCAGCGACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGaa                  <  1:1653558/68‑1 (MQ=255)
                                       gAGGCCCAGCGACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGaa                  <  1:701922/68‑1 (MQ=255)
                                                 gACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGa       <  1:3273610/69‑1 (MQ=255)
                                                 gACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGa       <  1:3199050/69‑1 (MQ=255)
                                                 gACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGa       <  1:3294481/69‑1 (MQ=255)
                                                 gACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGa       <  1:1010089/69‑1 (MQ=255)
                                                 gACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGa       <  1:2337681/69‑1 (MQ=255)
                                                 gACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGa       <  1:1182647/69‑1 (MQ=255)
                                                 gACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGa       <  1:2317794/69‑1 (MQ=255)
                                                 gACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGa       <  1:1707408/69‑1 (MQ=255)
                                                 gACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGa       <  1:2131405/69‑1 (MQ=255)
                                                 gACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGa       <  1:1786145/69‑1 (MQ=255)
                                                 gACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGAata                >  1:751751/1‑60 (MQ=255)
                                                 gACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGAATATCCGCCTGa       >  1:2044506/1‑69 (MQ=255)
                                                 gACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGAATATCCGCCTGa       >  1:985099/1‑69 (MQ=255)
                                                 gACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGAATATCCGCCTGa       >  1:956581/1‑69 (MQ=255)
                                                 gACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGAATATCCGCCTGa       >  1:137425/1‑69 (MQ=255)
                                                 gACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGAATATCCGCCTGa       >  1:1417889/1‑69 (MQ=255)
                                                 gACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGAATATCCGCCTGa       >  1:1452255/1‑69 (MQ=255)
                                                 gACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGAATATCCGCCTGa       >  1:717512/1‑69 (MQ=255)
                                                 gACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGAATATCCGCCTGa       >  1:1764025/1‑69 (MQ=255)
                                                 gACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGAATATCCGCCTGa       >  1:1806149/1‑69 (MQ=255)
                                                 gACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGAATATCCGCCTGa       >  1:2039060/1‑69 (MQ=255)
                                                 gACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGAATATCCGCCTGa       >  1:2849937/1‑69 (MQ=255)
                                                 gACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGAATATCCGCCTGa       >  1:590288/1‑69 (MQ=255)
                                                 gACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGAATATCCGCCTGa       >  1:488317/1‑69 (MQ=255)
                                                 gACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGAATATCCGCCTGa       >  1:426950/1‑69 (MQ=255)
                                                 gACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGAATATCCGCCTGa       >  1:2284396/1‑69 (MQ=255)
                                                 gACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGAATATCCGCCTGa       >  1:2307900/1‑69 (MQ=255)
                                                 gACAGCGGGATCACATCCAGCAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGAATATCCGCCTGa       >  1:2680290/1‑69 (MQ=255)
                                                 gACAGCGGGATCACATCCAGAAACAGCATTTCGCTGTCCGGCTTGTTACCCACCAGAATATCCGCCTGa       >  1:1261635/1‑69 (MQ=255)
                                                  aCAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGa       <  1:2221740/68‑1 (MQ=255)
                                                  aCAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGa       <  1:3136759/68‑1 (MQ=255)
                                                      tGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATcgc  >  1:360649/1‑69 (MQ=255)
                                                      tGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATcgc  >  1:3272479/1‑69 (MQ=255)
                                                      |                                                                    
ACGCGGAATCACTTTCTCCACCAGGCCGCCCATCGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGC  >  minE/1690389‑1690511

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: