Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 304,425 C→A 31.7% intergenic (+5/+26) ybaY → / ← ybaZ predicted outer membrane lipoprotein/predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE304,4250CA31.7% 48.1 / 23.0 41intergenic (+5/+26)ybaY/ybaZpredicted outer membrane lipoprotein/predicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (6/22);  new base A (8/5);  total (14/27)
Fisher's exact test for biased strand distribution p-value = 1.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCG  >  minE/304358‑304481
                                                                   |                                                        
ccTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGAc                                                          >  1:642791/1‑68 (MQ=255)
      aCTTCACCAACCCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCAACGCAATATTAAGCGAAGTTaa                                                     <  1:2888222/67‑1 (MQ=255)
      aCTTCACCAACCCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCAACGCAATATTAAGCGAAGTTaa                                                     <  1:2229039/67‑1 (MQ=255)
      aCTTCACCAACCCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCAACGCAATATTAAGCGAAGTTaa                                                     <  1:2101113/67‑1 (MQ=255)
      aCTTCACCAACCCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCAACGCAATATTAAGCGAAGTTaa                                                     <  1:1836926/67‑1 (MQ=255)
                     gTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGga                                      <  1:1102567/67‑1 (MQ=255)
                     gTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGga                                      <  1:1168919/67‑1 (MQ=255)
                     gTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGga                                      <  1:91973/67‑1 (MQ=255)
                     gTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGga                                      <  1:2275127/67‑1 (MQ=255)
                     gTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGga                                      <  1:1952668/67‑1 (MQ=255)
                     gTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGga                                      <  1:1785254/67‑1 (MQ=255)
                         aTCCGTCTTCGGCAGTTCCAGCTCCAACGCAATATTAAGCGAAGTTaa                                                     <  1:2490006/48‑1 (MQ=255)
                          tCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTaa                                                     >  1:350223/1‑47 (MQ=255)
                          tCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCg                                        >  1:563775/1‑60 (MQ=255)
                          tCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCg                                        >  1:2922158/1‑60 (MQ=255)
                             gTCTTCGGCAGTTCCAGCTCCAACGCAATATTAAGCGAAGt                                                        >  1:2647194/1‑41 (MQ=38)
                             gTCTTCGGCAGTTCCAGCTCCAACGCAATATTAAGCGAAGTTAACCCTCTCCGCCCGGAGAGGGTtagt                            >  1:1907846/1‑69 (MQ=39)
                             gTCTTCGGCAGTTCCAGCTCCAACGCAATATTAAGCGAAGTTAACCCTCTCCGCCCGGAGAGGGTtagt                            >  1:821594/1‑69 (MQ=39)
                             gTCTTCGGCAGTTCCAGCTCCAACGCAATATTAAGCGAAGTTAACCCTCTCCGCCCGGAGAGGGTtagt                            >  1:2931083/1‑69 (MQ=39)
                             gTCTTCGGCAGTTCCAGCTCCAACGCAATATTAAGCGAAGTTAACCCTCTCCGCCCGGAGAGGGTtagt                            >  1:2535327/1‑69 (MQ=39)
                             gTCTTCGGCAGTTCCAGCTCCAACGCAATATTAAGCGAAGTTAACCCTCTCCGCCCGGAGAGGGTtagt                            >  1:2509420/1‑69 (MQ=39)
                             gTCTTCGGCAGTTCCAGCTCCAACGCAATATTAAGCGAAGTTAACCCTCTCCGCCCGGAGAGGGTtagt                            >  1:239282/1‑69 (MQ=39)
                             gTCTTCGGCAGTTCCAGCTCCAACGCAATATTAAGCGAAGTTAACCCTCTCCGCCCGGAGAGGGTtagt                            >  1:1317365/1‑69 (MQ=39)
                              tCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagt                            <  1:2412281/68‑1 (MQ=255)
                              tCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagt                            <  1:787635/68‑1 (MQ=255)
                              tCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagt                            <  1:1193829/68‑1 (MQ=255)
                              tCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagt                            <  1:3010772/68‑1 (MQ=255)
                              tCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagt                            <  1:461188/68‑1 (MQ=255)
                              tCTTCGGCAGGTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagt                            <  1:3250199/68‑1 (MQ=255)
                                         tCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGAt                <  1:950496/69‑1 (MQ=255)
                                         tCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGAt                <  1:1012051/69‑1 (MQ=255)
                                         tCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGAt                <  1:1713562/69‑1 (MQ=255)
                                         tCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGAt                <  1:1616885/69‑1 (MQ=255)
                                         tCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGAt                <  1:2896846/69‑1 (MQ=255)
                                                  tACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCa       >  1:1240375/1‑69 (MQ=255)
                                                       aaTATTAAGCGACGTTAACCCTCTCCTCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCg  <  1:1525333/69‑1 (MQ=255)
                                                       aaTATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCg  <  1:388720/69‑1 (MQ=255)
                                                       aaTATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCg  <  1:191/69‑1 (MQ=255)
                                                       aaTATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCg  <  1:133256/69‑1 (MQ=255)
                                                       aaTATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCg  <  1:1190202/69‑1 (MQ=255)
                                                           ttAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAAc             >  1:3425133/1‑54 (MQ=255)
                                                           ttAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAAc             >  1:3271880/1‑54 (MQ=255)
                                                                   |                                                        
CCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCG  >  minE/304358‑304481

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: