Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 860,344 A→C 57.6% T519T (ACT→ACG trpE ← component I of anthranilate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE860,3440AC57.6% 34.8 / 60.3 66T519T (ACT→ACGtrpEcomponent I of anthranilate synthase
Reads supporting (aligned to +/- strand):  ref base A (9/19);  new base C (18/20);  total (27/39)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

TGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTAC  >  minE/860280‑860402
                                                                |                                                          
tGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGtct                                                        <  1:2829279/69‑1 (MQ=255)
tGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGtct                                                        <  1:3202866/69‑1 (MQ=255)
tGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGtct                                                        <  1:2070492/69‑1 (MQ=255)
tGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGtct                                                        <  1:523122/69‑1 (MQ=255)
tGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGtct                                                        <  1:1955262/69‑1 (MQ=255)
tGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGtct                                                        <  1:1273756/69‑1 (MQ=255)
           gTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCatga                                             >  1:3213245/1‑69 (MQ=255)
           gTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCatga                                             >  1:272295/1‑69 (MQ=255)
           gTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCatga                                             >  1:2024916/1‑69 (MQ=255)
             tGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATTAGAACGTCTCCTGTGc                                                 >  1:1729004/1‑63 (MQ=255)
             tGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATTAGAACGTCTCCTGTGCa                                                >  1:2498542/1‑64 (MQ=255)
             tGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATTAGAACGTCTCCTGTGCa                                                >  1:2347816/1‑64 (MQ=255)
             tGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATTAGAACGTCTCCTGTGCa                                                >  1:2123767/1‑64 (MQ=255)
             tGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATTAGAACGTCTCCTGTGCa                                                >  1:2777144/1‑64 (MQ=255)
             tGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATTAGAACGTCTCCTGTGCa                                                >  1:2262597/1‑64 (MQ=255)
             tGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATTAGAACGTCTCCTGTGCa                                                >  1:1750978/1‑64 (MQ=255)
             tGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATTAGAACGTCTCCTGTGCa                                                >  1:1577276/1‑64 (MQ=255)
             tGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATTAGAACGTCTCCTGTGCa                                                >  1:1682718/1‑64 (MQ=255)
             tGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATTAGAACGTCTCCTGTGCa                                                >  1:161617/1‑64 (MQ=255)
             tGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATTAGAACGTCTCCTGTGCa                                                >  1:3113416/1‑64 (MQ=255)
             tGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATTAGAACGTCTCCTGTGCa                                                >  1:1274510/1‑64 (MQ=255)
             tGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATTAGAACGTCTCCTGTGCa                                                >  1:379127/1‑64 (MQ=255)
             tGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATTAGAACGTCTCCTGTGCa                                                >  1:1216900/1‑64 (MQ=255)
             tGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATTAGAACGTCTCCTGTGCa                                                >  1:712608/1‑64 (MQ=255)
             tGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATTAGAACGTCTCCTGTGCa                                                >  1:760115/1‑64 (MQ=255)
             tGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATTAGAACGTCTCCTGAGCa                                                >  1:498992/1‑64 (MQ=255)
             tGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATTAGAACGTCCCCTGTGCa                                                >  1:115433/1‑64 (MQ=255)
              gTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCatgatg                                           >  1:760473/1‑68 (MQ=255)
              gTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATgc                                          >  1:808009/1‑69 (MQ=255)
                 cGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGcgcg                                       >  1:1058352/1‑69 (MQ=255)
                 cGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGcgcg                                       >  1:597961/1‑69 (MQ=255)
                 cGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGcgcg                                       >  1:1448499/1‑69 (MQ=255)
                  gTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATgcgc                                        <  1:797916/67‑1 (MQ=255)
                       agagTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGc                                  <  1:3211059/68‑1 (MQ=255)
                       agagTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGc                                  <  1:1138000/68‑1 (MQ=255)
                       agagTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGc                                  <  1:1505209/68‑1 (MQ=255)
                        gagTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCaa                                <  1:1117111/69‑1 (MQ=255)
                           tCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGtct                                                        <  1:1510488/42‑1 (MQ=255)
                                              gAATGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCTGTACAGCGCGGGc          <  1:1261639/69‑1 (MQ=255)
                                              gAATGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGc          <  1:2079218/69‑1 (MQ=255)
                                              gAATGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGc          <  1:2708357/69‑1 (MQ=255)
                                              gAATGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGc          <  1:23537/69‑1 (MQ=255)
                                              gAATGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGc          <  1:78283/69‑1 (MQ=255)
                                              gAATGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGc          <  1:1095178/69‑1 (MQ=255)
                                              gAATGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGc          <  1:2447646/69‑1 (MQ=255)
                                              gAATGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGc          <  1:688358/69‑1 (MQ=255)
                                              gAATGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGc          <  1:2027311/69‑1 (MQ=255)
                                              gAATGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGc          <  1:534621/69‑1 (MQ=255)
                                              gAATGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGc          <  1:118840/69‑1 (MQ=255)
                                              gAATGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGc          <  1:470647/69‑1 (MQ=255)
                                              gAATGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGc          <  1:328834/69‑1 (MQ=255)
                                              gAATGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGc          <  1:2625131/69‑1 (MQ=255)
                                              gAATGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGc          <  1:2251513/69‑1 (MQ=255)
                                              gAATGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGc          <  1:3003568/69‑1 (MQ=255)
                                              gAATGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGc          <  1:1915957/69‑1 (MQ=255)
                                              gAATGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGc          <  1:2769890/69‑1 (MQ=255)
                                               aaTGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGc          <  1:1201076/68‑1 (MQ=255)
                                               aaTGTCAGCCATTAGAACGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGc          <  1:3061788/68‑1 (MQ=255)
                                                   tCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGt     <  1:354856/69‑1 (MQ=255)
                                                   tCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGt     <  1:2665921/69‑1 (MQ=255)
                                                   tCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGt     <  1:571651/69‑1 (MQ=255)
                                                   tCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGt     <  1:590950/69‑1 (MQ=255)
                                                   tCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGt     <  1:3162311/69‑1 (MQ=255)
                                                   tCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGt     <  1:2902089/69‑1 (MQ=255)
                                                   tCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGt     <  1:1935637/69‑1 (MQ=255)
                                                      gCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTAc  >  1:1885557/1‑69 (MQ=255)
                                                                |                                                          
TGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTAC  >  minE/860280‑860402

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: