Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,533,542 G→A 35.4% P51S (CCT→TCT)  prmB ← N5‑glutamine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,533,5420GA35.4% 50.8 / 43.3 48P51S (CCT→TCT) prmBN5‑glutamine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (21/10);  new base A (11/6);  total (32/16)
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

TTTTTCGCTGGAGGTCAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCG  >  minE/1533474‑1533607
                                                                    |                                                                 
tttttCGCTGGAGGTCAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGc                                                                  <  1:1144480/70‑1 (MQ=255)
tttttCGCTGGAGGTCAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGc                                                                  <  1:4103914/70‑1 (MQ=255)
tttttCGCTGGAGGTCAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGc                                                                  >  1:3302717/1‑70 (MQ=255)
  tttCGCTGGAGGTCAGACGCGCGGTGCGCATTTCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGACaa                                                                <  1:713360/70‑1 (MQ=255)
  tttCGCTGGAGGTCAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGACaa                                                                <  1:3491441/70‑1 (MQ=255)
               cAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACgg                                                  >  1:1716070/1‑71 (MQ=255)
               cAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACgg                                                  >  1:625869/1‑71 (MQ=255)
               cAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGACAACACCAGCTg                                                       >  1:146520/1‑66 (MQ=255)
               cAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGACAACACCAGCTg                                                       >  1:3136695/1‑66 (MQ=255)
               cAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGACAACACCAGCTGTACgg                                                  >  1:406251/1‑71 (MQ=255)
               cAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGACAACACCAGCTGTACg                                                   >  1:2391633/1‑70 (MQ=255)
               cAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGACAACACCAGCTGTACg                                                   >  1:2515951/1‑70 (MQ=255)
                     cgcgGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCAt                                            <  1:1589016/71‑1 (MQ=255)
                              cATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAaca                                                              >  1:1730242/1‑44 (MQ=255)
                              cATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGt                                   >  1:2374656/1‑71 (MQ=255)
                              cATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGt                                   >  1:1511101/1‑71 (MQ=255)
                              cATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGt                                   >  1:618988/1‑71 (MQ=255)
                                 atCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAaca                                                              >  1:1291472/1‑41 (MQ=255)
                                 atCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAaca                                                              >  1:2059314/1‑41 (MQ=255)
                                 atCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAaca                                                              >  1:347182/1‑41 (MQ=255)
                                 atCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAaca                                                              >  1:741374/1‑41 (MQ=255)
                                    ttCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGTTATCgg                             >  1:526607/1‑71 (MQ=255)
                                    ttCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGTTATCgg                             >  1:1957459/1‑71 (MQ=255)
                                     tCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGtt                                  <  1:3405011/65‑1 (MQ=255)
                                      ccGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGTTATc                               <  1:2120385/67‑1 (MQ=255)
                                      ccGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGTTATc                               <  1:2584765/67‑1 (MQ=255)
                                      ccGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGTTATc                               <  1:2738131/67‑1 (MQ=255)
                                      ccGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGTTATc                               <  1:1837882/67‑1 (MQ=255)
                                      ccGGAATATCCAGCGGCAGGTAGAGCGAAGACAACACCAGCTGTACGGCTTCATCCCACGGGTTATCgg                             <  1:737961/69‑1 (MQ=255)
                                      ccGGAATATCCAGCGGCAGGTAGAGCGAAGACAACACCAGCTGTACGGCTTCATCCCACGGGTTATCgg                             <  1:1245170/69‑1 (MQ=255)
                                      ccGGAATATCCAGCGGCAGGTAGAGCGAAGACAACACCAGCTGTACGGCTTCATCCCACGGGTTATCgg                             <  1:749709/69‑1 (MQ=255)
                                              tCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGTTATCGGTACCGTGAcc                   >  1:71184/1‑71 (MQ=255)
                                              tCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGTTATCGGTACCGTGAcc                   >  1:1619163/1‑71 (MQ=255)
                                              tCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGTTATCGGTACCGTGAcc                   >  1:1592452/1‑71 (MQ=255)
                                                cAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCAt                                            >  1:2291991/1‑44 (MQ=255)
                                                cAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGTTATCGGTACCGTGACCg                  >  1:221454/1‑70 (MQ=255)
                                                cAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGTTATCGGTACCGTGACCg                  <  1:1192604/70‑1 (MQ=255)
                                                    ggCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGTTATCGGTACCGTGACCGTAc               >  1:1962272/1‑69 (MQ=255)
                                                         gaagagCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGTTATCGGTACCGTGACCGTAc               <  1:780222/62‑1 (MQ=255)
                                                          tAGAGCGAAGACAACACCAGCTGTACGGCTTCATCCCACGGGTTATCGGTACCGTGACCGTACCAGATAtt       <  1:2655619/71‑1 (MQ=255)
                                                           agagCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGTTATCGGTACCGTGACCGTACCAGATAttt      >  1:51418/1‑71 (MQ=255)
                                                           agagCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGTTATCGGTACCGTGACCGTACCAGATAttt      >  1:1147796/1‑71 (MQ=255)
                                                           agagCGAAGACAACACCAGCTGTACGGCTTCATcc                                          >  1:2571183/1‑35 (MQ=255)
                                                           agagCGAAGACAACACCAGCTGTACGGCTTCATCCCACGGGTTATCGGTa                           >  1:1476817/1‑50 (MQ=255)
                                                           agagCGAAGACAACACCAGCTGTACGGCTTCATCCCACGGGTTATCGGTACCGTGACCGTACCAGATAttt      >  1:1371595/1‑71 (MQ=255)
                                                           agagCGAAGACAACACCAGCTGTACGGCTTCATCCCACGGGTTATCGGTACCGTGACCGTACCAGATAtt       >  1:2246957/1‑70 (MQ=255)
                                                               cGAAGACAACACCAGCTGTACGGCTTCATCCCACGGGTTATCGGTACCGTGACCGTACCAGATATTTGCcg  >  1:337362/1‑71 (MQ=255)
                                                               cGAAGACAACACCAGCTGTACGGCTTCATCCCACGGGTTATCGGTACCGTGACCGTACCAGATATTTGCcg  >  1:1220147/1‑71 (MQ=255)
                                                                    |                                                                 
TTTTTCGCTGGAGGTCAGACGCGCGGTGCGCATATCTTCCGGAATATCCAGCGGCAGGTAGAGCGAAGGCAACACCAGCTGTACGGCTTCATCCCACGGGTTATCGGTACCGTGACCGTACCAGATATTTGCCG  >  minE/1533474‑1533607

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