Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 615,101 T→C 100% G250G (GGT→GGC serC → 3‑phosphoserine/phosphohydroxythreonine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE615,1010TC76.4% 22.3 / 2.4 21G250G (GGT→GGCserC3‑phosphoserine/phosphohydroxythreonine aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (4/1);  new base C (0/16);  total (4/17)
Fisher's exact test for biased strand distribution p-value = 8.35e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTATCTGGTCTGGTCTTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAG  >  minE/615037‑615161
                                                                |                                                            
ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg                                                          >  1:2079359/1‑69 (MQ=255)
ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg                                                          >  1:2927475/1‑69 (MQ=255)
ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg                                                          >  1:3232092/1‑69 (MQ=255)
ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg                                                          >  1:3372302/1‑69 (MQ=255)
ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg                                                          >  1:1859209/1‑69 (MQ=255)
ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg                                                          >  1:1784204/1‑69 (MQ=255)
ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg                                                          >  1:1621876/1‑69 (MQ=255)
ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg                                                          >  1:1470328/1‑69 (MQ=255)
ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg                                                          >  1:581050/1‑69 (MQ=255)
   tCCTCAACGATAACGATTCCATGTTCAACACGCCGCCGACATTTGCCTGGTACTTATctggcctggt                                                         <  1:601711/67‑1 (MQ=255)
   tCCTCAACGATAACGATTCCATGTTCAACACGCCGCCGACATTTGCCTGGTACTTATctggcctggt                                                         <  1:2080678/67‑1 (MQ=255)
   tCCTCAACGATAACGATTCCATGTTCAACACGCCGCCGACATTTGCCTGGTACTTATctggcctggt                                                         <  1:169106/67‑1 (MQ=255)
   tCCTCAACGATAACGATTCCATGTTCAACACGCCGCCGACATTTGCCTGGTACTTATctggcctggt                                                         <  1:2652213/67‑1 (MQ=255)
   tCCTCAACGATAACGATTCCATGTTCAACACGCCGCCGACATTTGCCTGGTACTTATctggcctggt                                                         <  1:268536/67‑1 (MQ=255)
   tCCTCAACGATAACGATTCCATGTTCAACACGCCGCCGACATTTGCCTGGTACTTATctggcctggt                                                         <  1:272242/67‑1 (MQ=255)
    ccTCAACGATAACGATTCCATGTTCAACACGCCGCCGACATTTGCCTGGTACTTATctggcctggt                                                         <  1:2123783/66‑1 (MQ=255)
                   cTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTATCTGGTCTGGTCTTTa                                                    <  1:2952560/56‑1 (MQ=255)
                                      ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg                      <  1:511090/67‑1 (MQ=255)
                                      ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg                      <  1:383505/67‑1 (MQ=255)
                                      ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg                      <  1:716339/67‑1 (MQ=255)
                                      ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg                      <  1:1186500/67‑1 (MQ=255)
                                      ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg                      <  1:3231803/67‑1 (MQ=255)
                                      ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg                      <  1:2939064/67‑1 (MQ=255)
                                      ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg                      <  1:576087/67‑1 (MQ=255)
                                      ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg                      <  1:2922007/67‑1 (MQ=255)
                                      ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg                      <  1:2765107/67‑1 (MQ=255)
                                      ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg                      <  1:2321859/67‑1 (MQ=255)
                                      ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg                      <  1:2136972/67‑1 (MQ=255)
                                      ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg                      <  1:2012570/67‑1 (MQ=255)
                                      ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg                      <  1:164962/67‑1 (MQ=255)
                                      ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg                      <  1:1649539/67‑1 (MQ=255)
                                      ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg                      <  1:1592519/67‑1 (MQ=255)
                                      ccgACATTTGCCTGGTACTTATCTGGCCTGGGCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg                      <  1:451959/67‑1 (MQ=255)
                                                         tatctGGTCTGGTCTTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATg                 >  1:201560/1‑53 (MQ=255)
                                                         tatctGGTCTGGTCTTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAatcaatca   >  1:2363037/1‑67 (MQ=255)
                                                         tatctGGTCTGGTCTTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAg  >  1:2506641/1‑68 (MQ=255)
                                                         tatctGGTCTGGTCTTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAg  >  1:1719596/1‑68 (MQ=255)
                                                                |                                                            
CCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTATCTGGTCTGGTCTTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAG  >  minE/615037‑615161

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: