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,7680CG98.1% 175.4 / NA 53V365L (GTA→CTA) hemGprotoporphyrinogen oxidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (29/23);  minor base T (1/0);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCA  >  CP000730/1963632‑1963898
                                                                                                                                        |                                                                                                                                  
cTTTTCGGCAAACGATTGAAAATAAAACAATCACGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATAc                                                                                                                                >  1:285832/1‑141 (MQ=255)
  tttCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGc                                                                                                                              <  1:68163/141‑1 (MQ=255)
      ggCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGa                                                                                                                          >  1:416035/1‑141 (MQ=255)
         aaaCGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATc                                                                                                                       <  1:1195477/141‑1 (MQ=255)
             gATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCaaaa                                                                                                                   >  1:1011381/1‑141 (MQ=255)
                   cAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttac                                                                                                             <  1:383451/141‑1 (MQ=255)
                   cAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttac                                                                                                             <  1:87305/141‑1 (MQ=255)
                   cAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttac                                                                                                             <  1:1009773/141‑1 (MQ=255)
                   cAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttac                                                                                                             <  1:1193947/141‑1 (MQ=255)
                   cAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttac                                                                                                             >  1:963129/1‑141 (MQ=255)
                   cAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttac                                                                                                             <  1:127270/141‑1 (MQ=255)
                      ttGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttacctt                                                                                                          >  1:378216/1‑141 (MQ=255)
                          aaaTTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTCGATACGCACGAATCAAAACCTTACCTTCCgg                                                                                                      >  1:376134/1‑141 (MQ=255)
                          aaaTTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCgg                                                                                                      >  1:183881/1‑141 (MQ=255)
                                       ccTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGc                                                                                         >  1:336044/1‑141 (MQ=255)
                                        cTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGcc                                                                                        >  1:472182/1‑141 (MQ=255)
                                             aaTGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGACCAtttt                                                                                   >  1:983666/1‑141 (MQ=255)
                                              aTGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCAttttt                                                                                  <  1:851866/141‑1 (MQ=255)
                                                  catcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATCCGCACGCATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGa                                                                              <  1:1122721/141‑1 (MQ=255)
                                                    tcacTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTATCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATg                                                                            <  1:176020/141‑1 (MQ=255)
                                                         tGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGCATCAAAACCGGACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCaa                                                                       >  1:682021/1‑141 (MQ=255)
                                                          gACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAg                                                                      >  1:480863/1‑141 (MQ=255)
                                                           aCTTAAATCTCTACGTACAATCGATACTAAGTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTCACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGt                                                                     <  1:704013/141‑1 (MQ=255)
                                                            cTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTa                                                                    <  1:265496/141‑1 (MQ=255)
                                                                aaTCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATg                                                                >  1:100460/1‑141 (MQ=255)
                                                                    tctACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTAAATGCGGt                                                            <  1:751708/141‑1 (MQ=255)
                                                                        cGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATg                                                        >  1:752235/1‑141 (MQ=255)
                                                                         gTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGt                                                       >  1:999005/1‑141 (MQ=255)
                                                                         gTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGt                                                       >  1:150953/1‑141 (MQ=255)
                                                                              aTCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTa                                                  >  1:1185718/1‑141 (MQ=255)
                                                                                 gATACTAATTCATTATCTGTATGATCATCAGCCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCa                                               <  1:298350/141‑1 (MQ=255)
                                                                                 gATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCa                                               >  1:949268/1‑141 (MQ=255)
                                                                                    aCTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTc                                            >  1:1058433/1‑141 (MQ=255)
                                                                                    aCTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATACGGTAATCTCTATATCACTc                                            >  1:1058422/1‑141 (MQ=255)
                                                                                        aTTCCTTATCTGTATGATCATCAACCACAGTATCCCCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCTGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTtc                                        <  1:166649/141‑1 (MQ=255)
                                                                                            cttctCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGc                                    <  1:528751/137‑1 (MQ=255)
                                                                                            aTTATCTGTATGATCATCAACCACAGTACCACCCGGTTTACCTGGATACGCACGAACCAAAACCTTACCTCCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGc                                    <  1:28330/141‑1 (MQ=255)
                                                                                               aTCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAAt                                 >  1:204371/1‑141 (MQ=255)
                                                                                               aTCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGAAATCACCCATGTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAAt                                 >  1:867832/1‑141 (MQ=255)
                                                                                                    tATGATCATCAACCACAGTAGCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTCCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGa                            >  1:937073/1‑141 (MQ=255)
                                                                                                     aTGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGaa                           <  1:76938/141‑1 (MQ=255)
                                                                                                      tGATCATCAACCACAGAAAACCAAGGGGTACCTATATACGCACGAATCAAAACCTTACCTTCCGGTGAAGTAAATGGCCAAGTTATCGCGGTCCAAGTACATGCGGTAAAGTCTGTATCACTCGTTCTCGCAATCACGAAg                          >  1:286487/1‑141 (MQ=255)
                                                                                                        atcatcAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCTGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGcc                        <  1:766678/141‑1 (MQ=255)
                                                                                                                ccACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcat                >  1:859766/1‑141 (MQ=255)
                                                                                                                ccACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcat                >  1:674261/1‑141 (MQ=255)
                                                                                                                ccACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcac                >  1:1184728/1‑140 (MQ=255)
                                                                                                                 cacaGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACCTGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAGGCCGGTGCcatc               >  1:246959/1‑141 (MQ=255)
                                                                                                                  acaGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTGCTCCCAATTACGAAGCCAGTACcatct              >  1:351978/1‑140 (MQ=255)
                                                                                                                    aGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATCGCCATTTTTTCGATGTCCAAGTACGTGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATg            <  1:772813/141‑1 (MQ=255)
                                                                                                                    aGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTCCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATg            <  1:160967/141‑1 (MQ=255)
                                                                                                                          cccTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTAtt      <  1:938236/140‑1 (MQ=255)
                                                                                                                             cggTTTACCTAGATACGCCCGAATCAAACCCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTc   <  1:32646/140‑1 (MQ=255)
                                                                                                                              ggTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGCTGTCCAAGTACATGCGGTAACGTCAGTCTCCTGCAGTTTCGCGATTCGGAAGCCAGTACCATCATGGGTATTTTCa  >  1:563509/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                  
CTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCA  >  CP000730/1963632‑1963898

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: