Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,271,926 T→C 100% D38G (GAC→GGC)  yrbA ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,271,9260TC100.0% 107.9 / NA 30D38G (GAC→GGC) yrbApredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/10);  total (20/10)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCATGGGATCTTCAATCATCAGTTAATAAAC  >  REL606/3271798‑3272065
                                                                                                                                |                                                                                                                                           
gTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCaa                                                                                                                        >  2:639019/1‑150 (MQ=255)
           cgcgcCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAg                                                                                                             <  2:149894/150‑1 (MQ=255)
            gcgcCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGt                                                                                                            <  1:604535/150‑1 (MQ=255)
                            ggTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGa                                                                                            >  2:223410/1‑150 (MQ=255)
                            ggTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTGTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCAAGCGCCGGa                                                                                            >  1:207990/1‑150 (MQ=255)
                              tataCGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGaaa                                                                                          <  2:191436/150‑1 (MQ=255)
                                     ttggATCGACACAGCATGAATGCGGTTATCCGCAATCTATTCCATCAGCGGCCCATAGTCCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAAAAACTCACCCTCGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGaa                                                                                   <  1:96798/147‑2 (MQ=255)
                                        gATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTc                                                                                >  2:646726/1‑150 (MQ=255)
                                                   cATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGgaga                                                                     >  2:225665/1‑150 (MQ=255)
                                                      gAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGagagag                                                                  <  2:92743/150‑1 (MQ=255)
                                                      gAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGagagag                                                                  >  1:453121/1‑150 (MQ=255)
                                                              ataagagacagatatTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtcatc                                                          <  2:489569/139‑1 (MQ=255)
                                                                        tatatTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGctct                                                >  2:366388/1‑150 (MQ=255)
                                                                         atatTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtcatcctgtctcttat                                               >  1:488183/1‑139 (MQ=255)
                                                                                 aagagacagcATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCat                                       <  2:667004/141‑1 (MQ=255)
                                                                                      ggACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTAttt                                  >  1:173637/1‑150 (MQ=255)
                                                                                       gACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTAtttt                                 >  1:373188/1‑150 (MQ=255)
                                                                                        aCCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCTGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTc                                >  1:138328/1‑150 (MQ=255)
                                                                                        aCCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTc                                >  2:468504/1‑150 (MQ=255)
                                                                                        aCCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACCTCCTGGAGGGAGAAAGCGTTCATCAACACGCTCTGAATTTCATTATTTTc                                >  1:381244/1‑150 (MQ=255)
                                                                                        aCCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAAAGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTc                                >  1:683030/1‑150 (MQ=255)
                                                                                         ccATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTcc                               >  2:85713/1‑150 (MQ=255)
                                                                                         ccATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTcc                               >  2:388202/1‑150 (MQ=255)
                                                                                         ccATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTcc                               >  2:629470/1‑150 (MQ=255)
                                                                                         ccATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTcc                               >  1:205873/1‑150 (MQ=255)
                                                                                         ccATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGACCGTTCCTCAACACGCCCTGAATTTCCTTATTTTcc                               >  1:507695/1‑150 (MQ=255)
                                                                                          cATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCatctgtctctt                              >  1:665040/1‑141 (MQ=255)
                                                                                                                  cccGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCATGGGATCTTCAATCATCAGTtaat      <  2:509143/150‑1 (MQ=255)
                                                                                                                     cgcTCATGCCGCCAAACCACTCAACCACCGCAATAAACTGAAAGTGGCTGCCATCGGCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCATGGGATCTTCAATCATCAGTTAATaaa   <  2:42440/148‑1 (MQ=255)
                                                                                                                      aCTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCATGGGATCTTCAATCATCAGTTAATAAAc  <  2:585357/150‑1 (MQ=255)
                                                                                                                                |                                                                                                                                           
GTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCATGGGATCTTCAATCATCAGTTAATAAAC  >  REL606/3271798‑3272065

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