Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,662,691 C→T 39.7% L46L (TTG→TTA rtcR ← sigma 54‑dependent transcriptional regulator of rtcBA expression

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,662,6910CT39.7% 46.6 / 61.5 93L46L (TTG→TTArtcRsigma 54‑dependent transcriptional regulator of rtcBA expression
Reads supporting (aligned to +/- strand):  ref base C (19/37);  new base T (11/26);  total (30/63)
Fisher's exact test for biased strand distribution p-value = 8.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGGCGAAACGCTGGCGATATCGCGTTTAAGCGTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGT  >  minE/2662623‑2662739
                                                                    |                                                
tCTGGCGAAACGCTGGCGATATCGCGTTTAAGCGTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGcaac                                                  <  1:2550573/69‑1 (MQ=255)
                    atCGCGTTTAAGCGTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGa                               >  1:489741/1‑68 (MQ=255)
                           ttAAGCGTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGAt                       <  1:246244/69‑1 (MQ=255)
                           ttAAGCGTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGAt                       <  1:687130/69‑1 (MQ=255)
                           ttAAGCGTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGAt                       <  1:835668/69‑1 (MQ=255)
                            tAAGCGTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGAtt                      >  1:3151558/1‑69 (MQ=255)
                            tAAGCGTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGAt                       <  1:3076120/68‑1 (MQ=255)
                            tAAGCGTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGAt                       <  1:650300/68‑1 (MQ=255)
                                 gTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTcc                                            >  1:406770/1‑42 (MQ=255)
                                 gTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTcc                                            >  1:699461/1‑42 (MQ=255)
                                 gTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTcc                                            >  1:2776108/1‑42 (MQ=255)
                                 gTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTcc                                            >  1:1483896/1‑42 (MQ=255)
                                 gTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTcc                                            >  1:1636627/1‑42 (MQ=255)
                                 gTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTcc                                            >  1:2471795/1‑42 (MQ=255)
                                 gTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTcc                                            >  1:2310094/1‑42 (MQ=255)
                                 gTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTcc                                            >  1:1290021/1‑42 (MQ=255)
                                 gTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTcc                                            >  1:3082361/1‑42 (MQ=255)
                                  tttCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGAt                       <  1:3342976/62‑1 (MQ=255)
                                        aTAGCGAGGGCGAGCGGGCGTCGTGCAATAATTCCAGccgatcg                                   <  1:55602/44‑1 (MQ=37)
                                        aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGccgatcg                                   <  1:1909686/44‑1 (MQ=38)
                                        aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGccgatcg                                   <  1:1703929/44‑1 (MQ=38)
                                        aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGccgatcg                                   <  1:1667511/44‑1 (MQ=38)
                                        aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGccgatcg                                   <  1:3280043/44‑1 (MQ=38)
                                        aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGccgatcg                                   <  1:1615559/44‑1 (MQ=38)
                                        aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGccgatcg                                   <  1:1470944/44‑1 (MQ=38)
                                        aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGccgatcg                                   <  1:1454524/44‑1 (MQ=38)
                                        aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGccgatcg                                   <  1:2794569/44‑1 (MQ=38)
                                        aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa          <  1:346974/69‑1 (MQ=255)
                                        aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa          <  1:2322287/69‑1 (MQ=255)
                                        aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa          <  1:1163004/69‑1 (MQ=255)
                                        aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa          <  1:1367879/69‑1 (MQ=255)
                                        aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa          <  1:2867891/69‑1 (MQ=255)
                                        aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa          <  1:2665142/69‑1 (MQ=255)
                                        aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa          <  1:708057/69‑1 (MQ=255)
                                        aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa          <  1:2692919/69‑1 (MQ=255)
                                        aTAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa          <  1:1797527/69‑1 (MQ=255)
                                        aTAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa          <  1:1593613/69‑1 (MQ=255)
                                        aTAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa          <  1:42392/69‑1 (MQ=255)
                                         tAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGccgatcg                                   <  1:1161679/43‑1 (MQ=38)
                                         tAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa          <  1:2989013/68‑1 (MQ=255)
                                         tAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa          <  1:3046492/68‑1 (MQ=255)
                                         tAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa          <  1:1561075/68‑1 (MQ=255)
                                          agcgagcgCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa          <  1:3343003/67‑1 (MQ=255)
                                          agcgagcgCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa          <  1:1955225/67‑1 (MQ=255)
                                          agcgagcgCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa          >  1:6111/1‑67 (MQ=255)
                                           gcgagcgcgagcGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCATCGATTCTTGCTGGCATa          >  1:970175/1‑66 (MQ=255)
                                           gcgagcgcgagcGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGAt                       <  1:2569264/53‑1 (MQ=255)
                                           gcgagcgcgagcGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa          >  1:1188771/1‑66 (MQ=255)
                                           gcgagcgcgagcGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa          >  1:1926787/1‑66 (MQ=255)
                                           gcgagcgcgagcGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa          >  1:98355/1‑66 (MQ=255)
                                           gcgagcgcgagcGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa          >  1:1290517/1‑66 (MQ=255)
                                           gcgagcgcgagcGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa          >  1:1378662/1‑66 (MQ=255)
                                           gcgagcgcgagcGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa          >  1:774344/1‑66 (MQ=255)
                                             gagcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACag     <  1:212182/69‑1 (MQ=255)
                                             gagcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACag     <  1:111227/69‑1 (MQ=255)
                                             gagcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACag     <  1:116077/69‑1 (MQ=255)
                                             gagcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACag     <  1:1628432/69‑1 (MQ=255)
                                             gagcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACag     <  1:2674322/69‑1 (MQ=255)
                                             gagcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACag     <  1:2524031/69‑1 (MQ=255)
                                             gagcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACag     <  1:2471419/69‑1 (MQ=255)
                                             cagcgcgtgcgGGCGTCGTGCAACTATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACag     <  1:313677/68‑1 (MQ=255)
                                              agcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACag     <  1:3347099/68‑1 (MQ=255)
                                               gcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACag     <  1:1507318/67‑1 (MQ=255)
                                                cgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCtt                   <  1:1431412/52‑1 (MQ=255)
                                                cgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCtt                   <  1:164516/52‑1 (MQ=255)
                                                cgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt  <  1:825495/69‑1 (MQ=255)
                                                cgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt  <  1:348728/69‑1 (MQ=255)
                                                cgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt  <  1:984866/69‑1 (MQ=255)
                                                cgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt  <  1:3279480/69‑1 (MQ=255)
                                                cgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt  <  1:2082057/69‑1 (MQ=255)
                                                cgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt  <  1:3010766/69‑1 (MQ=255)
                                                cgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt  <  1:2106219/69‑1 (MQ=255)
                                                cgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt  <  1:128661/69‑1 (MQ=255)
                                                cgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt  <  1:2698397/69‑1 (MQ=255)
                                                cgcgagcgGGCGTCGCGCAACAATTCCAGTCGATCGATGACCAACGATTCCTGCTGGCATAAACAGAgt  <  1:910220/69‑1 (MQ=255)
                                                  cgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt  <  1:1405917/67‑1 (MQ=255)
                                                    agcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt  <  1:3360311/65‑1 (MQ=255)
                                                     gcgGGCGTCGTGCAATAATTCCGGCCGATCGATAAct                             >  1:3123315/1‑36 (MQ=19)
                                                     gcgGGCGTCGTGCAATAATTCCAGCCGATCGATAAct                             >  1:2768351/1‑36 (MQ=25)
                                                     gcgGGCGTCGTGCAATAATTCCAGCCGATCGATAAct                             >  1:2591250/1‑36 (MQ=25)
                                                     gcgGGCGTCGTGCAATAATTCCAGCCGATCGATAAct                             >  1:2398931/1‑36 (MQ=25)
                                                     gcgGGCGTCGTGCAATAATTCCAGCCGATCGATAAct                             >  1:2123060/1‑36 (MQ=25)
                                                     gcgGGCGTCGTGCAATAATTCCAGCCGATCGATAAct                             >  1:2794144/1‑36 (MQ=25)
                                                     gcgGGCGTCGTGCAATAATTCCAGCCGATCGATAAct                             >  1:684072/1‑36 (MQ=25)
                                                     gcgGGCGTCGTGCAATAATTCCAGCCGATCGATAAct                             >  1:1628723/1‑36 (MQ=25)
                                                     gcgGGCGTCGTGCAATAATTCCAGCCGATCGATAAct                             >  1:141048/1‑36 (MQ=25)
                                                     gcgGGCGTCGTGCAATAATTCCAGCCGATCGATAAct                             >  1:1364824/1‑36 (MQ=25)
                                                     gcgGGCGTCGTGCAATAATTCCAGCCGATCGATAAct                             >  1:1361706/1‑36 (MQ=25)
                                                     gcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGAt                       <  1:938953/43‑1 (MQ=255)
                                                     gcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACag     <  1:1150052/61‑1 (MQ=255)
                                                      cgGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa          <  1:3172890/55‑1 (MQ=255)
                                                       gggCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa          <  1:2411379/54‑1 (MQ=255)
                                                        ggCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa          <  1:1295052/53‑1 (MQ=255)
                                                                gCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa          <  1:1603394/45‑1 (MQ=255)
                                                                    |                                                
TCTGGCGAAACGCTGGCGATATCGCGTTTAAGCGTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGT  >  minE/2662623‑2662739

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: