Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,268,951 T→G 100% I129S (ATC→AGC)  atoC → fused response regulator of ato operon, in two‑component system with AtoS: response regulator/sigma54 interaction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,268,9510TG93.5% 24.0 / ‑6.2 10I129S (ATC→AGC) atoCfused response regulator of ato operon, in two‑component system with AtoS: response regulator/sigma54 interaction protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (3/6);  total (3/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCCTATGCGGAAGTGGAAACCGCCGTCGAAGCGCTACGCTGCGGAGCCTTCGACTATGTTATTAAACCGTTTGATCTCGATGAGTTGAATTTAATCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGATCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATGGGGGCACATTCTCACCAACAGCCCGGCGATGATGGACATCTGCAAAGACACCGCCAAAAT  >  REL606/2268815‑2269050
                                                                                                                                        |                                                                                                   
gCCTATGCGGAAGTGGATACCGCCGTCGAAGCGCGACGCTGCGGAGCCTTCGACTATGTTATTAACCCGTTTGATCTCGATGAGTTGAATTTAATCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCAc                                                                                        <  1:436149/150‑1 (MQ=255)
gCCTATGCGGAAGTGGAAACCGCCGTCGAAGCGCTACGCTGCGGAGCCTTCGACTATGTTATTAAACCGTTTGATCTCGATGAGTTGAATTTAATCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCAc                                                                                        <  2:592707/150‑1 (MQ=255)
            gTGGAAACCGCCGTCGAAGCGCTACGCTGCGGAGCCTTCGACTATGTTATTAAACCGTTTGATCTCGATGAGTTGAATTTAATCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGAGc                                                                            >  1:491082/1‑150 (MQ=255)
                 aaCCGCCGTCGAAGCGCTACGCTGCGGAGCCTTCGACTATGTTATTAAACCGTTTGATCTCGATGAGTTGAATTTAATCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAg                                                                       <  1:362649/150‑1 (MQ=255)
                                 cTACGCTGCGGAGCCTTCGACTATGTTATTAAACCGTTTGATCTCGATGAGTTGAATTTAATCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATGGGGGcac                                                       >  2:244180/1‑150 (MQ=255)
                                             gCCTTCGACTATGTTATTAAACCGTTTGATCTCGATGAGTTGAATTTAATCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCATCATCTGCACCAGGCACTGAGCACCAGCTGGCAATGGGGGCACATTCTCACCAAc                                           <  2:23852/150‑1 (MQ=255)
                                                               aaaCCGTTTGATCTCGATGAGTTGAATTTAATCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATGGGGGCACATTCTCACCAACAGCCCGGCGATGATGGAc                         <  2:479416/150‑1 (MQ=255)
                                                                       ttatgtcGAGGCGTTGAATTTAGTCGTTCAGGCGGCTTTACAACTCCAGTCAATGAAAAACGAGGTCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATGGGGGCACATTCTCACCAACAGCCCGGCGATGATGGACATCTGCaa                 <  1:155240/145‑1 (MQ=255)
                                                                          tctcGATGAGTTGAATTTAATCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATGGGGGCACATTCTCACCAACAGCCCGGCGATGATGGACATCTGCAAAGa              <  1:127460/150‑1 (MQ=255)
                                                                                      gAATTTAATCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATGGGGGCACATTCTCACCAACAGCCCGGCGATGATGGACATCTGCAAAGACACCGCCAAAAt  >  2:139439/1‑150 (MQ=255)
                                                                                                                                        |                                                                                                   
GCCTATGCGGAAGTGGAAACCGCCGTCGAAGCGCTACGCTGCGGAGCCTTCGACTATGTTATTAAACCGTTTGATCTCGATGAGTTGAATTTAATCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGATCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATGGGGGCACATTCTCACCAACAGCCCGGCGATGATGGACATCTGCAAAGACACCGCCAAAAT  >  REL606/2268815‑2269050

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: