Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,302,032 C→T 53.8% Q13Q (CAG→CAA metL ← fused aspartokinase II and homoserine dehydrogenase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,302,0320CT53.8% 10.2 / 89.0 104Q13Q (CAG→CAAmetLfused aspartokinase II and homoserine dehydrogenase II
Reads supporting (aligned to +/- strand):  ref base C (23/25);  new base T (28/28);  total (51/53)
Fisher's exact test for biased strand distribution p-value = 8.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 109 total aligned reads displayed.

ATGCCCGCGACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGTTTGCAGCCCGGAAGCCA  >  minE/2301970‑2302100
                                                              |                                                                    
aTGCCCGCGACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTgacga                                                                >  1:830753/1‑69 (MQ=255)
aTGCCCGCGACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTgacga                                                                >  1:78836/1‑69 (MQ=255)
  gCCCGCGACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTgacgac                                                               <  1:2826057/68‑1 (MQ=255)
  gCCCGCGACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTgacgac                                                               <  1:2116342/68‑1 (MQ=255)
  gCCCGCGACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTgacgac                                                               <  1:2122621/68‑1 (MQ=255)
  gCCCGCGACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTgacgac                                                               <  1:538909/68‑1 (MQ=255)
  gCCCGCGACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTgacgac                                                               <  1:1688715/68‑1 (MQ=255)
  gCCCGCGACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTgacgac                                                               <  1:560447/68‑1 (MQ=255)
  gCCCGCGACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTgacgac                                                               <  1:1246425/68‑1 (MQ=255)
  gCCCGCGACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTgacgac                                                               <  1:3251414/68‑1 (MQ=255)
  gCCCGCGACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTgacgac                                                               <  1:1772975/68‑1 (MQ=255)
  gCCCGCGACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTgacgac                                                               <  1:1194654/68‑1 (MQ=255)
  gCCCGCGACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTgacgac                                                               <  1:836276/68‑1 (MQ=255)
  gCCCGCGACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTgacgac                                                               <  1:2750428/68‑1 (MQ=255)
  gCCCGCGACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTgacgac                                                               <  1:1057031/68‑1 (MQ=255)
       cgACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTgacgac                                                               <  1:2500091/63‑1 (MQ=255)
        gACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTGACGACCTTTCGc                                                        <  1:917807/69‑1 (MQ=255)
        gACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTGACGACCTTTCGc                                                        <  1:1682597/69‑1 (MQ=255)
        gACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTGACGACCTTTCGc                                                        <  1:1209047/69‑1 (MQ=255)
        gACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTGACGACCTTTCGc                                                        <  1:854667/69‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:1010394/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:2662673/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:2546229/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:2465791/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:263020/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:2398157/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:2923441/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:2932117/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:3038883/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:309965/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:3159320/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:3316324/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:3326593/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:389449/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:441149/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:464198/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:499931/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:545389/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:787880/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:1903679/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:156347/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:1364775/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:1259487/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:1171445/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:1781030/56‑1 (MQ=255)
             gCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:2042386/56‑1 (MQ=255)
              cAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:2304630/55‑1 (MQ=255)
                      acTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTGACGAcc                                                              >  1:1124117/1‑49 (MQ=255)
                      acTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTGACGAcc                                                              >  1:623753/1‑49 (MQ=255)
                      acTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTGACGAcc                                                              >  1:607880/1‑49 (MQ=255)
                      acTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTGACGAcc                                                              >  1:2227844/1‑49 (MQ=255)
                      acTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTGACGAcc                                                              >  1:3298203/1‑49 (MQ=255)
                      acTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTGACGAcc                                                              >  1:2243511/1‑49 (MQ=255)
                               cAGCCAGACTACTGCCACCAAATTTATGCAGCTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCa                                 <  1:2471729/69‑1 (MQ=255)
                                aGCCAGACTACTGCCACCAAATTTATGCAGCTGACGACCTTTCGCCcct                                                    >  1:1193380/1‑49 (MQ=255)
                                aGCCAGACTACTGCCACCAAATTTATGCAGCTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCa                                 >  1:1621388/1‑68 (MQ=255)
                                aGCCAGACTACTGCCACCAAATTTATGCAGCTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCAt                                <  1:2892065/69‑1 (MQ=255)
                                aGCCAGACTACTGCCACCAAATTTATGCAGCTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCAt                                >  1:2926467/1‑69 (MQ=255)
                                aGCCAGACTACTGCCACCAAATTTATGCAGCTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCAt                                <  1:2028112/69‑1 (MQ=255)
                                aGCCAGACTACTGCCACCAAATTTATGCAGCTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCAt                                <  1:2143334/69‑1 (MQ=255)
                                aGCCAGACTACTGCCACCAAATTTATGCAGCTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCAt                                <  1:107564/69‑1 (MQ=255)
                                 gCCAGACTACTGCCACCAAATTTATGCAGTTgacga                                                                <  1:1470695/36‑1 (MQ=255)
                                       ctactGCCACCAAATTTATGCAGCTGACGACCTTTCGCCCCTGCCTGCGCAATCacac                                    >  1:2693060/1‑58 (MQ=255)
                                       ctactGCCACCAAATTTATGCAGCTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTAcc                         >  1:1991537/1‑69 (MQ=255)
                                       ctactGCCACCAAATTTATGCAGCTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTAcc                         >  1:2893133/1‑69 (MQ=255)
                                       ctactGCCACCAAATTTATGCAGCTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTAcc                         >  1:702137/1‑69 (MQ=255)
                                       ctactGCCACCAAATTTATGCAGCTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTAcc                         >  1:1899973/1‑69 (MQ=255)
                                       ctactGCCACCAAATTTATGCAGCTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTAcc                         >  1:2540280/1‑69 (MQ=255)
                                       ctactGCCACCAAATTTATGCAGCTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTAcc                         >  1:1451225/1‑69 (MQ=255)
                                             ccaccaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGt                   >  1:1715056/1‑69 (MQ=255)
                                             ccaccaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGt                   >  1:625479/1‑69 (MQ=255)
                                             ccaccaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGt                   >  1:897566/1‑69 (MQ=255)
                                             ccaccaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGt                   >  1:59915/1‑69 (MQ=255)
                                             ccaccaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGt                   >  1:1784663/1‑69 (MQ=255)
                                             ccaccaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGt                   >  1:3272346/1‑69 (MQ=255)
                                             ccaccaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCCTGt                   >  1:762773/1‑69 (MQ=255)
                                               accaAATTTATGCGGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGtt                  >  1:2239565/1‑68 (MQ=255)
                                               accaAATTTATGCAGTTGACGTCCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGtt                  >  1:2811132/1‑68 (MQ=255)
                                               accaAATTTATGCAGTTGACGCCCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGtt                  >  1:578597/1‑68 (MQ=255)
                                               accaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTAccc                        >  1:2289382/1‑62 (MQ=255)
                                               accaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGtt                  >  1:615947/1‑68 (MQ=255)
                                               accaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGtt                  >  1:1049964/1‑68 (MQ=255)
                                               accaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGtt                  >  1:2401573/1‑68 (MQ=255)
                                               accaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGtt                  >  1:2031477/1‑68 (MQ=255)
                                               accaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGtt                  >  1:805686/1‑68 (MQ=255)
                                               accaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGtt                  >  1:2705202/1‑68 (MQ=255)
                                               accaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGtt                  >  1:1221562/1‑68 (MQ=255)
                                               accaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGtt                  >  1:301381/1‑68 (MQ=255)
                                               accaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGtt                  >  1:3066551/1‑68 (MQ=255)
                                               accaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGtt                  >  1:1780429/1‑68 (MQ=255)
                                               accaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGtt                  >  1:513720/1‑68 (MQ=255)
                                               accaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGtt                  >  1:1399635/1‑68 (MQ=255)
                                               accaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGtt                  >  1:1401641/1‑68 (MQ=255)
                                               accaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGtt                  >  1:3221856/1‑68 (MQ=255)
                                               accaAATTTATGCAGTTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGtt                  >  1:361228/1‑68 (MQ=255)
                                                              cTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGTTTGCAGCCCGGAAGCCa  <  1:1947773/69‑1 (MQ=255)
                                                              cTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGTTTGCAGCCCGGAAGCCa  <  1:2712287/69‑1 (MQ=255)
                                                              cTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGTTTGCAGCCCGGAAGCCa  <  1:1381360/69‑1 (MQ=255)
                                                              cTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGTTTGCAGCCCGGAAGCCa  <  1:821423/69‑1 (MQ=255)
                                                              cTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGTTTGCAGCCCGGAAGCCa  <  1:2124144/69‑1 (MQ=255)
                                                              |                                                                    
ATGCCCGCGACACGCAAATAACACTTCACATCAGCCAGACTACTGCCACCAAATTTATGCAGCTGACGACCTTTCGCCCCTGCCTGCGCAATCACACTCATTTTTACCCCTTGTTTGCAGCCCGGAAGCCA  >  minE/2301970‑2302100

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: