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% 111.5 / NA 29V365L (GTA→CTA) hemGprotoporphyrinogen oxidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (18/11);  total (18/11)

GCATACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCTTTTTCATC  >  CP000730/1963627‑1963912
                                                                                                                                             |                                                                                                                                                
gCATACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATATGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCa                                                                                                                                           <  1:180477/149‑1 (MQ=255)
    aCTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGaa                                                                                                                                       >  1:11957/1‑149 (MQ=255)
      ttttCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATc                                                                                                                                     >  1:221734/1‑149 (MQ=255)
         tCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCaaa                                                                                                                                  >  1:149995/1‑149 (MQ=255)
           ggCAAACGATTGACAATTGTAAATTCAGGATCACCTTTGAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAac                                                                                                                                >  2:306161/1‑149 (MQ=255)
           ggCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAac                                                                                                                                >  2:332407/1‑149 (MQ=255)
                     tGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTcc                                                                                                                      <  1:61799/149‑1 (MQ=255)
                           ttGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTgta                                                                                                                <  2:6432/149‑1 (MQ=255)
                                            ccTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCAtttttt                                                                                               >  1:97020/1‑149 (MQ=255)
                                              tttAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAACTGGCCATTTTTTCg                                                                                             >  1:230948/1‑149 (MQ=255)
                                                  aaTGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACGGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGt                                                                                         >  1:208521/1‑149 (MQ=255)
                                                    tGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTcc                                                                                       <  2:130381/149‑1 (MQ=255)
                                                     gTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCa                                                                                      >  1:139347/1‑149 (MQ=255)
                                                      tcatcaTTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCaa                                                                                     <  1:34578/149‑1 (MQ=255)
                                                      tcatcaTTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCaa                                                                                     <  1:38811/149‑1 (MQ=255)
                                                       catcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAg                                                                                    <  1:156860/149‑1 (MQ=255)
                                                                aCTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCgg                                                                           <  2:22109/149‑1 (MQ=255)
                                                                              gTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATcc                                                             >  1:275135/1‑148 (MQ=255)
                                                                                          cTAATTCCTAATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCa                                                 <  1:120046/149‑1 (MQ=255)
                                                                                            aaTTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAAt                                               >  2:189750/1‑149 (MQ=255)
                                                                                            aaTTCATCATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAAt                                               >  2:131724/1‑149 (MQ=255)
                                                                                               tCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTAc                                            >  1:72832/1‑149 (MQ=255)
                                                                                                     tCTGTATGATCATGAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGTTGTCTAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTNCGAAGcc                                      >  2:255531/1‑149 (MQ=255)
                                                                                                                atcaACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcatcat                           <  2:221499/149‑1 (MQ=255)
                                                                                                                  caACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATgg                         >  2:284867/1‑149 (MQ=255)
                                                                                                                              cACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATg             >  1:247318/1‑149 (MQ=255)
                                                                                                                                    gTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGACTGTATCAATCGTTATAGAACTTAAGAAGCCAGTACCATCATGGGTATTTTCAATGTCtttt       >  1:31223/1‑149 (MQ=255)
                                                                                                                                        aCCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCTTTTtcat   <  2:55788/149‑1 (MQ=255)
                                                                                                                                         ccTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTTTTTTTTTTTTTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCTTTTtcatc  >  2:148959/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                
GCATACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCTTTTTCATC  >  CP000730/1963627‑1963912

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