Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 660,873 G→A 16.6% V76V (GTG→GTA ycbF → predicted periplasmic pilini chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE660,8730GA16.6% 56.5 / 4.3 24V76V (GTG→GTAycbFpredicted periplasmic pilini chaperone
Reads supporting (aligned to +/- strand):  ref base G (5/15);  new base A (2/2);  total (7/17)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTG  >  minE/660810‑660942
                                                               |                                                                     
cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAc                                                                >  1:2082974/1‑71 (MQ=255)
        ttaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAg                                                                   >  1:2473976/1‑60 (MQ=255)
                 aaaTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTACCAGCACCATTACTGGCCGCTCCg                                               >  1:2053594/1‑71 (MQ=255)
                   aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGc                                              >  1:832407/1‑70 (MQ=255)
                      aaCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTACCAGCAcc                                                               <  1:2361143/50‑1 (MQ=255)
                      aaCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTACCAGCACCATTACTGGCCGCTCCGCCACt                                          >  1:414500/1‑71 (MQ=255)
                         cGCCCAACGGGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                                                       <  1:1488591/39‑1 (MQ=255)
                         cGCCCAACGCGTGGGGCCGGGGTGGAAGCGGCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc                                             <  1:2065874/65‑1 (MQ=255)
                                    ttGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcc                                                     <  1:4093854/46‑1 (MQ=255)
                                    ttGGGCCGGGTGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGcc                           <  1:992161/71‑1 (MQ=255)
                                     tGGGCCGGGGGGGAAGCGTCGACAGTACCAGCACCATTACt                                                         <  1:2462314/41‑1 (MQ=255)
                                        gCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGc                                              <  1:2823946/49‑1 (MQ=255)
                                          cGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGc                                                   <  1:3662519/42‑1 (MQ=255)
                                          cGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTc                                                 <  1:3681943/44‑1 (MQ=255)
                                          cGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtct                                    <  1:3814014/57‑1 (MQ=255)
                                             gggggAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGa                   <  1:3650069/71‑1 (MQ=255)
                                              ggggAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGa                   <  1:3701992/70‑1 (MQ=255)
                                               gggAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGa                   <  1:3374172/69‑1 (MQ=255)
                                                 gAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt                                         <  1:1463219/45‑1 (MQ=255)
                                                    gCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTAc                     <  1:3493556/62‑1 (MQ=255)
                                                     cGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGc              <  1:3497934/68‑1 (MQ=255)
                                                        cGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTa                                        >  1:995139/1‑39 (MQ=255)
                                                           cAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCa             >  1:1704236/1‑63 (MQ=255)
                                                            aGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACg                    <  1:930404/55‑1 (MQ=255)
                                                              tGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTg  <  1:2121298/71‑1 (MQ=255)
                                                               |                                                                     
CTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTG  >  minE/660810‑660942

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: