Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 663,712 C→G 34.9% V196V (GTG→GTC ycbX ← predicted 2Fe‑2S cluster‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE663,7120CG34.9% 55.2 / 49.9 63V196V (GTG→GTCycbXpredicted 2Fe‑2S cluster‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (22/19);  new base G (10/12);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

GGCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACC  >  minE/663660‑663772
                                                    |                                                            
ggCTGACGGTGGTGAAAATACAGCGGCTACAGGGCTTAACCACATCAAACACGACATCACCAATGCGa                                               >  1:1423506/1‑68 (MQ=255)
ggCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCAATGCGa                                               >  1:3320546/1‑68 (MQ=255)
ggCTGACGGTGGTGAAAATACAGCGGCTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGa                                               >  1:1384511/1‑68 (MQ=255)
ggCTGACGGTGGTGAAAATACAGCGGCTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGa                                               >  1:1489647/1‑68 (MQ=255)
ggCTGACGGTGGTGAAAATACAGCGGCTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGa                                               >  1:3121254/1‑68 (MQ=255)
ggCTGACGGTGGTGAAAATACAGCGGCTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGa                                               >  1:2461120/1‑68 (MQ=255)
ggCTGACGGTGGTGAAAATACAGCGGCTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGa                                               >  1:1576594/1‑68 (MQ=255)
ggCTGACGGTGGTGAAAATACAGCGGCTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGa                                               >  1:2837333/1‑68 (MQ=255)
ggCTGACGGTGGTGAAAATACAGCGGCTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGa                                               >  1:368589/1‑68 (MQ=255)
ggCTGACGGTGGTGAAAATACAGCGGCTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGa                                               >  1:1678059/1‑68 (MQ=255)
ggCTGACGGTGATGAAAATACAGCGGCTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGa                                               >  1:146536/1‑68 (MQ=255)
 gCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCAATGCGAAt                                             >  1:1421852/1‑69 (MQ=255)
 gCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCAATGCGAAt                                             >  1:1444002/1‑69 (MQ=255)
 gCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCAATGCGAAt                                             >  1:2914205/1‑69 (MQ=255)
 gCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCAATGCGAAt                                             >  1:2617200/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:2573006/50‑1 (MQ=255)
             gAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCaa                                                    <  1:188971/50‑1 (MQ=255)
             gAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCaa                                                    <  1:1301560/50‑1 (MQ=255)
             gAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCaa                                                    <  1:2403902/50‑1 (MQ=255)
             gAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCaa                                                    <  1:1813054/50‑1 (MQ=255)
                          cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg                    <  1:1040871/69‑1 (MQ=255)
                          cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg                    <  1:234976/69‑1 (MQ=255)
                          cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg                    <  1:299482/69‑1 (MQ=255)
                          cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg                    <  1:3004588/69‑1 (MQ=255)
                          cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg                    <  1:3173732/69‑1 (MQ=255)
                          cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg                    <  1:543608/69‑1 (MQ=255)
                          cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg                    <  1:774075/69‑1 (MQ=255)
                          cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg                    <  1:1307981/69‑1 (MQ=255)
                          cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg                    <  1:2169564/69‑1 (MQ=255)
                          cTACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg                    <  1:172113/69‑1 (MQ=255)
                           tACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg                    <  1:408590/68‑1 (MQ=255)
                           tACAAGGCTTAACCACATCAAACACGACATCACCAATGCGAATCACTTTCCAGCTATCTTCTTCCCAg                    <  1:490881/68‑1 (MQ=255)
                                   ttAACCACATCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGAcgcg           <  1:1642901/69‑1 (MQ=255)
                                   ttAACCACATCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGAcgcg           <  1:2887079/69‑1 (MQ=255)
                                         acaTCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     >  1:449323/1‑69 (MQ=255)
                                         acaTCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     >  1:81254/1‑69 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     <  1:560531/66‑1 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     >  1:1186824/1‑66 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     >  1:3022516/1‑66 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     >  1:575240/1‑66 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     <  1:504776/66‑1 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     >  1:1153882/1‑66 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     <  1:1246982/66‑1 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     <  1:1005316/66‑1 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     >  1:939435/1‑66 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     >  1:1327328/1‑66 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     >  1:1458195/1‑66 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     >  1:3000930/1‑66 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     >  1:1494329/1‑66 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     >  1:1562184/1‑66 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     >  1:2758839/1‑66 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     >  1:179959/1‑66 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     <  1:1839561/66‑1 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     >  1:2026184/1‑66 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     >  1:2037025/1‑66 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     >  1:2374986/1‑66 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGaa     <  1:2167546/66‑1 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAacc  <  1:1425045/69‑1 (MQ=255)
                                            tCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAacc  <  1:2908759/69‑1 (MQ=255)
                                                    |                                                            
GGCTGACGGTGGTGAAAATACAGCGGCTACAAGGTTTAACCACATCAAACACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACC  >  minE/663660‑663772

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: