Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,244,320 C→T 100% C113Y (TGC→TAC)  ldcA ← L,D‑carboxypeptidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,244,3200CT81.8% 70.5 / 10.8 33C113Y (TGC→TAC) ldcAL,D‑carboxypeptidase A
Reads supporting (aligned to +/- strand):  ref base C (1/5);  new base T (16/11);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 8.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGCCAAAGTTTGCTACCAGCATGGGGCCGCTAAAGGTGATGACATTGCCATGCGCCAGAAGACCGCACTGAATGGCGGTAAAATCGCTATGTCCGCAAATGAGCAACGGGTCATGTTGTTGGCGAGCCACCAGT  >  W3110S.gb/1244255‑1244388
                                                                 |                                                                    
cgcCAAAGTTTGCTACCAGCATGGGGCCGCTAAAGGTGATGACATTGCCATGCGCCAGAAGACCGCAc                                                                    <  1:2227155/68‑1 (MQ=255)
                 aGCATGGGGCCGCTAAAGGTGATGACATTGCCATGCGCCAGAAGACCGTACTGAATGGCGGTAAAATc                                                   >  1:3681157/1‑68 (MQ=255)
                 aGCATGGGGCCGCTAAAGGTGATGACATTGCCATGCGCCAGAAGACCGTACTGAATGGCGGTAAAATc                                                   >  1:3037819/1‑68 (MQ=255)
                  gCATGGGGCCGCTAAAGGTGATGACATTGCCATGCGCCAGAAGACCGTACTGAATGGCGGTAAAATCGc                                                 <  1:2781524/69‑1 (MQ=255)
                  gCATGGGGCCGCTAAAGGTGATGACATTGCCATGCGCCAGAAGACCGTACTGAATGGCGGTAAAATCGc                                                 <  1:248485/69‑1 (MQ=255)
                  gCATGGGGCCGCTAAAGGTGATGACATTGCCATGCGCCAGAAGACCGTACTGAATGGCGGTAAAATCGc                                                 <  1:2411636/69‑1 (MQ=255)
                  gCATGGGGCCGCTAAAGGTGATGACATTGCCATGCGCCAGAAGACCGTACTGAATGGCGGTAAAATCGc                                                 <  1:2247809/69‑1 (MQ=255)
                  gCATGGGGCCGCTAAAGGTGATGACATTGCCATGCGCCAGAAGACCGCACTGAATGGCGGTAAAATCGc                                                 >  1:1000290/1‑69 (MQ=255)
                     tGGGGCCGCTAAAGGTGATGACATTGCCATGCGCCAGAAGACCGTACTGAATg                                                              <  1:677587/53‑1 (MQ=255)
                     tGGGGCCGCTAAAGGTGATGACATTGCCATGCGCCAGAAGACCGTACTGAATGGCGGTAAAATCGCTAt                                              <  1:977659/69‑1 (MQ=255)
                     tGGGGCCGCTAAAGGTGATGACATTGCCATGCGCCAGAAGACCGTACTGAATGGCGGTAAAATCGCTAt                                              <  1:2005062/69‑1 (MQ=255)
                       gggCCGCTAAAGGTGATGACATTGCCATGCGCCAGAAGACCGTACTGAATGGCGGTAAAATCGCTATGt                                            >  1:1066013/1‑69 (MQ=255)
                       gggCCGCTAAAGGTGATGACATTGCCATGCGCCAGAAGACCGCACTGAATGGCGGTAAAATCGc                                                 <  1:3372646/64‑1 (MQ=255)
                         gCCGCTAAAGGTGATGACATTGCCATGCGCCAGAAGACCGTACTGAATGGCGGTAAAATCGCTAt                                              <  1:3201687/65‑1 (MQ=255)
                           cGCTAAAGGTGATGACATTGCCATGCGCCAGAAGACCGCACTg                                                                  <  1:2441442/43‑1 (MQ=255)
                                              gCCATGCGCCAGAAGACCGTACTGAATGGCGGTAAAATCGc                                                 >  1:845990/1‑41 (MQ=255)
                                              gCCATGCGCCAGAAGACCGTACTGAATGGCGGTAAAATCGc                                                 >  1:1206295/1‑41 (MQ=255)
                                              gCCATGCGCCAGAAGACCGTACTGAATGGCGGTAAAATCGCTATGTCCGCAAATGAGCa                               >  1:252841/1‑59 (MQ=255)
                                              gCCATGCGCCAGAAGACCGTACTGAATGGCGGTAAAATCGCTATGTCCGCAAATGAGCAACGGGTCAtg                     >  1:350855/1‑69 (MQ=255)
                                              gCCATGCGCCAGAAGACCGTACTGAATGGCGGTAAAATCGCTATGTCCGCAAATGAGCAACGGGTCAtg                     >  1:1581599/1‑69 (MQ=255)
                                              gCCATGCGCCAGAAGACCGTACTGAATGGCGGTAAAATCGCTATGTCCGCAAATGAGCAACGGGTCAtg                     >  1:2515316/1‑69 (MQ=255)
                                              gCCATGCGCCAGAAGACCGTACTGAATGGCGGTAAAATCGCTATGTCCGCAAATGAGCAACGGGTCAtg                     >  1:1218882/1‑69 (MQ=255)
                                              gCCATGCGCCAGAAGACCGTACTGAATGGCGGTAAAATCGCTATGTCCGCAAATGAGCAACGGGTCAtg                     >  1:787466/1‑69 (MQ=255)
                                                 aTGCGCCAGAAGACCGTACTGAATGGCGGTAAAATCGCTATGTCCGCAAATGAGCAACGGGTCAtgtt                   <  1:3010140/68‑1 (MQ=255)
                                                 aTGCGCCAGAAGACCGTACTGAATGGCGGTAAAATCGCTATGTCCGCAAATGAGCAACGGGTCAtgtt                   <  1:1389571/68‑1 (MQ=255)
                                                    cgcCAGAAGACCGTACTGAATGGCGGTAAAATCGCTATGt                                            >  1:2771395/1‑40 (MQ=255)
                                                    cgcCAGAAGACCGTACTGAATGGCGGTAAAATCGCTATGt                                            >  1:2592135/1‑40 (MQ=255)
                                                    cgcCAGAAGACCGTACTGAATGGCGGTAAAATCGCTATGt                                            >  1:1844052/1‑40 (MQ=255)
                                                        agaagaCCGTACTGAATGGCGGTAAAATCGCTATGTCCGCAAATGAGc                                >  1:228641/1‑48 (MQ=255)
                                                        agaagaCCGTACTGAATGGCGGTAAAATCGCTATGTCCGCAAATGAGc                                >  1:3238325/1‑48 (MQ=255)
                                                              ccGTACTGAATGGCGGTAAAATCGCTATGTCCGCAAATGAGCAACGGGTCATGTTGTTGGCGAGccacc     <  1:1420046/69‑1 (MQ=255)
                                                              ccGCACTGAATGGCGGTAAAATCGCTATGTCCGCAAATGAGCAACGGGTCATGTTGTTGGCGAGccacc     <  1:1651725/69‑1 (MQ=255)
                                                              ccGCACTGAATGGCGGTAAAATCGCTATGTCCGCAAATGAGCAACGGGTCATGTTGTTGGCGAGccacc     <  1:1149813/69‑1 (MQ=255)
                                                                 cACTGAATGGCGGTAAAATCGCTATGTCCGCAAATGAGCAACGGGTGATGTTGTTGGCGAGCCACCAGt  <  1:3385118/69‑1 (MQ=255)
                                                                 |                                                                    
CGCCAAAGTTTGCTACCAGCATGGGGCCGCTAAAGGTGATGACATTGCCATGCGCCAGAAGACCGCACTGAATGGCGGTAAAATCGCTATGTCCGCAAATGAGCAACGGGTCATGTTGTTGGCGAGCCACCAGT  >  W3110S.gb/1244255‑1244388

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