Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 373,918 T→A 18.7% V256V (GTT→GTA fepE → regulator of length of O‑antigen component of lipopolysaccharide chains

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE373,9180TA18.7% 90.5 / 14.2 48V256V (GTT→GTAfepEregulator of length of O‑antigen component of lipopolysaccharide chains
Reads supporting (aligned to +/- strand):  ref base T (19/20);  new base A (5/4);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

GCAAACATTCAGCGCCTCAATTATTCACTCGACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGTTAAAGATGACCCCGATTTTTCTATTTCTCTCGGTGCAGACGGTA  >  minE/373853‑373982
                                                                 |                                                                
gCAAACATTCAGCGCCTCAATTATTCACTCGACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTTTa                                                                <  1:3214134/68‑1 (MQ=255)
gCAAACATTCAGCGCCTCAATTATTCACTCGACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTTTa                                                                <  1:1773977/68‑1 (MQ=255)
 cAAACATTCAGCGCCTCAATTATTCACTTGACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTATAc                                                               >  1:972565/1‑68 (MQ=255)
 cAAACATTCAGCGCCTCAATTATTCACTTGACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTATAc                                                               >  1:813982/1‑68 (MQ=255)
 cAAACATTCAGCGCCTCAATTATTCACTTGACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTATAc                                                               >  1:1576180/1‑68 (MQ=255)
 cAAACATTCAGCGCCTCAATTATTCACTTGACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTATAc                                                               >  1:1801563/1‑68 (MQ=255)
 cAAACATTCAGCGCCTCAATTATTCACTTGACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTATAc                                                               >  1:302781/1‑68 (MQ=255)
         cAGCGCCTCAATTATTCACTCGACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTTTACAGTAATg                                                        >  1:1324746/1‑67 (MQ=255)
         cAGCGCCTCAATTATTCACTCGACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTTTACAGTAATg                                                        >  1:274603/1‑67 (MQ=255)
                     tattCACTCGACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGt                                              >  1:2054102/1‑65 (MQ=255)
                     tattCACTCGACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGt                                              >  1:1053824/1‑65 (MQ=255)
                     tattCACTCGACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGt                                              >  1:1205209/1‑65 (MQ=255)
                         cACTCGACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGTTAAAGAt                                       <  1:858066/68‑1 (MQ=255)
                         cACTCGACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGTTAAAGAt                                       <  1:1915345/68‑1 (MQ=255)
                         cACTCGACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGTTAAAGAt                                       <  1:588534/68‑1 (MQ=255)
                            ttgACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTATACAGTAATGGTCAGGCCGTTAAAGATGAc                                    <  1:2323616/66‑1 (MQ=255)
                            ttgACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTATACAGTAATGGTCAGGCCGTTAAAGATGAc                                    <  1:315675/66‑1 (MQ=255)
                            ttgACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTATACAGTAATGGTCAGGCCGTTAAAGATGAc                                    <  1:3363414/66‑1 (MQ=255)
                            ttgACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTATACAGTAATGGTCAGGCCGTTAAAGATGAc                                    <  1:3417613/66‑1 (MQ=255)
                                    gCCAACGCGGCAGGAATTAAAAAGCCCGTTTACAGTAATGGt                                                      >  1:195092/1‑42 (MQ=255)
                                    gCCAACGCGGCAGGAATTAAAAAGCCCGTTTACAGTAATGGt                                                      >  1:830492/1‑42 (MQ=255)
                                    gCCAACGCGGCAGGAATTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGTTAAAGATGACCCCGAtttt                           >  1:1905283/1‑69 (MQ=255)
                                    gCCAACGCGGCAGGAATTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGTTAAAGATGACCCCGAtttt                           >  1:772675/1‑69 (MQ=255)
                                    gCCAACGCGGCAGGAATTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGTTAAAGATGACCCCGAtttt                           >  1:1692251/1‑69 (MQ=255)
                                    gCCAACGCGGCAGGAATTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGTTAAAGATGACCCCGAtttt                           >  1:476217/1‑69 (MQ=255)
                                            ggCAGGAATTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGTTAAAGAt                                       <  1:1974727/49‑1 (MQ=255)
                                                   aTTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGTTaaa                                          >  1:1914911/1‑39 (MQ=255)
                                                   aTTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGTTaaa                                          >  1:3099233/1‑39 (MQ=255)
                                                   aTTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGTTaaa                                          >  1:3147766/1‑39 (MQ=255)
                                                   aTTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGTTaaa                                          >  1:478772/1‑39 (MQ=255)
                                                   aTTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGTTaaa                                          >  1:875030/1‑39 (MQ=255)
                                                   aTTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGTTaaa                                          >  1:359173/1‑39 (MQ=255)
                                                   aTTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGTTaaa                                          >  1:100665/1‑39 (MQ=255)
                                                   aTTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGTTAAAGATGACCCCGATTTTTCTATTTCTCTCgg             >  1:935991/1‑68 (MQ=255)
                                                           gCCCGTTTACAGTAATGGTCAGGCCGTTAAAGATGACCCCGATTTTTCTATTTCTCTCGGTGCAGACgg    <  1:1043761/69‑1 (MQ=255)
                                                           gCCCGTTTACAGTAATGGTCAGGCCGTTAAAGATGACCCCGATTTTTCTATTTCTCTCGGTGCAGACgg    <  1:837667/69‑1 (MQ=255)
                                                            cccGTTTACAGTAATGGTCAGGCCGTTAAAGATGACCCCGATTTTTCTATTTCTCTCGGTGCAGACGGt   <  1:2187079/69‑1 (MQ=255)
                                                            cccGTTTACAGTAATGGTCAGGCCGTTAAAGAAGACCCCGATTTTTCTATTTCTCTCGGTGCAGACGGt   <  1:3163521/69‑1 (MQ=255)
                                                             ccGTTTACAGTAATGGTCAGGCCGTTAAAGATGACCCCGATTTTTCTATTTCTCTCGGTTCAGACGGTa  <  1:3166771/69‑1 (MQ=255)
                                                             ccGTTTACAGTAATGGTCAGGCCGTTAAAGATGACCCCGATTTTTCTATTTCTCTCGGTGCAGACGGTa  <  1:3323474/69‑1 (MQ=255)
                                                             ccGTTTACAGTAATGGTCAGGCCGTTAAAGATGACCCCGATTTTTCTATTTCTCTCGGTGCAGACGGTa  <  1:2961762/69‑1 (MQ=255)
                                                             ccGTTTACAGTAATGGTCAGGCCGTTAAAGATGACCCCGATTTTTCTATTTCTCTCGGTGCAGACGGTa  <  1:1726374/69‑1 (MQ=255)
                                                             ccGTTTACAGTAATGGTCAGGCCGTTAAAGATGACCCCGATTTTTCTATTTCTCTCGGTGCAGACGGTa  <  1:1717251/69‑1 (MQ=255)
                                                             ccGTTTACAGTAATGGTCAGGCCGTTAAAGATGACCCCGATTTTTCTATTTCTCTCGGTGCAGACGGTa  <  1:171206/69‑1 (MQ=255)
                                                             ccGTTTACAGTAATGGTCAGGCCGTTAAAGATGACCCCGATTTTTCTATTTCTCTCGGTGCAGACGGTa  <  1:152732/69‑1 (MQ=255)
                                                             ccGTTTACAGTAATGGTCAGGCCGTTAAAGATGACCCCGATTTTTCTATTTCTCTCGGTGCAGACGGTa  <  1:1238488/69‑1 (MQ=255)
                                                             ccGTTTACAGTAATGATCAGGCCGTTAAAGATGACCCCGATTTTTCTATTTCTCTCGGTGCAGACGGTa  <  1:3107018/69‑1 (MQ=255)
                                                              cGTTTACAGTAATGGTCAGGCCGTTAAAGATGACCCCGATTTTTCTATTTCTCTCGGTGCAGACGGTa  <  1:1157805/68‑1 (MQ=255)
                                                                 |                                                                
GCAAACATTCAGCGCCTCAATTATTCACTCGACATTGCCAACGCGGCAGGAATTAAAAAGCCCGTTTACAGTAATGGTCAGGCCGTTAAAGATGACCCCGATTTTTCTATTTCTCTCGGTGCAGACGGTA  >  minE/373853‑373982

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: