Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,091,167 T→C 52.1% E178E (GAA→GAG sufB ← component of SufBCD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,091,1670TC52.1% 1.5 / 67.1 48E178E (GAA→GAGsufBcomponent of SufBCD complex
Reads supporting (aligned to +/- strand):  ref base T (12/11);  new base C (15/10);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 7.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

CCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGAT  >  minE/1091100‑1091237
                                                                   |                                                                      
ccgcATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTcc                                                                     <  1:2071346/71‑1 (MQ=255)
  gcATTAAGCGCGTCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCgg                                                                   <  1:2060041/71‑1 (MQ=255)
  gcATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCgg                                                                   <  1:2625242/71‑1 (MQ=255)
  gcATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCgg                                                                   <  1:1359093/71‑1 (MQ=255)
       aaGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGAt                                                              >  1:676362/1‑71 (MQ=255)
       aaGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGAt                                                              >  1:225160/1‑71 (MQ=255)
          cgcgGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGt                                                           <  1:2702694/71‑1 (MQ=255)
          cgcgGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGt                                                           <  1:1617315/71‑1 (MQ=255)
               cAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATc                                                      <  1:2078114/71‑1 (MQ=255)
               cAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATc                                                      <  1:1978844/71‑1 (MQ=255)
                  agaagTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGcc                                                   >  1:3530986/1‑71 (MQ=255)
                        tGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCAcc                                              >  1:2514481/1‑70 (MQ=255)
                         gTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCg                                                            >  1:2676701/1‑55 (MQ=255)
                         gTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCgg                                                                   >  1:1744903/1‑48 (MQ=255)
                          tCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGAt                                                              <  1:1047754/52‑1 (MQ=255)
                            aTTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAgg                                         <  1:1419729/71‑1 (MQ=255)
                            aTTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGAt                                                       <  1:3503460/57‑1 (MQ=255)
                            aTTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAg                                          <  1:2451826/70‑1 (MQ=255)
                             ttCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCa                                                >  1:1770516/1‑63 (MQ=255)
                             ttCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCa                                                >  1:1812430/1‑63 (MQ=255)
                             ttCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGa                                        <  1:1128103/71‑1 (MQ=255)
                                  cGGCACCGTGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAGCAgg                                   >  1:2015725/1‑70 (MQ=255)
                                  cGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGa                                   <  1:1068632/71‑1 (MQ=255)
                                     caccacGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGaaaa                                <  1:2514225/71‑1 (MQ=255)
                                         gtggTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCaa                                            <  1:3203599/53‑1 (MQ=255)
                                         acGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaataat                            >  1:494505/1‑71 (MQ=255)
                                          cGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAAtt                           >  1:1030784/1‑71 (MQ=255)
                                           ggTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCAcc                                              >  1:1490792/1‑51 (MQ=255)
                                           ggTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTc                          >  1:1011953/1‑71 (MQ=255)
                                             tGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAAc                                      <  1:1328921/57‑1 (MQ=255)
                                             tGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTccc                        <  1:3258934/71‑1 (MQ=255)
                                                cGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAAc                                      <  1:2494855/54‑1 (MQ=255)
                                                 gagaTATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCt                    >  1:3232024/1‑71 (MQ=255)
                                                 gagaTATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCt                    >  1:535695/1‑71 (MQ=255)
                                                   gaTATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTcc                  >  1:930958/1‑71 (MQ=255)
                                                   gaTATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAAc                                      >  1:2012162/1‑51 (MQ=255)
                                                   gaTATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAAc                                      >  1:366854/1‑51 (MQ=255)
                                                   gaTATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAAc                                      >  1:1761208/1‑51 (MQ=255)
                                                   gaTATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTcc                  >  1:974280/1‑71 (MQ=255)
                                                    atatTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCg                 >  1:231193/1‑71 (MQ=255)
                                                       tttACGCACCAGCTCCGGGTGATCGTGGATCGCCTCa                                                >  1:3258570/1‑37 (MQ=255)
                                                         tACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGt            >  1:43620/1‑71 (MQ=255)
                                                          aCGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAg                                          >  1:3016528/1‑40 (MQ=255)
                                                          aCGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAg                                          >  1:3304539/1‑40 (MQ=255)
                                                          aCGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAg                                          >  1:836133/1‑40 (MQ=255)
                                                            gCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGtttt         <  1:344388/71‑1 (MQ=255)
                                                               ccAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAAtt                           <  1:2719848/50‑1 (MQ=255)
                                                                cAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTcgc     >  1:11333/1‑71 (MQ=255)
                                                                   ttCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGAt  <  1:3595851/71‑1 (MQ=255)
                                                                   ttCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGAt  <  1:1907374/71‑1 (MQ=255)
                                                                   |                                                                      
CCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGAT  >  minE/1091100‑1091237

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: