Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,801,022 C→T 49.2% D160D (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,0220CT49.2% ‑0.1 / 64.5 74D160D (GAC→GATyjbNtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base C (16/21);  new base T (14/22);  total (31/43)
Fisher's exact test for biased strand distribution p-value = 8.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATG  >  minE/2800954‑2801088
                                                                    |                                                                  
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:3104722/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                                >  1:1731908/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                                >  1:2159797/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                                >  1:197300/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                                >  1:1903123/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                                >  1:1849667/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                                >  1:2343794/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                                >  1:2410543/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                                >  1:174787/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                                >  1:468749/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                                >  1:1590364/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                                >  1:1031997/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                                >  1:858890/1‑50 (MQ=255)
                                       acggtgaAAACGCGTATTGGCATCGATGATCAGGACAGCTATGAATTtct                                                >  1:662025/1‑50 (MQ=255)
                                               aaCGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCg                      >  1:733554/1‑68 (MQ=255)
                                               aaCGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCg                      >  1:1280408/1‑68 (MQ=255)
                                               aaCGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCg                      >  1:413749/1‑68 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt              <  1:1451975/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt              <  1:554086/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt              <  1:2922405/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt              <  1:799194/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt              <  1:2764009/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt              <  1:3084209/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt              <  1:1002546/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt              <  1:1140945/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt              <  1:3173854/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt              <  1:876120/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt              <  1:1256892/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt              <  1:3217546/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt              <  1:1286253/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt              <  1:1701023/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt              <  1:1708239/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt              <  1:1739827/69‑2 (MQ=255)
                                                      aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt              <  1:423762/69‑2 (MQ=255)
                                                      aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc              <  1:3097998/69‑1 (MQ=255)
                                                      aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc              <  1:922419/69‑1 (MQ=255)
                                                      aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc              <  1:2738232/69‑1 (MQ=255)
                                                      aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc              <  1:2104685/69‑1 (MQ=255)
                                                      aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc              <  1:748613/69‑1 (MQ=255)
                                                      aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc              <  1:1217470/69‑1 (MQ=255)
                                                      aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc              <  1:1363416/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:2360573/1‑67 (MQ=255)
                                                            aTCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgt          >  1:1294927/1‑67 (MQ=255)
                                                            aTCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgt          >  1:188080/1‑67 (MQ=255)
                                                            aTCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgt          >  1:1129143/1‑67 (MQ=255)
                                                            aTCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgt          >  1:2797894/1‑67 (MQ=255)
                                                            aTCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgt          >  1:2776741/1‑67 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga      <  1:3308995/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga      <  1:1372923/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga      <  1:133577/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga      <  1:168429/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga      <  1:3410258/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga      <  1:2278925/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga      <  1:2009620/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga      <  1:87341/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga      <  1:824595/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga      <  1:2060544/69‑1 (MQ=255)
                                                              cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga      <  1:2187629/69‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg             <  1:1091847/61‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg             <  1:1150491/61‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg             <  1:894739/61‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg             <  1:2359893/61‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg             <  1:842812/61‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg             <  1:49336/61‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg             <  1:2536293/61‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg             <  1:1760393/61‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg             <  1:276932/61‑1 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg             >  1:2238824/1‑61 (MQ=255)
                                                               gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg             <  1:2868832/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)
                                                                  gacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATg  <  1:505303/69‑1 (MQ=255)
                                                                  gacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATg  <  1:2016211/69‑1 (MQ=255)
                                                                  gacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATg  <  1:88368/69‑1 (MQ=255)
                                                                  gacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATg  <  1:2285798/69‑1 (MQ=255)
                                                                  gacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATg  <  1:999701/69‑1 (MQ=255)
                                                                    |                                                                  
TGCGTGAAAGCGATGCGCGATGTGGTGTCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATG  >  minE/2800954‑2801088

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: