Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,879,026 T→A 100% D380V (GAT→GTT)  yoaE ← fused predicted membrane protein/conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,879,0260TA100.0% 54.5 / NA 17D380V (GAT→GTT) yoaEfused predicted membrane protein/conserved protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (9/8);  total (9/8)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGAATTCACCCGCAATGGCTTCCAGCACATCCAGCGGCGTGACCAAGCCTTGTACCACACCAAACTCGTTGGTCACGATAACAAAGCTCCCGCGAGCGCGACGCAGCACGCCCAGCAGGTTGATCGGATCGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTACGATTCCGATGATTTCATCCAGTTCACCGCGACATACCGGGAAC  >  REL606/1878898‑1879166
                                                                                                                                |                                                                                                                                            
ggAAGTCACCCGCCACGGCTTCCAGCACAGCCAGCGGCGTGACCAAGCCTTGTACCACACCAAACTCGTTGGTCACGATAACAAAGCTCCCGCGAGCGCGACGCAGCACGCCCAGCAGGTTGATCGGAACGAGGGTTTCCGGGACGATaa                                                                                                                         <  1:462764/150‑1 (MQ=255)
                   ttCCAGCACATCCAGCGGCGTGACCAAGCCTTGTACCACACCAAACTCGTTGGTCACGATAACAAAGCTCCCGCGAGCGCGACGCAGCACGCCCAGCAGGTTGATCGGAACGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGc                                                                                                      <  2:535979/150‑1 (MQ=255)
                                         aCCAAGCCTTGTACCACACCAAACTCGTTGGTCACGATAACAAAGCTCCCGCGAGCGCGACCCAGCACGCCCAGCAGGTTGATCGGAACGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCCATCGCCGCCACATCAACGccc                                                                                >  1:303821/1‑150 (MQ=255)
                                              gCCTTGTACCACACCAAACTCGTTGGTCACGATAACAAAGCTCCCGCGAGCGCGACGCAGCACGCCCAGCAGGTTGATCGGAACGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTc                                                                           >  2:413109/1‑150 (MQ=255)
                                              gCCTTGTACCACACCAAACTCGTTGGTCACGATAACAAAGCTCCCGCGAGCGCGACGCAGCACGCCCAGCAGGTTGATCGGAACGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTc                                                                           >  2:545550/1‑150 (MQ=255)
                                              gCCTTGTACCACACCAAACTCGTTGGTCACGATAACAAAGCTCCCGCGAGCGCGACGCAGCACGCCCAGCAGGTTGATCGGAACGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTc                                                                           >  2:395548/1‑150 (MQ=255)
                                              gCCTTGTACCACACCAAACTCGTTGGTCACGATAACAAAGCTCCCGCGAGCGCGACGCAGCACGCCCAGCAGGTTGATCGGAACGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTc                                                                           >  2:333487/1‑150 (MQ=255)
                                              gCCTTGTACCACACCAAACTCGTTGGTCACGATAACAAAGCTCCCGCGAGCGCGACGAAGCACGCCCAGCAGGTTGAACGGCACGAGGGTTTCCCGGACGATAATCGCCGGAGACGCCGAACCACTCGCCGCCACATCAACGCCCTCTTc                                                                           >  1:140332/1‑150 (MQ=255)
                                               ccTTGTACCACACCAAACTCGTTGGTCACGATAACAAAGCTCCCGCGAGCGCGACGCAGCACGCCCAGCAGGTTGATCGGAACGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTcc                                                                          <  1:543972/150‑1 (MQ=255)
                                               ccTTGCACCACACCAACCGCGCAGGTCACGATAACCAAGCTCCCGCGAGCGCGACGCAGCACGCCCAGCAGGTTGATCGGAACGGGGGTTTCCGGGACGATAATCGCCGGGGACGCCGAAGCAATCGCCGCCACATCAACGCCCGCTTcc                                                                          <  1:621997/150‑1 (MQ=255)
                                                 ttGTACCACACCAAACTCGTTGGTCACGATAACAAAGCTCCCGCGAGCGCGACGCAGCACGCCCAGCAGGTTGATCGGAACGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAg                                                                        >  2:122093/1‑150 (MQ=255)
                                                  tGTACCAAACCAAACTCGTTGGTCACGATAACAAAGCTCCCGCGAGCGCGACGCAGCACGCCCAGCAGGTTGATCGGAACGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAgc                                                                       <  2:140734/150‑1 (MQ=255)
                                                                cTCGTTGGTCACGATAACAAAGCTCCCGCGAGCGCGACGCAGCACGCCCAGCAGGTTGATCGGAACGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGTTc                                                         >  1:538337/1‑150 (MQ=255)
                                                                      ggTCACGATAACAAAGCTCCCGCGAGCGCGACGCAGCACGCCCAGCAGGTTGATCGGAACGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGc                                                   <  2:304681/150‑1 (MQ=255)
                                                                           cGATAACAAAGCTCCCGCGAGCGCGACGCAGCACGCCCAGCAGGTTGATCGGAACGAGGGTTTCCGCGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTa                                              >  2:433169/1‑150 (MQ=255)
                                                                                                    aCGCAGCACGCCCAGCAGGTTGATCGGAACGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTACGATTCCGATGATTTCATCCAGTTc                     <  2:261284/150‑1 (MQ=255)
                                                                                                            cGCCCAGCAGGTTGATCGGAACGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTACGATTCCGATGATTTCATCCAGTTCACCGCGAc             <  1:394402/150‑1 (MQ=255)
                                                                                                                       ttGATCGGAACGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTACGATTCCGATGATTTCATCCAGTTCACCGCGACATACCGGGAAc  <  1:121738/150‑1 (MQ=255)
                                                                                                                                |                                                                                                                                            
GGAATTCACCCGCAATGGCTTCCAGCACATCCAGCGGCGTGACCAAGCCTTGTACCACACCAAACTCGTTGGTCACGATAACAAAGCTCCCGCGAGCGCGACGCAGCACGCCCAGCAGGTTGATCGGATCGAGGGTTTCCGGGACGATAATCGCCGGAGACGCCGAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTACGATTCCGATGATTTCATCCAGTTCACCGCGACATACCGGGAAC  >  REL606/1878898‑1879166

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