Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,315,996:1 +T 100% coding (594/738 nt) hisA → N‑(5'‑phospho‑L‑ribosyl‑formimino)‑5‑amino‑1‑ (5'‑phosphoribosyl)‑4‑imidazolecarboxamide isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,315,9961.T90.9% 32.7 / ‑2.1 11coding (594/738 nt)hisAN‑(5'‑phospho‑L‑ribosyl‑formimino)‑5‑amino‑1‑ (5'‑phosphoribosyl)‑4‑imidazolecarboxamide isomerase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (10/0);  total (11/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

AGGCTCTAACGTCTCTTTATATGAAGAAGTGTGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCGCGGCG  >  minE/1315942‑1316066
                                                       |                                                                      
aGGCTCTAACGTCTCTTTATATGAAGAAGTGTGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCcggcg                                                        >  1:401062/1‑71 (MQ=255)
      tAACGTCTCTTTATATGAAGAAGTGTGTGCCAGATATCCGCAGGTGGCATTTTCAGTCCTCCGGCGGTAt                                                    >  1:574210/1‑70 (MQ=255)
      tAACGTCTCTTTATATGAAGAAGTGTGCGCCAGATATCCGCAGGTGGCATTTTCAGTCCTCCGGCGGTAtt                                                   >  1:319481/1‑71 (MQ=255)
      tAACGTCTCTTTATATGAAGAAGTGTGCGCCAGATATCCGCAGGTGGCATTTTCAGTCCTCCGGCGGTAtt                                                   >  1:125975/1‑71 (MQ=255)
      tAACGTCTCTTTATATGAAGAAGTGTGCGCCAGATATCCGCAGGTGGCATTTTCAGTCCTCCGGCGGTAtt                                                   >  1:395163/1‑71 (MQ=255)
      tAACGTCTCTTTATATGAAGAAGTGTGCGCCAGATATCCGCAGGTGGCATTTTCAGTCCTCCGGCGGTAtt                                                   >  1:390433/1‑71 (MQ=255)
      tAACGTCTCTTTATATGAAGAAGTGTGCGCCAGATATCCGCAGGTGGCATTTTCAGTCCTCCGGCGGTAtt                                                   >  1:627017/1‑71 (MQ=255)
      tAACGTCTCTTTATATGAAGAAGTGTGCGCCAGATATCCGCAGGTGGCATTTTCAGTCCTCCGGCGGTAtt                                                   >  1:372090/1‑71 (MQ=255)
      tAACGTCTCTTTATATGAAGAAGTGTGCGCCAGATATCCGCAGGTGGCATTTTCAGTCCTCCGGCGGTAtt                                                   >  1:370744/1‑71 (MQ=255)
      tAACGTCTCTTTATATGAAGAAGTGTGCGCCAGATATCCGCAGGTGGCATTTTCAGTCCTCCGGCGGTAt                                                    >  1:372834/1‑70 (MQ=255)
      tAACGTCTCTTTATATGAAGAAGTGTGCGCCAGATATCCGCAGGAGGCATTTTCAGTCCTCCGGCGGTAtt                                                   >  1:253659/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGTACTGGTGTGCgcggcg  >  1:643664/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGgt            >  1:188727/1‑61 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:449407/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:451814/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:574252/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:462263/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:474971/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:478491/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:441760/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:505788/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:517025/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:527534/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:538439/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:540514/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:427963/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:542022/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:426045/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:644760/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:106612/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:91179/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:87494/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:85157/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:82864/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:82730/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:75622/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:68256/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:653875/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:646615/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:555216/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:641558/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:637835/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:634550/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:611825/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:59485/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:5834/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:579731/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:97178/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:568786/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:238449/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:298848/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:425175/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:295624/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:295569/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:286513/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:285455/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:283357/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:259664/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:301944/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:212601/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:212092/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:184587/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:183570/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:181157/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:178965/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:169304/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:159604/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:348068/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:117810/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:398230/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:386090/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:362439/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:361982/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:357717/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:142911/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:351613/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:350668/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:347370/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:346755/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:344438/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:335953/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:335617/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:330110/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:308933/1‑71 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggc   >  1:297875/1‑70 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggc   >  1:158255/1‑70 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggc   >  1:291876/1‑70 (MQ=255)
                                                      aTTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGAGGCGGCCCTGCGTGGCACTGGTGTGCgcggc   >  1:406675/1‑70 (MQ=255)
                                                      aTTTCAGTCATCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg  >  1:470238/1‑71 (MQ=255)
                                                       |                                                                      
AGGCTCTAACGTCTCTTTATATGAAGAAGTGTGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCGCGGCG  >  minE/1315942‑1316066

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: