Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,889,046 G→A P395S (CCC→TCC)  gudX ← predicted glucarate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,889,0460GA85.7% 233.4 / 28.1 91P395S (CCC→TCC) gudXpredicted glucarate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (13/0);  new base A (27/51);  total (40/51)
Fisher's exact test for biased strand distribution p-value = 6.28e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

TTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCC  >  minE/1888987‑1889113
                                                           |                                                                   
tttATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGACGCATCAGGaa                                                          >  1:1913389/1‑71 (MQ=255)
tttATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGACGCATCAGGaa                                                          >  1:652130/1‑71 (MQ=255)
  tataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCAt                                                                >  1:3803117/1‑63 (MQ=255)
  tataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCAt                                                                >  1:4158740/1‑63 (MQ=255)
  tataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCAt                                                                >  1:3341124/1‑63 (MQ=255)
  tataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCAt                                                                >  1:3283461/1‑63 (MQ=255)
  tataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCAt                                                                >  1:4027301/1‑63 (MQ=255)
  tataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCAt                                                                >  1:3839268/1‑63 (MQ=255)
  tataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCAt                                                                >  1:1031946/1‑63 (MQ=255)
  tataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCAt                                                                >  1:1344005/1‑63 (MQ=255)
  tataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCAt                                                                >  1:3612081/1‑63 (MQ=255)
  tataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCAt                                                                >  1:1266762/1‑63 (MQ=255)
  tataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCAt                                                                >  1:760144/1‑63 (MQ=255)
  tataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCAt                                                                >  1:3607541/1‑63 (MQ=255)
  tataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCAt                                                                >  1:868296/1‑63 (MQ=255)
  tataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGACGCATCAGGAACg                                                        >  1:2689616/1‑71 (MQ=255)
  tataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGACGCATCAGGAACg                                                        >  1:2098791/1‑71 (MQ=255)
  tataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGACGCATCAGGAACg                                                        >  1:483239/1‑71 (MQ=255)
  tataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGACGCATCAGGAACg                                                        >  1:1592141/1‑71 (MQ=255)
  tataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGACGCATCAGGAACg                                                        >  1:533497/1‑71 (MQ=255)
  tataGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGCACGCAACAgg                                                            >  1:1417736/1‑67 (MQ=255)
            aTGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGACGCATCAGGACCGGCAAttttt                                               >  1:1293256/1‑70 (MQ=255)
            aTGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTc                                              >  1:3912591/1‑71 (MQ=255)
            aTGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTc                                              >  1:619771/1‑71 (MQ=255)
               tgCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCg                                            <  1:412012/70‑1 (MQ=255)
               tgCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCg                                            <  1:356644/70‑1 (MQ=255)
               tgCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCg                                            <  1:2820897/70‑1 (MQ=255)
               tgCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCg                                            <  1:3202700/70‑1 (MQ=255)
               tgCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCg                                            <  1:1087348/70‑1 (MQ=255)
                     ttGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTaa                                     >  1:3954479/1‑71 (MQ=255)
                                  cagtccagtTCCATGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtt                        >  1:3951951/1‑71 (MQ=255)
                                  cagtccagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCg                            >  1:4362958/1‑67 (MQ=255)
                                  cagtccagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtt                        >  1:2818919/1‑71 (MQ=255)
                                  cagtccagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtt                        <  1:2899785/71‑1 (MQ=255)
                                  cagtccagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtt                        >  1:3566956/1‑71 (MQ=255)
                                  cagtccagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtt                        >  1:4746510/1‑71 (MQ=255)
                                  cagtccagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtt                        >  1:32672/1‑71 (MQ=255)
                                  cagtccagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtt                        >  1:1139762/1‑71 (MQ=255)
                                  cagtccagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtt                        >  1:845749/1‑71 (MQ=255)
                                  cagtccagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGGGAttt                       >  1:3225839/1‑71 (MQ=255)
                                     tccagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTg                     <  1:3190011/71‑1 (MQ=255)
                                     tccagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTg                     <  1:805980/71‑1 (MQ=255)
                                     tccagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTg                     <  1:2171256/71‑1 (MQ=255)
                                     tccagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTg                     <  1:919376/71‑1 (MQ=255)
                                     tccagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTg                     <  1:2946236/71‑1 (MQ=255)
                                     tccagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTg                     <  1:1640246/71‑1 (MQ=255)
                                     tccagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTg                     <  1:1653265/71‑1 (MQ=255)
                                      ccagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTg                     <  1:3773811/70‑1 (MQ=255)
                                       cagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCCGCGGATTTTggg                   >  1:4256665/1‑71 (MQ=255)
                                       cagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtttt                      <  1:4716013/68‑1 (MQ=255)
                                       cagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtttt                      <  1:4437802/68‑1 (MQ=255)
                                       cagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtttt                      <  1:4431958/68‑1 (MQ=255)
                                       cagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtttt                      <  1:954317/68‑1 (MQ=255)
                                       cagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtttt                      <  1:1790154/68‑1 (MQ=255)
                                       cagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtttt                      <  1:1319566/68‑1 (MQ=255)
                                       cagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtttt                      <  1:1202272/68‑1 (MQ=255)
                                       cagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtttt                      <  1:3151425/68‑1 (MQ=255)
                                       cagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtttt                      <  1:3274238/68‑1 (MQ=255)
                                       cagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtttt                      <  1:89441/68‑1 (MQ=255)
                                       cagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtttt                      <  1:2974499/68‑1 (MQ=255)
                                       cagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtttt                      <  1:806512/68‑1 (MQ=255)
                                       cagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTggg                   >  1:2024135/1‑71 (MQ=255)
                                       cagtTCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTggg                   >  1:3256489/1‑71 (MQ=255)
                                          ttCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGAtttt                      <  1:2275652/65‑1 (MQ=255)
                                          ttCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCa                <  1:3046498/71‑1 (MQ=255)
                                          ttCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCa                <  1:4651277/71‑1 (MQ=255)
                                          ttCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCa                <  1:4594402/71‑1 (MQ=255)
                                          ttCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCa                <  1:524919/71‑1 (MQ=255)
                                          ttCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCa                <  1:1778374/71‑1 (MQ=255)
                                          ttCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCa                <  1:247298/71‑1 (MQ=255)
                                          ttCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCa                <  1:2409934/71‑1 (MQ=255)
                                          ttCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCa                <  1:686953/71‑1 (MQ=255)
                                          ttCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCa                <  1:3754457/71‑1 (MQ=255)
                                          ttCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCa                <  1:2356485/71‑1 (MQ=255)
                                           tCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTg                     <  1:3475880/65‑1 (MQ=255)
                                           tCCACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCa                <  1:4404684/70‑1 (MQ=255)
                                            ccACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAgg              >  1:3053284/1‑71 (MQ=255)
                                            ccACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAgg              >  1:1480595/1‑71 (MQ=255)
                                            ccACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAgg              >  1:3066442/1‑71 (MQ=255)
                                             cACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGc             <  1:4637685/71‑1 (MQ=255)
                                             cACGCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGc             <  1:4037983/71‑1 (MQ=255)
                                                gCCCAGACCGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTg                     <  1:3599683/60‑1 (MQ=255)
                                                       ccGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCAGCGGATTTTGGGTCAGGCGACAATCGccc  <  1:1832078/71‑1 (MQ=255)
                                                        cGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGccc  <  1:969392/71‑1 (MQ=255)
                                                        cGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGccc  <  1:3271558/71‑1 (MQ=255)
                                                        cGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGccc  <  1:3155874/71‑1 (MQ=255)
                                                        cGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGccc  <  1:2610174/71‑1 (MQ=255)
                                                        cGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGccc  <  1:2400201/71‑1 (MQ=255)
                                                        cGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGccc  <  1:1821701/71‑1 (MQ=255)
                                                        cGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGccc  <  1:1497951/71‑1 (MQ=255)
                                                        cGGACGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGccc  <  1:1461261/71‑1 (MQ=255)
                                                           |                                                                   
TTTATAGGCCTCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCC  >  minE/1888987‑1889113

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: