Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,091,167 T→C 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,1670TC97.2% 127.3 / ‑3.0 36coding (534/1488 nt)sufBcomponent of SufBCD complex
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (15/20);  minor base . (1/0);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGC  >  minE/1091100‑1091234
                                                                   |                                                                   
ccgcATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTcc                                                                  <  1:2480419/71‑1 (MQ=255)
            cgGCAAAGAAGTTGTCATTCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTg                                                       <  1:2066330/69‑1 (MQ=255)
              gCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGAt                                                    <  1:2688863/71‑1 (MQ=255)
                     agTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCa                                             >  1:722758/1‑71 (MQ=255)
                      gTTGTCATTCCCCGGCACCACGGTGCCGAGGTATTTACGCACCAGCTCCGGGTGATCGt                                                        >  1:366913/1‑59 (MQ=255)
                      gTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGt                                                        >  1:801265/1‑59 (MQ=255)
                      gTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTgg                                                      <  1:2941519/61‑1 (MQ=255)
                      gTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCAc                                            <  1:2542999/71‑1 (MQ=255)
                         gTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCa                                          >  1:1643969/1‑70 (MQ=255)
                         gTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCa                                          >  1:1647526/1‑70 (MQ=255)
                           cATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAG‑‑CTGGGTGATCGTGGATc                                                   >  1:929244/1‑57 (MQ=255)
                            aTTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAgg                                      <  1:740279/71‑1 (MQ=255)
                             ttCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGa                                     <  1:1112898/71‑1 (MQ=255)
                             ttCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGa                                     <  1:551904/71‑1 (MQ=255)
                               ccccGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGaa                                    >  1:275844/1‑70 (MQ=255)
                               ccccGGCACCACGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGaa                                    >  1:1896878/1‑70 (MQ=255)
                                  cGGCACCACGGTGCCGAGATATTTACGCACCAGCGCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAg                                 >  1:2682070/1‑70 (MQ=255)
                                     caccacGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGaaaa                             <  1:1584480/71‑1 (MQ=255)
                                       ccacGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaata                           <  1:2398778/71‑1 (MQ=255)
                                        cacGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCa                                          <  1:649791/55‑1 (MQ=255)
                                        cacGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCa                                          <  1:1628929/55‑1 (MQ=255)
                                          cGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGAt                                                    >  1:1830218/1‑43 (MQ=255)
                                          cGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATc                                                   <  1:534803/44‑1 (MQ=255)
                                          cGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaataat                         >  1:3010105/1‑70 (MQ=255)
                                          cGGTGCCGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAAtt                        <  1:17789/71‑1 (MQ=255)
                                               ccGAGATATTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTg                   >  1:1076889/1‑71 (MQ=255)
                                                      atTTACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGc             >  1:799574/1‑70 (MQ=255)
                                                       tttACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaat                            <  1:1815033/54‑1 (MQ=255)
                                                        ttACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAg          >  1:2984787/1‑71 (MQ=255)
                                                        ttACGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAg          >  1:731722/1‑71 (MQ=255)
                                                           cGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGttt       <  1:1443742/71‑1 (MQ=255)
                                                           cGCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGttt       <  1:1184032/71‑1 (MQ=255)
                                                            gCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCg              <  1:2623944/63‑1 (MQ=255)
                                                            gCACCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCg              <  1:130236/63‑1 (MQ=255)
                                                              aCCAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTc    >  1:2799985/1‑71 (MQ=255)
                                                                cAGCTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTcgc  <  1:1397656/71‑1 (MQ=255)
                                                                   |                                                                   
CCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGC  >  minE/1091100‑1091234

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: