Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,450 G→T 100% L705L (CTG→CTT thrA → bifunctional aspartokinase I/homeserine dehydrogenase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,4500GT87.5% 37.4 / 1.2 16L705L (CTG→CTTthrAbifunctional aspartokinase I/homeserine dehydrogenase I
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (8/6);  total (8/8)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.17e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGAT  >  REL606/2318‑2591
                                                                                                                                    |                                                                                                                                             
tgAAATGGGTTATACCGCACCGTATCCGCGAGAGGATCTTGCTGGTTTGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACGGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTTACGc                                                                                                                              <  1:426117/149‑1 (MQ=255)
ggAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGc                                                                                                                              <  1:109650/150‑1 (MQ=255)
ggAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGc                                                                                                                              <  2:243134/150‑1 (MQ=255)
      gggTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAggg                                                                                                                        >  1:347379/1‑150 (MQ=255)
            taCCGAACCGGATCCGGGAGATGATCTTTCTGGTATGGAGGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGt                                                                                                                  <  1:127095/150‑1 (MQ=255)
                                    tCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTc                                                                                          >  1:609098/1‑150 (MQ=255)
                                      ttttTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTcac                                                                                        <  1:432321/150‑1 (MQ=255)
                                          gggTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCt                                                                                    <  1:614134/149‑1 (MQ=255)
                                           ggTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTc                                                                                   >  2:480575/1‑150 (MQ=255)
                                           ggTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTc                                                                                   >  2:615936/1‑150 (MQ=255)
                                           ggTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTc                                                                                   >  2:39729/1‑150 (MQ=255)
                                                         gcgTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGc                                                                     >  2:8415/1‑150 (MQ=255)
                                                             aaGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCcgcg                                                                 >  1:164261/1‑150 (MQ=255)
                                                                                gTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGtg                                              >  2:65243/1‑150 (MQ=255)
                                                                                                              atatTGAAATTGAACCTGTGCTTCCCGCAGAGTGTAACGCTGAGGGTGCTGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGCTCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCa                <  1:39608/150‑1 (MQ=255)
                                                                                                                            ccTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGAt  <  2:348375/150‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                             
GGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGAT  >  REL606/2318‑2591

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