Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,322,814 A→G 39.4% I40I (ATT→ATC galF ← predicted subunit with GalU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,322,8140AG39.4% 52.4 / 79.5 114I40I (ATT→ATCgalFpredicted subunit with GalU
Reads supporting (aligned to +/- strand):  ref base A (17/52);  new base G (31/14);  total (48/66)
Fisher's exact test for biased strand distribution p-value = 5.35e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 114 total aligned reads displayed.

GCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGT  >  minE/1322749‑1322873
                                                                  |                                                           
gcgTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATg                                                          <  1:3266192/69‑1 (MQ=255)
gcgTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATg                                                          <  1:896299/69‑1 (MQ=255)
gcgTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATg                                                          <  1:332532/69‑1 (MQ=255)
gcgTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATg                                                          <  1:1354303/69‑1 (MQ=255)
gcgTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATg                                                          <  1:432248/69‑1 (MQ=255)
gcgTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATg                                                          <  1:369281/69‑1 (MQ=255)
gcgTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTAGTCAACAATg                                                          <  1:1093737/69‑1 (MQ=255)
         gaagcgTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACGATGTACTGAATc                                                 >  1:586643/4‑69 (MQ=255)
         gaagcgTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACGATGTACTGAATc                                                 >  1:1505847/4‑69 (MQ=255)
         gaagcgTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACGATGTACTGAATc                                                 >  1:1811701/4‑69 (MQ=255)
         gaagcgTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACGATGTACTGAATc                                                 >  1:1829900/4‑69 (MQ=255)
         gaagcgTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACGATGTACTGAATc                                                 >  1:1835454/4‑69 (MQ=255)
         gaagcgTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACGATGTACTGAATc                                                 >  1:769044/4‑69 (MQ=255)
         gaagcgTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACGATGTACTGAATc                                                 >  1:1957152/4‑69 (MQ=255)
         gaagcgTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACGATGTACTGAATc                                                 >  1:2132600/4‑69 (MQ=255)
         gaagcgTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACGATGTACTGAATc                                                 >  1:886721/4‑69 (MQ=255)
         gaagcgTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACGATGTACTGAATc                                                 >  1:2470214/4‑69 (MQ=255)
         gaagcgTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACGATGTACTGAATc                                                 >  1:2830712/4‑69 (MQ=255)
         gaagcgAGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACGATGTACTGAATc                                                 >  1:435584/4‑69 (MQ=255)
          aCGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAAt                                                  <  1:3335036/67‑1 (MQ=255)
          aCGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAAt                                                  <  1:145675/67‑1 (MQ=255)
                          gaggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAAATGTACTGAATc                                                 >  1:3001551/1‑53 (MQ=255)
                          gaggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATc                                                 >  1:600160/1‑52 (MQ=255)
                          gaggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATc                                                 <  1:53094/52‑1 (MQ=255)
                          gaggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATc                                                 >  1:3051187/1‑52 (MQ=255)
                          gaggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATc                                                 >  1:1505599/1‑52 (MQ=255)
                          gaggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATc                                                 >  1:2512067/1‑52 (MQ=255)
                          gaggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATc                                                 >  1:2299732/1‑52 (MQ=255)
                          gaggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATc                                                 <  1:2259220/52‑1 (MQ=255)
                          gaggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATc                                                 >  1:1217571/1‑52 (MQ=255)
                          gaggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATc                                                 >  1:1235418/1‑52 (MQ=255)
                          gaggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATc                                                 >  1:1925587/1‑52 (MQ=255)
                          gaggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTCAATc                                                 >  1:2822262/1‑52 (MQ=255)
                          gaggaTTTCTTTGATCCCTGAAGCCACAATCTCGTCAACAATGTACTGatt                                                  >  1:11382/1‑49 (MQ=255)
                           aggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATc                                                 <  1:1506860/51‑1 (MQ=255)
                               tttCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGt                            >  1:3018617/1‑68 (MQ=255)
                               tttCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGt                            >  1:3151088/1‑68 (MQ=255)
                                ttCTTTGATCCCTGCAGCCACAATCTCGTCAACAATg                                                          >  1:659967/1‑37 (MQ=255)
                                 tCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGc                         <  1:281044/69‑1 (MQ=255)
                                 tCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGc                         <  1:2717993/69‑1 (MQ=255)
                                 tCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGc                         <  1:1369868/69‑1 (MQ=255)
                                 tCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGc                         <  1:2357085/69‑1 (MQ=255)
                                 tCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCCACGATTGGTAGc                         <  1:500841/69‑1 (MQ=255)
                                 tCTTTGATCCCTGCAGCCACAATCTCATCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGc                         <  1:955442/69‑1 (MQ=255)
                                       aTCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCAtct                     <  1:932872/67‑1 (MQ=255)
                                       aTCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCAtct                     <  1:1043551/67‑1 (MQ=255)
                                       aTCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCAtct                     <  1:1906946/67‑1 (MQ=255)
                                       aTCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATCGGCTTGTCGACGATTGGTAGCAtct                     <  1:334774/67‑1 (MQ=255)
                                       aTCCCTGCAGCCACAATCTCATCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCAtct                     <  1:128591/67‑1 (MQ=255)
                                                  cacaATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGc        <  1:3142137/69‑1 (MQ=255)
                                                  cacaATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGc        <  1:1522053/69‑1 (MQ=255)
                                                  cacaATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGc        <  1:623961/69‑1 (MQ=255)
                                                  cacaATCTCGTCAACAATGTACTGAATCATTGGATTGTCGACGATTGGTAGCATCTCTTTGGGTATCGc        <  1:8075/69‑1 (MQ=255)
                                                     aaTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCtt     <  1:3048215/69‑1 (MQ=255)
                                                     aaTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCtt     <  1:3052536/69‑1 (MQ=255)
                                                     aaTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCtt     <  1:855889/69‑1 (MQ=255)
                                                     aaTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCtt     <  1:417349/69‑1 (MQ=255)
                                                     aaTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCtt     <  1:2637155/69‑1 (MQ=255)
                                                     aaTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCtt     <  1:2858214/69‑1 (MQ=255)
                                                     aaTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCtt     <  1:1268559/69‑1 (MQ=255)
                                                     aaTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCtt     <  1:162715/69‑1 (MQ=255)
                                                     aaTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCtt     <  1:181540/69‑1 (MQ=255)
                                                     aaTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCtt     <  1:1234184/69‑1 (MQ=255)
                                                     aaTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCtt     <  1:118408/69‑1 (MQ=255)
                                                     aaTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCtt     <  1:2132728/69‑1 (MQ=255)
                                                      aTCTCGTCAACGATGTACTGAATCATTGGCTTGTCGACGATCGGTAGCATCTCTTTGGGAATCGCCTTa    <  1:796206/69‑1 (MQ=255)
                                                      aTCTCGTCAACGATGTACTGAATCATTGGCTTGTCGACGATCGGTAGCATCTCTTTGGGAATCGCCTTa    <  1:1265464/69‑1 (MQ=255)
                                                      aTCTCGTCAACGATGTACTGAATCATTGGCTTGTCGACGATCGGTAGCATCTCTTTGGGAATCGCCTTa    <  1:38032/69‑1 (MQ=255)
                                                      aTCTCGTCAACGATGTACTGAATCATTGGCTTGTCGACGATCGGTAGCATCTCTTTGGGAATCGCCTTa    <  1:2178689/69‑1 (MQ=255)
                                                      aTCTCGTCAACGATGTACTGAATCATTGGCTTGTCGACGATCGGTAGCATCTCTTTGGGAATCGCCTTa    <  1:2396689/69‑1 (MQ=255)
                                                      aTCTCGTCAACGATGTACTGAATCATTGGCTTGTCGACGATCGGTAGCATCTCTTTGGGAATCGCCTTa    <  1:2771730/69‑1 (MQ=255)
                                                      aTCTCGTCAACGATGTACTGAATCATTGGCTTGTCGACGATCGGTAGCATCTCTTTGGGAATCGCCTTa    <  1:247057/69‑1 (MQ=255)
                                                      aTCTCGTCAACGATGTACTGAATCATTGGCTTGTCGACGATCGGTAGCATCTCTTTGGGAATCGCCTTa    <  1:2416301/69‑1 (MQ=255)
                                                      aTCTCGTCAACGATGTACTGAATCATTGGCTTGTCGACGATCGGTAGCATCTCTTTGGGAATCGCCTTa    <  1:3194939/69‑1 (MQ=255)
                                                         tcGTCAACAATGTACTTAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTggg              <  1:1305189/56‑1 (MQ=255)
                                                         tcGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTTGTAGCATCTCTTTggg              <  1:3276158/56‑1 (MQ=255)
                                                         tcGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTggg              <  1:1456931/56‑1 (MQ=255)
                                                         tcGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTggg              <  1:2340839/56‑1 (MQ=255)
                                                         tcGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTggg              <  1:423569/56‑1 (MQ=255)
                                                         tcGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTggg              <  1:18846/56‑1 (MQ=255)
                                                         tcGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTggg              <  1:100249/56‑1 (MQ=255)
                                                         tcGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTggg              <  1:1683604/56‑1 (MQ=255)
                                                         tcGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTggg              <  1:3032487/56‑1 (MQ=255)
                                                         tcGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGt  >  1:3341716/1‑68 (MQ=255)
                                                         tcGTCAACAATGTACTGAATCATTGGCTTGTCGACGAATGGTAGCATCTCTTTggg              <  1:1759437/56‑1 (MQ=255)
                                                            tCAACGATGTACTGAATCATTGGGTTGTCGACGATCgg                             >  1:3250244/1‑38 (MQ=25)
                                                            tCAACGATGTACTGAATCATTGGCTTGTCGACGATCgg                             >  1:841078/1‑38 (MQ=38)
                                                            tCAACGATGTACTGAATCATTGGCTTGTCGACGATCgg                             >  1:1079180/1‑38 (MQ=38)
                                                            tCAACGATGTACTGAATCATTGGCTTGTCGACGATCgg                             >  1:887270/1‑38 (MQ=38)
                                                            tCAACGATGTACTGAATCATTGGCTTGTCGACGATCgg                             >  1:108491/1‑38 (MQ=38)
                                                            tCAACGATGTACTGAATCATTGGCTTGTCGACGATCgg                             >  1:2970916/1‑38 (MQ=38)
                                                            tCAACGATGTACTGAATCATTGGCTTGTCGACGATCgg                             >  1:103793/1‑38 (MQ=38)
                                                            tCAACGATGTACTGAATCATTGGCTTGTCGACGATCgg                             >  1:739077/1‑38 (MQ=38)
                                                            tCAACGATGTACTGAATCATTGGCTTGTCGACGATCgg                             >  1:1634350/1‑38 (MQ=38)
                                                            tCAACGATGTACTGAATCATTGGCTTGTCGACGATCgg                             >  1:1901731/1‑38 (MQ=38)
                                                            tCAACGATGTACTGAATCATTGGCTTGTCGACGATCgg                             >  1:1905889/1‑38 (MQ=38)
                                                            tCAACGATGTACTGAATCATTGGCTTGTCGACGATCgg                             >  1:2216766/1‑38 (MQ=38)
                                                            tCAACGATGTACTGAATCATTGGCTTGTCGACGATCgg                             >  1:3231518/1‑38 (MQ=38)
                                                            tCAACGATGTACTGAATCATTGGCTTGTCGACGATCgg                             >  1:3181401/1‑38 (MQ=38)
                                                            tCAACGATGTACTGAATCATTGGCTTGTCGACGATCgg                             >  1:2854158/1‑38 (MQ=38)
                                                                  |                                                           
GCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGT  >  minE/1322749‑1322873

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: