Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,363,569 A→C 25.3% A353A (GCA→GCC metG → methionyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,363,5690AC25.3% 94.7 / 33.6 71A353A (GCA→GCCmetGmethionyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (24/29);  new base C (11/7);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATGTCCAAGTCTCGCGGCACCTTTATTAAAGCCAGCACCTGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATTGATGATATCGATCTCAACCTGG  >  minE/1363508‑1363636
                                                             |                                                                   
aaaTGTCCAAGTCTCGCGGCACCTTTATTAAAGCCAGCACCTGGCTGAATCATTTTGACGCCGACAGcc                                                              >  1:1959627/3‑69 (MQ=255)
aaaTGTCCAAGTCTCGCGGCACCTTTATTAAAGCCAGCACCTGGCTGAATCATTTTGACGCCGACAGcc                                                              >  1:2582510/3‑69 (MQ=255)
aaaTGTCCAAGTCTCGCGGCACCTTTATTAAAGCCAGCACCTGGCTGAATCATTTTGACGCCGACAGcc                                                              >  1:1300954/3‑69 (MQ=255)
aaaTGTCCAAGTCTCGCGGCACCTTTATTAAAGCCAGCACCTGGCTGAATCATTTTGACGCCGACAGcc                                                              >  1:1347536/3‑69 (MQ=255)
aaaTGTCCAAGTCTCGCGGCACCTTTATTAAAGCCAGCACCTGGCTGAATCATTTTGACGCCGACAGcc                                                              >  1:2505751/3‑69 (MQ=255)
aaaTGTCCAAGTCTCGCGGCACCTTTATTAAAGCCAGCACCTGGCTGAATCATTTTGACGCCGACAGcc                                                              >  1:236707/3‑69 (MQ=255)
aaaTGTCCAAGTCTCGCGGCACCTTTATTAAAGCCAGCACCTGGCTGAATCATTTTGACGCCGACAGcc                                                              >  1:3405702/3‑69 (MQ=255)
aaaTGTCCAAGTCTCGCGGCACCTTTATTAAAGCCAGCACCTGGCTGAATCATTTTGACGCCGACAGc                                                               >  1:2371615/3‑68 (MQ=255)
aaaTGTCCAAGTCTCGCGGCACCTTTATTAAAGCCAGCACCTGGCTGAATCATTTTGACGCCGACAGc                                                               >  1:1790585/3‑68 (MQ=255)
aaaTGTCCAAGTCTCGCGGCACCTTTATTAAAGCCAGCACCTGGCTGAATCATTTTGACGCCGACAGc                                                               >  1:1696437/3‑68 (MQ=255)
aaaTGTCCAAGTCTCGCAGCACCTTTATTAAAGCCAGCACCTGGCTGAATCATTTTGACGCCGACAGcc                                                              >  1:2040622/3‑69 (MQ=255)
aGATGTCCAAGTCTCGCGGCACCTTTATTAAAGCCAGCACCTGGCTGAATCATTTTGACGCAGACAGcc                                                              <  1:440392/69‑1 (MQ=255)
aGATGTCCAAGTCTCGCGGCACCTTTATTAAAGCCAGCACCTGGCTGAATCATTTTGACGCAGACAGcc                                                              <  1:816966/69‑1 (MQ=255)
aGATGTCCAAGTCTCGCGGCACCTTTATTAAAGCCAGCACCTGGCTGAATCATTTTGACGCAGACAGcc                                                              <  1:2803512/69‑1 (MQ=255)
     tCCAAGTCTCGCGGCACCTTTATTAAAGCCAGCACCTGGCTGAATCATTTTGACGCAGACAGCCTGCGt                                                         <  1:1751602/69‑1 (MQ=255)
     tCCAAGTCTCGCGGCACCTTTATTAAAGCCAGCACCTGGCTGAATCATTTTGACGCAGACAGCCTGCGt                                                         <  1:492309/69‑1 (MQ=255)
                                      aCCTGGCTGAATCATTTTGACGCAGAGAGCCTGCGTtactactac                                                <  1:1788614/45‑1 (MQ=255)
                                      aCCTGGCTGAATCATTTTGACGCAGACAGCCTGCGTtactactac                                                >  1:2497810/1‑45 (MQ=255)
                                      aCCTGGCTGAATCATTTTGACGCAGACAGCCTGCGTtactactac                                                <  1:387958/45‑1 (MQ=255)
                                      aCCTGGCTGAATCATTTTGACGCAGACAGCCTGCGTtactactac                                                >  1:2989953/1‑45 (MQ=255)
                                      aCCTGGCTGAATCATTTTGACGCAGACAGCCTGCGTtactactac                                                >  1:3091139/1‑45 (MQ=255)
                                      aCCTGGCTGAATCATTTTGACGCAGACAGCCTGCGTtactactac                                                >  1:3018303/1‑45 (MQ=255)
                                      aCCTGGCTGAATCATTTTGACGCAGACAGCCTGCGTtactactac                                                <  1:1788455/45‑1 (MQ=255)
                                      aCCTGGCTGAATCATTTTGACGCAGACAGCCTGCGTtactactac                                                <  1:733294/45‑1 (MQ=255)
                                      aCCTGGCTGAATCATTTTGACGCAGACAGCCTGCGTtactactac                                                <  1:797293/45‑1 (MQ=255)
                                      aCCTGGCTGAATCATTTTGACGCAGACAGCCTGCGTtactactac                                                <  1:137118/45‑1 (MQ=255)
                                      aCCTGGCTGAATCATTTTGACGCAGACAGCCTGCGTtactactac                                                <  1:2376738/45‑1 (MQ=255)
                                      aCCTCGCTGAATCATTTTGACGCAGACAGCCTGCGTtactactac                                                <  1:2156542/45‑1 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGATACTCTCTTCGCGCATtgat                     <  1:2067024/69‑1 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCAt                         >  1:2180013/1‑65 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     >  1:2949547/1‑69 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     >  1:278289/1‑69 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     >  1:2920472/1‑69 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     >  1:985034/1‑69 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     <  1:3095406/69‑1 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     <  1:427771/69‑1 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     <  1:443707/69‑1 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     >  1:480946/1‑69 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     >  1:868116/1‑69 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     >  1:932088/1‑69 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     >  1:977979/1‑69 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     >  1:1054434/1‑69 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     <  1:1660486/69‑1 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     <  1:2043516/69‑1 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     >  1:208346/1‑69 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     >  1:1872226/1‑69 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     >  1:2168578/1‑69 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     >  1:1832235/1‑69 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     <  1:1831774/69‑1 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     >  1:2183868/1‑69 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     >  1:1549429/1‑69 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     <  1:1519261/69‑1 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     >  1:2568113/1‑69 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     >  1:2655664/1‑69 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATtgat                     >  1:2746314/1‑69 (MQ=255)
                                         tGGCTGAATCATTTTGACGCAGACAGCCTGCGTAACTACTACACTGCGaaa                                       >  1:1126227/1‑51 (MQ=255)
                                          ggCTGAATCATTTTGACGCCGACAGCCTGCGCTACTACTACACCGCCAAACTGTCTTCGCGCATtgatg                    <  1:2744765/69‑1 (MQ=255)
                                          ggCTGAATCATTTTGACGCCGACAGCCTGCGCTACTACTACACCGCCAAACTGTCTTCGCGCATtgatg                    <  1:727663/69‑1 (MQ=255)
                                          ggCTGAATCATTTTGACGCCGACAGCCTGCGCTACTACTACACCGCCAAACTGTCTTCGCGCATtgatg                    <  1:926494/69‑1 (MQ=255)
                                          ggCTGAATCATTTTGACGCCGACAGCCTGCGCTACTACTACACCGCCAAACTGTCTTCGCGCATtgatg                    <  1:2272733/69‑1 (MQ=255)
                                          ggCTGAATCATTTTGACGCCGACAGCCTGCGCTACTACTACACCGCCAAACTGTCTTCGCGCATtgatg                    <  1:2436936/69‑1 (MQ=255)
                                          ggCTGAATCATTTTGACGCCGACAGCCTGCGCTACTACTACACCGCCAAACTGTCTTCGCGCATtgatg                    <  1:2535844/69‑1 (MQ=255)
                                          ggCTGAATCATTTTGACGCCGACAGCCTGCGCTACTACTACACCGCCAAACTGTCTTCGCGCATtgatg                    <  1:2806270/69‑1 (MQ=255)
                                                aTCATTTTGACGCAGAAAGCCTGCGTtactactac                                                <  1:3214143/35‑1 (MQ=255)
                                                            cAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATTGATGATATCGATCTCAACCTgg  <  1:2119294/69‑1 (MQ=255)
                                                            cAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATTGATGATATCGATCTCAACCTgg  <  1:3210351/69‑1 (MQ=255)
                                                            cAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATTGATGATATCGATCTCAACCTgg  <  1:3143572/69‑1 (MQ=255)
                                                            cAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATTGATGATATCGATCTCAACCTgg  <  1:2313909/69‑1 (MQ=255)
                                                            cAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATTGATGATATCGATCTCAACCTgg  <  1:1447670/69‑1 (MQ=255)
                                                            cAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATTGATGATATCGATCTCAACCTgg  <  1:1256509/69‑1 (MQ=255)
                                                            cAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATTGATGATATCGATCTCAACCTgg  <  1:1227216/69‑1 (MQ=255)
                                                             |                                                                   
AGATGTCCAAGTCTCGCGGCACCTTTATTAAAGCCAGCACCTGGCTGAATCATTTTGACGCAGACAGCCTGCGTTACTACTACACTGCGAAACTCTCTTCGCGCATTGATGATATCGATCTCAACCTGG  >  minE/1363508‑1363636

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: