Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 663,694 T→C 51.2% K202K (AAA→AAG ycbX ← predicted 2Fe‑2S cluster‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE663,6940TC51.2% ‑2.8 / 35.6 43K202K (AAA→AAGycbXpredicted 2Fe‑2S cluster‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (5/16);  new base C (10/12);  total (15/28)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

CTGCCGGATGTTTTTGCCCTTTTTCTGGGCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAG  >  minE/663633‑663754
                                                             |                                                            
cTGCCGGATGTTTTTGCCCTTTTTCTGGGCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACca                                                       <  1:2615275/69‑1 (MQ=255)
cTGCCGGATGTTTTTGCCCTTTTTCTGGGCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACca                                                       <  1:3011569/69‑1 (MQ=255)
     ggATGTTTTTGCCCTTTTTCTGGGCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCa                                                  <  1:2376064/69‑1 (MQ=255)
     ggATGTTTTTGCCCTTTTTCTGGGCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCa                                                  <  1:1882323/69‑1 (MQ=255)
     ggATGTTTTTGCCCTTTTTCTGGGCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCa                                                  <  1:277896/69‑1 (MQ=255)
     ggATGTTTTTGCCCTTTTTCTGGGCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCa                                                  <  1:1311071/69‑1 (MQ=255)
     ggATGTTTTTGCCCTTTTTCTGGGCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCa                                                  <  1:530751/69‑1 (MQ=255)
                           ggCTGACGGTGGTGAAAATACAGCGGCTACAGGGCTTAACCACATCAAACACGACATCACCAATGCGa                             >  1:1423506/1‑68 (MQ=255)
                           ggCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCAATGCGa                             >  1:3320546/1‑68 (MQ=255)
                           ggCTGACGGTGGTGAAAATACAGCGGCTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGa                             >  1:1384511/1‑68 (MQ=255)
                           ggCTGACGGTGGTGAAAATACAGCGGCTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGa                             >  1:2461120/1‑68 (MQ=255)
                           ggCTGACGGTGGTGAAAATACAGCGGCTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGa                             >  1:1489647/1‑68 (MQ=255)
                           ggCTGACGGTGGTGAAAATACAGCGGCTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGa                             >  1:1576594/1‑68 (MQ=255)
                           ggCTGACGGTGGTGAAAATACAGCGGCTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGa                             >  1:1678059/1‑68 (MQ=255)
                           ggCTGACGGTGGTGAAAATACAGCGGCTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGa                             >  1:3121254/1‑68 (MQ=255)
                           ggCTGACGGTGGTGAAAATACAGCGGCTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGa                             >  1:368589/1‑68 (MQ=255)
                           ggCTGACGGTGGTGAAAATACAGCGGCTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGa                             >  1:2837333/1‑68 (MQ=255)
                           ggCTGACGGTGATGAAAATACAGCGGCTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGa                             >  1:146536/1‑68 (MQ=255)
                            gCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCAATGCGAAt                           >  1:2617200/1‑69 (MQ=255)
                            gCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCAATGCGAAt                           >  1:2914205/1‑69 (MQ=255)
                            gCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCAATGCGAAt                           >  1:1444002/1‑69 (MQ=255)
                            gCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCAATGCGAAt                           >  1:1421852/1‑69 (MQ=255)
                             cTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCAATGCGAAtg                          <  1:1405869/69‑2 (MQ=255)
                             cTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCAATGCGAATc                          <  1:2737241/69‑1 (MQ=255)
                                    tggtgAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCAATGCGAATCACTTTcc                   <  1:933206/69‑1 (MQ=255)
                                    tggtgAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCAATGCGAATCACTTTcc                   <  1:1277827/69‑1 (MQ=255)
                                        gAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCaa                                  <  1:2403902/50‑1 (MQ=255)
                                        gAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCaa                                  <  1:2573006/50‑1 (MQ=255)
                                        gAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCaa                                  <  1:1301560/50‑1 (MQ=255)
                                        gAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCaa                                  <  1:188971/50‑1 (MQ=255)
                                        gAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCaa                                  <  1:1813054/50‑1 (MQ=255)
                                                     cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg  <  1:1040871/69‑1 (MQ=255)
                                                     cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg  <  1:774075/69‑1 (MQ=255)
                                                     cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg  <  1:543608/69‑1 (MQ=255)
                                                     cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg  <  1:1307981/69‑1 (MQ=255)
                                                     cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg  <  1:172113/69‑1 (MQ=255)
                                                     cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg  <  1:3173732/69‑1 (MQ=255)
                                                     cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg  <  1:3004588/69‑1 (MQ=255)
                                                     cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg  <  1:299482/69‑1 (MQ=255)
                                                     cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg  <  1:2169564/69‑1 (MQ=255)
                                                     cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg  <  1:234976/69‑1 (MQ=255)
                                                      tACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg  <  1:408590/68‑1 (MQ=255)
                                                      tACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg  <  1:490881/68‑1 (MQ=255)
                                                             |                                                            
CTGCCGGATGTTTTTGCCCTTTTTCTGGGCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAG  >  minE/663633‑663754

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: