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

TATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTAC  >  USA300TCH1516_ALE/1513803‑1514064
                                                                                                                               |                                                                                                                                      
gactgcgcGCTGTATTAGTAGCCATTGGTATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGat                                                                                                                           <  1:138521/134‑1 (MQ=255)
   gtatCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGattat                                                                                                                        <  1:598299/141‑1 (MQ=255)
         cTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAGCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGc                                                                                                                  <  1:533294/141‑1 (MQ=255)
            tattaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCat                                                                                                               <  1:658621/141‑1 (MQ=255)
            tattaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCat                                                                                                               >  1:464082/1‑141 (MQ=255)
            tattaGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCat                                                                                                               <  1:44920/141‑1 (MQ=255)
               tagtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCattat                                                                                                            >  2:544438/1‑141 (MQ=255)
                  tagaaaatGGCATAGTGTTTGTCTCTTGGATAAAGCGTAAAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTg                                                                                                         <  2:68806/134‑1 (MQ=255)
                           gCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCAc                                                                                                >  2:206599/1‑141 (MQ=255)
                             aTGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACtt                                                                                              <  2:6746/141‑1 (MQ=255)
                             aTGGGGTTTGTCTATTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCAcct                                                                                              >  2:89715/1‑139 (MQ=255)
                                   tttGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTAGTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCtt                                                                                        >  1:514614/1‑141 (MQ=255)
                                   tttGGCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCtt                                                                                        >  2:134017/1‑141 (MQ=255)
                                    ttGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCttt                                                                                       >  1:164054/1‑140 (MQ=255)
                                     tGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGt                                                                                      >  1:695017/1‑141 (MQ=255)
                                       tctcTTGGACAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTATGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTAt                                                                                    <  1:58215/141‑1 (MQ=255)
                                           ttGAATAAAGCATAGAAAAAATACTGCAATAAAAATACCACATCGCATAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTtat                                                                                <  2:547779/141‑1 (MQ=255)
                                                tAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTg                                                                           <  1:382655/141‑1 (MQ=255)
                                                  aaGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTACAAGGTTGTGGAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGaa                                                                         >  1:548398/1‑141 (MQ=255)
                                                   aGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAAt                                                                        >  1:603170/1‑141 (MQ=255)
                                                      gTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTaa                                                                     >  2:359165/1‑141 (MQ=255)
                                                      gTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTaa                                                                     >  1:643979/1‑141 (MQ=255)
                                                            aaaaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAAGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAAtattat                                                               >  2:672641/1‑141 (MQ=255)
                                                             aaaTACGTCGATAAAAGTAACACATCTCGTATAAAAACCAGCATCTATATCACTCATGATTTGTTTTATTGTATTAGGATTATTAATGTCATGATTTGGAGTATCACTTATCCGTGGGTTTATAATTGAATTAATATtata                                                              <  2:144707/141‑1 (MQ=255)
                                                              aaTACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGTATTATTAATGCCATTATTTGTATTATCACTTATCCTTGTATTTTTAATTGAATTACTATtatat                                                             >  1:42382/1‑141 (MQ=255)
                                                                    gCAATAAAGGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTGTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTaaa                                                       >  2:361288/1‑141 (MQ=255)
                                                                          aaaTTATCACATCGAGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAtcgt                                                 <  2:371062/141‑1 (MQ=255)
                                                                          aaaGTACCACATCGCGTAAAAAAACCAGCATATATATCATTTATAATATATTTTATTGTATTAGGATTATTAATGCCATTATTAGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAtcgt                                                 <  1:586052/141‑1 (MQ=255)
                                                                           aaGTACCACATCGAGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGAtcgtc                                                <  2:307346/141‑1 (MQ=255)
                                                                                ccACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTATAGATCGTCGTGCt                                           >  1:215340/1‑141 (MQ=255)
                                                                                ccACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCt                                           <  1:89715/141‑1 (MQ=255)
                                                                                     tCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACa                                      <  1:394196/141‑1 (MQ=255)
                                                                                     ccgcgTAAAAATATCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACa                                      <  1:581238/140‑1 (MQ=255)
                                                                                      cgcgTAAAAAAACCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTACACaa                                     >  1:329363/1‑141 (MQ=255)
                                                                                           aaaaaaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACaataatt                                >  1:625650/1‑140 (MQ=255)
                                                                                            aaaaaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGCTTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGc                               >  2:343543/1‑141 (MQ=255)
                                                                                              aaaaCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCAc                             <  2:695017/141‑1 (MQ=255)
                                                                                                 aCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTa                          >  1:689955/1‑141 (MQ=255)
                                                                                                 aCCAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTa                          >  1:298481/1‑141 (MQ=255)
                                                                                                   cAGCATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTaaa                        <  1:277730/141‑1 (MQ=255)
                                                                                                      cATCTATATCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATAGTCGTGCTAAACAATAATGCACTTAAAGtt                     >  2:57220/1‑141 (MQ=255)
                                                                                                           atatCACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGaa                >  1:155166/1‑141 (MQ=255)
                                                                                                              tCACTCATAATATGTTTTATTGTATTAGGATTATTAATGTCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTg             <  2:646974/141‑1 (MQ=255)
                                                                                                               cACTCATAATATGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGa            <  1:735654/141‑1 (MQ=255)
                                                                                                                 tttataataTGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACg          <  2:622696/138‑1 (MQ=255)
                                                                                                                     taataTGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAAt      >  2:24665/1‑141 (MQ=255)
                                                                                                                       ttatGTTTTATTGTATTATTATTATTAATGCCATTATTTGGATTAGCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAAttt    <  1:168064/140‑1 (MQ=255)
                                                                                                                        tatGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTa   >  1:430973/1‑141 (MQ=255)
                                                                                                                         atGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTAc  >  1:557528/1‑141 (MQ=255)
                                                                                                                         atGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAATTTTGAACTGACGAAATTTAc  >  2:40810/1‑141 (MQ=255)
                                                                                                                         atGTTTTATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACGTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTAc  <  2:607708/141‑1 (MQ=255)
                                                                                                                               |                                                                                                                                      
TATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTATCACTTATCCTTGTATTTATAATTGAATTAATATTATATATTAAAGATCGTCGTGCTAAACAATAATGCACTTAAAGTTTTGAACTGACGAAATTTAC  >  USA300TCH1516_ALE/1513803‑1514064

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: