Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 232,901 C→T 10.0% intergenic (‑140/+30) hemB ← / ← ykiB porphobilinogen synthase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported232,9010CT10.0% 115.8 / 4.1 40intergenic (‑140/+30)hemB/ykiBporphobilinogen synthase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (18/18);  new base T (2/2);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTT  >  Exported/232831‑232970
                                                                      |                                                                     
gTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAAc                                                                       >  1:843355/1‑71 (MQ=255)
   gTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACtt                                                                     >  1:1746118/1‑70 (MQ=255)
     cgcgATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTta                                                                   >  1:835463/1‑70 (MQ=255)
               aCAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGa                                                         >  1:1191031/1‑70 (MQ=255)
                   gaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTatat                                                                <  1:1241264/59‑1 (MQ=255)
                         aaTGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATAttttt                                               >  1:1370422/1‑70 (MQ=255)
                             gAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTtat                                                                  <  1:2308883/47‑1 (MQ=255)
                              aaGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAatat                                                   >  1:1139609/1‑61 (MQ=255)
                              aaGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCtt                                         >  1:3444513/1‑71 (MQ=255)
                                gagagGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCttta                                       >  1:1461308/1‑71 (MQ=255)
                                  gagGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCtttatt                                     >  1:246809/1‑71 (MQ=255)
                                   agGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCtttattt                                    <  1:1921078/71‑1 (MQ=255)
                                   agGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCtttattt                                    <  1:1105051/71‑1 (MQ=255)
                                    ggCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCtttattta                                   >  1:314191/1‑71 (MQ=255)
                                    ggCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATACTTTTAATCtttattt                                    <  1:900364/70‑1 (MQ=255)
                                        aGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGc                               >  1:2711472/1‑71 (MQ=255)
                                          gACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTg                             <  1:2844478/71‑1 (MQ=255)
                                           aCATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGt                            >  1:3258840/1‑71 (MQ=255)
                                            cATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGtt                           >  1:1614413/1‑71 (MQ=255)
                                            cATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGtt                           <  1:3153533/71‑1 (MQ=255)
                                             aTCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAatat                                                   >  1:3398298/1‑46 (MQ=255)
                                               cTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATAttt                                                 >  1:3369744/1‑46 (MQ=255)
                                                  tttAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGttat                     >  1:3445995/1‑71 (MQ=255)
                                                      aTAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGttatta                   <  1:1526346/69‑1 (MQ=255)
                                                       tAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCtttatt                                     <  1:1135481/50‑1 (MQ=255)
                                                         gAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAttttt              >  1:3036537/1‑71 (MQ=255)
                                                         gAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAttttt              <  1:2739729/71‑1 (MQ=255)
                                                         gAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAttttt              <  1:1582520/71‑1 (MQ=255)
                                                           aaCAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGttatt                    >  1:3066173/1‑63 (MQ=255)
                                                           aaCAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGttat                     >  1:2232436/1‑62 (MQ=255)
                                                             cAGCAATAATTTTATATTCACTGAAAATATTTTTa                                              <  1:2783965/35‑1 (MQ=255)
                                                             cAGCAATAATTTTATATTCACTGAAAATATTTTTa                                              <  1:2662108/35‑1 (MQ=255)
                                                              aGCAATAACTTTATATTCACTGTAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAg         <  1:3617960/71‑1 (MQ=255)
                                                               gCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGt        <  1:3024746/71‑1 (MQ=255)
                                                                cAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGtt       >  1:2444217/1‑71 (MQ=255)
                                                                cAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGtt       <  1:806604/71‑1 (MQ=255)
                                                                 aaTAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGtt       <  1:1269262/70‑1 (MQ=255)
                                                                   tAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTgg    <  1:785825/71‑1 (MQ=255)
                                                                     aCTTTATATTCACTGAAAATATTTTTAATCTttatttat                                  >  1:1723124/1‑39 (MQ=255)
                                                                     aCTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTgg    <  1:2454692/69‑1 (MQ=255)
                                                                     aCTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGtt  <  1:1362992/71‑1 (MQ=255)
                                                                      |                                                                     
GTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTT  >  Exported/232831‑232970

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