Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,251,077 T→C 29.8% A72A (GCA→GCG ruvC ← component of RuvABC resolvasome, endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,251,0770TC29.8% 81.2 / 39.7 87A72A (GCA→GCGruvCcomponent of RuvABC resolvasome, endonuclease
Reads supporting (aligned to +/- strand):  ref base T (30/31);  new base C (11/15);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACCGCCGCCACAATCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTGGAACTGGGTGATGATTTCCGTCAC  >  minE/1251013‑1251139
                                                                |                                                              
tCACCGCCGCCACAATCGCCACGCCGCGCGCCTGGCCCAGTTTCATGGCTGAGTCAGCGTTCTTTGCCa                                                            <  1:840289/69‑1 (MQ=255)
tCACCGCCGCCACAATCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCa                                                            <  1:1088867/69‑1 (MQ=255)
tCACCGCCGCCACAATCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCa                                                            <  1:842427/69‑1 (MQ=255)
tCACCGCCGCCACAATCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCa                                                            <  1:3032446/69‑1 (MQ=255)
tCACCGCCGCCACAATCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCa                                                            <  1:2039394/69‑1 (MQ=255)
tCACCGCCGCCACAATCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCa                                                            <  1:678258/69‑1 (MQ=255)
tCACCGCCGCCACAATCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGAGTTCTTTGCCa                                                            <  1:1999646/69‑1 (MQ=255)
tCACCGCCGCCACAATCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAACGTTCTTTGCCa                                                            <  1:42426/69‑1 (MQ=255)
 cACCGCCGCCACAATCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCa                                                            >  1:1150983/1‑68 (MQ=255)
   ccgccgCCACAATCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCa                                                            <  1:2217043/66‑1 (MQ=255)
           acaATCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACtt                                                   >  1:784654/1‑67 (MQ=255)
           acaATCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACtt                                                   >  1:2145472/1‑67 (MQ=255)
           acaATCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTg                                                  >  1:1953025/1‑68 (MQ=255)
           acaATCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTg                                                  >  1:1974571/1‑68 (MQ=255)
           acaATCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTg                                                  >  1:373592/1‑68 (MQ=255)
           acaATCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTg                                                  >  1:3053494/1‑68 (MQ=255)
           acaATCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTg                                                  >  1:808838/1‑68 (MQ=255)
           acaATCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTg                                                  >  1:2972533/1‑68 (MQ=255)
           acaATAGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTg                                                  >  1:3377251/1‑68 (MQ=255)
             aaTCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTc                                               <  1:331657/69‑1 (MQ=255)
             aaTCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTc                                               <  1:2499261/69‑1 (MQ=255)
             aaTCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTc                                               <  1:1883307/69‑1 (MQ=255)
              aTCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTc                                               >  1:1827556/1‑68 (MQ=255)
              aTCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCa                                              >  1:1082699/1‑69 (MQ=255)
              aTCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCa                                              >  1:798615/1‑69 (MQ=255)
              aTCGCCACGCCGCGCGCCTGCCCCAGTTTCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTc                                               >  1:386459/1‑68 (MQ=255)
              aTCGCCACGCCGCGCGCCTGCCCCAGTTTCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTc                                               >  1:177884/1‑68 (MQ=255)
              aTCGCCACGCCGCGCGCCTGCCCCAGTTTCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCa                                              >  1:3260444/1‑69 (MQ=255)
              aTCGCCACGCCGCGCGCCTGCCCCAGTTTCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCa                                              >  1:2949666/1‑69 (MQ=255)
              aTCGCCACGCCGCGCGCCTGCCCCAGTTTCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCa                                              >  1:1239713/1‑69 (MQ=255)
              aTCGCCACGCCGCGCGCCTGCCCCAGTTTCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCa                                              >  1:81966/1‑69 (MQ=255)
              aTCGCCACGCCGCGCGCCTGCCCCAGTTTCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCa                                              >  1:145862/1‑69 (MQ=255)
              aTCGCCACGCCGCGCGCCTGCCCCAGTTTCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCa                                              >  1:658588/1‑69 (MQ=255)
              aTCGCCACGCCGCGCGCCTGCCCCAGTTTCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCa                                              >  1:616535/1‑69 (MQ=255)
              aTCGCCACGCCGCGCGCCTGCCCCAGTTTCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCa                                              >  1:2098883/1‑69 (MQ=255)
              aTCGCCACGCCGCGCGCCTGCCCCAGTTTCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCa                                              >  1:2086793/1‑69 (MQ=255)
                   cACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATgg                                          >  1:3000105/1‑68 (MQ=255)
                   cACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATgg                                          >  1:2977027/1‑68 (MQ=255)
                   cACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATgg                                          >  1:470549/1‑68 (MQ=255)
                   cACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATgg                                          >  1:2667372/1‑68 (MQ=255)
                   cACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATgg                                          >  1:948964/1‑68 (MQ=255)
                   cACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATg                                           >  1:474589/1‑67 (MQ=255)
                       ccgcGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGAGTTGTTCAATGGTGaa                                      >  1:1296165/1‑68 (MQ=255)
                       ccgcGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGaa                                      >  1:2191210/1‑68 (MQ=255)
                       ccgcGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGaa                                      >  1:1201137/1‑68 (MQ=255)
                       ccgcGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGaa                                      >  1:2199024/1‑68 (MQ=255)
                       ccgcGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGaa                                      >  1:885634/1‑68 (MQ=255)
                       ccgcGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGaa                                      >  1:887554/1‑68 (MQ=255)
                           gcgcCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAaata                                   >  1:2537255/1‑67 (MQ=255)
                           gcgcCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAaata                                   >  1:817470/1‑67 (MQ=255)
                           gcgcCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAaata                                   >  1:1784721/1‑67 (MQ=255)
                           gcgcCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAaata                                   >  1:2168732/1‑67 (MQ=255)
                           gcgcCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAaata                                   >  1:2346462/1‑67 (MQ=255)
                                 gccccAGTTTCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCaa                                             <  1:497222/49‑1 (MQ=255)
                                   cccAGTTTCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCaa                                             <  1:16888/49‑1 (MQ=255)
                                    ccAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTgg                         <  1:3270871/68‑1 (MQ=255)
                                    ccAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTgg                         <  1:410884/68‑1 (MQ=255)
                                    ccAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTgg                         <  1:1507718/68‑1 (MQ=255)
                                    ccAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTgg                         <  1:1801997/68‑1 (MQ=255)
                                         ttCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCaa                                             <  1:220693/43‑1 (MQ=38)
                                         ttCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCAATg                                           <  1:3253014/45‑1 (MQ=255)
                                          tCAGGGCTGAGTCAGCGTTCTTCGCCCTAAAGACCTGTTCAATg                                           <  1:3072755/44‑1 (MQ=38)
                                          tCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCAATg                                           <  1:2211311/44‑1 (MQ=38)
                                          tCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCAATg                                           <  1:1500478/44‑1 (MQ=38)
                                          tCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCAATg                                           <  1:2358009/44‑1 (MQ=38)
                                          tCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCAATg                                           <  1:3005390/44‑1 (MQ=38)
                                          tCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCAATg                                           <  1:2962993/44‑1 (MQ=38)
                                          tCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCAATg                                           <  1:2333987/44‑1 (MQ=38)
                                          tCAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCAATg                                           <  1:2240489/44‑1 (MQ=38)
                                           cAGGGCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCAATg                                           <  1:3326591/43‑1 (MQ=38)
                                              ggCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCAATg                                           <  1:2677905/40‑1 (MQ=38)
                                               gCTGAGTCAGCGTTCTTCGCCATAAAGACCTGTTCaa                                             <  1:2789511/37‑1 (MQ=38)
                                                  gAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTGGAACTGGGTGATGAtt          <  1:1681497/69‑1 (MQ=255)
                                                  gAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTGGAACTGGGTGATGAtt          <  1:2426511/69‑1 (MQ=255)
                                                  gAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTGGAACTGGGTGATGAtt          <  1:758639/69‑1 (MQ=255)
                                                       aGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTGGAACTGGGtgat              <  1:1190157/60‑1 (MQ=255)
                                                       aGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTGGAACTGGGtgat              <  1:778089/60‑1 (MQ=255)
                                                       aGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTGGAACTGGGtgat              <  1:611052/60‑1 (MQ=255)
                                                       aGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTGGAACTGGGtgat              <  1:3271125/60‑1 (MQ=255)
                                                       aGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTGGAACTGGGtgat              <  1:3149686/60‑1 (MQ=255)
                                                       aGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTGGAACTGGGtgat              <  1:2918797/60‑1 (MQ=255)
                                                          gTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTGGAACTGGGTGATGATTTCCGTCAc  <  1:1518814/69‑1 (MQ=255)
                                                          gTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTGGAACTGGGTGATGATTTCCGTCAc  <  1:1374081/69‑1 (MQ=255)
                                                          gTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTGGAACTGGGTGATGATTTCCGTCAc  <  1:1734260/69‑1 (MQ=255)
                                                          gTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTGGAACTGGGTGATGATTTCCGTCAc  <  1:1371446/69‑1 (MQ=255)
                                                          gTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTGGAACTGGGTGATGATTTCCGTCAc  <  1:2071545/69‑1 (MQ=255)
                                                          gTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTGGAACTGGGTGATGATTTCCGTCAc  <  1:2211864/69‑1 (MQ=255)
                                                                |                                                              
TCACCGCCGCCACAATCGCCACGCCGCGCGCCTGGCCCAGTTTCAGGGCTGAGTCAGCGTTCTTTGCCATAAAGACTTGTTCAATGGCGAAATAATCAGGCTGGAACTGGGTGATGATTTCCGTCAC  >  minE/1251013‑1251139

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