Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,851,260 C→A 100% A115S (GCG→TCG)  acsA_1 ← Acetyl‑coenzyme A synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,851,2600CA100.0% 74.3 / NA 28A115S (GCG→TCG) acsA_1Acetyl‑coenzyme A synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/12);  total (16/12)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AATGCCTTATTAGTAATTAACACTTTAGCTTCACTGTTCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGCAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGATAACCGTTGCATATCTTTATAAGTATACGATTCTTTTC  >  USA300TCH1516_ALE/1851136‑1851394
                                                                                                                            |                                                                                                                                      
aaTGCCTTATTAGTAATTAACACTTTAGCTTCACTGTTCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGgt                                                                                                                        <  2:204393/141‑1 (MQ=255)
aaTGCCTTAGTAGTAATTAACACTTTAGCTTCACTGTTCTCTAATCTATCCGCAACTGCCGTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGgt                                                                                                                        <  2:378109/141‑1 (MQ=255)
 aTGCCTGATTAGTAATTAACACTTTAGCTTCACTGTTCTCTAATCTATCCTCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGgtt                                                                                                                       >  2:291302/1‑140 (MQ=255)
    ccTAATTAGTAATTAACACTTTAGCTTCACTGTTCTCTAATCTATCCGCAACTGCATTTTCCATAAATGCTTCAAATAACGGCCCTACAATTGCACCAATGATTAAATGACCTAACAACGAAAAATAGAGTTCAGGTGTAc                                                                                                                    <  2:45984/141‑1 (MQ=255)
     cTTATTAGTAATTAACACTTTAGCTTCACTGTTCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTAcg                                                                                                                   <  1:394705/141‑1 (MQ=255)
     cTTATTAGTAACTAACACTTTAGCTTCACTGTTCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTAcg                                                                                                                   <  2:225072/141‑1 (MQ=255)
               aTTAACACTTTAGCTTCACTGTTCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTATCAACGAAAAATATAGCCCAGGTGCACGCAACataaaa                                                                                                         >  2:31139/1‑140 (MQ=255)
                           gCTTCACTGTTCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTg                                                                                             >  1:153242/1‑141 (MQ=255)
                            cTTCACTGTTCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAAAGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACCTAAATATAAATACACTGt                                                                                            >  2:161826/1‑141 (MQ=255)
                             ttCACTGTTCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTc                                                                                           <  2:346467/141‑1 (MQ=255)
                                    ttCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATTCTCTGTCACCTTTg                                                                                    <  2:126344/141‑1 (MQ=255)
                                      ctctAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTTTc                                                                                  >  2:177353/1‑141 (MQ=255)
                                       tctAATCTATCCGCAACTGCCATTTCCATAAATGCTTCAAATAACGGCCCAACAACTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTAAGGTGTACGCCACAAAAATATAAATACTCTGCCACCttttcca                                                                                 >  2:123329/1‑137 (MQ=255)
                                           aTCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCCCCTTTGTCAACtt                                                                             >  1:23119/1‑141 (MQ=255)
                                           aTCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACtt                                                                             >  1:26603/1‑141 (MQ=255)
                                                tCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGCCATAAATATAAATACTCTGTCACCTTTGTCAACTTCTTCa                                                                        >  2:294820/1‑141 (MQ=255)
                                                         gCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGAc                                                               >  2:202298/1‑141 (MQ=255)
                                                                ccATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACAt                                                        >  2:251230/1‑141 (MQ=255)
                                                                            cAAATAACGGCCCAACAATTGCACGAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATTCTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTttat                                            <  2:251232/141‑1 (MQ=255)
                                                                                tAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCTACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGa                                        <  1:364013/141‑1 (MQ=255)
                                                                                aaaCGGCCCAACAATTGCACCAATTGTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGa                                        <  2:73096/140‑1 (MQ=255)
                                                                                     gCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGATAAcc                                   >  1:61272/1‑141 (MQ=255)
                                                                                       ccAACAATTGCACCAATTTTTAAAACACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGATACCCGt                                 >  1:4498/1‑141 (MQ=255)
                                                                                                gCACCAATTTCTTAAACACCTAACAACAAAAAACATCGTTCACGTGTACGCGACATAAATATAAATACTCTGTCCATTTTGGCATCTTGTGCATTTTCAGACAAAACATTCGCTTCTTTTTTAGATAACCGTTGCATATCt                        <  2:258044/141‑1 (MQ=255)
                                                                                                             aaaCACCTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAAGATAAATACACTGTCACCTTGGTCAACTTCTGCATGTTCAGACAAAACATTCGCTACATTATTAGATAACGGTTGCATATCTGTATAAGTATACg           >  2:296528/1‑141 (MQ=255)
                                                                                                               acacCTAACAACGAAAAATATAGTTCAGGTGTAAGCGAAGTAAATATAAATACACTGTCACATTTGTTTAATACTAAATGTAAAGACACAACAGGCGCCGCGCTATTAGATAACCCTTTCATATCTTTATAAGTATACGAt         >  2:190160/1‑141 (MQ=255)
                                                                                                                   cTAACAACGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGATAACCGTTGCATATCTTTATAAGTATACGATTCtt     >  1:162387/1‑141 (MQ=255)
                                                                                                                      acaacGAAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGATAACCGTTGCATATCTTTATAAGTATACGATTCTTTTc  <  1:16272/141‑1 (MQ=255)
                                                                                                                            |                                                                                                                                      
AATGCCTTATTAGTAATTAACACTTTAGCTTCACTGTTCTCTAATCTATCCGCAACTGCCTTTTCCATAAATGCTTCAAATAACGGCCCAACAATTGCACCAATTTTTAAAACACCTAACAACGCAAAATATAGTTCAGGTGTACGCGACATAAATATAAATACTCTGTCACCTTTGTCAACTTCTGCATGTTCAGACAAAACATTCGCTGCTTTATTAGATAACCGTTGCATATCTTTATAAGTATACGATTCTTTTC  >  USA300TCH1516_ALE/1851136‑1851394

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