Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,662,700 T→C 32.5% R43R (CGA→CGG 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,7000TC32.5% 66.4 / 44.0 83R43R (CGA→CGGrtcRsigma 54‑dependent transcriptional regulator of rtcBA expression
Reads supporting (aligned to +/- strand):  ref base T (10/46);  new base C (11/16);  total (21/62)
Fisher's exact test for biased strand distribution p-value = 3.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

ATCGCGTTTAAGCGTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTG  >  minE/2662643‑2662765
                                                         |                                                                 
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:650300/68‑1 (MQ=255)
        tAAGCGTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGAt                                                 <  1:3076120/68‑1 (MQ=255)
              tttCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGAt                                                 <  1:3342976/62‑1 (MQ=255)
                    aTAGCGAGGGCGAGCGGGCGTCGTGCAATAATTCCAGccgatcg                                                             <  1:55602/44‑1 (MQ=37)
                    aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGccgatcg                                                             <  1:1667511/44‑1 (MQ=38)
                    aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGccgatcg                                                             <  1:2794569/44‑1 (MQ=38)
                    aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGccgatcg                                                             <  1:1909686/44‑1 (MQ=38)
                    aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGccgatcg                                                             <  1:1703929/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:3280043/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:708057/69‑1 (MQ=255)
                    aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa                                    <  1:2867891/69‑1 (MQ=255)
                    aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa                                    <  1:2665142/69‑1 (MQ=255)
                    aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa                                    <  1:2692919/69‑1 (MQ=255)
                    aTAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa                                    <  1:1367879/69‑1 (MQ=255)
                    aTAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa                                    <  1:42392/69‑1 (MQ=255)
                    aTAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa                                    <  1:1797527/69‑1 (MQ=255)
                    aTAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa                                    <  1:1593613/69‑1 (MQ=255)
                     tAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGccgatcg                                                             <  1:1161679/43‑1 (MQ=38)
                     tAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa                                    <  1:3046492/68‑1 (MQ=255)
                     tAGCGAGCGCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa                                    <  1:2989013/68‑1 (MQ=255)
                     tAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa                                    <  1:1561075/68‑1 (MQ=255)
                      agcgagcgCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa                                    <  1:1955225/67‑1 (MQ=255)
                      agcgagcgCGAGCGGGCGTCGTGCAATAATTCCAGCCGATCGATAACTAACGATTCTTGCTGGCATa                                    <  1:3343003/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:1378662/1‑66 (MQ=255)
                       gcgagcgcgagcGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa                                    >  1:1188771/1‑66 (MQ=255)
                       gcgagcgcgagcGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa                                    >  1:1290517/1‑66 (MQ=255)
                       gcgagcgcgagcGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa                                    >  1:1926787/1‑66 (MQ=255)
                       gcgagcgcgagcGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa                                    >  1:98355/1‑66 (MQ=255)
                       gcgagcgcgagcGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATa                                    >  1:774344/1‑66 (MQ=255)
                         gagcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACag                               <  1:212182/69‑1 (MQ=255)
                         gagcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACag                               <  1:2471419/69‑1 (MQ=255)
                         gagcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACag                               <  1:2524031/69‑1 (MQ=255)
                         gagcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACag                               <  1:116077/69‑1 (MQ=255)
                         gagcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACag                               <  1:111227/69‑1 (MQ=255)
                         gagcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACag                               <  1:2674322/69‑1 (MQ=255)
                         gagcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACag                               <  1:1628432/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:2698397/69‑1 (MQ=255)
                            cgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt                            <  1:984866/69‑1 (MQ=255)
                            cgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt                            <  1:128661/69‑1 (MQ=255)
                            cgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt                            <  1:3279480/69‑1 (MQ=255)
                            cgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt                            <  1:825495/69‑1 (MQ=255)
                            cgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt                            <  1:3010766/69‑1 (MQ=255)
                            cgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt                            <  1:2082057/69‑1 (MQ=255)
                            cgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt                            <  1:2106219/69‑1 (MQ=255)
                            cgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt                            <  1:348728/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:2591250/1‑36 (MQ=25)
                                 gcgGGCGTCGTGCAATAATTCCAGCCGATCGATAAct                                                       >  1:684072/1‑36 (MQ=25)
                                 gcgGGCGTCGTGCAATAATTCCAGCCGATCGATAAct                                                       >  1:1361706/1‑36 (MQ=25)
                                 gcgGGCGTCGTGCAATAATTCCAGCCGATCGATAAct                                                       >  1:2794144/1‑36 (MQ=25)
                                 gcgGGCGTCGTGCAATAATTCCAGCCGATCGATAAct                                                       >  1:2768351/1‑36 (MQ=25)
                                 gcgGGCGTCGTGCAATAATTCCAGCCGATCGATAAct                                                       >  1:1364824/1‑36 (MQ=25)
                                 gcgGGCGTCGTGCAATAATTCCAGCCGATCGATAAct                                                       >  1:2398931/1‑36 (MQ=25)
                                 gcgGGCGTCGTGCAATAATTCCAGCCGATCGATAAct                                                       >  1:141048/1‑36 (MQ=25)
                                 gcgGGCGTCGTGCAATAATTCCAGCCGATCGATAAct                                                       >  1:2123060/1‑36 (MQ=25)
                                 gcgGGCGTCGTGCAATAATTCCAGCCGATCGATAAct                                                       >  1:1628723/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)
                                                     ccAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCt   <  1:2127079/69‑1 (MQ=255)
                                                     ccAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCt   <  1:829567/69‑1 (MQ=255)
                                                     ccAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCt   <  1:689392/69‑1 (MQ=255)
                                                     ccAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCt   <  1:2128582/69‑1 (MQ=255)
                                                     ccAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCt   <  1:2565415/69‑1 (MQ=255)
                                                     ccAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCt   <  1:2838133/69‑1 (MQ=255)
                                                     ccAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCt   <  1:314475/69‑1 (MQ=255)
                                                       aGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTg  <  1:376439/68‑1 (MQ=255)
                                                        gTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAgt                            <  1:534375/41‑1 (MQ=255)
                                                         |                                                                 
ATCGCGTTTAAGCGTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTG  >  minE/2662643‑2662765

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: