Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 588,490:1 +AA 100% intergenic (‑217/‑52) ybdK ← / → insJ‑1 gamma‑glutamyl:cysteine ligase/IS150 hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606588,4901.A100.0% 72.9 / NA 22intergenic (‑217/‑52)ybdK/insJ‑1gamma‑glutamyl:cysteine ligase/IS150 hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (12/10);  total (12/10)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*REL606588,4902.A100.0% 72.1 / NA 22intergenic (‑217/‑52)ybdK/insJ‑1gamma‑glutamyl:cysteine ligase/IS150 hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (12/10);  total (12/10)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTCTTTTTTTTGGTGATGGTGAACTGATGGTGCCTGAAGCAATTTGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑AAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTC  >  REL606/588352‑588627
                                                                                                                                           ||                                                                                                                                         
ttCTTTTTTTTGGTGATGGTTAACTGATGGTGCCTGAAGCAATTTGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaagttgtac                                                                                                                                <  1:222932/150‑10 (MQ=255)
ttCTTTTTTTTGGTGATGGTGAACTGATGGTGCCTGAAGCAATTTGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑aaaagttgtac                                                                                                                                <  2:608947/150‑10 (MQ=255)
      tttttGGTGATGGTGAACTGATGGTGCCTGAAGCAATTTGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCAcc                                                                                                                            >  2:55864/1‑150 (MQ=255)
                     aaCTGATGGTGCCTGAAGCAATTTGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGa                                                                                                             >  2:266160/1‑150 (MQ=255)
                                    aaGCAATTTGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGcc                                                                                              >  2:580496/1‑150 (MQ=255)
                                     aGCAATTTGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGccc                                                                                             <  2:653835/150‑1 (MQ=255)
                                     aGCAATTTGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGccc                                                                                             <  2:154589/150‑1 (MQ=255)
                                           ttGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAata                                                                                       <  2:630417/150‑1 (MQ=255)
                                            tGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAatat                                                                                      <  1:607318/150‑1 (MQ=255)
                                                            gtgACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATCTGTctcttatacacatc                                                                      >  2:342964/1‑137 (MQ=37)
                                                            gtgACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATAc                                                                      >  2:262268/1‑150 (MQ=255)
                                                            gtgACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTCACGATTTAGCCATGAAAAAAAAAACGGGTGCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATAAACCGGAATAGCCCCTAATTTGTAGAATCCAAAATaa                                                                      >  1:82312/1‑149 (MQ=21)
                                                                  aGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtttt                                                                >  2:522061/1‑150 (MQ=255)
                                                                  aGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtttt                                                                >  2:138470/1‑150 (MQ=255)
                                                                  aGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAAAGCCCCTAATATGTCAAAGCCAAAATACCCtttt                                                                >  2:606232/1‑150 (MQ=255)
                                                                   gTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTg                                                               >  1:399497/1‑150 (MQ=255)
                                                                      aTGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGaaa                                                            <  1:578798/150‑1 (MQ=255)
                                                                       tGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGaaaa                                                           >  1:272447/1‑150 (MQ=255)
                                                                                                     gACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTcac                             <  2:400651/150‑1 (MQ=255)
                                                                                                        gtgtGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAAc                          <  1:309666/147‑1 (MQ=255)
                                                                                                        gATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAAc                          <  1:34594/150‑1 (MQ=255)
                                                                                                            gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACtgat                      >  2:180706/1‑150 (MQ=255)
                                                                                                                  tGCAGAAAAAAACGGGTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGt                <  1:520563/150‑1 (MQ=255)
                                                                                                                                ggTTCAGCTTCAAAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCAtctc  <  1:55700/150‑1 (MQ=255)
                                                                                                                                           ||                                                                                                                                         
TTCTTTTTTTTGGTGATGGTGAACTGATGGTGCCTGAAGCAATTTGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATGCAGAAAAAAACGGGTTCAGCTTC‑‑AAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTC  >  REL606/588352‑588627

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