Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 28,911 T→C 33.3% I139I (ATT→ATC carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE28,9110TC33.3% 56.8 / 40.0 75I139I (ATT→ATCcarBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base T (33/17);  new base C (21/4);  total (54/21)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.40e-01

TGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCTGACG  >  minE/28845‑28978
                                                                  |                                                                   
tGCCACTAACGATACGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTg                                                                    <  1:2271783/68‑1 (MQ=255)
      tGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGa                                                             >  1:2289096/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGCCGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:1186057/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:2957700/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:2274541/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:2535721/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:2625017/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:2781537/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:2834490/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:2935317/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:2741196/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:3021128/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:316107/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:3190657/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:3274641/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:445573/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:657651/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:993821/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:1855122/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:1124829/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:1363758/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:1536973/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:1720787/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCAGGaa                                                            >  1:1481257/1‑69 (MQ=255)
       gCCGATGCGATTGATAAAGCAGAAGACAGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:1082443/1‑69 (MQ=255)
       gCCGATACGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGaa                                                            >  1:2235/1‑69 (MQ=255)
            tGCGATTGATAAAGCCGAAGACCGCCGCCGTTTCGACGTAGCGATGAAGAAAATCGGTCTGGAAACtgc                                                       <  1:3105160/69‑1 (MQ=255)
            tGCGATTGATAAAGCCGAAGACCGCCGCCGTTTCGACGTAGCGATGAAGAAAATCGGTCTGGAAACtgc                                                       <  1:3098739/69‑1 (MQ=255)
            tGCGATTGATAAAGCCGAAGACCGCCGCCGTTTCGACGTAGCGATGAAGAAAATCGGTCTGGAAACtgc                                                       <  1:3052764/69‑1 (MQ=255)
            tGCGATTGATAAAGCCGAAGACCGCCGCCGTTTCGACGTAGCGATGAAGAAAATCGGTCTGGAAACtgc                                                       <  1:3023660/69‑1 (MQ=255)
                            gaagaCCGCCGTCGTTTCGCCGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGca                                       >  1:1891921/1‑69 (MQ=255)
                               gaCCGCCGTCGTTTCGACGTAGCGATGATGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGcacac                                    >  1:2452215/1‑69 (MQ=255)
                               gaCCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGcacac                                    >  1:2707994/1‑69 (MQ=255)
                               gaCCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGcacac                                    >  1:538553/1‑69 (MQ=255)
                               gaCCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGcacac                                    >  1:2899471/1‑69 (MQ=255)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:3299193/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:1215641/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:1262207/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:31166/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:1984240/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:3047085/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:1163105/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:1849401/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:1464513/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:1802774/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:2880289/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:2854361/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:770544/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:2831336/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:1624198/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:2775678/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:2720360/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:1656523/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:2544170/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGaa                                                            >  1:2482673/1‑40 (MQ=38)
                                    ccgccgTTTCGACGTAGCGATGAAGAAAATCGGTCTGGa                                                             >  1:3117627/1‑39 (MQ=37)
                                         gtTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCgg                                              <  1:670679/49‑1 (MQ=255)
                                         gtTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCgg                                              <  1:835345/49‑1 (MQ=255)
                                         gtTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCgg                                              <  1:404553/49‑1 (MQ=255)
                                         gtTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCgg                                              <  1:950284/49‑1 (MQ=255)
                                         gtTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCgg                                              <  1:210307/49‑1 (MQ=255)
                                         gtTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCgg                                              <  1:1203445/49‑1 (MQ=255)
                                         gtTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCgg                                              <  1:2209339/49‑1 (MQ=255)
                                         gtTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCgg                                              <  1:2908339/49‑1 (MQ=255)
                                         gtTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCgg                                              <  1:1737967/49‑1 (MQ=255)
                                         gtTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCgg                                              <  1:2325127/49‑1 (MQ=255)
                                         ctttCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCgg                                              <  1:3324306/48‑1 (MQ=255)
                                                    gCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCg               >  1:348272/1‑69 (MQ=255)
                                                    gCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCg               >  1:2839256/1‑69 (MQ=255)
                                                    gCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCg               >  1:826836/1‑69 (MQ=255)
                                                       aTGAAGAAAATTGGGCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGtt            <  1:3339003/69‑1 (MQ=255)
                                                             aaaaTTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCt      <  1:1999153/69‑1 (MQ=255)
                                                             aaaaTTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCt      <  1:657610/69‑1 (MQ=255)
                                                             aaaaTTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCt      <  1:1962675/69‑1 (MQ=255)
                                                             aaaaTTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCt      <  1:1035517/69‑1 (MQ=255)
                                                                 ttGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCTGACg  <  1:1954962/69‑1 (MQ=255)
                                                                  |                                                                   
TGCCACTGCCGATGCGATTGATAAAGCAGAAGACCGCCGTCGTTTCGACGTAGCGATGAAGAAAATTGGTCTGGAAACCGCGCGTTCCGGTATCGCACACACGATGGAAGAAGCGCTGGCGGTTGCCGCTGACG  >  minE/28845‑28978

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: