Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,343,489 G→A 100% L542L (CTG→TTG)  yciR ← modulator of Rnase II stability

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,343,4890GA91.7% 63.4 / ‑1.3 24L542L (CTG→TTG) yciRmodulator of Rnase II stability
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (6/16);  total (6/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAGTTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGAAAATGGTTTGATCGGCGAGCTGACGTGCAGAAATATTTAC  >  REL606/1343342‑1343624
                                                                                                                                                   |                                                                                                                                       
gAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAAtt                                                                                                                                       <  2:464073/150‑1 (MQ=255)
 aaTATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAGTtc                                                                                                                                      <  1:418696/150‑1 (MQ=255)
 aaTATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATtc                                                                                                                                      <  2:461241/150‑1 (MQ=255)
                       tCAAGTTTGATGGCATCGATCGGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAGTTCATCATTCTCAATCAGACAActc                                                                                                                <  1:441841/150‑1 (MQ=255)
                             tttgtGGCATCCATCGGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGGACCTAGTTGGCTAAATTGTTGAATCACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCGGTc                                                                                                          <  1:210978/146‑1 (MQ=255)
                             ttGATGGCATCGATCGGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTc                                                                                                          <  2:496311/150‑1 (MQ=255)
                                       cgatcgGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACAt                                                                                                <  1:299002/150‑1 (MQ=255)
                                        gttcgGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGAACCTAGTTGGCTAAATTGTTGAATAACGGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATc                                                                                               <  1:17542/148‑1 (MQ=255)
                                        gatcgGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAACAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATc                                                                                               <  2:638959/150‑1 (MQ=255)
                                         atcgGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCt                                                                                              >  1:554001/1‑150 (MQ=255)
                                               aaGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAggg                                                                                        <  2:31774/150‑1 (MQ=255)
                                                taagagacAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGc                                                                                       <  2:398067/143‑1 (MQ=255)
                                                agacactcAGATGCGAAAATTAAGAGAAACCGGTACCAAAATCAGACAGTTGTACTTGCGCACCTAGTTTGCTAAAATGTTGAATAACAGACGGTGGCAATTCATCATTCTCAATCAGACAAATCTCTGTCAGTTCAACATCTATAGGGc                                                                                       <  1:498471/143‑1 (MQ=255)
                                                aGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGc                                                                                       <  2:25447/150‑1 (MQ=255)
                                                aGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGc                                                                                       <  2:199524/150‑1 (MQ=255)
                                                aGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGc                                                                                       <  1:277119/150‑1 (MQ=255)
                                                          ttGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGcctgtctctta                                                                             >  1:396917/1‑140 (MQ=255)
                                                                                               aGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGaaa                                        <  1:61255/150‑1 (MQ=255)
                                                                                                   gTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGAAAATgg                                    <  2:555608/150‑1 (MQ=255)
                                                                                                          ggcACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGAAAATGGTTTGATc                             <  1:580628/149‑1 (MQ=255)
                                                                                                                      gCTAAATTGTTGAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACGACTCCCTGTCAGTTCCACATCTATAGGGCAGTCTTCAACATTGAGCTCCTGGAGAACCTTTTTCAGGGCGGTGAAAATGGTTTGATCGGCGAGCTGACg                 >  1:112404/1‑150 (MQ=255)
                                                                                                                              gttgAATAACAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGAAAATGGTTTGATCGGCGAGCTGACGTGCAGaaa         >  2:346674/1‑150 (MQ=255)
                                                                                                                                     aaCAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGAAAATGGTTTGATCGGCGAGCTGACGTGCAGAAATATTTAc  >  2:37780/1‑150 (MQ=255)
                                                                                                                                     aaCAGACAGTGCCAATTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGAAAATGGTTTGATCGGCGAGCTGACGTGCAGAAATATTTAc  >  2:3376/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                       
GAATATCTCGAACAAAAACCTGGTCAAGTTTGATGGCATCGATCGGAAAGCGCGCCAGTTGCGAAAGTGAAGAGTAGCCGGTACCAAAATCATCCAGATGTACTTGCGCACCTAGTTGGCTAAATTGTTGAATAACAGACAGTGCCAGTTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGAAAATGGTTTGATCGGCGAGCTGACGTGCAGAAATATTTAC  >  REL606/1343342‑1343624

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