Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,718,421 C→T 14.1% I77I (ATC→ATT argD → bifunctional acetylornithine aminotransferase and succinyldiaminopimelate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,718,4210CT14.1% 138.4 / 9.2 77I77I (ATC→ATTargDbifunctional acetylornithine aminotransferase and succinyldiaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (36/30);  new base T (4/7);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.30e-01

GCCATCCTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTAAACTGATTGAGGCAACGTTTGCCG  >  minE/2718366‑2718488
                                                       |                                                                   
gCCATCCTGCGCTGGTGAACGCGTTAAAACCCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTc                                                         >  1:2948255/1‑68 (MQ=255)
gCCATCCTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTc                                                         >  1:1566569/1‑68 (MQ=255)
gCCATCCTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTc                                                         >  1:2374051/1‑68 (MQ=255)
gCCATCCTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTc                                                         >  1:2373033/1‑68 (MQ=255)
gCCATCCTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTc                                                         >  1:2825567/1‑68 (MQ=255)
gCCATCCTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTc                                                         >  1:2221468/1‑68 (MQ=255)
gCCATCCTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTc                                                         >  1:2960953/1‑68 (MQ=255)
gCCATCCTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTc                                                         >  1:1465338/1‑68 (MQ=255)
gCCATCCTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTc                                                         >  1:1325698/1‑68 (MQ=255)
gCCATCCTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTc                                                         >  1:1277408/1‑68 (MQ=255)
gCCATCCTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTc                                                         >  1:1070099/1‑68 (MQ=255)
gCCATCCTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTc                                                         >  1:1148129/1‑68 (MQ=255)
  cATCCTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCAcc                                                      >  1:1805389/1‑69 (MQ=255)
   aTCCTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCAcc                                                      <  1:1209825/68‑1 (MQ=255)
   aTCCTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCAcc                                                      <  1:917659/68‑1 (MQ=255)
   aTCCTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCAcc                                                      <  1:2271846/68‑1 (MQ=255)
   aTCCTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCAcc                                                      <  1:6601/68‑1 (MQ=255)
   aTCCTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCAcc                                                      <  1:3094786/68‑1 (MQ=255)
      cTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATg                                                  <  1:1062395/69‑1 (MQ=255)
      cTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATg                                                  <  1:8431/69‑1 (MQ=255)
      cTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATg                                                  <  1:1734448/69‑1 (MQ=255)
      cTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATg                                                  <  1:1244781/69‑1 (MQ=255)
              gTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGgcg                                          >  1:3025723/1‑69 (MQ=255)
              gTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGgcg                                          >  1:494688/1‑69 (MQ=255)
              gTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGgcg                                          >  1:2572626/1‑69 (MQ=255)
              gTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGgcg                                          >  1:1221854/1‑69 (MQ=255)
              gTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGATTTCTCCAATGAACCGgcg                                          >  1:3343060/1‑69 (MQ=255)
                   cgcgTTAAAAACCCAGGGCGAAACTCTGTGGCATATCATTAACGTTTTCACCAATGAACCGGCGCTGc                                      <  1:404475/68‑1 (MQ=255)
                   cgcgTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGGCGCTGc                                      <  1:350566/68‑1 (MQ=255)
                   cgcgTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGGCGCTGc                                      <  1:967821/68‑1 (MQ=255)
                   cgcgTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGGCGCTGc                                      <  1:667101/68‑1 (MQ=255)
                   cgcgTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGGCGCTGc                                      <  1:2926691/68‑1 (MQ=255)
                   cgcgTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGGCGCTGc                                      <  1:476762/68‑1 (MQ=255)
                   cgcgTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGGCGCTGc                                      <  1:2813294/68‑1 (MQ=255)
                   cgcgTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGGCGCTGc                                      <  1:167991/68‑1 (MQ=255)
                   cgcgTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGGCGCTGc                                      <  1:1904002/68‑1 (MQ=255)
                   cgcgTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGGCGCTGc                                      <  1:2146885/68‑1 (MQ=255)
                         aaaaaCCCAGGGCGAAACTCTGTGGCATATTAGTAACGTTTTCACCAATGagc                                               <  1:3281197/53‑3 (MQ=255)
                         aaaaaCCCAGGGCGAAACTCTGTGGCATATTAGTAACGTTTTCACCAATGagc                                               <  1:254462/53‑3 (MQ=255)
                         aaaaaCCCAGGGCGAAACTCTGTGGCATATTAGTAACGTTTTCACCAATGagc                                               <  1:1178497/53‑3 (MQ=255)
                         aaaaaCCCAGGGCGAAACTCTGTGGCATATTAGTAACGTTTTCACCAATGagc                                               <  1:736334/53‑3 (MQ=255)
                         aaaaaCCCAGGGCGAAACTCTGTGGCATATTAGTAACGTTTTCACCAATGagc                                               <  1:2082626/53‑3 (MQ=255)
                         aaaaaCCCAGGGCGAAACTCTGTGGCATATTAGTAACGTTTTCACCAATGagc                                               <  1:2803872/53‑3 (MQ=255)
                         aaaaaCCCAGGGCGAAACTCTGTGGCATATTAGTAACGTTTTCACCAATGagc                                               <  1:868080/53‑3 (MQ=255)
                              cccAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTa                          <  1:1046618/69‑1 (MQ=255)
                              cccAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTa                          <  1:2822668/69‑1 (MQ=255)
                              cccAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTa                          <  1:886264/69‑1 (MQ=255)
                              cccAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTa                          <  1:1075711/69‑1 (MQ=255)
                              cccAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTa                          <  1:2372508/69‑1 (MQ=255)
                                  gggCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTAAACt                      <  1:2002820/69‑1 (MQ=255)
                                    gCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATg                                                  <  1:811396/39‑1 (MQ=255)
                                                   atatTAGTAACGTTTTCACCAATGAGCCGGCGCTGCGTCTTGGCCGTAAACTGATTGAGGCAACGTTTg     >  1:2358861/1‑69 (MQ=255)
                                                   atatTAGTAACGTTTTCACCAATGAGCCGGCGCTGCGTCTTGGCCGTAAACTGATTGAGGCAACGTTTg     >  1:2639360/1‑69 (MQ=255)
                                                   atatTAGTAACGTTTTCACCAATGAGCCGGCGCTGCGTCTTGGCCGTAAACTGATTGAGGCAACGTTTg     >  1:2556404/1‑69 (MQ=255)
                                                   atatTAGTAACGTTTTCACCAATGAGCCGGCGCTGCGTCTTGGCCGTAAACTGATTGAGGCAACGTTTg     >  1:1205086/1‑69 (MQ=255)
                                                     atCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTa                          <  1:2505724/46‑1 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGCACCGGCGCTGCGTCTTGGGCGTAAACTGATTGAGGCAACGt        >  1:228378/1‑63 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCg                            >  1:3326950/1‑43 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTa                          <  1:1581969/45‑1 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTa                          <  1:2493435/45‑1 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTa                          <  1:3214019/45‑1 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTAAACTGATTGAGGCaa           >  1:1080894/1‑60 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTAAACTGATTGAGGCAACGtt       >  1:327944/1‑64 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTAAACTGATTGAGGCAACGtt       >  1:522381/1‑64 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTAAACTGATTGAGGCAACGtt       >  1:2624512/1‑64 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTAAACTGATTGAGGCAACGtt       >  1:2924757/1‑64 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTAAACTGATTGAGGCAACGtt       >  1:3279688/1‑64 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTAAACTGATTGAGGCAACGtt       >  1:2971753/1‑64 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTAAACTGATTGAGGCAACGtt       >  1:294912/1‑64 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTAAACTGATTGAGGCAACGtt       >  1:1076164/1‑64 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTAAACTGATTGAGGCAACGTTTGCCg  >  1:1158919/1‑69 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTAAACTGATTGAGGCAACGTTTGCCg  >  1:23776/1‑69 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTAAACTGATTGAGGCAACGTTTGCCg  >  1:2361191/1‑69 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTAAACTGATTGAGGCAACGTTTGCCg  >  1:578908/1‑69 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTAAACTGATTGAGGCAACGTTTGCCg  >  1:577009/1‑69 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTAAACTGATTGAGGCAACGTTTGCCg  >  1:541389/1‑69 (MQ=255)
                                                      tCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTAAACTGATTGAGGCAACGTTTGCCg  >  1:1350397/1‑69 (MQ=255)
                                                       |                                                                   
GCCATCCTGCGCTGGTGAACGCGTTAAAAACCCAGGGCGAAACTCTGTGGCATATCAGTAACGTTTTCACCAATGAACCGGCGCTGCGTCTTGGGCGTAAACTGATTGAGGCAACGTTTGCCG  >  minE/2718366‑2718488

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: