Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 651,069 T→C 45.5% N4S (AAT→AGT)  citC ← citrate lyase synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb651,0690TC45.5% ‑0.9 / 27.8 22N4S (AAT→AGT) citCcitrate lyase synthetase
Reads supporting (aligned to +/- strand):  ref base T (5/7);  new base C (7/3);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

TTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATG  >  W3110S.gb/650999‑651138
                                                                      |                                                                     
ttCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCAt                                                                       >  1:112099/1‑71 (MQ=255)
           aTTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCACTGCCGAACata                                                            <  1:1651816/71‑1 (MQ=255)
           aTTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCACTGCCGAACat                                                             >  1:1480848/1‑70 (MQ=255)
              gggCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACataa                                                           >  1:53143/1‑69 (MQ=255)
                 cGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATaataaa                                                       >  1:1268708/1‑70 (MQ=255)
                              ttttttATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACataat                                                          <  1:647054/54‑1 (MQ=255)
                              ttttttATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACataat                                                          <  1:1606677/54‑1 (MQ=255)
                                ttttATTTTCTGAACGTTTTACGCGGGTGAAAATATCACTGCCGAACATAATAAATAGTATCCTGAAGGTg                                       >  1:1066608/1‑71 (MQ=255)
                                ttttATTTTCTGAACGTTTTACGCGGGTGAAAATATCACTGCCGAACATAATAAATAGTATCCTGAAGGTg                                       >  1:1355677/1‑71 (MQ=255)
                                       ttCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgttg                                >  1:1529855/1‑71 (MQ=255)
                                        tCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCAtgttg                                <  1:1014173/70‑1 (MQ=255)
                                            aaCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATc                           >  1:17704/1‑71 (MQ=255)
                                            aaCGTTTTACGCGGGTGAAAATATCACTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATc                           >  1:1396091/1‑71 (MQ=255)
                                             aCGTTTTACGCGGGTGAAAATATCACTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATc                           >  1:123680/1‑70 (MQ=255)
                                             aCGTTTTACGCGGGTGAAAATATCACTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCg                          <  1:641924/71‑1 (MQ=255)
                                                       cgGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACg                <  1:1119521/71‑1 (MQ=255)
                                                          gTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGc                    <  1:573676/64‑1 (MQ=255)
                                                           tGAAAATATCACTGCCGAACATAATAAATAGTATCCTGAAGGTg                                       <  1:1722206/44‑1 (MQ=255)
                                                            gAAAATATCACTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAAtgt           >  1:1633545/1‑71 (MQ=255)
                                                             aaaaTATCATTGCCGAACATAATAAATAGTATCCTg                                             <  1:1663121/36‑1 (MQ=255)
                                                               aaTATCACTGCCGAACATAATAAATAGTATCCTGAAg                                          >  1:782210/1‑37 (MQ=255)
                                                                atatCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTa                             >  1:474275/1‑49 (MQ=255)
                                                                     aTTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATg  <  1:962338/71‑1 (MQ=255)
                                                                      |                                                                     
TTCATGCAGGAATTGGGCGATTTCCGCCATTTTTTTATTTTCTGAACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATG  >  W3110S.gb/650999‑651138

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