Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,785,788 A→C 32.9% V79V (GTA→GTC yjbC → 23S rRNA pseudouridine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,785,7880AC32.9% 58.0 / 54.8 73V79V (GTA→GTCyjbC23S rRNA pseudouridine synthase
Reads supporting (aligned to +/- strand):  ref base A (24/25);  new base C (21/3);  total (45/28)
Fisher's exact test for biased strand distribution p-value = 1.83e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

TTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCAAACGCGTGT  >  minE/2785723‑2785852
                                                                 |                                                                
ttgattgaACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGc                                                               <  1:1584315/69‑1 (MQ=255)
  gattgaACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTa                                                                  >  1:426872/1‑64 (MQ=255)
  gattgaACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTa                                                                  >  1:1521659/1‑64 (MQ=255)
  gattgaACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGcac                                                             <  1:2227778/69‑1 (MQ=255)
  gattgaACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGcac                                                             <  1:2323389/69‑1 (MQ=255)
  gattgaACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGcac                                                             <  1:1767847/69‑1 (MQ=255)
  gattgaACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGcac                                                             <  1:2562653/69‑1 (MQ=255)
      gaACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGcaccac                                                          >  1:3351117/1‑68 (MQ=255)
      gaACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGcaccac                                                          >  1:1214044/1‑68 (MQ=255)
      gaACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGcaccac                                                          >  1:942451/1‑68 (MQ=255)
      gaACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGcaccac                                                          >  1:2608818/1‑68 (MQ=255)
         ccTCGTGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa                                                      <  1:1739771/69‑1 (MQ=255)
         ccTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa                                                      <  1:159812/69‑1 (MQ=255)
         ccTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa                                                      <  1:2616796/69‑1 (MQ=255)
         ccTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa                                                      <  1:1687823/69‑1 (MQ=255)
         ccTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa                                                      <  1:2681996/69‑1 (MQ=255)
         ccTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa                                                      <  1:2830762/69‑1 (MQ=255)
         ccTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa                                                      <  1:3017630/69‑1 (MQ=255)
          cTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa                                                      <  1:1290272/68‑1 (MQ=255)
                            tGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGcac                                                             <  1:2862656/43‑1 (MQ=255)
                            tGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGATGGCGAGCGCGATAAc                                    <  1:3275101/68‑1 (MQ=255)
                                   tATCGCCCTGAACAAGCCCGTTGGTATTGTAAGcaccac                                                          >  1:2393826/1‑39 (MQ=255)
                                   tATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa                                                      <  1:1663428/43‑1 (MQ=255)
                                   tATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTc                              >  1:55721/1‑67 (MQ=255)
                                   tATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTc                              >  1:370266/1‑67 (MQ=255)
                                   tATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTc                              >  1:1917950/1‑67 (MQ=255)
                                    aTCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa                                                      >  1:606010/1‑42 (MQ=255)
                                    aTCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa                                                      >  1:3430934/1‑42 (MQ=255)
                                    aTCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa                                                      >  1:847955/1‑42 (MQ=255)
                                    aTCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa                                                      >  1:1747717/1‑42 (MQ=255)
                                    aTCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa                                                      >  1:1445491/1‑42 (MQ=255)
                                    aTCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa                                                      >  1:127114/1‑42 (MQ=255)
                                    aTCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa                                                      >  1:438857/1‑42 (MQ=255)
                                    aTCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa                                                      >  1:1233715/1‑42 (MQ=255)
                                    aTCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa                                                      >  1:2094155/1‑42 (MQ=255)
                                    aTCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGAt                           >  1:1213961/1‑69 (MQ=255)
                                    aTCGACCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGaa                                                      >  1:2905320/1‑42 (MQ=255)
                                          cTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGtt                     <  1:2446178/69‑1 (MQ=255)
                                          cTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGtt                     <  1:2714635/69‑1 (MQ=255)
                                             aaCAAGCCGGTTGGTATTGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTa                    <  1:1690489/67‑1 (MQ=255)
                                               cAAGCCCGTTGGTATTGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAAc                  <  1:1899039/67‑1 (MQ=255)
                                                    ccGTTGGTATTGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGct           <  1:349813/69‑2 (MQ=255)
                                                    ccGTTGGTATTGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCa           <  1:1988097/69‑1 (MQ=255)
                                                    ccGTTGGTATTGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCa           <  1:595345/69‑1 (MQ=255)
                                                    ccGTTGGTATTGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCa           <  1:2163236/69‑1 (MQ=255)
                                                    ccGTTGGTATTGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCa           <  1:498957/69‑1 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCaa          >  1:87665/1‑69 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCaa          >  1:369363/1‑69 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCaa          >  1:785975/1‑69 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCaa          >  1:2468416/1‑69 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCaa          >  1:961727/1‑69 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCaa          >  1:3039264/1‑69 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCaa          >  1:1832416/1‑69 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCaa          >  1:2677096/1‑69 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCaa          >  1:2673105/1‑69 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCaa          >  1:1091253/1‑69 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCaa          >  1:2450322/1‑69 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCaa          >  1:2362168/1‑69 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCaa          >  1:1836584/1‑69 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCaa          >  1:1904312/1‑69 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCaa          >  1:2109430/1‑69 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCaa          >  1:1571065/1‑69 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCaa          >  1:1977943/1‑69 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCa           >  1:1739530/1‑68 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCa           >  1:73886/1‑68 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCa           >  1:578133/1‑68 (MQ=255)
                                                     cGTCGGTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCa           >  1:424997/1‑68 (MQ=255)
                                                     cGTCGGTATTGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCaa          >  1:3281411/1‑69 (MQ=255)
                                                         ggTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCAAAcgc      <  1:592724/69‑1 (MQ=255)
                                                         ggTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCAAAcgc      <  1:478619/69‑1 (MQ=255)
                                                         ggTATTGTCAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCAAAcgc      <  1:144632/69‑1 (MQ=255)
                                                             ttGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCAAACGCgtgt  >  1:1134124/1‑69 (MQ=255)
                                                             ttGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCAAACGCgtgt  >  1:542488/1‑69 (MQ=255)
                                                             ttGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCAAACGCgtgt  >  1:2311684/1‑69 (MQ=255)
                                                             ttGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCAAACGCgtgt  >  1:3369731/1‑69 (MQ=255)
                                                                 |                                                                
TTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGATGGCGAGCGCGATAACATTGTCGATTTCGTTAACCACAGCAAACGCGTGT  >  minE/2785723‑2785852

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: