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

AAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCTT  >  USA300TCH1516_ALE/1963641‑1963905
                                                                                                                               |                                                                                                                                         
aaaCTATTGACAATTGTAAATTCAGGATCACCTTTAAATGCCATCACTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGCATGATCAGCAACCACACTATCACCTGGTTTACCTAGATACGCACGAATc                                                                                                                              <  2:482285/141‑1 (MQ=255)
   cGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCaaa                                                                                                                           >  1:507619/1‑141 (MQ=255)
      ttGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGACTCAAAccc                                                                                                                        >  2:343766/1‑141 (MQ=255)
       tGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAacct                                                                                                                       >  1:175153/1‑141 (MQ=255)
        gACAATTTTAACTTCAGGATCACCTTTAAATTTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACCGTATCACCTGGTTTACCTGGATACGCACGAATCAAAacctt                                                                                                                      <  2:330181/141‑1 (MQ=255)
            aTTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTTTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttacct                                                                                                                  <  2:112762/141‑1 (MQ=255)
            aTTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttacct                                                                                                                  <  2:112038/141‑1 (MQ=255)
            aTTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttacct                                                                                                                  <  2:622135/141‑1 (MQ=255)
            aTTGTAAATTCAGGATCACCTTTAAATGTCAGCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttacct                                                                                                                  <  1:244258/141‑1 (MQ=255)
                   aTTCACGAACACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGgtg                                                                                                           <  2:283854/141‑1 (MQ=255)
                             aCCTTTAAATGTCATCATTTGACTTAAATCTCTTCGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGATTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTCGTAAATgg                                                                                                 >  2:440391/1‑141 (MQ=255)
                             aCCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATgg                                                                                                 <  2:265606/141‑1 (MQ=255)
                              ccTTTAAATGTCATCATTTGAATTAAATCTCTACGTACAATCGATACTAATTCATAATCTGTATGATCATCAACCACAGTATCACCTGGGTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGc                                                                                                >  2:630612/1‑141 (MQ=255)
                                 ttAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCAt                                                                                             >  2:81460/1‑141 (MQ=255)
                                 ttAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCAt                                                                                             >  1:75350/1‑141 (MQ=255)
                                 ttAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCAt                                                                                             >  2:582659/1‑141 (MQ=255)
                                        tcatcaTTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCg                                                                                      >  1:628993/1‑141 (MQ=255)
                                              tttGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTcc                                                                                >  1:327632/1‑141 (MQ=255)
                                                   cTTAAACCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTa                                                                           <  2:266912/141‑1 (MQ=255)
                                                     tAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACa                                                                         <  1:74537/141‑1 (MQ=255)
                                                     tAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACCCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACa                                                                         <  1:197084/141‑1 (MQ=255)
                                                     tAAATCTCTACGTACAATCCATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACa                                                                         <  2:419526/141‑1 (MQ=255)
                                                       aaTCTCTACGGACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATg                                                                       <  1:647773/141‑1 (MQ=255)
                                                            ctACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTa                                                                  <  2:70736/141‑1 (MQ=255)
                                                            ctACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAAACTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCACGTACATGCGGTa                                                                  >  2:416774/1‑141 (MQ=255)
                                                             tACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTaa                                                                 >  1:75426/1‑141 (MQ=255)
                                                                gTACAACCGATTCTAATTTATTATCTGCATGATCATCAACCACAGGATCACCTTGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGt                                                              <  1:300543/141‑1 (MQ=255)
                                                                cttcAATCGATACTAATTCCTTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGt                                                              <  2:250750/138‑1 (MQ=255)
                                                                      tCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCCCCTAGATTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTAt                                                        >  1:456955/1‑141 (MQ=255)
                                                                      tCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTAt                                                        >  1:53920/1‑141 (MQ=255)
                                                                            cTAATTCCTTCTCTGTATGATCATCAACCACAGTAACACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCg                                                  <  1:504108/141‑1 (MQ=255)
                                                                                  cATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTTTTTGTCTTATTACATGCGGTAATGTCTGTATCACTAGTTCTCg                                            >  2:428600/1‑141 (MQ=255)
                                                                                     tATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCaa                                         <  2:507619/141‑1 (MQ=255)
                                                                                      aTCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATAGCCAAGTACATGCGGTAATGTCTGTATCACTCGTtctctccat                                        >  2:310209/1‑136 (MQ=255)
                                                                                       tCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAAtt                                       <  2:137360/141‑1 (MQ=255)
                                                                                       ccTGTACGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCCCTAATGTCTGTATCACTCGTTCTCGCAAtt                                       <  2:369207/140‑1 (MQ=255)
                                                                                         tGTATGATCATCAACCACAGGATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTAc                                     <  1:119518/141‑1 (MQ=255)
                                                                                         tGTATGATCAACAACCACATTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTAc                                     <  2:158533/141‑1 (MQ=255)
                                                                                             ggagcttcaACCACAGTATCCCCTGGTTTGCCTAGATACGCATGAATCAAAACCTTACCGTCCGTTGTAGTAAATGGCCATCTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAg                                 <  2:483526/135‑1 (MQ=255)
                                                                                               atcatcAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGcc                               >  1:655986/1‑141 (MQ=255)
                                                                                                tcatcaACCACAGTATCACCTGGTTTACCTGGATACGCACGAATGAAAACCTTACCTTCCGGTGTAGTAACTGGCCATTTTTGCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCa                              <  2:434612/141‑1 (MQ=255)
                                                                                                tcagcaACCACGGTATCACATGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCa                              <  2:162965/141‑1 (MQ=255)
                                                                                                  atcaACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGGTCTCGCAATTACGAAAccgtt                            >  2:165224/1‑138 (MQ=255)
                                                                                                    ccaCCACAGTATCCCCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTAc                          <  2:304918/139‑1 (MQ=255)
                                                                                                    caACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTAc                          >  2:103465/1‑141 (MQ=255)
                                                                                                     aaCCCCAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTAGATGAACAACCACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACTAAGCCAGTAcc                         >  1:453971/1‑141 (MQ=255)
                                                                                                        cacaGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGTTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcatc                      >  2:272783/1‑141 (MQ=255)
                                                                                                         ccaGTGTCCCGTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTGCCGGTGTAGTAAATGGCCATTTTTGCGATGTCCAAGTGCATGCGGGAATGTCGGTATCACTCGTGCTCGCAATAACGATGCCAGTAAcataa                     >  2:408009/2‑139 (MQ=255)
                                                                                                         acaGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCCATTACGCAACCACTACcagct                     >  2:601659/1‑138 (MQ=255)
                                                                                                           aGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATg                   <  2:374058/141‑1 (MQ=255)
                                                                                                           aGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATg                   <  1:552213/141‑1 (MQ=255)
                                                                                                             tATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATggg                 >  1:120903/1‑141 (MQ=255)
                                                                                                              aTCACCTGGTTTACCTAGATAAGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGAGGTCCAAGTACATGCGGTAATGCCTGTATAACTCGTTCTCGCAATTACGAAGCCAGTACCATCAAggta                >  2:5306/1‑139 (MQ=255)
                                                                                                                cACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTAt              >  1:482503/1‑141 (MQ=255)
                                                                                                                cACCCGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTAt              <  2:665981/141‑1 (MQ=255)
                                                                                                                    tGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGTTGTGGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTc          >  1:628981/1‑141 (MQ=255)
                                                                                                                    tGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAACTGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTCCGAAGCCAGTGACATCATGGGTATATTc          >  1:349474/1‑141 (MQ=255)
                                                                                                                    tGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTc          >  1:173589/1‑141 (MQ=255)
                                                                                                                    tGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGACCGAGTCATTCGTTATCGCAATTACGGAGCCAGTACCATCATGGGTATTTTc          >  1:651764/1‑141 (MQ=255)
                                                                                                                      gTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCaa        >  1:547342/1‑141 (MQ=255)
                                                                                                                        ttcccTAGATACGCACGAATCAAAACCTTACCTTCCGTTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCGCTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATg      <  2:376683/138‑1 (MQ=255)
                                                                                                                        ttACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATg      >  1:196998/1‑141 (MQ=255)
                                                                                                                            cTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCtt  <  1:440391/141‑1 (MQ=255)
                                                                                                                               |                                                                                                                                         
AAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCTT  >  USA300TCH1516_ALE/1963641‑1963905

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