Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,963,768 C→G 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% 102.6 / NA 31V365L (GTA→CTA) hemYProtoporphyrinogen oxidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (15/16);  total (15/16)

TTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTC  >  USA300TCH1516_ALE/1963635‑1963903
                                                                                                                                     |                                                                                                                                       
ttCGGCAAACGATTGACAAGTGTAACTTCAGGACCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACTATTGATACTAATTCATTATCTGTATGATCATCAACCACTGTATCACCTGGTTCACCTAGATACGCa                                                                                                                                  <  1:181160/141‑1 (MQ=255)
 tCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTTTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCAc                                                                                                                                 <  1:463723/141‑1 (MQ=255)
                cAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttac                                                                                                                  >  1:44199/1‑141 (MQ=255)
                    tGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTc                                                                                                              <  1:24087/141‑1 (MQ=255)
                            cAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGACACGCACGAATCAAAACCTTACCTTCCGGTgtag                                                                                                      <  1:445883/141‑1 (MQ=255)
                             aGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTgtagt                                                                                                     >  1:476686/1‑141 (MQ=255)
                                aTCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCGAAACCTTACTTTCCGGTGTAGTaaa                                                                                                  <  1:251177/141‑1 (MQ=255)
                                    ccTTTAAATGTCATCACTTCACTTAAATCTCTACGTACAATCGATACTAATTCATGATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGc                                                                                              <  1:201486/141‑1 (MQ=255)
                                     cTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTATTTCATTATTTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGAACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGcc                                                                                             <  1:30098/141‑1 (MQ=255)
                                     cTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGcc                                                                                             <  1:213431/141‑1 (MQ=255)
                                     cTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGcc                                                                                             <  1:277165/141‑1 (MQ=255)
                                            tGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCAtttttt                                                                                      >  1:269341/1‑141 (MQ=255)
                                              tcatcaTTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCg                                                                                    <  1:421034/141‑1 (MQ=255)
                                               catcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTGGa                                                                                   >  1:387015/1‑141 (MQ=255)
                                                      tGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCaa                                                                            >  1:175813/1‑141 (MQ=255)
                                                          ttAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTAc                                                                        <  1:7482/141‑1 (MQ=255)
                                                                ctctACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCgg                                                                  >  1:245987/1‑141 (MQ=255)
                                                                 tctACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAACCCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGt                                                                 <  1:476142/141‑1 (MQ=255)
                                                                     cGTACAATCGATACTTATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTGCCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGAGATCCAAGTACATGCGGTAATg                                                             >  1:370368/1‑141 (MQ=255)
                                                                      gTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGt                                                            >  1:427256/1‑141 (MQ=255)
                                                                        aCAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCt                                                          >  1:2470/1‑141 (MQ=255)
                                                                           aTCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTa                                                       >  1:72109/1‑141 (MQ=255)
                                                                               aTACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTCCATGCGGTAATTTCTGTATCAc                                                   >  1:486200/1‑141 (MQ=255)
                                                                                tattAATCCATTATCTGTATGATCATCAACCACAGCATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACt                                                  <  1:177466/138‑1 (MQ=255)
                                                                                       tCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCGGTATCACTCGTtctc                                           >  1:164240/1‑141 (MQ=255)
                                                                                                             ccACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcat                     >  1:61454/1‑141 (MQ=255)
                                                                                                               acaGTATCCCCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcatca                   <  1:108988/141‑1 (MQ=255)
                                                                                                               acaGTATCACCTTGTTTACCTAGATACTCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcatca                   <  1:484893/141‑1 (MQ=255)
                                                                                                                     tCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTa             >  1:74694/1‑141 (MQ=255)
                                                                                                                            gTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCaa      <  1:495784/141‑1 (MQ=255)
                                                                                                                                aCCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTc  >  1:271923/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                                       
TTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTC  >  USA300TCH1516_ALE/1963635‑1963903

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: