Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,438,531 T→G D110A (GAC→GCC)  fbp ← fructose‑1,6‑bisphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,438,5310TG96.8% 118.1 / ‑4.0 31D110A (GAC→GCC) fbpfructose‑1,6‑bisphosphatase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (13/17);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAA  >  REL606/4438461‑4438603
                                                                      |                                                                        
gcggcgGTAGATGGAGAAAATGGTACCGACAGAGACGTTAAGATCGCTGTTGGACGAGCCATCCAGGAGGGCCATc                                                                     >  2:1310528/1‑76 (MQ=255)
 cggcggTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCa                                                                    <  2:172599/76‑1 (MQ=255)
 cggcggTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCa                                                                    <  2:607995/76‑1 (MQ=255)
     ggTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGcac                                                                <  1:1311214/76‑1 (MQ=255)
       tAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGcacca                                                              >  2:83300/1‑76 (MQ=255)
       tAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGcacca                                                              >  2:326712/1‑76 (MQ=255)
       tAGAGGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGcacca                                                              <  1:1666195/76‑1 (MQ=255)
                    tGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGTGCCATCAGCACCACGTACTTTGCgtg                                                 >  1:178919/1‑76 (MQ=255)
                    tGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCgtg                                                 >  2:444990/1‑76 (MQ=255)
                         ccGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTcac                                            <  1:1069102/76‑1 (MQ=255)
                            aCAGAAACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGc                                         <  2:1823558/76‑1 (MQ=255)
                                  aCGTTAACATCGGTGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCaa                                   <  2:190417/76‑1 (MQ=255)
                                  aCGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCaa                                   >  1:1378142/1‑76 (MQ=255)
                                    gTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAg                                 <  1:409832/76‑1 (MQ=255)
                                     ttAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAga                                <  1:665978/76‑1 (MQ=255)
                                       aaCATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAgacg                              >  2:616375/1‑76 (MQ=255)
                                       aaCATCGATGTCGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAgacg                              >  2:525640/1‑76 (MQ=255)
                                        aCATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAgacga                             >  2:1046763/1‑76 (MQ=255)
                                         cATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAgacgac                            <  1:1694044/76‑1 (MQ=255)
                                          aTCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACa                           <  2:1056966/76‑1 (MQ=255)
                                          aTCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACa                           <  2:1572039/76‑1 (MQ=255)
                                            cGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAAt                         >  2:1567841/1‑76 (MQ=255)
                                                 ttGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAt                    <  1:1111475/76‑1 (MQ=255)
                                                  tGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAtc                   >  2:1782923/1‑76 (MQ=255)
                                                   ggACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAtct                  >  2:116965/1‑76 (MQ=255)
                                                         gCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAtcttcttct            >  1:888984/1‑76 (MQ=255)
                                                          ccATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATcttcttctt           <  1:427901/76‑1 (MQ=255)
                                                          ccATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATcttcttctt           <  1:219716/76‑1 (MQ=255)
                                                           cATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCttcttcttc          <  1:1530123/76‑1 (MQ=255)
                                                              ccAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCaga       >  2:420624/1‑76 (MQ=255)
                                                                 ggggggCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGc    >  2:1110881/1‑76 (MQ=255)
                                                                 ggggggCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGc    <  2:815824/76‑1 (MQ=255)
                                                                  ggggGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCa   <  1:934713/76‑1 (MQ=255)
                                                                   ggggCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCaa  >  1:187148/1‑76 (MQ=255)
                                                                      |                                                                        
GCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAA  >  REL606/4438461‑4438603

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