Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 615,055 G→A 45.8% G235D (GGC→GAC) ‡ serC → 3‑phosphoserine/phosphohydroxythreonine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE615,0550GA45.8% 24.7 / 101.3 123G235D (GGC→GAC) ‡serC3‑phosphoserine/phosphohydroxythreonine aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (30/36);  new base A (28/28);  total (58/65)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

Only 100 of 123 total aligned reads displayed.

TGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTATCTGGTCTGGT  >  minE/614992‑615106
                                                               |                                                   
tGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCa                                                >  1:1826383/1‑69 (MQ=255)
tGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCa                                                >  1:1016786/1‑69 (MQ=255)
tGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCa                                                >  1:2320464/1‑69 (MQ=255)
tGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCa                                                >  1:2973246/1‑69 (MQ=255)
   tGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATg                                              <  1:2495653/68‑1 (MQ=255)
   tGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATg                                              <  1:2012771/68‑1 (MQ=255)
   tGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATg                                              <  1:1704575/68‑1 (MQ=255)
       cAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGttt                                           <  1:1436857/67‑1 (MQ=255)
       cAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGttt                                           <  1:2266509/67‑1 (MQ=255)
       cAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGttt                                           <  1:1662265/67‑1 (MQ=255)
       cAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGttt                                           <  1:3315192/67‑1 (MQ=255)
       cAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTaa                                         <  1:2531606/69‑1 (MQ=255)
       cAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTaa                                         <  1:3391691/69‑1 (MQ=255)
       cAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTaa                                         <  1:595993/69‑1 (MQ=255)
       cAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTaa                                         <  1:599385/69‑1 (MQ=255)
       cAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTaa                                         <  1:3224386/69‑1 (MQ=255)
       cAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTaa                                         <  1:77056/69‑1 (MQ=255)
       cAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTaa                                         <  1:1128620/69‑1 (MQ=255)
        aaaGCGAATATCGCCTGCCCATCGATTCTGGATTACTCCATCCTCAACGATAACGATTCCAt                                               <  1:2128909/62‑1 (MQ=255)
        aaaGCGAATATCGCCTGCCCATCGATTCTGGATTACTCCATCCTCAACGATAACGATTCCAt                                               <  1:2191608/62‑1 (MQ=255)
        aaaGCGAATATCGCCTGCCCATCGATTCTGGATTACTCCATCCTCAACGATAACGATTCCAt                                               <  1:2770392/62‑1 (MQ=255)
        aaaGCGAATATCGCCTGCCCATCGATTCTGGATTACTCCATCCTCAACGATAACGATTCCAt                                               <  1:2042774/62‑1 (MQ=255)
        aaaGCGAATATCGCCTGCCCATCGATTCTGGATTACTCCATCCTCAACGATAACGATTCCAt                                               <  1:2283802/62‑1 (MQ=255)
        aaaGCGAATATCGCCTGCCCATCGATTCTGGATTACTCCATCCTCAACGATAACGATTCCAt                                               <  1:272270/62‑1 (MQ=255)
        aaaGCGAATATCGCCTGCCCATCGATTCTGGATTACTCCATCCTCAACGATAACGATTCCAt                                               <  1:1929735/62‑1 (MQ=255)
        aaaGCGAATATCGCCTGCCCATCGATTCTGGATTACTCCATCCTCAACGATAACGATTCCAt                                               <  1:434736/62‑1 (MQ=255)
        aaaGCGAATATCGCCTGCCCATCGATTCTGGATTACTCCATCCTCAACGATAACGATTCCAt                                               <  1:621242/62‑1 (MQ=255)
        aaaGCGAATATCGCCTGCCCATCGATTCTGGATTACTCCATCCTCAACGATAACGATTCCAt                                               <  1:698796/62‑1 (MQ=255)
        aaaGCGAATATCGCCTGCCCATCGATTCTGGATTACTCCATCCTCAACGATAACGATTCCAt                                               <  1:157780/62‑1 (MQ=255)
        aaaGCGAATATCGCCTGCCCATCGATTCTGGATTACTCCATCCTCAACGATAACGATTCCAt                                               <  1:1738739/62‑1 (MQ=255)
        aaaGCGAATATCGCCTGCCCATCGATTCTGGATTACTCCATCCTCAACGATAACGATTCCAt                                               <  1:1804135/62‑1 (MQ=255)
        aaaGCGAATATCGCCTGCCCATCGATTCTGGATTACTCCATCCTCAACAATAACGATTCCAt                                               <  1:2996315/62‑1 (MQ=255)
         aaGCGAATATCGCCTGCCCATCGATTCTGGATTACTCCATCCTCAACGATAACGATTCCAt                                               <  1:548238/61‑1 (MQ=255)
          aGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAacac                                      <  1:863835/69‑1 (MQ=255)
          aGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAacac                                      <  1:272812/69‑1 (MQ=255)
              aaTATCGGCTGGCCATCGATTCTGGATTACTCCATCCTCAACGATAACGATTCCAt                                               <  1:2132655/56‑1 (MQ=255)
                  tCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCAGCCGa                              >  1:1963873/1‑69 (MQ=255)
                    gcgTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGCCa                            <  1:1461563/69‑1 (MQ=255)
                    gcgTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACa                            <  1:261932/69‑1 (MQ=255)
                    gcgTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACa                            <  1:1227620/69‑1 (MQ=255)
                    gcgTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACa                            <  1:975698/69‑1 (MQ=255)
                    gcgTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACa                            <  1:69972/69‑1 (MQ=255)
                    gcgTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACa                            <  1:145142/69‑1 (MQ=255)
                    gcgTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACa                            <  1:153080/69‑1 (MQ=255)
                    gcgTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACa                            <  1:3186298/69‑1 (MQ=255)
                    gcgTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACa                            <  1:1992178/69‑1 (MQ=255)
                          ccGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGcc                      >  1:1552367/1‑69 (MQ=255)
                          ccGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGcc                      >  1:3077426/1‑69 (MQ=255)
                          ccGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGcc                      >  1:755896/1‑69 (MQ=255)
                          ccGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGcc                      >  1:3381704/1‑69 (MQ=255)
                          ccGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGcc                      >  1:2443991/1‑69 (MQ=255)
                          ccGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGcc                      >  1:1802825/1‑69 (MQ=255)
                           cGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCt                     <  1:1057227/69‑1 (MQ=255)
                               gATTCTGGTTTACTCCATCCTCAACGATAACGATTCCAt                                               <  1:2472758/39‑1 (MQ=25)
                               gATTCTGGATTACTCCATCCTTAACGATAACGATTCCATGTTCAac                                        >  1:2997454/1‑46 (MQ=25)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCaa                                         >  1:761690/1‑45 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:1088381/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:1776450/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:597674/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:1749551/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:1890171/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:693954/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:1399831/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:1389405/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:2082691/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:316965/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:1269048/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:3022708/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:2159136/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:217016/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:2962818/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:953750/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:2926321/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:2812168/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:2202095/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:2351231/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:2412733/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:967627/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:2483330/1‑46 (MQ=255)
                               gATTCTGGATTACTCCATCCTCAACGATAACGATTCCATGTTCAac                                        >  1:2638824/1‑46 (MQ=255)
                                  tCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGtatct              <  1:2066864/69‑1 (MQ=255)
                                   cTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGAc                             >  1:1032027/1‑53 (MQ=255)
                                   cTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGAc                             >  1:757916/1‑53 (MQ=255)
                                   cTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGAc                             >  1:642930/1‑53 (MQ=255)
                                   cTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGAc                             >  1:495457/1‑53 (MQ=255)
                                   cTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGAc                             >  1:47696/1‑53 (MQ=255)
                                   cTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGAc                             >  1:1986862/1‑53 (MQ=255)
                                   cTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGAc                             >  1:3064436/1‑53 (MQ=255)
                                   cTGGATTATTCCATCCTCAACGATAACGGCTCAATGTTTAACACGCCGCCGAc                             >  1:3284511/1‑53 (MQ=255)
                                     ggATTATTCCATCCTCAACGATAACGTCTCCATGttt                                           <  1:2047288/37‑1 (MQ=255)
                                             ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg   >  1:1859209/1‑69 (MQ=255)
                                             ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg   >  1:1621876/1‑69 (MQ=255)
                                             ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg   >  1:1470328/1‑69 (MQ=255)
                                             ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg   >  1:3372302/1‑69 (MQ=255)
                                             ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg   >  1:2927475/1‑69 (MQ=255)
                                                tCCTCAACGATAACGATTCCATGTTCAACACGCCGCCGACATTTGCCTGGTACTTATctggcctggt  <  1:601711/67‑1 (MQ=255)
                                                tCCTCAACGATAACGATTCCATGTTCAACACGCCGCCGACATTTGCCTGGTACTTATctggcctggt  <  1:2080678/67‑1 (MQ=255)
                                                tCCTCAACGATAACGATTCCATGTTCAACACGCCGCCGACATTTGCCTGGTACTTATctggcctggt  <  1:272242/67‑1 (MQ=255)
                                                tCCTCAACGATAACGATTCCATGTTCAACACGCCGCCGACATTTGCCTGGTACTTATctggcctggt  <  1:2652213/67‑1 (MQ=255)
                                                 ccTCAACGATAACGATTCCATGTTCAACACGCCGCCGACATTTGCCTGGTACTTATctggcctggt  <  1:2123783/66‑1 (MQ=255)
                                                               |                                                   
TGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTATCTGGTCTGGT  >  minE/614992‑615106

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: