Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 880,747 A→C 67.7% intergenic (+66/‑128) 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,7470AC67.7% 34.7 / 18.3 34intergenic (+66/‑128)yciM/pyrFconserved hypothetical protein/orotidine‑5'‑phosphate decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (8/3);  new base C (12/11);  total (20/14)
Fisher's exact test for biased strand distribution p-value = 2.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

ATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCAC  >  minE/880681‑880807
                                                                  |                                                            
aTTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAac                                                             <  1:2625722/68‑1 (MQ=255)
aTTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAac                                                             <  1:1758095/68‑1 (MQ=255)
   tttAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAacgacg                                                         >  1:2505746/1‑69 (MQ=255)
   tttAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAacgacg                                                         >  1:2361791/1‑69 (MQ=255)
   tttAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAacgacg                                                         >  1:2680813/1‑69 (MQ=255)
           aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa                                                 <  1:1763014/69‑1 (MQ=38)
           aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa                                                 <  1:2587617/69‑1 (MQ=38)
           aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa                                                 <  1:2026931/69‑1 (MQ=38)
           aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa                                                 <  1:2993575/69‑1 (MQ=38)
           aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa                                                 <  1:3076454/69‑1 (MQ=38)
           aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa                                                 <  1:1436205/69‑1 (MQ=38)
           aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa                                                 <  1:3135649/69‑1 (MQ=38)
           aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa                                                 <  1:3211392/69‑1 (MQ=38)
           aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa                                                 <  1:3222331/69‑1 (MQ=38)
           aaaTCCGACTTTAGTTACAACATACTAATTATTGATGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa                                                 <  1:3417501/69‑1 (MQ=38)
                           acaacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAacgacg                                                         >  1:126230/1‑45 (MQ=255)
                           acaacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAg                                             >  1:1141794/1‑57 (MQ=255)
                           acaacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAg                                             >  1:2740941/1‑57 (MQ=255)
                           acaacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAg                                             >  1:2554712/1‑57 (MQ=255)
                           acaacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAg                                             >  1:695882/1‑57 (MQ=255)
                           acaacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAg                                             >  1:1039924/1‑57 (MQ=255)
                           acaacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAg                                             >  1:2054759/1‑57 (MQ=255)
                           acaacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAg                                             >  1:1279016/1‑57 (MQ=255)
                           acaacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAg                                             >  1:1685367/1‑57 (MQ=255)
                                             aTGTTCTATTGTGCTCCGGCACCGACGGCGCGGaa                                                 <  1:633800/35‑1 (MQ=25)
                                                          cTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCAc  <  1:2901481/69‑1 (MQ=255)
                                                          cTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCAc  <  1:1801186/69‑1 (MQ=255)
                                                          cTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCAc  <  1:1432217/69‑1 (MQ=255)
                                                              ggCACCGACGGCGCGGAAAAGCCTGTCAGAGGAGAAATCGCAACTGTTAAtt              >  1:1201482/1‑52 (MQ=255)
                                                              ggCACCGACGGCGCGGAAAAGCCTGCCAGAGGAGAAATCGCACCTGTTAAtt              >  1:2380115/1‑52 (MQ=25)
                                                              ggCACCGACGGCGCGGAAAAGCCTGCCAGAGGAGAAATCGCAACTGtt                  >  1:3088495/1‑48 (MQ=25)
                                                              ggCACCGACGGCGCGGAAAAGCCTGCCAGAGGAGAAATCGCAACTGTTAAtt              >  1:294435/1‑52 (MQ=37)
                                                              ggCACCGACGGCGCGGAAAAGCCTGCCAGAGGAGAAATCGCAACTGTTAAtt              >  1:1251027/1‑52 (MQ=37)
                                                              ggCACCGACGGCGCGGAAAAGCCTGCCAGAGGAGAAATCGCAACTGTTAAtt              >  1:550660/1‑52 (MQ=37)
                                                              ggCACCGACGGCGCGGAAAAGCCTGCCAGAGGAGAAATCGCAACTGTTAAtt              >  1:627587/1‑52 (MQ=37)
                                                              ggCACCGACGGCGCGGAAAAGCCTGCCAGAGGAGAAATCGCAACTGTTAAtt              >  1:2167438/1‑52 (MQ=37)
                                                              ggCACCGACGGCGCGGAAAAGCCTGCCAGAGGAGAAATCGCAACTGTTAAtt              >  1:674426/1‑52 (MQ=37)
                                                              ggCACCGACGGCGCGGAAAAGCCTGCCAGAGGAGAAATCGCAACTGTTAAtt              >  1:2418404/1‑52 (MQ=37)
                                                              ggCACCGACGGCGCGGAAAAGCCTGCCAGAGGAGAAATCGCAACTGTTAAtt              >  1:967431/1‑52 (MQ=37)
                                                              ggCACCGACGGCGCGGAAAAGCCTGCCAGAGGAGAAATCGCAAATGTTAAtt              >  1:1870687/1‑52 (MQ=25)
                                                                  |                                                            
ATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCAC  >  minE/880681‑880807

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: