Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 750,056 T→A 9.7% D1055V (GAT→GTT)  mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE750,0560TA 9.7% 197.1 / 12.8 72D1055V (GAT→GTT) mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base T (25/40);  new base A (3/4);  total (28/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

TCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGC  >  minE/749991‑750123
                                                                 |                                                                   
tCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTgg                                                                 >  1:3830519/1‑70 (MQ=255)
tCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTgg                                                                 >  1:158345/1‑70 (MQ=255)
 ccTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTgg                                                                 >  1:2206455/1‑69 (MQ=255)
  cTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGgcag                                                              <  1:579851/71‑1 (MQ=255)
     aTCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGcagcag                                                           >  1:2493071/1‑71 (MQ=255)
       ccccAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCATCAGGGTACGCGCCGGATCTGGCAGCAGGc                                                         >  1:2833906/1‑71 (MQ=255)
       ccccAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGAAGGGTACGCGCCGGATCTGGCAGCAGGc                                                         >  1:4220747/1‑71 (MQ=255)
               tCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcg                                                  <  1:1086207/70‑1 (MQ=255)
                cTGCGCTTGCTGGCGCAGACGGGCAATATCCAGTAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgat                                                >  1:4314337/1‑71 (MQ=255)
                cTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgat                                                >  1:3497374/1‑71 (MQ=255)
                cTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgat                                                >  1:1477742/1‑71 (MQ=255)
                    gcTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgatcgat                                            <  1:4028773/71‑1 (MQ=255)
                     cTTGCTGGCGCAGCCGGGCAATATCCAGCATGGTACGCGCCGATCTGGCAGCAGGCCGAAGcgatcgat                                            >  1:1619424/1‑69 (MQ=255)
                     cTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgatcgat                                            >  1:4290838/1‑70 (MQ=255)
                     cTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgatcgat                                            >  1:4378522/1‑70 (MQ=255)
                      atGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGAACTGGCAGCAGGCCGAAGcgatcgat                                            >  1:2040698/2‑69 (MQ=255)
                        gCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgat                                                <  1:2038703/63‑1 (MQ=255)
                        gCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgat                                                <  1:4004288/63‑1 (MQ=255)
                        gCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgat                                                <  1:4126279/63‑1 (MQ=255)
                        gCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgat                                                <  1:1219196/63‑1 (MQ=255)
                        gCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCAGAAGcgat                                                <  1:3560254/63‑1 (MQ=255)
                         cTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGc                                        <  1:1783273/70‑1 (MQ=255)
                         cTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGc                                        <  1:664650/70‑1 (MQ=255)
                         cTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGc                                        <  1:1368901/70‑1 (MQ=255)
                            gcgcAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTc                                      <  1:4608564/69‑1 (MQ=255)
                            gcgcAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTc                                      <  1:583222/69‑1 (MQ=255)
                            gcgcAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTc                                      <  1:4722884/69‑1 (MQ=255)
                            gcgcAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGa                                    >  1:4222683/1‑71 (MQ=255)
                              gcAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTg                                                                  >  1:3786614/1‑39 (MQ=255)
                              gcAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTg                                                                  >  1:3424846/1‑39 (MQ=255)
                              gcAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACt                                  <  1:1887862/71‑1 (MQ=255)
                              gcAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACt                                  <  1:2514594/71‑1 (MQ=255)
                                 gACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTg                               <  1:4720836/71‑1 (MQ=255)
                                 gACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTg                               <  1:1459396/71‑1 (MQ=255)
                                  aCGGGCAATATCCAGCAGGGTACGCGCCGGAACTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGa                              >  1:3134455/1‑71 (MQ=255)
                                   cGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGAt                             >  1:1429354/1‑71 (MQ=255)
                                    gggCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgatc                                               >  1:1332964/1‑52 (MQ=255)
                                      gCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGAtctc                          <  1:898081/71‑1 (MQ=255)
                                         atatCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGc                                                   <  1:4365370/43‑1 (MQ=255)
                                         atatCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGc                                                   <  1:3350401/43‑1 (MQ=255)
                                         atatCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGc                                                   <  1:2359014/43‑1 (MQ=255)
                                         atatCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGc                                                   <  1:466646/43‑1 (MQ=255)
                                          tatCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAg                                                    <  1:777271/41‑1 (MQ=255)
                                          tatCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAg                                                    <  1:98425/41‑1 (MQ=255)
                                          tatCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAg                                                    <  1:2515413/41‑1 (MQ=255)
                                          tatCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAg                                                    <  1:3546940/41‑1 (MQ=255)
                                          tatCCAGCAGGGTACGCGCCGGAACTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTc                       <  1:165040/70‑1 (MQ=255)
                                          tatCCAGCAGGGTACGCGCCGGAACTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTc                       <  1:2984516/70‑1 (MQ=255)
                                          tatCCAGCAGGGTACGCGCCGGAACTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTc                       <  1:4249874/70‑1 (MQ=255)
                                          tatCCAGCAGGGTACGCGCCGGAACTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTc                       <  1:3547823/70‑1 (MQ=255)
                                            tCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAg                    >  1:32534/1‑71 (MQ=255)
                                            tCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAg                    >  1:1933175/1‑71 (MQ=255)
                                            tCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAg                    >  1:2756166/1‑71 (MQ=255)
                                               agcagGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTc                 <  1:4291524/71‑1 (MQ=255)
                                               agcagGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTc                 <  1:3749908/71‑1 (MQ=255)
                                                 cagGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCg                                     <  1:3174666/49‑1 (MQ=255)
                                                 cagGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCg                                     <  1:704056/49‑1 (MQ=255)
                                                 cagGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTcc                      <  1:830527/64‑1 (MQ=255)
                                                 cagGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTcc                      <  1:3961581/64‑1 (MQ=255)
                                                 cagGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGt               <  1:3996345/71‑1 (MQ=255)
                                                   gggTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAg                    >  1:4163947/1‑64 (MQ=255)
                                                   gggTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAg                    >  1:211793/1‑64 (MQ=255)
                                                   gggTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAg                    >  1:2610990/1‑64 (MQ=255)
                                                      tACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTcc          >  1:2717532/1‑71 (MQ=255)
                                                      tACGCGCCGGAACTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCa                     >  1:3944606/1‑60 (MQ=255)
                                                        cgcgCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGtttt            <  1:57344/67‑1 (MQ=255)
                                                        cgcgCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGt        <  1:3691290/71‑1 (MQ=255)
                                                         gcgcCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGtt       <  1:4577981/71‑1 (MQ=255)
                                                              ggATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGgc  <  1:4222019/71‑1 (MQ=255)
                                                              ggATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGgc  >  1:445598/1‑71 (MQ=255)
                                                              ggATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGgc  <  1:874376/71‑1 (MQ=255)
                                                              ggATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGgc  <  1:4216589/71‑1 (MQ=255)
                                                                 |                                                                   
TCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGC  >  minE/749991‑750123

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: