Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 615,152 A→G 52.8% K267K (AAA→AAG serC → 3‑phosphoserine/phosphohydroxythreonine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE615,1520AG52.8% 9.2 / 77.8 91K267K (AAA→AAGserC3‑phosphoserine/phosphohydroxythreonine aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (9/34);  new base G (21/27);  total (30/61)
Fisher's exact test for biased strand distribution p-value = 2.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

TATCTGGTCTGGTCTTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCTATATGGGGTGATTGATAACAGCGATTTCTACCGCAATGA  >  minE/615094‑615217
                                                          |                                                                 
tatctGGTCTGGTCTTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAatcaatca                                                           >  1:2363037/1‑67 (MQ=255)
tatctGGTCTGGTCTTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAg                                                          >  1:1719596/1‑68 (MQ=255)
tatctGGTCTGGTCTTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAg                                                          >  1:2506641/1‑68 (MQ=255)
         tggtctTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAg                                                <  1:613084/69‑1 (MQ=255)
         tggtctTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAg                                                <  1:710005/69‑1 (MQ=255)
         tggtctTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAg                                                <  1:2533610/69‑1 (MQ=255)
         tggtctTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAg                                                <  1:3359934/69‑1 (MQ=255)
         tggtctTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAg                                                <  1:1382052/69‑1 (MQ=255)
         tggtctTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAg                                                <  1:2713493/69‑1 (MQ=255)
         tggtctTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAg                                                <  1:2727206/69‑1 (MQ=255)
         tggtctTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAg                                                <  1:1944513/69‑1 (MQ=255)
         tggtctTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAg                                                <  1:1910901/69‑1 (MQ=255)
         tggtctTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAg                                                <  1:2750113/69‑1 (MQ=255)
                  aaTGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCta                                        <  1:2585734/68‑1 (MQ=255)
                  aaTGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCta                                        <  1:2632714/68‑1 (MQ=255)
                  aaTGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCta                                        <  1:2027381/68‑1 (MQ=255)
                  aaTGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCta                                        <  1:2023596/68‑1 (MQ=255)
                  aaTGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCta                                        <  1:1692041/68‑1 (MQ=255)
                  aaTGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCta                                        <  1:1529405/68‑1 (MQ=255)
                  aaTGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCta                                        <  1:3202281/68‑1 (MQ=255)
                  aaTGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCta                                        <  1:3207136/68‑1 (MQ=255)
                  aaTGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCta                                        <  1:3279989/68‑1 (MQ=255)
                  aaTGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCta                                        <  1:3375483/68‑1 (MQ=255)
                  aaTGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCta                                        <  1:3394545/68‑1 (MQ=255)
                  aaTGGCTGAAAGAGAACGGCGGGGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCta                                        <  1:315573/68‑1 (MQ=255)
                      gCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCTATATgg                                   >  1:2911288/1‑69 (MQ=255)
                      gCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCTATATgg                                   >  1:3078506/1‑69 (MQ=255)
                      gCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCTATATgg                                   >  1:1994033/1‑69 (MQ=255)
                      gCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCTATATgg                                   >  1:3144157/1‑69 (MQ=255)
                      gCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCTATATgg                                   >  1:1067910/1‑69 (MQ=255)
                         gAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCGGAACTGCTATATGGGGt                                <  1:2702281/69‑1 (MQ=255)
                         gAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCTATATGGGGt                                <  1:2632136/69‑1 (MQ=255)
                         gAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCTATATGGGGt                                <  1:297780/69‑1 (MQ=255)
                         gAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCTATATGGGGt                                <  1:3224266/69‑1 (MQ=255)
                         gAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCTATATGGGGt                                <  1:2296259/69‑1 (MQ=255)
                         gAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCTATATGGGGt                                <  1:1162490/69‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:1136373/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:1757322/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:2947172/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:1592440/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:3019674/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:1545554/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:1501116/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:319995/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:1350975/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:1169954/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:3130137/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:379234/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:453140/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:460258/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:547955/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:1104059/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:622849/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:873791/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:919356/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:2833632/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:2267373/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:2185599/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:2647785/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:2019909/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:2460681/58‑1 (MQ=255)
                           aaGCGAACGGCGGTGTGGCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:1996429/58‑1 (MQ=255)
                                         gtAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCTATATGGGGTGATTGATAACAGCGAt                <  1:2706987/69‑1 (MQ=255)
                                         gtAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCTATATGGGGTGATTGATAACAGCGAt                >  1:1329035/1‑69 (MQ=255)
                                            gCTGCAATGGATAAGATCAATCGGCAAAAAGCAGAACTGCt                                         <  1:2983883/41‑1 (MQ=37)
                                            gCTGCAATGGATAAGATCAATCAGCAAAAAGCAGAACTGCtgt                                       >  1:780744/1‑41 (MQ=37)
                                            gCTGCAATGGATAAGATCAATCAGCAAAAAGCAGAACTGCtgt                                       >  1:2489990/1‑41 (MQ=37)
                                            gCTGCAATGGATAAGATCAATCAGCAAAAAGCAGAACTGCtgt                                       >  1:2514159/1‑41 (MQ=37)
                                            gCTGCAATGGATAAGATCAATCAGCAAAAAGCAGAACTGCtgt                                       >  1:2269432/1‑41 (MQ=37)
                                            gCTGCAATGGATAAGATCAATCAGCAAAAAGCAGAACTGCtgt                                       >  1:429285/1‑41 (MQ=37)
                                            gCTGCAATGGATAAGATCAATCAGCAAAAAGCAGAACTGCtgt                                       >  1:3422945/1‑41 (MQ=37)
                                            gCTGCAATGGATAAGATCAATCAGCAAAAAGCAGAACTGCtgt                                       >  1:1833439/1‑41 (MQ=37)
                                            gCTGCAATGGATAAGATCAATCAGCAAAAAGCAGAACTGCtgt                                       >  1:2985294/1‑41 (MQ=37)
                                            gCTGCAATGGATAAGATCAATCAGCAAAAAGCAGAACTGCtgt                                       >  1:2644190/1‑41 (MQ=37)
                                            gCTGCAATGGATAAGATCAATCAGCAAAAAGCAGAACTGCtgt                                       >  1:3130839/1‑41 (MQ=37)
                                            gCTGCAATGGATAAGATCAATAAGCAAAAAGCAGAACTGCtgt                                       >  1:1406881/1‑41 (MQ=37)
                                                   tGGATAAGATCAATCAGCAAAAAGCAGAACTGCTGTATGGGGTGATTGATAACAGCGATTTCTACCGc       >  1:816640/1‑68 (MQ=255)
                                                   tGGATAAGATCAATCAGCAAAAAGCAGAACTGCTGTATGGGGTGATTGATAACAGCGATTTCTACCGCa      >  1:566582/1‑69 (MQ=255)
                                                   tGGATAAGATCAATCAGCAAAAAGCAGAACTGCTGTATGGGGTGATTGATAACAGCGATTTCTACCGCa      >  1:2446284/1‑69 (MQ=255)
                                                   tGGATAAGATCAATCAGCAAAAAGCAGAACTGCTGTATGGGGTGATTGATAACAGCGATTTCTACCGCa      >  1:907189/1‑69 (MQ=255)
                                                   tGGATAAGATCAATCAGCAAAAAGCAGAACTGCTGTATGGGGTGATTGATAACAGCGATTTCTACCGCa      >  1:1080140/1‑69 (MQ=255)
                                                   tGGATAAGATCAATCAGCAAAAAGCAGAACTGCTGTATGGGGTGATTGATAACAGCGATTTCTACCGCa      >  1:635627/1‑69 (MQ=255)
                                                   tGGATAAGATCAATCAGCAAAAAGCAGAACTGCTGTATGGGGTGATTGATAACAGCGATTTCTACCGCa      >  1:555938/1‑69 (MQ=255)
                                                   tGGATAAGATCAATCAGCAAAAAGCAGAACTGCTGTATGGGGTGATTGATAACAGCGATTTCTACCGCa      >  1:2563299/1‑69 (MQ=255)
                                                   tGGATAAGATCAATCAGCAAAAAGCAGAACTGCTGTATGGGGTGATTGATAACAGCGATTTCTACCGCa      >  1:1171988/1‑69 (MQ=255)
                                                   tGGATAAGATCAATCAGCAAAAAGCAGAACTGCTGTATGGGGTGATTGATAACAGCGATTTCTACCGCa      >  1:1166947/1‑69 (MQ=255)
                                                       tAAAATCAATCAGCAAAAAGCAGAACTGCTATATGGGGTGATTGATAACAGCGATTTCTACCGCAATGa  <  1:2501745/69‑1 (MQ=255)
                                                       tAAAATCAATCAGCAAAAAGCAGAACTGCTATATGGGGTGATTGATAACAGCGATTTCTACCGCAATGa  <  1:2404681/69‑1 (MQ=255)
                                                       tAAAATCAATCAGCAAAAAGCAGAACTGCTATATGGGGTGATTGATAACAGCGATTTCTACCGCAATGa  <  1:1009190/69‑1 (MQ=255)
                                                       tAAAATCAATCAGCAAAAAGCAGAACTGCTATATGGGGTGATTGATAACAGCGATTTCTACCGCAATGa  <  1:2000090/69‑1 (MQ=255)
                                                       tAAAATCAATCAGCAAAAAGCAGAACTGCTATATGGGGTGATTGATAACAGCGATTTCTACCGCAATGa  <  1:256291/69‑1 (MQ=255)
                                                          |                                                                 
TATCTGGTCTGGTCTTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAGCAAAAAGCAGAACTGCTATATGGGGTGATTGATAACAGCGATTTCTACCGCAATGA  >  minE/615094‑615217

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: