Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,558,775 Δ1 bp 100% coding (663/1425 nt) hsdS ← specificity determinant for hsdM and hsdR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,558,7750T.100.0% 116.7 / NA 26coding (663/1425 nt)hsdSspecificity determinant for hsdM and hsdR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (12/14);  total (12/14)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AATTAATTCTATACCGTCTTTTATATCATCTCCAGGTTGCACCACACCATAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAGGCGTGCTTTGGTGCTGTCTACCTG  >  REL606/4558634‑4558903
                                                                                                                                             |                                                                                                                                
aattaattCTATACCGTCTTTTATATCATCTCCAGGTTGCACCACACCATAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttttttatt                                                                                                                          >  1:101843/1‑147 (MQ=255)
     attCTATACCGTCTTTTATATCATCTCCAGGTTGCACCACACCATAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATaaa                                                                                                                    >  1:474332/1‑150 (MQ=255)
     attCTATACCGTCTTTTATATCATCTCCAGGTTGCACCACACCATAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATaaa                                                                                                                    >  2:609335/1‑150 (MQ=255)
     attCTATACCGTCTTTTATATCATCTCCAGGTTGCACCACACCATAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATaaa                                                                                                                    >  2:220589/1‑150 (MQ=255)
               gTCTTTTATATCATCTCCAGGTTGCACCACACCATAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTAtct                                                                                                          <  1:391374/150‑1 (MQ=255)
                 cTTTTATATCATCTCCAGGTTGCACCACACCATAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTAtcttc                                                                                                        >  2:221067/1‑150 (MQ=255)
                            tctcCAGGTTGCACCACACCATAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAAcc                                                                                             >  1:141904/1‑150 (MQ=255)
                                  tataagagacagccaTAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAg                                                                                       <  2:54669/138‑1 (MQ=255)
                                  ggTTTCACCACACCATAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAg                                                                                       <  1:501593/150‑1 (MQ=255)
                                  ggTTGCACCACACCATAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAg                                                                                       <  1:607557/150‑1 (MQ=255)
                                  ggTTGCACCACACCATAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAg                                                                                       >  2:210858/1‑150 (MQ=255)
                                  ggTTGCACCACACCATAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAg                                                                                       <  1:385831/150‑1 (MQ=255)
                                  ggTTGCACCACACCATAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAg                                                                                       <  2:102141/150‑1 (MQ=255)
                                  ggTTGCACCACACCATAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAg                                                                                       <  2:633046/150‑1 (MQ=255)
                                  ggTTGCACCACACCATAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAg                                                                                       <  1:143492/150‑1 (MQ=255)
                                      agacagcaccaTAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTAc                                                                                   <  2:30066/144‑1 (MQ=255)
                                         ccacaccaTAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGc                                                                                >  1:371036/1‑150 (MQ=255)
                                            caccaTAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTAcctgtct                                                                             >  1:29976/1‑144 (MQ=255)
                                            caccaTAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGc                                                                             >  1:91144/1‑150 (MQ=255)
                                              ccaTAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAgctgtctcttata                                                                           >  1:54507/1‑138 (MQ=255)
                                                ataacataaTGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGcc                                                                         <  2:372111/150‑1 (MQ=255)
                                                   acataaTGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAac                                                                      >  2:579686/1‑150 (MQ=255)
                                                                                                 cTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAGGc                        >  2:86643/1‑150 (MQ=255)
                                                                                                                 tAAGGATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAGGCGTGCTTTGGTGCTGTc        <  1:411929/150‑1 (MQ=255)
                                                                                                                  aaGGATTTTGTAGGTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTATTTGGTCAACCTATCAGTTACTGCCGCAGCCAACACCGCCTGTCGAAAGCGTTTCAGGATTTGTGGGATTTGCTCAAGGCGTGCTTTGGTGCTGTCt       <  1:285641/150‑1 (MQ=255)
                                                                                                                     gATTTTGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAGGCGTGCTTTGGTGCTGTCTAcc    <  1:435852/150‑1 (MQ=255)
                                                                                                                       ttttGTAGTTATCTAATTCAACTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAGGCGTGCTTTGGTGCTGTCTACCTg  <  2:333782/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                
AATTAATTCTATACCGTCTTTTATATCATCTCCAGGTTGCACCACACCATAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAGGCGTGCTTTGGTGCTGTCTACCTG  >  REL606/4558634‑4558903

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