Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,076,941 C→T 17.8% intergenic (+57/+88) ebgC → / ← uxaA cryptic beta‑D‑galactosidase, beta subunit/altronate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,076,9410CT17.8% 57.9 / 6.3 28intergenic (+57/+88)ebgC/uxaAcryptic beta‑D‑galactosidase, beta subunit/altronate hydrolase
Reads supporting (aligned to +/- strand):  ref base C (13/10);  new base T (2/3);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

AACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTC  >  Exported/2076876‑2076986
                                                                 |                                             
aacaaataacaaCTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGtt                                          >  1:696637/1‑71 (MQ=255)
 acaaataacaaCTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGATGTGtta                                         >  1:674936/1‑70 (MQ=255)
 acaaataacaaCTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGATGTGtta                                         >  1:4490368/1‑70 (MQ=255)
             tACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTGATAAATCCTTAg                             >  1:996873/1‑71 (MQ=255)
             tACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTGATAAATCCTTAg                             >  1:4125962/1‑71 (MQ=255)
               cggcggCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGATGTGtt                                          <  1:6115/56‑1 (MQ=255)
                  cggcAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTGATAAATCCTTAGCAGGt                        <  1:2556103/71‑1 (MQ=255)
                      aaaaGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTGATAAATCCTTAGCa                           >  1:1609770/1‑64 (MQ=255)
                          ggAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTGATAAATCCTTAGCAGGTATGGaa                  >  1:4520091/1‑69 (MQ=255)
                             gTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTGATAAATCCTTAGCAGGTATGGAAAAGCa             <  1:1292725/71‑1 (MQ=255)
                              tttGCCGCCACCGCTACCCTACTCATTTTCGGAGATgtgt                                           <  1:80979/40‑1 (MQ=255)
                              tttGCCGCCACCGCTACCCTACTCATTTTCGGAGATgtgt                                           <  1:4641401/40‑1 (MQ=255)
                              tttGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTGATAAATCCTTAGCAGGTATGGAAAAGCaa            <  1:2835140/71‑1 (MQ=255)
                                tGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTGATAAATCCTTAGCAGGTATGGAAAAGCAAAc          >  1:1938776/1‑71 (MQ=255)
                                tGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGGTGATaa                                     >  1:1157794/1‑44 (MQ=255)
                                        ccGCTACCCTACTCATTTTCGGAGACGTGTTGATAAATCCTTAGCAGGTATGGAAAAGCa             <  1:3774690/60‑1 (MQ=255)
                                        ccGCTACCCTACTCATTTTCGGAGACGTGTTGATAAATCCTTAGCAGGTATGGAAAAGCa             <  1:2510973/60‑1 (MQ=255)
                                        ccGCTACCCTACTCATTTTCGGAGACGTGTTGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTc  <  1:81953/71‑1 (MQ=255)
                                                  aCTCATTTTCGGAGACGTGTTGATAAATCCTTAGCAGGTATGGa                   <  1:3717421/44‑1 (MQ=255)
                                                  aCTCATTTTCGGAGACGTGTTGATAAATCCTTAGCAGGTATGGa                   <  1:433317/44‑1 (MQ=255)
                                                  aCTCATTTTCGGAGACGTGTTGATAAATCCTTAGCAGGTATGGa                   <  1:2299471/44‑1 (MQ=255)
                                                  aCTCATTTTCGGAGACGTGTTGATAAATCCTTAGCAGGTATGGa                   <  1:1486249/44‑1 (MQ=255)
                                                  aCTCATTTTCGGAGACGTGTTGATAAATCCTTAGCAGGTATGGAAAAGCa             >  1:1460200/1‑50 (MQ=255)
                                                  aCTCATTTTCGGAGACGTGTTGATAAATCCTTAGCAGGTATGGAAAAGCa             >  1:2584456/1‑50 (MQ=255)
                                                      aTTTTCGGAGACGTGTTGATAAATCCTTAGCAGGTATGGaaaa                >  1:754749/1‑43 (MQ=255)
                                                      aTTTTCGGAGACGTGTTGATAAATCCTTAGCAGGTATGGaaa                 >  1:995356/1‑42 (MQ=255)
                                                         ttCGGAGACGTGTTGATAAATCCTTAGCAGGTATGGAAAAGCAAACgg        >  1:2963432/1‑48 (MQ=255)
                                                         ttCGGAGACGTGTTGATAAATCCTTAGCAGGTATGGAAAAGCAAACgg        >  1:1765498/1‑48 (MQ=255)
                                                                 |                                             
AACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACGTGTTGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTC  >  Exported/2076876‑2076986

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: