Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 1,963,768 C→G V365L (GTA→CTA)  hemG ← protoporphyrinogen oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,963,7680CG100.0% 102.2 / NA 31V365L (GTA→CTA) hemGprotoporphyrinogen oxidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (15/16);  total (15/16)

TTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTC  >  CP000730/1963635‑1963903
                                                                                                                                     |                                                                                                                                       
ttCGGCAAACGATTGACAAGTGTAACTTCAGGACCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACTATTGATACTAATTCATTATCTGTATGATCATCAACCACTGTATCACCTGGTTCACCTAGATACGCa                                                                                                                                  <  1:179804/141‑1 (MQ=255)
 tCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTTTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCAc                                                                                                                                 <  1:460684/141‑1 (MQ=255)
                cAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttac                                                                                                                  >  1:43880/1‑141 (MQ=255)
                    tGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTc                                                                                                              <  1:23902/141‑1 (MQ=255)
                            cAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGACACGCACGAATCAAAACCTTACCTTCCGGTgtag                                                                                                      <  1:442899/141‑1 (MQ=255)
                             aGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTgtagt                                                                                                     >  1:473557/1‑141 (MQ=255)
                                aTCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCGAAACCTTACTTTCCGGTGTAGTaaa                                                                                                  <  1:249296/141‑1 (MQ=255)
                                    ccTTTAAATGTCATCACTTCACTTAAATCTCTACGTACAATCGATACTAATTCATGATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGc                                                                                              <  1:200013/141‑1 (MQ=255)
                                     cTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTATTTCATTATTTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGAACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGcc                                                                                             <  1:29877/141‑1 (MQ=255)
                                     cTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGcc                                                                                             <  1:211881/141‑1 (MQ=255)
                                     cTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGcc                                                                                             <  1:275134/141‑1 (MQ=255)
                                            tGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCAtttttt                                                                                      >  1:267370/1‑141 (MQ=255)
                                              tcatcaTTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCg                                                                                    <  1:418163/141‑1 (MQ=255)
                                               catcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTGGa                                                                                   >  1:384264/1‑141 (MQ=255)
                                                      tGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCaa                                                                            >  1:174507/1‑141 (MQ=255)
                                                       tacttacacttctACGTACAAGCGATACTAATCCATGATCTGTATGACCATCAACCACAGTATCCCCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAg                                                                           <  1:169154/131‑1 (MQ=255)
                                                          ttAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTAc                                                                        <  1:7429/141‑1 (MQ=255)
                                                                ctctACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCgg                                                                  >  1:244130/1‑141 (MQ=255)
                                                                 tctACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAACCCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGt                                                                 <  1:473014/141‑1 (MQ=255)
                                                                     cGTACAATCGATACTTATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTGCCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGAGATCCAAGTACATGCGGTAATg                                                             >  1:367701/1‑141 (MQ=255)
                                                                      gTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGt                                                            >  1:424360/1‑141 (MQ=255)
                                                                        aCAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCt                                                          >  1:2442/1‑141 (MQ=255)
                                                                           aTCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTa                                                       >  1:71611/1‑141 (MQ=255)
                                                                               aTACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTCCATGCGGTAATTTCTGTATCAc                                                   >  1:482990/1‑141 (MQ=255)
                                                                                       tCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCGGTATCACTCGTtctc                                           >  1:163010/1‑141 (MQ=255)
                                                                                                             ccACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcat                     >  1:61028/1‑141 (MQ=255)
                                                                                                               acaGTATCCCCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcatca                   <  1:108221/141‑1 (MQ=255)
                                                                                                               acaGTATCACCTTGTTTACCTAGATACTCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcatca                   <  1:481686/141‑1 (MQ=255)
                                                                                                                     tCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTa             >  1:74190/1‑141 (MQ=255)
                                                                                                                            gTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCaa      <  1:492471/141‑1 (MQ=255)
                                                                                                                                aCCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTc  >  1:269946/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                                       
TTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTC  >  CP000730/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: