Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 880,753 G→A 51.1% intergenic (+72/‑122) yciM → / → pyrF conserved hypothetical protein/orotidine‑5'‑phosphate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE880,7530GA51.1% ‑6.2 / 66.4 45intergenic (+72/‑122)yciM/pyrFconserved hypothetical protein/orotidine‑5'‑phosphate decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (15/7);  new base A (11/12);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

TTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAATGCTCGC  >  minE/880684‑880822
                                                                     |                                                                     
tttAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGgc                                                                      >  1:6679/1‑71 (MQ=255)
tttAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGgc                                                                      >  1:720363/1‑71 (MQ=255)
tttAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGgc                                                                      >  1:2532717/1‑71 (MQ=255)
tttAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGgc                                                                      >  1:3992169/1‑71 (MQ=255)
        aaaTCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGaaaa                                                              <  1:2650354/71‑1 (MQ=255)
           tCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGc                                                           <  1:686069/71‑1 (MQ=255)
           tCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGc                                                           <  1:891489/71‑1 (MQ=255)
                cTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGcc                                                      <  1:915358/71‑1 (MQ=255)
                cTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGcc                                                      <  1:137728/71‑1 (MQ=255)
                   tAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAg                                                    <  1:1181153/70‑1 (MQ=255)
                         caacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGgaga                                              <  1:1864385/70‑1 (MQ=255)
                         caacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGgaga                                              <  1:408269/70‑1 (MQ=255)
                          aacaTACTAATTATTAATGTTCCGTTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGaaa                                            >  1:4045672/1‑71 (MQ=255)
                              tACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGc                                        <  1:2730455/71‑1 (MQ=255)
                                cTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa                                      >  1:3902355/1‑71 (MQ=255)
                                cTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa                                      >  1:193302/1‑71 (MQ=255)
                                cTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa                                      >  1:1636740/1‑71 (MQ=255)
                                    ttattaATGTTCCATTGTGCTCCGGCAACGACGACGCAGaa                                                                <  1:547236/41‑1 (MQ=255)
                                    ttattaATGTTCCATTGTGCTCCGGCAACGACGACGCAGaa                                                                <  1:3135587/41‑1 (MQ=255)
                                    ttattaATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGgaga                                              <  1:2655695/59‑1 (MQ=255)
                                    ttattaATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGgaga                                              <  1:3866670/59‑1 (MQ=255)
                                    ttattaATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGc                                        <  1:3344099/65‑1 (MQ=255)
                                      attaATGTTCCATTGTGCTCCGGCAACGACGACGCAGaaaa                                                              >  1:914953/1‑41 (MQ=255)
                                      attaATGTTCCATTGTGCTCCGGCAACGACGACGCAGaaaa                                                              >  1:3653553/1‑41 (MQ=255)
                                      attaATGTTCCATTGTGCTCCGGCAACGACGACGCAGaaaa                                                              >  1:1929480/1‑41 (MQ=255)
                                      attaATGTTCCATTGTGCTCCGGCAACGACGACGCAGaaaa                                                              >  1:2060598/1‑41 (MQ=255)
                                        taATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGc                                        <  1:431435/61‑1 (MQ=255)
                                          aTGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa                                      >  1:685277/1‑61 (MQ=255)
                                          aTGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa                                      >  1:1090347/1‑61 (MQ=255)
                                          aTGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa                                      >  1:752537/1‑61 (MQ=255)
                                          aTGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa                                      >  1:542723/1‑61 (MQ=255)
                                           tGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAAtttt                           >  1:1178124/1‑71 (MQ=255)
                                           tGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGAGAAATCGCAACTGTTAAttttt                          >  1:50430/1‑71 (MQ=255)
                                             ttCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa                                      >  1:2101127/1‑58 (MQ=255)
                                             ttCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa                                      >  1:3008136/1‑58 (MQ=255)
                                               ccATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTAt                       >  1:2194931/1‑71 (MQ=255)
                                                 aTTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGcc                                                      <  1:746297/38‑1 (MQ=255)
                                                  ttGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGc                                        <  1:430226/51‑1 (MQ=255)
                                                  ttGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGt                                  <  1:1969786/57‑1 (MQ=255)
                                                   tgtgCTCCGGCAACGACGACGCAGAAAAAGCCTGcc                                                      <  1:3320028/36‑1 (MQ=255)
                                                    gtgCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCa                  >  1:924570/1‑71 (MQ=255)
                                                    gtgCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCa                  >  1:925021/1‑71 (MQ=255)
                                                    gtgCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCa                  >  1:4080046/1‑71 (MQ=255)
                                                    gtgCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCa                  >  1:3122385/1‑71 (MQ=255)
                                                    gtgCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCa                  >  1:1221110/1‑71 (MQ=255)
                                                                 gacgGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa                                      >  1:2281718/1‑38 (MQ=255)
                                                                 gacgGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACg                >  1:3263852/1‑60 (MQ=255)
                                                                 gacgGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACg                >  1:2522280/1‑60 (MQ=255)
                                                                 gacgGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAATGCt     >  1:3520444/1‑71 (MQ=255)
                                                                    ggCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAATGCTCGc  <  1:619196/71‑1 (MQ=255)
                                                                     |                                                                     
TTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAATGCTCGC  >  minE/880684‑880822

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: