Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,584,576 G→A 19.0% P21S (CCA→TCA)  hemN ← coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,584,5760GA19.0% 81.9 / 13.4 37P21S (CCA→TCA) hemNcoproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent
Reads supporting (aligned to +/- strand):  ref base G (14/15);  new base A (4/3);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

TAAAAACGCCTGTTCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTACAGACATGCTCACTCCT  >  W3110S.gb/3584511‑3584646
                                                                 |                                                                      
tAAAAACGCCTGTTCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCgg                                                                   <  1:5091943/71‑1 (MQ=255)
  aaaaCGCCTGTTCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGACCAGGaa                                                                 >  1:1447373/1‑71 (MQ=255)
        ccTGTTCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTta                                                           <  1:4118479/71‑1 (MQ=255)
            ttCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCgg                                                                   <  1:2294190/59‑1 (MQ=255)
              cGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATAttt                                                      <  1:1504330/70‑1 (MQ=255)
              cGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTc                                                     >  1:1135847/1‑71 (MQ=255)
              cGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTc                                                     >  1:2546811/1‑71 (MQ=255)
                      tcttctGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGTCCCGGAATAGTTATAtt                                                       <  1:5015730/61‑1 (MQ=255)
                      tcttctGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCgg                                                                   <  1:1139399/49‑1 (MQ=255)
                          ctGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAggg                                           <  1:3838584/69‑1 (MQ=255)
                          ctGAAAACTACAGCGCGGTCGGGTACGAGGTGTATCGTGGcac                                                                     >  1:4274991/1‑41 (MQ=255)
                            gAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGc                                          >  1:4618418/1‑68 (MQ=255)
                             aaaaCTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCa                                        <  1:2111899/69‑1 (MQ=255)
                             aaaaCTCCAGCGCGGTCGGGTACGAGGTGTATCGTGACCCGGAATAGTTATATTTCTGGATCAGGGCCAGa                                      <  1:1180121/71‑1 (MQ=255)
                               aaCTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATc                                    >  1:2494233/1‑71 (MQ=255)
                                 cTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAgg                                            >  1:1089639/1‑61 (MQ=255)
                                 cTCCAGCGCGGTCGGGTACGAGGTGTATCGTGACCCGGAATAGTTATATTTCTTGATCAGGGCCAGATcc                                   >  1:4217076/1‑70 (MQ=255)
                                     aGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTc                              <  1:1479891/71‑1 (MQ=255)
                                      gcgcGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGa                                                 >  1:2930030/1‑51 (MQ=255)
                                      gcgcGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGa                                                 >  1:4333759/1‑51 (MQ=255)
                                         cgGTCGGGTACGAGGTGTATCGTGACCCGGAATAGTTATATTTc                                                     >  1:314883/1‑44 (MQ=255)
                                            tggggTACGAGGCGTATCGTGGCCCGGAATAGTTATATTTCTGGa                                                 <  1:4817488/43‑1 (MQ=38)
                                            tCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGa                                                 <  1:1149233/45‑1 (MQ=255)
                                            tCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGa                                                 <  1:4020875/45‑1 (MQ=255)
                                            tCGGGTACGAGGCGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTc                              <  1:48250/64‑1 (MQ=255)
                                              gggTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTg                     >  1:3181044/1‑71 (MQ=255)
                                                   cGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGAt                                                <  1:1223484/39‑1 (MQ=255)
                                                   cGAGGTGTATCGTGACCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTa                   >  1:2694224/1‑68 (MQ=255)
                                                   cGAGGTGTATCGTGACCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTACAg                <  1:3193887/71‑1 (MQ=255)
                                                   cGAGGTGTATCGTGACCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTACAg                <  1:305623/71‑1 (MQ=255)
                                                      ggTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTACAGACa             >  1:3905686/1‑71 (MQ=255)
                                                        tgtATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTACAGACATg           >  1:2798011/1‑71 (MQ=255)
                                                        tgtATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTACAGACATg           <  1:1665028/71‑1 (MQ=255)
                                                          tATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTACAGACATGCt         >  1:5007539/1‑71 (MQ=255)
                                                          tATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTACAGACATGCt         <  1:2454862/71‑1 (MQ=255)
                                                              gTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCa                                 >  1:3353493/1‑43 (MQ=255)
                                                              gTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCa                                 >  1:4862915/1‑43 (MQ=255)
                                                                 gCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTACAGACATGCTCACTCCt  <  1:3029314/71‑1 (MQ=255)
                                                                 |                                                                      
TAAAAACGCCTGTTCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTACAGACATGCTCACTCCT  >  W3110S.gb/3584511‑3584646

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