Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,551,306 C→T 30.8% A26V (GCC→GTC)  pydX → NADP‑dependent dihydropyrimidine dehydrogenase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,551,3060CT30.8% 11.7 / 5.9 13A26V (GCC→GTC) pydXNADP‑dependent dihydropyrimidine dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base T (2/2);  total (6/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.81e-01

CGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGACCAGTTGGCCGAGCGCTTCAGCGACCTGGCCCCACCCCTCACCGCGCGTCAGGCCGCCCTGGAAAGTGCACGCTGCCTGTACTGCTACGACGCCCCCTGCGTCAACGCTTGCCCCAGCGAAATAGACATTCC  >  NC_002947/4551184‑4551408
                                                                                                                          |                                                                                                      
cGCAAAACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGACCAGTTGGCCGAGCGCTACAGCGACCTGGCCCCACCCCACAcc                                                                                           >  1:72391/1‑136 (MQ=255)
       ccTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGACCAGTTGGCCGAGCGCTTCAGCGACCTGGCCCCACCCCTCACCGCGCGTc                                                                                    <  2:20530/136‑1 (MQ=255)
               cTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGACCAGTTGGCCGAGCGCTTCAGCGACCTGGTCCCACCCCTCACCGCGCGTCAGgccgcc                                                                            >  2:202818/1‑136 (MQ=255)
               cTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGACCAGTTGGCCGAGCGCTTCAGCGACCTGGTCCCACCCCTCACCGCGCGTCAGgccgcc                                                                            >  2:41477/1‑136 (MQ=255)
               cTGCACCAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCCACCAGTTGGCCGAGCGCTTCAGCGACCTGGCCCCACCCCTCACCGCGCGTCAGGCCGccc                                                                           <  2:185332/136‑1 (MQ=255)
                           aaaTAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGTCCGGCGCCGACCAGTTGGCCGAGCGCTTCAGCGACCTGGCCCCACCCCTCACCGCGCGTCAGGCCGCCCTGGAAAGTGCa                                                                <  2:180609/136‑1 (MQ=255)
                                      atACAACCGTGACCAACGCCCTGAACCTTTTGCCGCGCCCCCGGGCCGGGGCCGACCAGTTGGCCGAGCGCTTCAGCGACCTGGCCCCACCCCTCACCGCGCGTCAGGCCGCCCTGGAAAGTGCACGCTGCCTGTa                                                     <  2:31881/135‑1 (MQ=255)
                                                  ccAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGACCAGTTGGCCGAGCGCTTCAGCGACCTGGCCCCCCCCCTCACCGCGCGTCAGGCCGCCCTGGAAAGTGCACGCTGCCTGTACTGCTACGACGc                                         >  1:217329/1‑136 (MQ=255)
                                                                cATTTGCCGCGCCCCCGGGCCGGCGCCGACCAGTTGGCCGAGCGCTTCAGCGACCTGGCCCCCCCCCTCACCGCGCGTCAGGCCGCCCTGGAAAGTGCACGCTGCCGGTACTGCTACGACGCCCCCTGCGTCAAcc                           >  1:216442/1‑135 (MQ=255)
                                                                           cccccGGGCCGGCGCCGACCAGTTGGCCGAGCGCTTCAGCGACCTGGCCCCACCCCTCACCGCGCGTCAGGCCGCCCTGGAAAGTGCACGCTGCCTGTACTGCTACGACGCCCCCTGCGTCATCGCTTGCCCCAGc                <  1:90656/136‑1 (MQ=255)
                                                                                      gcgcCGACCAGTTGGCCGAGCGCTTCAGCGACCTGGTCCCACCCCTCACCGCGCGTCAGGCCGCCCTGGAAAGTGCACGCTGCCTGTACTGCTACGACGCCCCCTGCGTCAACGCTTGCCCCAGCGAAATAGACAt     <  1:202818/136‑1 (MQ=255)
                                                                                      gcgcCGACCAGTTGGCCGAGCGCTTCAGCGACCTGGTCCCACCCCTCACCGCGCGTCAGGCCGCCCTGGAAAGTGCACGCTGCCTGTACTGCTACGACGCCCCCTGCGTCAACGCTTGCCCCAGCGAAATAGACAt     <  1:41477/136‑1 (MQ=255)
                                                                                         ccGACCAGTTGGCCGAGCGCTTCAGCGACCTGGCCCCACCCCTCACCGCGCGTCAGGCCGCCCTGGAAAGTGCACGCTGCCTGTACTGCTACGACGCCCCCTGCGTCAACGCTTGCCCCAGCGAAATAGACATTcc  >  2:51182/1‑136 (MQ=255)
                                                                                                                          |                                                                                                      
CGCAATACCTGTCGCCTGCACCAAAAAAAATAGAGAGGCTACAACCGTGACCAACGCCCTGAACCATTTGCCGCGCCCCCGGGCCGGCGCCGACCAGTTGGCCGAGCGCTTCAGCGACCTGGCCCCACCCCTCACCGCGCGTCAGGCCGCCCTGGAAAGTGCACGCTGCCTGTACTGCTACGACGCCCCCTGCGTCAACGCTTGCCCCAGCGAAATAGACATTCC  >  NC_002947/4551184‑4551408

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