Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 516,365 T→A 42.8% S189S (TCT→TCA bioB → biotin synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE516,3650TA42.8% 21.5 / 54.9 56S189S (TCT→TCAbioBbiotin synthase
Reads supporting (aligned to +/- strand):  ref base T (22/10);  new base A (19/5);  total (41/15)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCA  >  minE/516299‑516428
                                                                  |                                                               
ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:487252/2‑69 (MQ=255)
ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:2171489/2‑69 (MQ=255)
ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:2394844/2‑69 (MQ=255)
ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:2480938/2‑69 (MQ=255)
ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:2141110/2‑69 (MQ=255)
ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:2611706/2‑69 (MQ=255)
ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:2108761/2‑69 (MQ=255)
ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:2690297/2‑69 (MQ=255)
ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:2740157/2‑69 (MQ=255)
ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:2911204/2‑69 (MQ=255)
ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:1871857/2‑69 (MQ=255)
ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:3267887/2‑69 (MQ=255)
ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:1580494/2‑69 (MQ=255)
ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:729876/2‑69 (MQ=255)
ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:1428991/2‑69 (MQ=255)
ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:83141/2‑69 (MQ=255)
ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:133067/2‑69 (MQ=255)
ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:1153245/2‑69 (MQ=255)
ccGCACCTATCAGGAACGTCTTGATACGCTGGAAAAAGTGCGAGAAGCTGGGATTAAAGTCTGTTCAgg                                                               >  1:835675/2‑69 (MQ=255)
       tATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCAtt                                                        >  1:497114/1‑69 (MQ=255)
       tATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCAtt                                                        >  1:1902532/1‑69 (MQ=255)
       tATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCAtt                                                        >  1:525333/1‑69 (MQ=255)
       tATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCAtt                                                        >  1:115452/1‑69 (MQ=255)
       tATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCAtt                                                        >  1:894929/1‑69 (MQ=255)
        aTCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCAtt                                                        >  1:1116876/1‑68 (MQ=255)
               aCGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCtt                                                >  1:943425/1‑69 (MQ=255)
               aCGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCtt                                                >  1:3424375/1‑69 (MQ=255)
               aCGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCtt                                                >  1:1974029/1‑69 (MQ=255)
                                 aaaaGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACTGt                                    >  1:952649/1‑63 (MQ=255)
                                 aaaaGTGCGAGAAGCTGGGATTAAAGTCTGTTCAGGTGGCATTGTGGGCTTAGGCGAAACGGTGAAAg                               <  1:3379258/68‑1 (MQ=255)
                                 aaaaGTGCGAGAAGCTGGGATTAAAGTCTGTTCAGGTGGCATTGTGGGCTTAGGCGAAACGGTGAAAg                               <  1:105884/68‑1 (MQ=255)
                                 aaaaGTGCGAGAAGCTGGGATTAAAGTCTGTTCAGGTGGCATTGTGGGCTTAGGCGAAACGGTGAAAg                               <  1:3346759/68‑1 (MQ=255)
                                 aaaaGTGCGAGAAGCTGGGATTAAAGTCTGTTCAGGTGGCATTGTGGGCTTAGGCGAAACGGTGAAAg                               <  1:1974904/68‑1 (MQ=255)
                                 aaaaGTGCGAGAAGCTGGGATTAAAGTCTGTTCAGGTGGCATTGTGGGCTTAGGCGAAACGGTGAAAg                               <  1:2160617/68‑1 (MQ=255)
                                              gCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGatta                 >  1:2959342/1‑69 (MQ=255)
                                              gCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGatta                 >  1:1613775/1‑69 (MQ=255)
                                              gCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGatta                 >  1:711984/1‑69 (MQ=255)
                                              gCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGatta                 >  1:7056/1‑69 (MQ=255)
                                                     tCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATtgctgc          >  1:11453/1‑69 (MQ=255)
                                                             tGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCa  >  1:2109822/1‑69 (MQ=255)
                                                             tGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCa  >  1:3181885/1‑69 (MQ=255)
                                                             tGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCa  >  1:1707278/1‑69 (MQ=255)
                                                             tGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCa  >  1:268721/1‑69 (MQ=255)
                                                             tGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCa  >  1:2285675/1‑69 (MQ=255)
                                                             tGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCa  >  1:2645870/1‑69 (MQ=255)
                                                             tGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGAGCCGGATTATTGCTGCAACTGGCa  >  1:1975528/1‑69 (MQ=255)
                                                                 cTGGCGGCATTGTGGGCTTAGGCTAAACGGTAAAAGATCGCGCCGGattat                <  1:429035/51‑1 (MQ=255)
                                                                 cTGGCGGCATTGTGGGCTTAGGCTAAACGGTAAAAGATCGCGCCGGattat                <  1:2723332/51‑1 (MQ=255)
                                                                 cTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGattat                <  1:607195/51‑1 (MQ=255)
                                                                 cTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGattat                <  1:1495902/51‑1 (MQ=255)
                                                                 cTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGattat                <  1:1527601/51‑1 (MQ=255)
                                                                 cTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGattat                <  1:1488136/51‑1 (MQ=255)
                                                                 cTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGattat                <  1:811334/51‑1 (MQ=255)
                                                                 cTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGattat                <  1:1382334/51‑1 (MQ=255)
                                                                 cTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGattat                <  1:3005915/51‑1 (MQ=255)
                                                                 cTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGattat                <  1:2153026/51‑1 (MQ=255)
                                                                  |                                                               
ACGCACTTATCAGGAACGCCTCGATACGCTGGAAAAAGTGCGCGATGCCGGGATCAAAGTCTGTTCTGGCGGCATTGTGGGCTTAGGCGAAACGGTAAAAGATCGCGCCGGATTATTGCTGCAACTGGCA  >  minE/516299‑516428

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: