Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 474,558 C→T A202T (GCT→ACT)  tesB ← acyl‑CoA thioesterase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913474,5580CT93.8% 36.3 / ‑6.1 16A202T (GCT→ACT) tesBacyl‑CoA thioesterase II
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (5/10);  minor base G (1/0);  total (6/10)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACGGGCGATGGAACCACATGGAATGGTCAATGGTGGCAATCTGAATCCCCGGTTCGAGAAAACCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCAGCGATTGGGCGATTTGCGTTTCCGAAGGGAGGCCATCAGGCGCTGGCGCGGACGGC  >  NC_000913/474359‑474799
                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
tGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACGGGCGATGGAACCACATGGAATGGTCAATGGTGGAAATCTGAATCCCCGGTTCGAGAAAACCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGTGTAACCGAGCAGATACTGATGCACGCGCAGGTCATCCGGCACGCTACCAtt                                                                                                                                                                                                >  1:458715/1‑251 (MQ=255)
     nnnnnnnnncncgnnCACANAGCCACGTGCGCTGGACGCCGAGGNGCTCTCCACGCGATACAGCAGCCATTCATTCNANTTAAACGGGCGATTGAACCACATGGAATGGTCAATGGTGGCAATCTGANTCCCCGGTTCNAGANNACCGNTNCNGTGCGGCTGTAGAGCTACCGGCAGGAAGTNAAGAACAGAAGTGTAACCGAGCANNGANNNATNNAcgcg                                                                                                                                                                                                                        <  2:373655/211‑1 (MQ=255)
                   gcACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACGGGCGATGGAACCACATGGAATGGTCAATGGTGGCAATCTGAATCCCCGGTTCGAGAAAACCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGTGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGc                                                                                                                                                                                                       <  1:29210/225‑1 (MQ=255)
                    cacaAAGCCACGTGCGCTGGACGCCGCGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTCAACGGGCGATGGAACCACATGGCATGGTCAATGGTGGCAATCTGAATCCCCGGTTCGAGAAAACCGATGCCGTGCGGCTGTAGCGCTACCGGCAGGAAGTTAAGATCAGAAGTGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGAc                                                                                                                                                                            >  1:220611/1‑251 (MQ=255)
                                                                      gcCATTCATTCAAATTAAACGGGCGATGGAACCACATGGAATGGTCAATGGTGGCAATCTGAATCCCCGGTTCGAGAAAACCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGTGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTgcg                                                                                                                                                                                            <  1:164968/185‑1 (MQ=255)
                                                                                caAATTAAGCGGGCGATGGAACCACATGGAATGGTCAATGGTGGCAATCTGAATCCCCGGTACGAGAAACCCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGTGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGGACGCTACCATTTGCGCGGATCCACACCGGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACcggacgg                                                                                                                       <  1:500143/244‑1 (MQ=255)
                                                                                     tAAACNNGCNNNGGNNCAACATGGAATGGTCAATGGTGGNAATCTGAATCCCCGGTTCGAGAAAACCGANGNCNTGCGNNTGTNGCGCTACCGGNAGGAAGTTAAGATCAGAAGTGTAACCGAGCAGATACTGCGGAACGCGCNGNTCATCAGGCCCGTTACCATTTGCGCGGGTCCNCACCTGACGATGTGGTTCTGCGAcncatcnnggnannnnnnnnnn                                                                                                                                       >  2:377276/1‑202 (MQ=255)
                                                                                     tAAACGGGCGATGGAACCACATTGAATGGTCAATGGTGGCAATCTGAATCCCCGGTTCGAGAATACCGATGCCGTGCGTCTGTAGAGCTACCGGCATGAAGTTAAGATCAGAAGTGTAACCGATCAGATACTGATTAACGCGCAGGGCATCCGGCACGCTACCATTTGCGCGGATCCACACCTTACGATGTGGTTCTGCGAc                                                                                                                                                            <  1:377276/202‑1 (MQ=255)
                                                                                           ggCGATGGAACCACATGGAATGGTCAATAGTGGCAATCTGAATCCCCGGTTCGAGAAAACCGAGGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGTGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCt                                                                                                                                                                               <  1:497663/177‑1 (MQ=255)
                                                                                            gCGATGGAACCACATGGAATGGTCAATGGTGGCAATCTGAATCCCAGGTTCGAGAAAACCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGTGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCgn                                                                                                    <  1:57949/251‑2 (MQ=255)
                                                                                                                                       cccGGTTCGAGAAAACCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGCTCAGAAGTGTCACCGAGCAGATACTGATGAACGCGCAGGTCATCCGACACCCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGAt                                                                                                >  1:79242/1‑212 (MQ=255)
                                                                                                                                                           gCCGTGCGGCTGTAGAGCTACCGGCAGGAAGGTAAGATCAGAAGGGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGNGACNNNNNAg                                                                                                                                                >  2:299177/1‑144 (MQ=255)
                                                                                                                                                            ccGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGTGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCAGCGATTGGGCGATTTGCGtttn                                    <  1:299177/251‑2 (MQ=255)
                                                                                                                                                                                            aaGATCAGAAGTGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGAGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAg                                                                                                  >  1:167864/1‑157 (MQ=255)
                                                                                                                                                                                              gATCAGAAGTGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTGCTGCGACGTGACCTGTCAGTGGGTTATGAAACTCCACCGGACGGACTGCCAGCGTACGACCGCAGATGAATGTATCTTTCAGCACGGGCGGCAGCAGGAGCGCCAGCGATGGGGCGATTTGCGTTTCCGAAGGGAGGCCATCTGGCGCTGGCGCGGACGGc  <  1:453316/251‑1 (MQ=255)
                                                                                                                                                                                                 cAGAAGTGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATc                                                                                                      <  1:235783/148‑1 (MQ=255)
                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
TGGGTATAAAACTCACCGCGCACAAAGCCACGTGCGCTGGACGCCGAGGTGCTCTCCACGCTATACAGCAGCCATTCATTCAAATTAAACGGGCGATGGAACCACATGGAATGGTCAATGGTGGCAATCTGAATCCCCGGTTCGAGAAAACCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGCCAGCGATTGGGCGATTTGCGTTTCCGAAGGGAGGCCATCAGGCGCTGGCGCGGACGGC  >  NC_000913/474359‑474799

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