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,7680CG83.3% 101.2 / 15.1 36V365L (GTA→CTA) hemGprotoporphyrinogen oxidase
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base G (14/16);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

TACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTC  >  CP000730/1963630‑1963903
                                                                                                                                          |                                                                                                                                       
tACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACAt                                                                                                                                       <  1:330133/141‑1 (MQ=255)
tACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACAt                                                                                                                                       <  1:330257/141‑1 (MQ=255)
       cGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACg                                                                                                                                <  1:268680/141‑1 (MQ=255)
        ggCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGa                                                                                                                               <  1:153063/141‑1 (MQ=255)
         gCAAACGATTGACAATTGTAAATTCAGGATCACCGTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTCATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGATTACCTAGATACGCACGaa                                                                                                                              >  2:422515/1‑141 (MQ=255)
          cAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAAt                                                                                                                             >  2:216581/1‑141 (MQ=255)
          cAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAAt                                                                                                                             <  2:281878/141‑1 (MQ=255)
             aCGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATCATCTGTATGATCATCAACCACAGTATCACCTGCTTTACCTAGATACGCACGAATCaa                                                                                                                          >  2:327071/1‑141 (MQ=255)
                 ttGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAacc                                                                                                                      >  2:114314/1‑141 (MQ=255)
                    aCAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAacctta                                                                                                                   >  2:100029/1‑141 (MQ=255)
                    aCAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAacctta                                                                                                                   >  2:273996/1‑141 (MQ=255)
                               ttCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGgtgt                                                                                                        <  2:362377/141‑1 (MQ=255)
                                 cAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTgtag                                                                                                      >  2:151453/1‑141 (MQ=255)
                                 cAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTgtag                                                                                                      >  2:142220/1‑141 (MQ=255)
                                               aaTGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCAtttt                                                                                        >  2:221395/1‑141 (MQ=255)
                                                acgTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCAttttt                                                                                       <  2:414054/139‑1 (MQ=255)
                                                aTGTCATCATTTGACTTAAATCTCTACGTACAGTCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCAttttt                                                                                       <  2:196935/141‑1 (MQ=255)
                                                  gTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTc                                                                                     <  1:100029/141‑1 (MQ=255)
                                                    catcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGa                                                                                   >  1:82971/1‑141 (MQ=255)
                                                           tGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCaa                                                                            <  1:237207/141‑1 (MQ=255)
                                                                  aaTCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATg                                                                     <  2:6241/141‑1 (MQ=255)
                                                                       ctACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTa                                                                <  1:273996/141‑1 (MQ=255)
                                                                          cGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATg                                                             <  2:184349/141‑1 (MQ=255)
                                                                          cGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATg                                                             <  1:337428/141‑1 (MQ=255)
                                                                            tACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTc                                                           <  1:87083/141‑1 (MQ=255)
                                                                                     tACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACt                                                  >  2:154568/1‑141 (MQ=255)
                                                                                      aCTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCCTTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTc                                                 >  1:38685/1‑141 (MQ=255)
                                                                                          aTTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTtc                                             <  2:210195/141‑1 (MQ=255)
                                                                                                  tCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAAtt                                     <  1:116630/141‑1 (MQ=255)
                                                                                                  tCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAAtt                                     <  2:342144/141‑1 (MQ=255)
                                                                                                    tGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTAc                                   >  2:348388/1‑141 (MQ=255)
                                                                                                                      aGTATCACCTGGTTTACCTAGATAAGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATg                 >  1:351286/1‑141 (MQ=255)
                                                                                                                        tATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATCTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATggg               <  2:400833/141‑1 (MQ=255)
                                                                                                                           cACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTAt            <  1:225021/141‑1 (MQ=255)
                                                                                                                           cACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTAt            >  2:30440/1‑141 (MQ=255)
                                                                                                                                     aCCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTc  >  1:183652/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                       
TACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTC  >  CP000730/1963630‑1963903

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: