Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,279,749 Δ1 bp coding (633/1101 nt) trmA → tRNA (uracil‑5‑)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,279,7460A.98.5% 271.7 / ‑1.6 65coding (630/1101 nt)trmAtRNA (uracil‑5‑)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (50/14);  minor base C (1/0);  total (51/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTAAT  >  minE/2279679‑2279809
                                                                   |                                                               
ttACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaag                                                              >  1:2563630/1‑70 (MQ=255)
ttACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaag                                                              >  1:4193132/1‑70 (MQ=255)
 tACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaag                                                              <  1:527827/70‑2 (MQ=255)
 tACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaag                                                              <  1:2898387/70‑2 (MQ=255)
 tACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaag                                                              <  1:3930627/70‑2 (MQ=255)
 tACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaag                                                              <  1:3000736/70‑2 (MQ=255)
 tACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaag                                                              <  1:3607347/70‑2 (MQ=255)
  aCCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaag                                                              <  1:1705836/69‑2 (MQ=255)
                ggcgATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCt                                             >  1:324468/1‑71 (MQ=255)
                           ttCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCa                                  <  1:1075938/71‑1 (MQ=255)
                           ttCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCa                                  <  1:2237136/71‑1 (MQ=255)
                           ttCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCa                                  <  1:4278080/71‑1 (MQ=255)
                           ttCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCa                                  <  1:1666248/71‑1 (MQ=255)
                                        aTGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAAc                                >  1:602801/1‑60 (MQ=255)
                                        aTGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACtttt                      >  1:2403187/1‑70 (MQ=255)
                                        aTGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACtttt                      >  1:316493/1‑70 (MQ=255)
                                        aTGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACtttt                      >  1:3543723/1‑70 (MQ=255)
                                        aTGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACtttt                      >  1:3715283/1‑70 (MQ=255)
                                        aTGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACtttt                      >  1:391726/1‑70 (MQ=255)
                                        aTGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACtttt                      >  1:4261446/1‑70 (MQ=255)
                                        aTGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACtttt                      >  1:428613/1‑70 (MQ=255)
                                        aTGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACtttt                      >  1:740334/1‑70 (MQ=255)
                                        aTGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACtttt                      >  1:80027/1‑70 (MQ=255)
                                        aTGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACtttt                      >  1:1204639/1‑70 (MQ=255)
                                        aTGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACtttt                      >  1:2055585/1‑70 (MQ=255)
                                        aTGGGCGCTGGACGTAACCAAAGGCACAAAGCCGATTTACTGGAGCTGTACTGCGGCAACGCTAACtttt                      >  1:597341/1‑70 (MQ=255)
                                        aTGGGCGCTGGACGTAACCAAAGGCCCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACtttt                      >  1:85388/1‑70 (MQ=255)
                                        aTGGGCGCTGGACGTAACAAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACtttt                      >  1:2839047/1‑70 (MQ=255)
                                        aTGGGCGCCGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACtttt                      >  1:2059903/1‑70 (MQ=255)
                                            gcgcTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCAt                   >  1:1227092/1‑69 (MQ=255)
                                            gcgcTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCAt                   >  1:3090464/1‑69 (MQ=255)
                                            gcgcTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCAt                   >  1:2974113/1‑69 (MQ=255)
                                            gcgcTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCAt                   >  1:1311569/1‑69 (MQ=255)
                                            gcgcTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCAt                   >  1:124988/1‑69 (MQ=255)
                                             cgcTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTa                 >  1:4138578/1‑70 (MQ=255)
                                             cgcTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAg                >  1:2775622/1‑71 (MQ=255)
                                                   aCGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTgg          <  1:3238782/71‑1 (MQ=255)
                                                   aCGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTgg          <  1:3602396/71‑1 (MQ=255)
                                                   aCGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTgg          <  1:3308723/71‑1 (MQ=255)
                                                   aCGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTgg          <  1:1305167/71‑1 (MQ=255)
                                                   aCGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTgg          <  1:3249769/71‑1 (MQ=255)
                                                   aCGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTgg          <  1:575947/71‑1 (MQ=255)
                                                   aCGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTgg          <  1:1147524/71‑1 (MQ=255)
                                                   aCGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTgg          <  1:2255784/71‑1 (MQ=255)
                                                         ccAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGgcgc       >  1:662045/1‑68 (MQ=255)
                                                         ccAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGgcgc       >  1:3730732/1‑68 (MQ=255)
                                                         ccAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGcgcg      >  1:1221354/1‑69 (MQ=255)
                                                         ccAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGcgcg      >  1:1121769/1‑69 (MQ=255)
                                                         ccAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGt     >  1:2470584/1‑70 (MQ=255)
                                                         ccAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGt     >  1:1184551/1‑70 (MQ=255)
                                                         ccAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGt     >  1:1212956/1‑70 (MQ=255)
                                                         ccAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGt     >  1:709093/1‑70 (MQ=255)
                                                         ccAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGt     >  1:3921071/1‑70 (MQ=255)
                                                         ccAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGt     >  1:885347/1‑70 (MQ=255)
                                                         ccAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGt     >  1:154837/1‑70 (MQ=255)
                                                         ccAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGt     >  1:2505992/1‑70 (MQ=255)
                                                         ccAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGt     >  1:3378791/1‑70 (MQ=255)
                                                         ccAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGt     >  1:3288708/1‑70 (MQ=255)
                                                         ccAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGt     >  1:2691055/1‑70 (MQ=255)
                                                         ccAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGt     >  1:3086543/1‑70 (MQ=255)
                                                         ccAAAGACTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGt     >  1:3721/1‑70 (MQ=255)
                                                          cAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTaa   >  1:2364198/1‑71 (MQ=255)
                                                          cAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTaa   >  1:3894711/1‑71 (MQ=255)
                                                          cAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTaa   >  1:86637/1‑71 (MQ=255)
                                                          cAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTaa   >  1:2635917/1‑71 (MQ=255)
                                                          cAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTaa   >  1:1159697/1‑71 (MQ=255)
                                                          cAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTaa   >  1:1061676/1‑71 (MQ=255)
                                                          cAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTaa   >  1:3220731/1‑71 (MQ=255)
                                                          cAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTaa   >  1:3503936/1‑71 (MQ=255)
                                                          cAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTaa   >  1:4243374/1‑71 (MQ=255)
                                                          cAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTaa   >  1:3610765/1‑71 (MQ=255)
                                                           aaaGGCTCAAAGGCGATTTACTGGAGCTGGACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTAAt  <  1:2232793/71‑1 (MQ=255)
                                                           aaaGGCTCAAAGGCGATTTACTGGAGCTGCACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTAAt  <  1:2227862/71‑1 (MQ=255)
                                                                   |                                                               
TTACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTAAT  >  minE/2279679‑2279809

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: