Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,801,058 T→C 30.9% N172N (AAT→AAC yjbN → tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,801,0580TC30.9% 113.1 / 63.4 126N172N (AAT→AACyjbNtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base T (30/57);  new base C (17/22);  total (47/79)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

Only 100 of 128 total aligned reads displayed.

AACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGG  >  minE/2801001‑2801124
                                                         |                                                                  
aaCGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCg                                                          >  1:733554/1‑68 (MQ=255)
       aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt                                                  <  1:3217546/69‑2 (MQ=255)
       aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt                                                  <  1:1739827/69‑2 (MQ=255)
       aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt                                                  <  1:423762/69‑2 (MQ=255)
       aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt                                                  <  1:1286253/69‑2 (MQ=255)
       aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt                                                  <  1:1256892/69‑2 (MQ=255)
       aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt                                                  <  1:876120/69‑2 (MQ=255)
       aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt                                                  <  1:1708239/69‑2 (MQ=255)
       aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt                                                  <  1:1701023/69‑2 (MQ=255)
       aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt                                                  <  1:1140945/69‑2 (MQ=255)
       aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt                                                  <  1:3173854/69‑2 (MQ=255)
       aTTGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt                                                  <  1:3084209/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:1002546/69‑2 (MQ=255)
       aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc                                                  <  1:1217470/69‑1 (MQ=255)
       aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc                                                  <  1:748613/69‑1 (MQ=255)
       aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc                                                  <  1:1363416/69‑1 (MQ=255)
       aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc                                                  <  1:2104685/69‑1 (MQ=255)
       aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc                                                  <  1:3097998/69‑1 (MQ=255)
       aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc                                                  <  1:2738232/69‑1 (MQ=255)
       aTTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGc                                                  <  1:922419/69‑1 (MQ=255)
        ttGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCaaatgt                                                  <  1:156212/68‑4 (MQ=255)
        ttGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt                                                  <  1:1797131/68‑2 (MQ=255)
        ttGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt                                                  <  1:1137996/68‑2 (MQ=255)
        ttGGCATCGATGATCAGGACAGCTATGAATTTCTCTGCGATTTCATCAACACCGTTTCCGGCAAAggt                                                  <  1:1139508/68‑2 (MQ=255)
             aTCGACGAGCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgt                                              >  1:3408893/1‑67 (MQ=255)
             aTCGACGACCAGGAGAGCTATGAATTTCTCTGCGATTTCATCAATACAGTTTCCGGCAAAGGCGAgt                                              >  1:695437/1‑67 (MQ=255)
             aTCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgt                                              >  1:2797894/1‑67 (MQ=255)
             aTCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgt                                              >  1:2776741/1‑67 (MQ=255)
             aTCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgt                                              >  1:2360573/1‑67 (MQ=255)
             aTCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgt                                              >  1:188080/1‑67 (MQ=255)
             aTCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAgt                                              >  1:1294927/1‑67 (MQ=255)
               cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga                                          <  1:168429/69‑1 (MQ=255)
               cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga                                          <  1:2278925/69‑1 (MQ=255)
               cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga                                          <  1:3308995/69‑1 (MQ=255)
               cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga                                          <  1:2187629/69‑1 (MQ=255)
               cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga                                          <  1:2009620/69‑1 (MQ=255)
               cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga                                          <  1:3410258/69‑1 (MQ=255)
               cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga                                          <  1:824595/69‑1 (MQ=255)
               cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga                                          <  1:1372923/69‑1 (MQ=255)
               cgacgaCCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga                                          <  1:87341/69‑1 (MQ=255)
                gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg                                                 <  1:49336/61‑1 (MQ=255)
                gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg                                                 <  1:842812/61‑1 (MQ=255)
                gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg                                                 <  1:894739/61‑1 (MQ=255)
                gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg                                                 <  1:2536293/61‑1 (MQ=255)
                gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg                                                 <  1:1091847/61‑1 (MQ=255)
                gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCg                                                 <  1:2359893/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:1738683/1‑68 (MQ=255)
                gacgacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTga                                          <  1:2904635/68‑1 (MQ=255)
                   gacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATg                                      <  1:2285798/69‑1 (MQ=255)
                   gacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATg                                      <  1:505303/69‑1 (MQ=255)
                   gacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATg                                      <  1:2016211/69‑1 (MQ=255)
                         gACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTtcatc                                >  1:1112230/1‑69 (MQ=255)
                         gACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTtcatc                                >  1:1072018/1‑69 (MQ=255)
                         gACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTtcatc                                >  1:1068078/1‑69 (MQ=255)
                            aGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCatcatc                             <  1:102789/69‑1 (MQ=255)
                            aGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCatcatc                             <  1:3104083/69‑1 (MQ=255)
                            aGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCatcatc                             <  1:3333704/69‑1 (MQ=255)
                            aGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCatcatc                             <  1:289543/69‑1 (MQ=255)
                            aGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCatcatc                             <  1:146846/69‑1 (MQ=255)
                            aGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCatcatc                             <  1:1414521/69‑1 (MQ=255)
                              cTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCca                           <  1:2281281/69‑1 (MQ=255)
                                aTGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCcacg                         <  1:26016/69‑1 (MQ=255)
                                aTGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCcacg                         <  1:401480/69‑1 (MQ=255)
                                aTGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCcacg                         <  1:855322/69‑1 (MQ=255)
                                 tGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCcacgc                        >  1:360310/1‑69 (MQ=255)
                                 tGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCcacgc                        >  1:2227676/1‑69 (MQ=255)
                                     tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCcacgcacg                     >  1:1986737/1‑68 (MQ=255)
                                     tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCcacgcacg                     >  1:1987399/1‑68 (MQ=255)
                                     tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCcacgcacg                     >  1:1470069/1‑68 (MQ=255)
                                     tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCcacgcacg                     >  1:787228/1‑68 (MQ=255)
                                     tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCcacgcacg                     >  1:1032384/1‑68 (MQ=255)
                                     tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCcacgcacg                     >  1:2704023/1‑68 (MQ=255)
                                     tttCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGAGAGTGTGAGATGCTCATCATCcacgcacg                     >  1:1034067/1‑68 (MQ=255)
                                      ttCTCTGCGATTTCATCAACACCGTTTCCGGCAAAGGTGAGTGTga                                          >  1:1572060/1‑46 (MQ=255)
                                      ttCTCTGCGATTTCATCAACACCGTTTCCGGCAAAGGTGAGTGTGAGATg                                      >  1:583046/1‑50 (MQ=255)
                                      ttCTCTGCGATTTCATCAACACCGTTTCCGGCAAAGGTGAGTGTGAGATGt                                     >  1:1377911/1‑51 (MQ=255)
                                      ttCTCTGCGATTTCATCAACACCGTTTCCGGCAAAGGTGAGTGTGAGATGt                                     >  1:1324350/1‑51 (MQ=255)
                                      ttCTCTGCGATTTCATCAACACCGTTTCCGGCAAAGGTGAGTGTGAGATGt                                     >  1:2813025/1‑51 (MQ=255)
                                      ttCTCTGCGATTTCATCAACACCGTTTCCGGCAAAGGTGAGTGTGAGATGt                                     >  1:1851472/1‑51 (MQ=255)
                                      ttCTCTGCGATTTCATCAACACCGTTTCCGGCAAAGGTGAGTGTGAGATGt                                     >  1:945048/1‑51 (MQ=255)
                                      ttCTCTGCGATTTCATCAACACCGTTTCCGGCAAAGGTGAGTGTGAGATGt                                     >  1:1852350/1‑51 (MQ=255)
                                      ttCTCTGCGATTTCATCAACACCGTTTCCGGCAAAGGTGAGTGTGAGATGt                                     >  1:1985114/1‑51 (MQ=255)
                                      ttCTCTGCGATTTCATCAACACCGTTTCCGGCAAAGGTGAGTGTGAGATGt                                     >  1:3344697/1‑51 (MQ=255)
                                      ttCTCTGCGATTTCATCAACACCGTTTCCGGCAAAGGTGAGTGTGAGATGt                                     >  1:2127627/1‑51 (MQ=255)
                                      ttCTCTGCGATTTCATCAACACCGTTTCCGGCAAAGGTGAGTGTGAGATGt                                     >  1:3283197/1‑51 (MQ=255)
                                      ttCTCTGCGATTTCATCAACACCGTTTCCGGCAAAGGTGAGTGTGAGATGt                                     >  1:2385992/1‑51 (MQ=255)
                                      ttCTCTGCGATTTCATCAACACCGTTTCCGGCAAAGGTGAGTGTGAGATGt                                     >  1:3153843/1‑51 (MQ=255)
                                      ttCTCTGCGATTTCATCAACACCGTTTCCGGCAAAGGTGAGTGTGAGATGt                                     >  1:2790971/1‑51 (MQ=255)
                                               aTTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCt             <  1:1928970/66‑1 (MQ=255)
                                               aTTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCt             <  1:3311888/66‑1 (MQ=255)
                                               aTTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCt             <  1:249210/66‑1 (MQ=255)
                                                       aaTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCgggg  <  1:1251096/69‑1 (MQ=255)
                                                       aaTACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCgggg  <  1:2456766/69‑1 (MQ=255)
                                                         tACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCca                           <  1:2432461/42‑1 (MQ=255)
                                                         |                                                                  
AACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGG  >  minE/2801001‑2801124

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: