Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 910,820 C→T G77D (GGC→GAC)  poxB ← pyruvate oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913910,8200CT88.9% 16.2 / ‑3.2 9G77D (GGC→GAC) poxBpyruvate oxidase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (4/4);  total (4/5)
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.69e-01

CGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTC  >  NC_000913/910613‑911047
                                                                                                                                                                                                               |                                                                                                                                                                                                                                   
cgcATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGCGCTAGAAACCAGCCCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGCTGGGTTTCCTGGAAATAGCCGCTGCCACTTTCGCTGGAGGGCATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGTCGGGGCCGCACGATCCGgc                                                                                                                                                                                                                  >  1:186546/1‑227 (MQ=255)
       cAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGTCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCg                                                                                                                                                                                               >  1:303241/1‑239 (MQ=255)
                                cTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGgccgn                                                                                                                                                            <  1:203191/249‑2 (MQ=255)
                                                  cTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGTCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGcc                                                                                                                                                                           >  1:6396/1‑216 (MQ=255)
                                                                                                                           cTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACTGGATTGCGGTGGTAATAGAACAGGCAGTTGATTAAGTGCAGGTTGTCGGGGCCGCACGATCCGGCGCAGTCCGCCAGTTCTCCGCTAAGTTGTTCTTCAgc                                                                                                                                                                              <  1:269542/140‑1 (MQ=255)
                                                                                                                                gggAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGTCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACtn                                                                                                                              <  1:29377/183‑2 (MQ=255)
                                                                                                                                                aTCGCCAGTACCGGAACGTGATTTCGGTGGCAATCGAACAGGCCGTTGATTAAGGGCATGTTGTCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAgcgc                                                                                                                                                                            <  1:125092/121‑1 (MQ=255)
                                                                                                                                                                       gNNGTGGCANNCGAACAGGCCGTTGATTAAGTGCAGGTTGTCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGcc                                                                                                                                                                           <  2:6396/99‑1 (MQ=255)
                                                                                                                                                                                        ggCCGTTGATTAAGTGCAGGTTGTCGGGGCCGCACGATCCGGCGCAGACCGGCAGTTCTCCGCTAAGTTGTGCTTCAGAGCCAGCGGCCAAGGCAGCCACTTATTCGTGGCGGGGGGACATCCCCTCCATGGTGACCATGGGCGTAAATCTGTCACTAAGACCGTTCTGAGAGGCACCTGTGCCTCCCCAGCTGCGTTTCACCCCCGACGATTCGGGTGGTTTTGCGCTATAAggtggccaaggttgtttt  >  1:223611/1‑234 (MQ=255)
                                                                                                                                                                                                               |                                                                                                                                                                                                                                   
CGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCGGAATAGCTCTTGTGGGTGGGTTTCCTGGAAATAGCCGCTGCCAATTTCGCTGGAGGGAATATGAGCGGCAATCGCCAGTACCGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTC  >  NC_000913/910613‑911047

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