Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,101,091 T→C 100% E47G (GAA→GGA)  hslU ← ATP‑dependent protease ATP‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,101,0910TC100.0% 98.0 / NA 27E47G (GAA→GGA) hslUATP‑dependent protease ATP‑binding subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/15);  total (12/15)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCA  >  REL606/4100947‑4101220
                                                                                                                                                |                                                                                                                                 
tCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATgg                                                                                                                              <  2:483584/150‑1 (MQ=255)
   ttACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGgcgc                                                                                                                           <  2:479280/150‑1 (MQ=255)
         aCGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGctct                                                                                                                     >  2:419326/1‑150 (MQ=255)
         aCGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGctct                                                                                                                     >  2:559461/1‑150 (MQ=255)
         aCGGAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCTTTCGCCAGCTTAGCCAGCCGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCCTTAGGATAGTTAGCGGGGTCCCTCCTTGGCGCAGccct                                                                                                                     >  1:186167/1‑150 (MQ=255)
               ccctcgtCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCTCATTCGCCAGCGTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTg                                                                                                               <  1:147456/144‑1 (MQ=255)
                   cTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGTCGCAGCTCTTAGTCGAGCt                                                                                                           >  1:303078/1‑150 (MQ=255)
                      cGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCACCTCTTCGTTGCGCTGCa                                                                                                        >  1:268630/1‑150 (MQ=255)
                          gAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCg                                                                                                    >  1:520915/1‑150 (MQ=255)
                          ataagagacagCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCg                                                                                                    <  1:515917/140‑1 (MQ=255)
                                     cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCgctgtctcttat                                                                                         >  2:517397/1‑139 (MQ=255)
                                          aCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGc                                                                                    <  1:669754/150‑1 (MQ=255)
                                          aCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGc                                                                                    <  2:440955/150‑1 (MQ=255)
                                             ttGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCaga                                                                                 >  1:617232/1‑150 (MQ=255)
                                                            ttCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATGGCCACAGaa                                                                  >  1:171876/1‑150 (MQ=255)
                                                            ttCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGaa                                                                  <  2:353919/150‑1 (MQ=255)
                                                                      tAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGt                                                        >  2:537546/1‑150 (MQ=255)
                                                                        gCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTg                                                      >  2:499326/1‑150 (MQ=255)
                                                                           aGACGACGTGCGATGTCTGTTTTACCGGCACCGGTCGGGCCGGTCCTCAGGATATTTTTCGGGTGCCCTCCATGGCTCAGCTCTTCGTTGGGCTGCATGCGACGCCAGCGGTTACTCAGAGCAATCGCCAGAGAACGCTTGGCGGTGTcc                                                   <  1:85340/150‑1 (MQ=255)
                                                                           aGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTcc                                                   <  2:522415/150‑1 (MQ=255)
                                                                           aGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTcc                                                   <  2:125958/150‑1 (MQ=255)
                                                                                acgGTCTATTTCAGGTTTACCGCCTCCGGTCGGGCCGATCATCAGGATAGTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGcc                                              <  1:514585/150‑1 (MQ=255)
                                                                                     cGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTGCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCgatga                                         <  1:418136/150‑1 (MQ=255)
                                                                                             gTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTa                                 >  2:256215/1‑150 (MQ=255)
                                                                                                ttACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATcc                              <  2:533420/150‑1 (MQ=255)
                                                                                                ttACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATcc                              <  2:619041/150‑1 (MQ=255)
                                                                                                  aCCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAg                            <  2:339337/150‑1 (MQ=255)
                                                                                                  aCCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAg                            <  2:542266/150‑1 (MQ=255)
                                                                                                           ggTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGAc                   >  1:8812/1‑150 (MQ=255)
                                                                                                                            ggATATTTTTCGGGGTCACTCCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCa  >  1:273533/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                 
TCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCA  >  REL606/4100947‑4101220

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: