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

ATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTT  >  USA300TCH1516_ALE/1513804‑1514062
                                                                                                                              |                                                                                                                                    
ttgtatCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGatt                                                                                                                        <  2:806881/140‑1 (MQ=255)
atgtatCGCTTTAGTAGTAGCCATAGGCATGGGGTTTGTATCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGatt                                                                                                                        <  2:240548/141‑1 (MQ=255)
   tatCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTCGAattatt                                                                                                                     >  2:540926/1‑141 (MQ=255)
   taccacTTTATTACTAGCCATTGACATGCGTTTTGTCGCTGGGATAAAGCGTATAAAAAATACTGCAATAAAAGTACCAGATCGCGTAAAAAAACCAGCATCTATATCACTCACAATTTGTTTTATTGTATTAGGattatt                                                                                                                     <  1:143843/136‑1 (MQ=255)
       gCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATg                                                                                                                 >  1:434666/1‑141 (MQ=255)
        cTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGc                                                                                                                >  1:618404/1‑141 (MQ=255)
                gtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTAttt                                                                                                        >  1:367420/1‑141 (MQ=255)
                gtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTAttt                                                                                                        >  1:367607/1‑141 (MQ=255)
                     tagtGGCACGGGGTTTGTCTCCTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAACAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGAtt                                                                                                   <  2:458439/138‑1 (MQ=255)
                      aTTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTa                                                                                                  >  1:414812/1‑141 (MQ=255)
                       ttGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTAt                                                                                                 <  2:434666/141‑1 (MQ=255)
                         ggCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTTGGATTATTAATGCCATTATTTGGATTATCa                                                                                               >  1:121803/1‑141 (MQ=255)
                         ggCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCa                                                                                               >  2:104754/1‑141 (MQ=255)
                         ggCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCAATTATTTGGATTATcc                                                                                               >  2:104852/1‑140 (MQ=255)
                         ggCATGGGGTTTGTCTCTTGCATAAAGCGTACAAAAAATACTGCTATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATGAGTCATAATATGTTTTATTGTTTTAGGATTATTAATGCCATTATTTGGATTATcc                                                                                               >  2:498619/1‑140 (MQ=255)
                           cATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACt                                                                                             <  1:837192/141‑1 (MQ=255)
                            aTGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACtt                                                                                            >  2:887498/1‑141 (MQ=255)
                             tGGGGTTTGTCTCTGGGATCAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACttt                                                                                           >  1:500469/1‑140 (MQ=255)
                             tGGGGTTTGGCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTa                                                                                           <  1:877199/141‑1 (MQ=255)
                             aggggTTTGTCTATTGGATAAAACGTACAAAAAATACTGCAATAAAAATACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTa                                                                                           <  2:847549/140‑1 (MQ=255)
                               ggttttGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATc                                                                                         <  2:427820/138‑1 (MQ=255)
                                             gATAAAGCGTACAAAAATTACTGCAATAATAGTACCACATCGCGTTAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTGTTTATAAt                                                                           <  2:401969/141‑1 (MQ=255)
                                                     gTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTat                                                                   >  1:522315/1‑140 (MQ=255)
                                                     gTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTaa                                                                   >  1:521968/1‑141 (MQ=255)
                                                      tACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATAACATATACTTGTATTTCTAATTGAATTAat                                                                  >  2:267638/1‑141 (MQ=255)
                                                        cAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAatat                                                                <  2:858848/141‑1 (MQ=255)
                                                         aaaaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAAtatt                                                               >  2:397324/1‑141 (MQ=255)
                                                           aaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAAtattat                                                             >  2:510289/1‑141 (MQ=255)
                                                            aaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATtata                                                            >  1:57296/1‑141 (MQ=255)
                                                            aaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACGCATAATATGGGTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATtata                                                            >  1:718617/1‑141 (MQ=255)
                                                                   gCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATAAAAACCACTCAGAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTaaa                                                     >  2:598072/1‑141 (MQ=255)
                                                                   gCAATAAAAGTACCACATCGCGTAAAAAAACCAGAAACTAGAGCAGTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATAACTTATCTTTGTATTTATAATTGTATTTATATTCTATTTTaaa                                                     >  2:598175/1‑141 (MQ=255)
                                                                    cAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAg                                                    >  2:795833/1‑141 (MQ=255)
                                                                    cAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAg                                                    >  2:640187/1‑141 (MQ=255)
                                                                      aTAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATCTTAAAGAt                                                  >  2:234258/1‑141 (MQ=255)
                                                                         aaaGTACCACATGGCGTAAAAAAACCAGTATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAtcgt                                               <  1:739819/141‑1 (MQ=255)
                                                                         aaaGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAtcgt                                               <  1:139078/141‑1 (MQ=255)
                                                                         aaaGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAtcgt                                               <  2:502540/141‑1 (MQ=255)
                                                                         aaaGTACCACAGCGCGTAAAACAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAtcgt                                               <  1:316396/141‑1 (MQ=255)
                                                                            gTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTGTTTATAATTGAATTAATATTATATATTAAAGATcgtcgt                                            >  2:869336/1‑141 (MQ=255)
                                                                                           aaaaaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGc                             >  1:74995/1‑141 (MQ=255)
                                                                                           aaaaaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGc                             >  1:273407/1‑141 (MQ=255)
                                                                                           aaaaaaCCAGCATCTATATCACCCATAATATGTTTTATTATATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATAAAAGATCGTCGTGCTAAAAAAAAATga                             >  2:690901/1‑140 (MQ=255)
                                                                                                aCCAGCATCTCTATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTAGATATTAAAGATCGTCGTGCTAAACAATAATGCACTTa                        <  2:207405/141‑1 (MQ=255)
                                                                                                 ccAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTaa                       <  2:57296/141‑1 (MQ=255)
                                                                                                 ccAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTaa                       <  2:718617/141‑1 (MQ=255)
                                                                                                   aGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAg                     >  2:645061/1‑141 (MQ=255)
                                                                                                     cATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGtt                   <  1:780368/141‑1 (MQ=255)
                                                                                                      aTCTATATCACTCATAATATGTTTTACGGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGttt                  <  1:386621/141‑1 (MQ=255)
                                                                                                         tataTCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCGATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGa               >  2:643286/1‑141 (MQ=255)
                                                                                                               tcttatattaTGTTTTATTTTATTAGGTTTATTAATTCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGAc         <  2:121803/137‑1 (MQ=255)
                                                                                                                   ataataTGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGaaa     <  2:367420/141‑1 (MQ=255)
                                                                                                                   ataataTGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGaaa     <  2:367607/141‑1 (MQ=255)
                                                                                                                    taataTGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAAt    >  1:331901/1‑141 (MQ=255)
                                                                                                                     aatgTGTTTTATTGTATTAGTATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAAtt   <  2:385609/141‑1 (MQ=255)
                                                                                                                      atatGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAAttt  <  2:110421/141‑1 (MQ=255)
                                                                                                                              |                                                                                                                                    
ATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTT  >  USA300TCH1516_ALE/1513804‑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: ‑