Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,963,768 C→G 100% V365L (GTA→CTA)  hemY ← Protoporphyrinogen oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,963,7680CG100.0% 162.5 / NA 50V365L (GTA→CTA) hemYProtoporphyrinogen oxidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (24/26);  total (24/26)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCT  >  USA300TCH1516_ALE/1963631‑1963904
                                                                                                                                         |                                                                                                                                        
aCTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATa                                                                                                                                       >  1:158491/1‑141 (MQ=255)
aCTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATa                                                                                                                                       <  1:662630/141‑1 (MQ=255)
    ttCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCa                                                                                                                                   <  1:412367/141‑1 (MQ=255)
      cGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACg                                                                                                                                 <  1:660480/141‑1 (MQ=255)
           aaCGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGCCTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATTAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCa                                                                                                                            >  1:715155/1‑141 (MQ=255)
            aCGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCaa                                                                                                                           <  1:702494/141‑1 (MQ=255)
                   aCAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGAGACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGCATCAAAccctta                                                                                                                    >  1:127392/1‑141 (MQ=255)
                   aCAATTGTAAATTCAGGAGCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAacctta                                                                                                                    <  1:322816/141‑1 (MQ=255)
                       ttGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttacctt                                                                                                                >  1:596951/1‑141 (MQ=255)
                           aaaTTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCgg                                                                                                            >  1:162141/1‑141 (MQ=255)
                            aaTTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCggg                                                                                                           >  1:230836/1‑140 (MQ=255)
                               tCAGGATCACCTTTAAATGCCATCATTTGACTTAAATCTCTACGCACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTgta                                                                                                        >  1:405490/1‑141 (MQ=255)
                                   gATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGaaa                                                                                                    >  1:275045/1‑141 (MQ=255)
                                        ccTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGc                                                                                               <  1:26387/141‑1 (MQ=255)
                                        ccTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGc                                                                                               >  1:448814/1‑141 (MQ=255)
                                        ccTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATCCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGc                                                                                               <  1:738962/141‑1 (MQ=255)
                                           ttAAATGTCATCATTTGACTTAAATCTCTATGTCCACTCTACACTAACTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCAt                                                                                            <  1:590276/141‑1 (MQ=255)
                                             aaaTGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCAttt                                                                                          >  1:235985/1‑141 (MQ=255)
                                                 gTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTc                                                                                      <  1:122977/141‑1 (MQ=255)
                                                  tcatcaTTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCg                                                                                     >  1:584151/1‑141 (MQ=255)
                                                  tcatcaCTTGACCTAAACCACTACGTACAATCGATAGTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAACTGGCCATTTTTTCg                                                                                     <  1:486714/141‑1 (MQ=255)
                                                    atcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGAt                                                                                   >  1:575932/1‑141 (MQ=255)
                                                               tAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACa                                                                        <  1:343020/141‑1 (MQ=255)
                                                               tAAATCTCTACGTACAAGCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTCACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACa                                                                        <  1:407769/141‑1 (MQ=255)
                                                                      ctACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTa                                                                 >  1:409999/1‑141 (MQ=255)
                                                                            tcAATCGATCTTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTTCCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCt                                                           <  1:194693/140‑1 (MQ=255)
                                                                               aTCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTa                                                        >  1:421486/1‑141 (MQ=255)
                                                                                  gATACTAATTCATTATCTGTATGATCATCAACCACAGTACCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTTCATGCGGTTATGTCTGTATCa                                                     <  1:395395/141‑1 (MQ=255)
                                                                                       tAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGTGGTAATGTCTGTATCACTCGt                                                <  1:395244/141‑1 (MQ=255)
                                                                                         aTTCATTATCTGTATGCTCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTtc                                              <  1:28836/141‑1 (MQ=255)
                                                                                                aTCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAAt                                       >  1:318242/1‑141 (MQ=255)
                                                                                                aTCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAAt                                       >  1:318207/1‑141 (MQ=255)
                                                                                                 tcagTATGATCATCAACCACAGTATCACCTTGTTTACCTAGATCCGCACGAATCAAAACCTTTCCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTAAATTCCTTTTTTTTTGTATCACTCGTTCTCGCAAtt                                      <  1:368333/138‑1 (MQ=255)
                                                                                                  cTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTa                                     <  1:61408/141‑1 (MQ=255)
                                                                                                   tGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTAAATGCGGTAATGTCTGTATCACTTGTTCTCGCAATTAc                                    >  1:45718/1‑141 (MQ=255)
                                                                                                    gTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACg                                   <  1:134000/141‑1 (MQ=255)
                                                                                                      aTGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGaa                                 >  1:647231/1‑141 (MQ=255)
                                                                                                           catcaACCACAGTATCACCTGGTTTACCTAGATACGCACGAACCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATCTCCAACTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAg                            <  1:254453/141‑1 (MQ=255)
                                                                                                             tcaACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGCAGCCAGTa                          >  1:137980/1‑141 (MQ=255)
                                                                                                             tcaACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTa                          >  1:375761/1‑141 (MQ=255)
                                                                                                                       taccACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATggg                <  1:457496/138‑1 (MQ=255)
                                                                                                                       tATCACCTGGTTCACCTAGATACGCCCTAATCAAAACCTTACCTTTCGGTGTAGTAAATGGCCTTTTTTTCGTTTTCCAAGTACAGGGGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATggg                <  1:756062/141‑1 (MQ=255)
                                                                                                                       gaTCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATggg                <  1:180869/140‑1 (MQ=255)
                                                                                                                        aTCACCTGGTTTACCTAGATACGCACGAATCAAATCCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTATGAAGTTAGTACCATCATGGGt               <  1:563932/141‑1 (MQ=255)
                                                                                                                                 tttACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAAt      <  1:624395/141‑1 (MQ=255)
                                                                                                                                  ttACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCaaag     >  1:112131/1‑139 (MQ=255)
                                                                                                                                  ttACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATg     <  1:37000/141‑1 (MQ=255)
                                                                                                                                    aCCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTCTTTTCAATGTc   >  1:279979/1‑141 (MQ=255)
                                                                                                                                     ccTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCt  >  1:455710/1‑141 (MQ=255)
                                                                                                                                     ccTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCt  >  1:768060/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                        
ACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCT  >  USA300TCH1516_ALE/1963631‑1963904

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: