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% 118.5 / NA 37V365L (GTA→CTA) hemYProtoporphyrinogen oxidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (17/20);  total (17/20)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCT  >  USA300TCH1516_ALE/1963637‑1963904
                                                                                                                                   |                                                                                                                                         
cGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCCCg                                                                                                                                  >  1:692323/1‑141 (MQ=255)
cGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACg                                                                                                                                  <  1:735016/141‑1 (MQ=255)
cGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACg                                                                                                                                  <  1:704254/141‑1 (MQ=255)
cGGCAAACGATTGACAATTGTAAATTCAGGAGCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCCATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACg                                                                                                                                  <  1:331337/141‑1 (MQ=255)
      aCGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCaa                                                                                                                            >  1:673104/1‑141 (MQ=255)
        gATTGACAATTGTAAATTCAGGATCACCCTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCaaaa                                                                                                                          <  1:484173/141‑1 (MQ=255)
           tGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAacct                                                                                                                       >  1:716179/1‑141 (MQ=255)
                 ttGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttacctt                                                                                                                 <  1:513524/141‑1 (MQ=255)
                      aaTTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATGCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAAACTTACCTTCCGgt                                                                                                            >  1:787912/1‑141 (MQ=255)
                          cAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACTAAACTACCCATAGGTTTACCTCGCTACGAGAAAATCAAAAACTTACATTGCGGTgtag                                                                                                        >  1:185145/1‑141 (MQ=255)
                               tCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAAt                                                                                                   <  1:479146/141‑1 (MQ=255)
                                      tAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTAGGATCATCAACCACAGTATCACCTGGGTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGGGTAGAAAAGGGACAtt                                                                                            >  1:258331/1‑141 (MQ=255)
                                         ctGTCGTCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCAttttt                                                                                         <  1:707657/140‑1 (MQ=255)
                                         aTGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCAttttt                                                                                         <  1:268784/141‑1 (MQ=255)
                                           gTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTc                                                                                       <  1:448764/141‑1 (MQ=255)
                                           gTCAGCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTc                                                                                       <  1:227346/141‑1 (MQ=255)
                                                        ttAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCCACCAAAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTGTCGATGTCCAAGTAc                                                                          >  1:69224/1‑141 (MQ=255)
                                                          aaaTCTCTTCGTACAATCGATTCTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACAt                                                                        <  1:338045/141‑1 (MQ=255)
                                                           aaTCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATg                                                                       >  1:447900/1‑141 (MQ=255)
                                                                 caCGTACAATCGATACTAATTCATTATCTGCATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTaa                                                                 <  1:499039/140‑1 (MQ=255)
                                                                    gTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGt                                                              >  1:170243/1‑141 (MQ=255)
                                                                      aCAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTGCATGCGGTAATGTCt                                                            >  1:327235/1‑141 (MQ=255)
                                                                         aTCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTa                                                         >  1:145438/1‑141 (MQ=255)
                                                                            gATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCa                                                      >  1:692376/1‑141 (MQ=255)
                                                                                ttcattcattcTCTGTATGATCATCAACCACAGTATCACCTGGTTTCCCTAGATCCGCACGACTCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCg                                                 <  1:400788/138‑1 (MQ=255)
                                                                                gtAATTCATTCTCTGTATTATCATCAACCACAGTATCCCCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCg                                                 <  1:400836/140‑1 (MQ=255)
                                                                                cTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCg                                                 >  1:284880/1‑141 (MQ=255)
                                                                                        ttATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGACACGCACGAATCGAAACCTTACCTTCCGATGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATTCGGTAATGTCTGTATCACTCGTTCTCGCa                                         >  1:323522/1‑141 (MQ=255)
                                                                                         tATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCaa                                        <  1:481664/141‑1 (MQ=255)
                                                                                         tATCTGGATGATCATCAACCACAGTATCACCTGGTTTACCTAGACACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATTGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCaa                                        <  1:38567/141‑1 (MQ=255)
                                                                                                aTGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGaa                                 <  1:745909/141‑1 (MQ=255)
                                                                                                     catcaACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTAGCACTCGTTCTCGCAATTACGAGGCCAg                            >  1:342447/1‑141 (MQ=255)
                                                                                                     cagcaACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTTCATATCTTGCTCGCAATTACGAAGCCAg                            <  1:599728/138‑1 (MQ=255)
                                                                                                       tcaACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTAGCAATTACGAAGCCAGTa                          >  1:641512/1‑141 (MQ=255)
                                                                                                            cacaGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcatc                     <  1:782021/141‑1 (MQ=255)
                                                                                                                    cACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCATTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTAt             <  1:432364/141‑1 (MQ=255)
                                                                                                                               ccTAGATACGCACGAATCAAAACCTTACCTGCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCt  >  1:241728/1‑141 (MQ=255)
                                                                                                                                   |                                                                                                                                         
CGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCT  >  USA300TCH1516_ALE/1963637‑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: