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% 113.6 / ‑6.5 38L413F (TTA→TTT) USA300HOU_RS07375hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (21/16);  minor base C (0/1);  total (21/17)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTACA  >  USA300TCH1516_ALE/1513796‑1514065
                                                                                                                                      |                                                                                                                                       
ttatttGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAACAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTa                                                                                                                                   >  2:510132/1‑141 (MQ=255)
            atCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGAttatta                                                                                                                       <  1:767088/141‑1 (MQ=255)
                      tagtagCTATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattat                                                                                                             <  2:353495/141‑1 (MQ=255)
                      tagtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTTTATTAGGATTATTAATGCCatgat                                                                                                             >  2:250985/1‑141 (MQ=255)
                        gtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTCTTAGGCTTATTTATGCCATTAttt                                                                                                           >  2:361194/1‑141 (MQ=255)
                        gtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTAttt                                                                                                           <  1:550218/141‑1 (MQ=255)
                            ccATGGGCATGGGGTTTGTCTCTTGGATAAAACGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGAt                                                                                                       <  2:411030/141‑1 (MQ=255)
                             cATTGGGATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGAtt                                                                                                      <  2:316292/141‑1 (MQ=255)
                                tGGCATGGGGTCTGTCTATTAGATCAATCGTACAAAAAATAACACCATTATAGTTTCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATc                                                                                                   <  1:559567/141‑1 (MQ=255)
                                  gCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATTCCATTATTTGGATTATCAc                                                                                                 >  2:583519/1‑141 (MQ=255)
                                       gggTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATc                                                                                            >  2:315639/1‑141 (MQ=255)
                                          tttGTCTCGTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACGACATCGCGTAAAAAAACCAGCATCTATATGACTCATAATATGTTGTATTGTATTAGGATTATTAATGCCATTATTTGGATTATTACTTAGCCtt                                                                                         >  2:204614/1‑141 (MQ=255)
                                              tctcTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTAt                                                                                     <  1:2800/141‑1 (MQ=255)
                                               ctctTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTAGTAGGATTATTAATGCAATGAATTGGATTATCACTTATCCTTGTAtt                                                                                    >  2:123038/1‑141 (MQ=255)
                                                      aTAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAAtt                                                                             >  2:39431/1‑141 (MQ=255)
                                                       tAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTg                                                                            <  2:516327/141‑1 (MQ=255)
                                                         aaGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGaa                                                                          >  2:489150/1‑141 (MQ=255)
                                                                    aaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATtata                                                               >  2:770477/1‑141 (MQ=255)
                                                                          tGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTaa                                                         >  2:334659/1‑141 (MQ=255)
                                                                           gCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTaaa                                                        >  2:488437/1‑141 (MQ=255)
                                                                                      aCCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGc                                             >  1:762809/1‑141 (MQ=255)
                                                                                       ccACATCGCGTAAACAAACCAGCATCTATATCACTCATAATATGTTTCATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCt                                            <  1:202937/141‑1 (MQ=255)
                                                                                          caTCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTaaa                                         <  1:39431/141‑1 (MQ=255)
                                                                                           aacgcgTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAAc                                        <  2:675814/139‑1 (MQ=255)
                                                                                            tCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACa                                       <  2:55010/141‑1 (MQ=255)
                                                                                                    aaaaaCCAGCATCTTTATCACTCATAATGTGTTTTTTTGTGCGAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCa                               <  1:524046/141‑1 (MQ=255)
                                                                                                      aaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTACCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACt                             >  2:617676/1‑141 (MQ=255)
                                                                                                         ccAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATAACTTATCCTTGTATTTATAGTTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTaa                          <  1:370025/141‑1 (MQ=255)
                                                                                                              aTCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGttt                     <  1:55871/141‑1 (MQ=255)
                                                                                                               tCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATACACTTAAAGtttt                    >  2:563079/1‑141 (MQ=255)
                                                                                                                     tCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTg              >  2:442451/1‑141 (MQ=255)
                                                                                                                      cACTCATAATATTTTTTATTGTACTAGGATTATTAATGACATTATTTGGATTATCACTTATCCTTTTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGa             <  2:504512/141‑1 (MQ=255)
                                                                                                                        cTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACg           >  1:246755/1‑141 (MQ=255)
                                                                                                                          cATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGaa         <  1:419264/141‑1 (MQ=255)
                                                                                                                            taataTGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAAt       >  2:435972/1‑141 (MQ=255)
                                                                                                                             aataTGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAAtt      >  2:600521/1‑141 (MQ=255)
                                                                                                                             aataGGTTTTATTGTATTGGGATTATTAATGCCATTATTTGGATTATCAATTATTCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCTTGCTAAACAATAATTAACTTAAAGATTTAAACTGACGAAAtt      >  2:535901/1‑141 (MQ=255)
                                                                                                                                 tGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTACa  >  2:55105/1‑141 (MQ=255)
                                                                                                                                      |                                                                                                                                       
TTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTACA  >  USA300TCH1516_ALE/1513796‑1514065

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: ‑