Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,801,019 C→T 57.2% D159D (GAC→GAT yjbN → tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,801,0190CT57.2% 13.9 / 47.7 63D159D (GAC→GATyjbNtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base C (17/10);  new base T (14/22);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 7.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAG  >  minE/2800954‑2801085
                                                                 |                                                                  
tGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCAGCgacgac                                                                 <  1:717783/69‑1 (MQ=255)
               cgcgATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGaa                                                  >  1:28519/1‑69 (MQ=255)
               cgcgATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGaa                                                  >  1:2998319/1‑69 (MQ=255)
                        gtgtCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGc                                         >  1:2396639/1‑69 (MQ=255)
                        gtgtCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGc                                         >  1:1041136/1‑69 (MQ=255)
                        gtgtCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGc                                         >  1:1066307/1‑69 (MQ=255)
                        gtgtCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGc                                         >  1:3353789/1‑69 (MQ=255)
                                 ccGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTtcatc                                <  1:1981918/69‑1 (MQ=255)
                                 ccGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTtcatc                                <  1:697887/69‑1 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                             >  1:1590364/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                             >  1:2343794/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                             >  1:1731908/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                             >  1:174787/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                             >  1:3104722/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                             >  1:1849667/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                             >  1:1903123/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                             >  1:197300/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                             >  1:2159797/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                             >  1:2410543/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                             >  1:468749/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                             >  1:858890/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                             >  1:662025/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                             >  1:1031997/1‑50 (MQ=255)
                                               aaCGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCg                   >  1:413749/1‑68 (MQ=255)
                                               aaCGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCg                   >  1:1280408/1‑68 (MQ=255)
                                               aaCGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCg                   >  1:733554/1‑68 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:1002546/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:876120/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:2764009/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:2922405/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:799194/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:3084209/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:1140945/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:423762/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:1256892/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:3173854/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:1286253/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:1739827/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:3217546/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:1708239/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:1701023/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:554086/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:1451975/69‑2 (MQ=255)
                                                      aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc           <  1:922419/69‑1 (MQ=255)
                                                      aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc           <  1:2738232/69‑1 (MQ=255)
                                                      aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc           <  1:748613/69‑1 (MQ=255)
                                                      aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc           <  1:3097998/69‑1 (MQ=255)
                                                      aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc           <  1:1363416/69‑1 (MQ=255)
                                                      aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc           <  1:1217470/69‑1 (MQ=255)
                                                      aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc           <  1:2104685/69‑1 (MQ=255)
                                                       ttGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCaaatgt           <  1:156212/68‑4 (MQ=255)
                                                       ttGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:2306272/68‑2 (MQ=255)
                                                       ttGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:1797131/68‑2 (MQ=255)
                                                       ttGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:1139508/68‑2 (MQ=255)
                                                       ttGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt           <  1:1137996/68‑2 (MQ=255)
                                                       ttGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc           <  1:1065097/68‑1 (MQ=255)
                                                            aTCGACGAGCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgt       >  1:3408893/1‑67 (MQ=255)
                                                            aTCGACGACCAGGAGAGCTATGAATTTCTCTGCGATTTCATCAATACAGTTTCCGGCAAAGGCGAgt       >  1:695437/1‑67 (MQ=255)
                                                            aTCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgt       >  1:2776741/1‑67 (MQ=255)
                                                            aTCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgt       >  1:2797894/1‑67 (MQ=255)
                                                            aTCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgt       >  1:1129143/1‑67 (MQ=255)
                                                            aTCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgt       >  1:1294927/1‑67 (MQ=255)
                                                            aTCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgt       >  1:188080/1‑67 (MQ=255)
                                                            aTCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgt       >  1:2360573/1‑67 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga   <  1:824595/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga   <  1:133577/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga   <  1:87341/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga   <  1:1372923/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga   <  1:3410258/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga   <  1:168429/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga   <  1:3308995/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga   <  1:2278925/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga   <  1:2009620/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga   <  1:2187629/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga   <  1:2060544/69‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg          <  1:2359893/61‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg          <  1:894739/61‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg          <  1:1091847/61‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg          >  1:2238824/1‑61 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg          <  1:842812/61‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg          <  1:2536293/61‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg          <  1:276932/61‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg          <  1:1150491/61‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg          <  1:2868832/61‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg          <  1:49336/61‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg          <  1:1760393/61‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTgagat  >  1:1172555/1‑68 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTgag  >  1:2046462/1‑69 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga   <  1:2904635/68‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga   >  1:1738683/1‑68 (MQ=255)
                                                               gacgacCAGCACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg          <  1:3370928/61‑1 (MQ=255)
                                                                acgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga   <  1:2001828/67‑1 (MQ=255)
                                                                 |                                                                  
TGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAG  >  minE/2800954‑2801085

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