Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,067,926 G→A 100% A40V (GCA→GTA)  yghS ← predicted protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,067,9260GA83.4% 23.1 / ‑0.8 12A40V (GCA→GTA) yghSpredicted protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (7/3);  total (7/5)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTGTTGAGGGCTTTTCGTGCACATGTGCCGCTTTACTTCGCAGATAACGCCCAAAAGGTGCGCCAATAACAGGGAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTGCTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGAACCATCACATCCTACAATGGCAATAACACGCACTGGCGTGGAATTAATTATTGACATATTACGTTGATTCACGAAAAACCCGGCAGTAAATAATGTATTG  >  REL606/3067783‑3068072
                                                                                                                                               |                                                                                                                                                  
gtgtTGAGGGCTTTTCGTGCACATGTGCCGCTTTACTTCGCAGATAACGCCCAAAAGGTGCGCCAATAACAGGGAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCa                                                                                                                                              >  1:167541/1‑150 (MQ=255)
             ttCGTGCACATGTGCCGCTTTACTTCGCAGATAACGCCCAAAAGGTGCGCCAATAACAGGGAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAg                                                                                                                                 >  1:530786/1‑150 (MQ=255)
                           ccGCTTTACTTCGCAGATAACGCCCAAAAGGTGCGCCAATAACAGGGAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGCACGAGATAAATGTGTTCTGTTGGCAATCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAggg                                                                                                                   >  1:141787/1‑150 (MQ=255)
                                  aCTTCGCAGATAACGCCCAAAAGGTGCGCCAATAACAGGGAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAttt                                                                                                            >  1:669759/1‑150 (MQ=255)
                                           aTAACGCCCAAAAGGTGCGCCAATAACAGGGAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGAAcc                                                                                                   <  2:532323/150‑1 (MQ=255)
                                                    aaaaGGTGCGCCAATAACAGGGAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTGCTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGAACCATCACATcc                                                                                          <  2:249633/150‑1 (MQ=255)
                                                          tGCGCCAATAACAGGGAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGAACCATCACATCCTACAAt                                                                                    <  2:582077/150‑1 (MQ=255)
                                                            cgcCAATAACAGGGAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTGCTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGAACCATCACATCCTACAATgg                                                                                  <  2:513402/150‑1 (MQ=255)
                                                                 aTAACAGGGAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGAACCATCACATCCTACAATGGCAATa                                                                             >  1:235472/1‑150 (MQ=255)
                                                                                               aTTCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGAACCATCACATCCTACAATGGCAATAACACGCACTGGCGTGGAATTAATTATTGAc                                               <  2:168017/150‑1 (MQ=255)
                                                                                                                         aaTGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGAACCATCACATCCTACAATGGCAATAACACGCACTGGCGTGGAATTAATTATTGACATATTACGTTGATTCACGAAAAAccc                     >  1:443055/1‑150 (MQ=255)
                                                                                                                                            cTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGAACCATCACATCCTACAATGGCAATAACACGCACTGGCGTGGAATTAATTATTGACATATTACGTTGATTCACGAAAAACCCGGCAGTAAATAATGTATTg  >  1:265865/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GTGTTGAGGGCTTTTCGTGCACATGTGCCGCTTTACTTCGCAGATAACGCCCAAAAGGTGCGCCAATAACAGGGAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTGCTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGAACCATCACATCCTACAATGGCAATAACACGCACTGGCGTGGAATTAATTATTGACATATTACGTTGATTCACGAAAAACCCGGCAGTAAATAATGTATTG  >  REL606/3067783‑3068072

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