Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,513,930 A→T 100% L413F (TTA→TTT USA300HOU_RS07375 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,513,9300AT100.0% 64.7 / NA 22L413F (TTA→TTTUSA300HOU_RS07375hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (13/9);  total (13/9)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATTTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTT  >  USA300TCH1516_ALE/1513794‑1514062
                                                                                                                                        |                                                                                                                                    
atttatttGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTg                                                                                                                                  <  2:24841/141‑1 (MQ=255)
 tttatttGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGt                                                                                                                                 >  2:316214/1‑141 (MQ=255)
     tttGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGAGAAAGCGTACAAAAAATCATGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTa                                                                                                                             >  2:300941/1‑141 (MQ=255)
              atCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGAttatta                                                                                                                    >  1:26488/1‑141 (MQ=255)
                 gCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATg                                                                                                                 >  1:352333/1‑141 (MQ=255)
                     tattaGTAGCCATTGGGAGGGGGTTTGTATCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCat                                                                                                             <  2:33803/141‑1 (MQ=255)
                       ttaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCatta                                                                                                           >  1:163622/1‑141 (MQ=255)
                         agtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCCTAATATGTTTTATTGTATTAGGATTATTAATGCCattctt                                                                                                         >  1:133058/1‑141 (MQ=255)
                                 ttGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAATATTTATGGAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTAt                                                                                                 <  2:383878/141‑1 (MQ=255)
                                           gTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGCTTATCACTTATCCt                                                                                       <  1:189039/141‑1 (MQ=255)
                                                 ctctTGGATAAAACGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGAATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTAtt                                                                                 >  2:73406/1‑141 (MQ=255)
                                                     tGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATACCTTATCCTTGTATTtata                                                                             >  2:370927/1‑141 (MQ=255)
                                                             gCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTCTTTg                                                                                                       >  1:249889/1‑107 (MQ=255)
                                                             gCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTg                                                                                                       <  2:249889/107‑1 (MQ=255)
                                                             gCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAAtt                                                                     >  2:126245/1‑141 (MQ=255)
                                                                     aaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAAtattat                                                             <  1:313180/141‑1 (MQ=255)
                                                                           cTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTCAATTAATTTTCTATATTa                                                       >  2:71940/1‑141 (MQ=255)
                                                                           atGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATATATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTa                                                       <  1:267527/140‑1 (MQ=255)
                                                                                  aaaaGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTCTTTATAATTGAATTAATATTATATATTAAAGAtcg                                                <  1:356642/141‑1 (MQ=255)
                                                                                           acaTAGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTaa                                       <  2:204058/141‑1 (MQ=255)
                                                                                                                               aataTGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACTATAATGCACTTAAAGTTTTGAACTGACGAAAtt   >  2:259079/1‑141 (MQ=255)
                                                                                                                                atatGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAAttt  >  2:309768/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                    
ATTTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTT  >  USA300TCH1516_ALE/1513794‑1514062

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: