Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,950,925 G→A 56.9% H221Y (CAT→TAT)  visC ← predicted oxidoreductase, FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,950,9250GA56.9% 14.9 / 47.9 65H221Y (CAT→TAT) visCpredicted oxidoreductase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base G (18/10);  new base A (29/8);  total (47/18)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

TGCTGCATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCC  >  minE/1950870‑1950983
                                                       |                                                          
tgctgcATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGc                                               >  1:660793/1‑69 (MQ=255)
tgctgcATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGc                                               >  1:1931387/1‑69 (MQ=255)
tgctgcATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGc                                               >  1:1297122/1‑69 (MQ=255)
tgctgcATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGc                                               >  1:2592673/1‑69 (MQ=255)
tgctgcATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGc                                               >  1:1778805/1‑69 (MQ=255)
tgctgcATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGc                                               >  1:256203/1‑69 (MQ=255)
   tgcATCCGCAGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCgg                                             >  1:2129979/1‑68 (MQ=255)
   tgcATCCGCAGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCgg                                             >  1:3327377/1‑68 (MQ=255)
   tgcATCCGCAGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCgg                                             >  1:3162384/1‑68 (MQ=255)
   tgcATCCGCAGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCgg                                             >  1:2385263/1‑68 (MQ=255)
         cGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTaaaaa                                       >  1:2592633/1‑68 (MQ=255)
         cGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAg                                      >  1:2792062/1‑69 (MQ=255)
         cGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAg                                      <  1:2158018/69‑1 (MQ=255)
         cGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAg                                      >  1:2700004/1‑69 (MQ=255)
         cGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAg                                      >  1:2910252/1‑69 (MQ=255)
         cGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAg                                      <  1:3214317/69‑1 (MQ=255)
         cGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAg                                      <  1:960592/69‑1 (MQ=255)
         cGCTGCGCTTCCTCTGGCGACAGTGACAAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAg                                      >  1:2772445/1‑69 (MQ=255)
                cTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAg                                 <  1:1986146/67‑1 (MQ=255)
                   ccTCTGGCGACAGTGCCCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATg                             >  1:1695551/1‑68 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAgg                                     >  1:2982570/1‑60 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:1447845/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:1124164/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:3058192/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:3069702/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:3215479/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:3318781/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:3386075/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:489199/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:562994/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:591905/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:621228/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:708024/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:772615/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:843432/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:913404/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:1518259/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:2835007/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:1798558/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:2284973/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:245024/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:2173704/1‑69 (MQ=255)
                   ccTCTGGCGACAGTGACCAGACAATCGAGCAAAGATAAGGATCGCTAAGCGGTAAAAAGGCCAGAATGc                            >  1:2133958/1‑69 (MQ=255)
                   ccTCTGGCGAAAGTGACCAGACAATCGAGCAAAGATACGGATc                                                      >  1:2909037/1‑43 (MQ=38)
                           gACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCa                    <  1:2820516/69‑1 (MQ=255)
                           gACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCa                    <  1:1033427/69‑1 (MQ=255)
                           gACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCa                    <  1:635844/69‑1 (MQ=255)
                           gACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCa                    <  1:714310/69‑1 (MQ=255)
                           gACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCa                    <  1:882821/69‑1 (MQ=255)
                           gACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCa                    <  1:1892562/69‑1 (MQ=255)
                                 gACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGaaa              >  1:2239254/1‑69 (MQ=255)
                                            cGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTc                        >  1:705289/1‑48 (MQ=255)
                                            cGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCcgcg    >  1:1157367/1‑68 (MQ=255)
                                            cGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCgcgc   >  1:483010/1‑69 (MQ=255)
                                            cGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCgcgc   >  1:2543301/1‑69 (MQ=255)
                                            cGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCgcgc   >  1:1196229/1‑69 (MQ=255)
                                            cGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCgcgc   >  1:1432858/1‑69 (MQ=255)
                                             gAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGcc  <  1:2101461/69‑1 (MQ=255)
                                             gAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGcc  <  1:1154773/69‑1 (MQ=255)
                                             gAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGcc  <  1:2314226/69‑1 (MQ=255)
                                             gAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGcc  <  1:577235/69‑1 (MQ=255)
                                             gAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGcc  <  1:2381700/69‑1 (MQ=255)
                                             gAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGcc  <  1:182218/69‑1 (MQ=255)
                                             gAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGcc  <  1:3257558/69‑1 (MQ=255)
                                             gAGCAAAGATACGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGcc  <  1:178531/69‑1 (MQ=255)
                                                       |                                                          
TGCTGCATCCGCTGCGCTTCCTCTGGCGACAGTGACCAGACAATCGAGCAAAGATGCGGATCGCTAAGCGGTAAAAAGGCCAGAATGCCTTCGCCATGGAAAACCTGCCGCGCC  >  minE/1950870‑1950983

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: