Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,353,879 C→T 19.2% G179D (GGT→GAT)  yihF ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,353,8790CT19.2% 101.8 / 13.7 67G179D (GGT→GAT) yihFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (35/19);  new base T (10/3);  total (45/22)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

GGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTCCATACTCTGCAGGGGC  >  minE/2353816‑2353940
                                                               |                                                             
ggTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATc                                                           >  1:3276745/1‑68 (MQ=255)
ggTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATc                                                           >  1:1287216/1‑68 (MQ=255)
ggTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATc                                                           >  1:2460206/1‑68 (MQ=255)
ggTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATc                                                           >  1:2324525/1‑68 (MQ=255)
ggTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATc                                                           >  1:2763871/1‑68 (MQ=255)
         tGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTATCATCACCATTaaa                                                  >  1:2209233/1‑68 (MQ=255)
         tGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTATCATCACCATTaaa                                                  >  1:1459763/1‑68 (MQ=255)
         tGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTATCATCACCATTaaa                                                  >  1:149863/1‑68 (MQ=255)
         tGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTATCATCACCATTaaa                                                  >  1:2862945/1‑68 (MQ=255)
         tGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTATCATCACCATTaaa                                                  >  1:1988152/1‑68 (MQ=255)
         tGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTATCATCACCATTaaa                                                  >  1:2716734/1‑68 (MQ=255)
         tGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTATCATCACCATTaaa                                                  >  1:1808067/1‑68 (MQ=255)
         tGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTATCATCACCATTaaa                                                  >  1:716749/1‑68 (MQ=255)
         tGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTATCATCACCATTaaa                                                  >  1:728739/1‑68 (MQ=255)
         tGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTATCATCACCATTaaa                                                  >  1:1135725/1‑68 (MQ=255)
             gAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTa                                              <  1:2111205/68‑1 (MQ=255)
             gAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTa                                              <  1:2895000/68‑1 (MQ=255)
             gAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTa                                              <  1:3434297/68‑1 (MQ=255)
             gAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTa                                              <  1:3430557/68‑1 (MQ=255)
             gAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTa                                              <  1:2266583/68‑1 (MQ=255)
             gAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTa                                              <  1:2489389/68‑1 (MQ=255)
               aCAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAAc                                           >  1:1484175/1‑69 (MQ=255)
               aCAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAAc                                           >  1:140490/1‑69 (MQ=255)
               aCAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAAc                                           >  1:2548401/1‑69 (MQ=255)
                cAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACt                                          >  1:1088276/1‑69 (MQ=255)
                cAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACt                                          >  1:1722637/1‑69 (MQ=255)
                     tCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTa                                              <  1:523928/60‑1 (MQ=255)
                       ttCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTc                                   >  1:1523450/1‑69 (MQ=255)
                               tGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCaaaa                            >  1:2894656/1‑68 (MQ=255)
                               tGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCaaaa                            >  1:1257553/1‑68 (MQ=255)
                               tGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCaaaa                            >  1:569299/1‑68 (MQ=255)
                               tGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCaaaa                            >  1:2408086/1‑68 (MQ=255)
                               tGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCaaaa                            >  1:638780/1‑68 (MQ=255)
                                    tCAATATCCAGGTTAGACAATGAACTATCATCACCATTAAAAGTAAACTGGCCTTCACCAAAAGAAAAc                      <  1:1799489/69‑1 (MQ=255)
                                    tCAATATCCAGGTTAGACAATGAACTATCATCACCATTAAAAGTAAACTGGCCTTCACCAAAAGAAAAc                      <  1:2849461/69‑1 (MQ=255)
                                    tCAATATCCAGGTTAGACAATGAACTATCATCACCATTAAAAGTAAACTGGCCTTCACCAAAAGAAAAc                      <  1:1327696/69‑1 (MQ=255)
                                    tCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAAc                      <  1:1652279/69‑1 (MQ=255)
                                    tCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAAc                      <  1:3156604/69‑1 (MQ=255)
                                         atCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGc                                       >  1:30190/1‑47 (MQ=255)
                                         atCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTc                  >  1:961181/1‑68 (MQ=255)
                                         atCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTc                  >  1:43319/1‑68 (MQ=255)
                                         atCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTc                  >  1:2942704/1‑68 (MQ=255)
                                         atCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTc                  >  1:1134145/1‑68 (MQ=255)
                                         atCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTc                  >  1:237129/1‑68 (MQ=255)
                                         atCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTc                  >  1:1719204/1‑68 (MQ=255)
                                          tCCAGTTTAGACAATGAACTACCATCACCATTAAAAGTa                                              >  1:1395963/1‑39 (MQ=255)
                                          tCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTa                                              >  1:2664669/1‑39 (MQ=255)
                                          tCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTa                                              >  1:1318911/1‑39 (MQ=255)
                                          tCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTa                                              >  1:481147/1‑39 (MQ=255)
                                          tCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTa                                              >  1:1610933/1‑39 (MQ=255)
                                          tCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACtt                    >  1:115642/1‑65 (MQ=255)
                                          tCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACtt                    >  1:3409306/1‑65 (MQ=255)
                                          tCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACt                     >  1:1816481/1‑64 (MQ=255)
                                          tCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTcc                 >  1:1397082/1‑68 (MQ=255)
                                          tCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTCCa                <  1:202638/69‑1 (MQ=255)
                                          tCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTCCa                <  1:2947564/69‑1 (MQ=255)
                                          tCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTCCa                <  1:646161/69‑1 (MQ=255)
                                          tCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTCCa                >  1:3252870/1‑69 (MQ=255)
                                          tCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTCCa                >  1:2799432/1‑69 (MQ=255)
                                          tCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTCGccc                                     >  1:2973677/1‑48 (MQ=255)
                                            cAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTCCATa              <  1:520434/69‑1 (MQ=255)
                                                      aaTGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTCCAt               <  1:1320829/58‑1 (MQ=255)
                                                      aaTGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTCCAt               <  1:3296694/58‑1 (MQ=255)
                                                       aTGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTc                  <  1:2273837/54‑1 (MQ=255)
                                                        tGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTCCATACTCTGCAGGGGc  <  1:754495/69‑1 (MQ=255)
                                                         gAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTCCATACTCTGCAGGGGc  <  1:3366989/68‑1 (MQ=255)
                                                         gAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTCCATACTCTGCAGGGGc  <  1:1605849/68‑1 (MQ=255)
                                                               ccATCACCATTAAAAGTAAAGTGGCCCTCACCAAAAGAAAACtt                    <  1:2390555/44‑1 (MQ=255)
                                                               ccATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACtt                    <  1:1396340/44‑1 (MQ=255)
                                                               ccATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACtt                    >  1:1392307/1‑44 (MQ=255)
                                                               ccATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACtt                    >  1:3273506/1‑44 (MQ=255)
                                                               ccATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACtt                    <  1:1370853/44‑1 (MQ=255)
                                                               ccATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACtt                    <  1:53952/44‑1 (MQ=255)
                                                               ccATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACtt                    <  1:1462253/44‑1 (MQ=255)
                                                               ccATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACtt                    >  1:3154791/1‑44 (MQ=255)
                                                               ccATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACtt                    >  1:1236877/1‑44 (MQ=255)
                                                               ccATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACtt                    <  1:933989/44‑1 (MQ=255)
                                                               ccATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACtt                    >  1:264043/1‑44 (MQ=255)
                                                               |                                                             
GGTGATAATTGCAGAACAATATCTTCGACTTTGCCTTCAATATCCAGGTTAGACAATGAACTACCATCACCATTAAAAGTAAACTGGCCCTCACCAAAAGAAAACTTTCCATACTCTGCAGGGGC  >  minE/2353816‑2353940

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: