Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,809,462 T→C 42.5% I117I (ATT→ATC yjbQ → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,809,4620TC42.5% 24.1 / 66.4 87I117I (ATT→ATCyjbQconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (21/29);  new base C (13/24);  total (34/53)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

GCTGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTA  >  minE/2809396‑2809516
                                                                  |                                                      
gCTGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTg                                                      <  1:2272980/69‑1 (MQ=255)
gCTGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTg                                                      <  1:2965023/69‑1 (MQ=255)
gCTGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTg                                                      <  1:2004336/69‑1 (MQ=255)
    ggAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTg                                                      <  1:2701684/65‑1 (MQ=255)
                 tATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCa                                     >  1:2971754/1‑69 (MQ=255)
                 tATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCa                                     >  1:2134211/1‑69 (MQ=255)
                 tATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCa                                     >  1:2292517/1‑69 (MQ=255)
                 tATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCa                                     >  1:1977246/1‑69 (MQ=255)
                 tATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCa                                     >  1:2494225/1‑69 (MQ=255)
                 tATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCa                                     >  1:1879671/1‑69 (MQ=255)
                 tATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCa                                     >  1:3340382/1‑69 (MQ=255)
                 tATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCa                                     >  1:257611/1‑69 (MQ=255)
                 tATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCa                                     >  1:3261291/1‑69 (MQ=255)
                 tATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCa                                     >  1:3249753/1‑69 (MQ=255)
                 tATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCa                                     >  1:2911468/1‑69 (MQ=255)
                  aTTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGc                                      <  1:587521/67‑1 (MQ=255)
                  aTTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGc                                      <  1:3419436/67‑1 (MQ=255)
                                aaGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTTGGGGAACATCGCATCCACGGCGGATcg                       >  1:2195486/1‑68 (MQ=255)
                                aaGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCg                                       >  1:3195749/1‑52 (MQ=255)
                                aaGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATcgc                      <  1:1748403/69‑1 (MQ=255)
                                aaGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATcgc                      <  1:3431966/69‑1 (MQ=255)
                                aaGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATcgc                      <  1:2286573/69‑1 (MQ=255)
                                aaGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATcgc                      >  1:2945035/1‑69 (MQ=255)
                                aaGGGCGTATTCAGACCGGCACCTGGCAAGCCATTTGGCTGGGGGAACATCGCATCCACGGCGGATcgc                      >  1:1731470/1‑69 (MQ=255)
                                         ttaaGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGcatca             >  1:1414539/4‑69 (MQ=255)
                                         ttCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGa                          >  1:370622/1‑56 (MQ=255)
                                         ttCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGcatca             >  1:973314/1‑69 (MQ=255)
                                           cAGACCGGGACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:1976152/50‑1 (MQ=255)
                                           cAGACCGGCACCTGTCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:1260142/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:2638290/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:3161299/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:3183919/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:3195214/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:1335963/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:825765/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:773751/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:2745505/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:708830/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:1233016/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:1667712/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:2185894/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:3299407/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:2634402/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:331936/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:627282/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:253137/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:1917563/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:1953868/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:2340795/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:2066332/50‑1 (MQ=255)
                                           cAGACCGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCAcgg                              <  1:457841/50‑1 (MQ=255)
                                               ccGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGa                          >  1:682211/1‑50 (MQ=255)
                                               ccGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       <  1:3199682/69‑1 (MQ=255)
                                               ccGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       >  1:83575/1‑69 (MQ=255)
                                               ccGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       <  1:2738483/69‑1 (MQ=255)
                                               ccGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       <  1:463439/69‑1 (MQ=255)
                                               ccGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       <  1:1855536/69‑1 (MQ=255)
                                               ccGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       <  1:2995221/69‑1 (MQ=255)
                                               ccGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       <  1:1289805/69‑1 (MQ=255)
                                               ccGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       <  1:1781860/69‑1 (MQ=255)
                                               ccGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       <  1:2742603/69‑1 (MQ=255)
                                               ccGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       <  1:2243681/69‑1 (MQ=255)
                                               ccGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       >  1:2564115/1‑69 (MQ=255)
                                               ccGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCg                  >  1:2040062/1‑58 (MQ=255)
                                               ccGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       >  1:263749/1‑69 (MQ=255)
                                               ccGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       >  1:48934/1‑69 (MQ=255)
                                               ccGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       >  1:3152299/1‑69 (MQ=255)
                                               ccGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       >  1:2377861/1‑69 (MQ=255)
                                               ccGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       >  1:1812612/1‑69 (MQ=255)
                                               ccGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       >  1:2599687/1‑69 (MQ=255)
                                               ccGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       >  1:2770650/1‑69 (MQ=255)
                                               ccGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       >  1:1486858/1‑69 (MQ=255)
                                               ccGGCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       >  1:1252390/1‑69 (MQ=255)
                                               ccGGCACCTGGAAAGGCATCTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       >  1:3294755/1‑69 (MQ=255)
                                               caggCAGCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGa       >  1:2166913/3‑69 (MQ=255)
                                               caggCACCTGGCAAGGCATCTGGCTGGGGGAACATCGCATCCACGGCGGa                          >  1:2825379/3‑50 (MQ=255)
                                                cGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGac      <  1:814736/69‑1 (MQ=255)
                                                    aCCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTa  <  1:539246/69‑1 (MQ=255)
                                                    aCCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTa  <  1:2459341/69‑1 (MQ=255)
                                                    aCCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTa  <  1:67876/69‑1 (MQ=255)
                                                    aCCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTa  <  1:2761415/69‑1 (MQ=255)
                                                    aCCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTa  <  1:708373/69‑1 (MQ=255)
                                                    aCCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTa  <  1:1528380/69‑1 (MQ=255)
                                                    aCCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTa  <  1:76859/69‑1 (MQ=255)
                                                    aCCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTa  <  1:3196477/69‑1 (MQ=255)
                                                    aCCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTa  <  1:3077599/69‑1 (MQ=255)
                                                     ccTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTa  <  1:1346503/68‑1 (MQ=255)
                                                                  |                                                      
GCTGGGAACATCGCTTGTATTGCCGGTGCATAAAGGGCGTATTCAGACCGGCACCTGGCAAGGCATTTGGCTGGGGGAACATCGCATCCACGGCGGATCGCGTCGCATCATCGCGACACTA  >  minE/2809396‑2809516

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: