Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,156,704 Δ1 bp coding (1255/3093 nt) ebgA → cryptic beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,156,7040C.96.3% 118.1 / ‑1.0 27coding (1255/3093 nt)ebgAcryptic beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (13/13);  minor base T (0/1);  total (13/14)
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.97e-01

CGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGAT  >  REL606/3156636‑3156771
                                                                    |                                                                   
cGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGAGTGGGaa                                                             <  2:485781/76‑1 (MQ=255)
 ggAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGAGTGGGaaa                                                            >  1:221376/1‑76 (MQ=255)
   aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGAGTGGGaaaaa                                                          <  2:686281/76‑1 (MQ=255)
   aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGAGTGGGaaaaa                                                          <  2:1246688/76‑1 (MQ=255)
      ccGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGAGTGGGAAAAAGTc                                                       <  2:574959/76‑1 (MQ=255)
        gACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCGTGGGAAAAAGTCTa                                                     >  2:807744/1‑76 (MQ=255)
        gACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGAGTGGGGAAAAGTCTa                                                     >  2:807727/1‑76 (MQ=255)
         aCGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGAGTGGGAAAAAGTCTAc                                                    <  2:1060870/76‑1 (MQ=255)
           gTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGAGTGGGAAAAAGTCTACGt                                                  >  1:965057/1‑76 (MQ=255)
             cGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGAGTGGGAAAAAGTCTACGTCg                                                <  2:453234/76‑1 (MQ=255)
             cGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGAGTGGGAAAAAGTCTACGTCg                                                >  1:575409/1‑76 (MQ=255)
             cGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGAGTGGGAAAAAGTCTACGTCg                                                <  2:1022966/76‑1 (MQ=255)
               aaTCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGAGTGGGAAAAAGTCTACGTCGAg                                              >  2:1356049/1‑76 (MQ=255)
                         cTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTc                                    >  2:1249102/1‑76 (MQ=255)
                          tttGCTAATTTCGGCGATATTAGCCGTATTACCGACGATCCTAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCg                                   <  1:1105300/76‑1 (MQ=255)
                             gCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCa                                <  2:80910/76‑1 (MQ=255)
                                aaTGTCGGCGATATTAGCCGTATTACCGACGATCCGAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCatat                             >  2:320550/1‑76 (MQ=255)
                                     cGGCGATATTAGCCGTATTACCGACGATCCGAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCAcg                        >  2:829522/1‑76 (MQ=255)
                                         gATATTAGGCGTATTACCGACGATCCGAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCa                    >  1:54253/1‑76 (MQ=255)
                                          atatTAGCCGTATTACCGACGATCCGAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAg                   <  1:377083/76‑1 (MQ=255)
                                             ttAGCCGTATTACCGACGATCCGAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGaaa                >  2:292097/1‑76 (MQ=255)
                                               aGCCGTATTACCGACGATCCGAGTGGGACAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGaaaac              >  1:142555/1‑75 (MQ=255)
                                               aGCCGTATTACCGACGATCCGAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGaaaaa              <  1:804615/76‑1 (MQ=255)
                                                     aTTACCGACGATCCGAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATcc        >  2:1135436/1‑76 (MQ=255)
                                                       tACCGACGATCCGAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGt      <  1:978573/76‑1 (MQ=255)
                                                          cgacgaTCCGAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGa   >  2:226701/1‑76 (MQ=255)
                                                           gacgaTCCGAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGat  <  1:987924/76‑1 (MQ=255)
                                                           gacgaTCCGAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGat  <  2:636334/76‑1 (MQ=255)
                                                                    |                                                                   
CGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGAT  >  REL606/3156636‑3156771

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: