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.2 / NA 37L413F (TTA→TTT) USA300HOU_RS07375hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (23/14);  total (23/14)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TAATTTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTACAA  >  USA300TCH1516_ALE/1513792‑1514066
                                                                                                                                          |                                                                                                                                        
tAATTTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTAt                                                                                                                                        >  1:353904/1‑141 (MQ=255)
 aaTTTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTAtt                                                                                                                                       <  1:115051/141‑1 (MQ=255)
               tatCGATTTATAAGTAGACAATGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGattatt                                                                                                                         <  2:131107/141‑1 (MQ=255)
                    cTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGc                                                                                                                    >  2:215344/1‑141 (MQ=255)
                         ttaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAACAAATACTGCAATAAGGGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATCATATGTTCTATTGTATTAGTCTTATTACTGCCatca                                                                                                               >  2:12935/1‑139 (MQ=255)
                         ttaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCatta                                                                                                               <  1:300211/141‑1 (MQ=255)
                          tagtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattat                                                                                                              >  1:598894/1‑141 (MQ=255)
                          tagtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattat                                                                                                              >  1:477719/1‑141 (MQ=255)
                             tagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTg                                                                                                           <  2:374085/141‑1 (MQ=255)
                             tagCCATTAGAATGGGGTTTGTATCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTg                                                                                                           <  2:373334/141‑1 (MQ=255)
                                 cATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGAtt                                                                                                       >  1:339159/1‑141 (MQ=255)
                                 cATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGAtt                                                                                                       >  1:171593/1‑141 (MQ=255)
                                     ggCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCa                                                                                                   >  2:87662/1‑141 (MQ=255)
                                          ggggTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCATTAACAGTACCACATCGCGTAAAAAAACAAGCATCTATATAATTCATTATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTAt                                                                                              >  2:422058/1‑141 (MQ=255)
                                          ggggTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAGGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTAt                                                                                              >  2:271959/1‑141 (MQ=255)
                                            ggTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATcc                                                                                            >  1:405553/1‑141 (MQ=255)
                                             ttttGTATCTTGTATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCt                                                                                           <  1:13794/140‑1 (MQ=255)
                                             gTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACGCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCt                                                                                           <  1:143194/141‑1 (MQ=255)
                                               ttGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAAGTCCATTATTTGGATTATCACTTATCCttt                                                                                         >  2:537152/1‑140 (MQ=255)
                                                    tctTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAGAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTAttt                                                                                    >  1:151005/1‑141 (MQ=255)
                                                           tAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATAACTTATCCTTGTATTTATAATTg                                                                             >  1:579800/1‑141 (MQ=255)
                                                            aaaGCGTACAAAAAATACTGCAAGAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGa                                                                            <  1:112810/141‑1 (MQ=255)
                                                                      aaaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAAtatta                                                                  <  1:215344/141‑1 (MQ=255)
                                                                        aaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATtata                                                                <  1:275024/141‑1 (MQ=255)
                                                                                 aaTAAAAGTACCACATCGCGTAAAAAACCCAGCATATATATAACTCCTAATATGTTTTATCGTATTCCCAACACTACGCTCTTTATTTGGATTATCACTTATCCTTGCATTTATCATCTACATAATATTATATATTAAAGa                                                       <  2:353904/141‑1 (MQ=255)
                                                                                    aaaaGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCCTTAATTCGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAtcg                                                    >  1:148218/1‑141 (MQ=255)
                                                                                    aaaaGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAATAtcg                                                    >  2:624677/1‑141 (MQ=255)
                                                                                       aGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAtcgtcg                                                 >  1:25218/1‑141 (MQ=255)
                                                                                          aCCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGc                                              >  2:372903/1‑141 (MQ=255)
                                                                                          aCCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTAGTGc                                              >  2:305308/1‑141 (MQ=255)
                                                                                              caTCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCAATTATTTGGATTATCACTTATCCTTATATTTATAATTGAATTAATATTATATATTAAAGATCTTCTTTCTaaa                                          >  2:108644/1‑141 (MQ=255)
                                                                                                 cgcgTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACaa                                       >  2:271211/1‑141 (MQ=255)
                                                                                                      aaaaaaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGCATTAATATTATATATTAAAGTTCGTCGTTCTAAACAATACTg                                  >  2:315652/1‑141 (MQ=255)
                                                                                                                 aaTCTATATCAATAATAATATTTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGtt                       <  2:179252/140‑1 (MQ=255)
                                                                                                                               ataataTGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGaaa         <  1:87662/141‑1 (MQ=255)
                                                                                                                                    atGTTTTATTGTATTAGGATTATTAATGCCATGATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTAc    <  1:412311/141‑1 (MQ=255)
                                                                                                                                      gTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTACaa  >  1:228269/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                        
TAATTTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTACAA  >  USA300TCH1516_ALE/1513792‑1514066

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