Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,460,012 T→C V320A (GTT→GCT)  paaH → 3‑hydroxyadipyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,460,0120TC97.1% 118.5 / NA 35V320A (GTT→GCT) paaH3‑hydroxyadipyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (17/17);  total (18/17)
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.00e+00

GGCAAGAAAAGTGGGCTGGGCGTGTACGACTGGCGCGCGGAACGTGAGGCAGTTGTTGGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGTTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCC  >  NC_000913/1459876‑1460135
                                                                                                                                        |                                                                                                                           
ggCAAGAAAAGTGGGCTGGGCGTGTACGACTGGCGCGCGGAACGTGAGGCAGTTGTTGGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGTAAATAAAGGTGACGGTTTCACGGAAATTGACGATGcttt                                                                                                                          >  2:202423/1‑140 (MQ=255)
       aaaGTGGGCTGGGCGTGTACGACTGGCGCGCGGAACGTGAGGCAGTTGTTGGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTattgatt                                                                                                                   <  2:476383/140‑1 (MQ=255)
          gTGGGCTGGGCGTGTACGACTGGCGCGCGGAACGTGAGGCAGTTGTTGGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGCAATTGACGAGGCTTTATTGATTgag                                                                                                                >  1:286838/1‑140 (MQ=255)
           tgggctgggcGTGTACGACTGGCGCGCGGAACGTGAGGCAGTTGTTGGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTgaga                                                                                                               >  1:164454/1‑140 (MQ=255)
             ggctgggcGTGTACGACTGGCGCGCGGAACGTGAGGCAGTTGTTGGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGaca                                                                                                             >  1:434781/1‑140 (MQ=255)
             ggctgggcGTGTACGACTGGCGCGCGGAACGTGAGCCAGTTGTTGGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGaca                                                                                                             <  2:284340/140‑1 (MQ=255)
                  ggcGTGTACGACTGGCGCGCGGCACGTGAGGCAGTGGTTGGCCTGGAATCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAgg                                                                                                        <  2:290961/140‑1 (MQ=255)
                            aCTGGCGCGCGGAACGTGAGGCAGTTGTTGGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGca                                                                                              <  1:343599/140‑1 (MQ=255)
                            aCTGGCGCGCGGAACGTGAGGCAGTTGTTGGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGca                                                                                              >  1:598050/1‑140 (MQ=255)
                            aCTGGCGCGCGGAACATGAGGCAGTTGTTCGCCTGGATGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGca                                                                                              <  1:232938/140‑1 (MQ=255)
                                gcgcgcGGAACGTGAGGCAGTTGTTGGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGTTTTTTTGATTGAAACACGAGGCGAGACGGCACAgg                                                                                          >  2:261748/1‑140 (MQ=255)
                                gcgcgcGGAACGTGAGGCAGTTGTTGGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAgg                                                                                          >  1:618830/1‑140 (MQ=255)
                                gcgcgcGGAACGAGAGGCAGTTGTTGGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCAGGGAAATTGACGATGGTTTTGTGATTGATACACATGGCGAGACGGCACAgg                                                                                          >  2:318215/1‑140 (MQ=255)
                                 cgcgcgGAACGTGAGGCAGTTGTTGGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGgc                                                                                         <  2:370451/140‑1 (MQ=255)
                                                   gttgttGGCCTGGAAGCGGTAATCGACAGTTCTAGGCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGc                                                                       <  2:286838/140‑1 (MQ=255)
                                                    ttgttgGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGGACCGGCGCTTGCAATACGACTGGCa                                                                      >  1:217899/1‑140 (MQ=255)
                                                    ttgttgGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTTGCAATACGACTGGCa                                                                      >  1:247166/1‑140 (MQ=255)
                                                    ttgttgGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCa                                                                      >  2:611165/1‑140 (MQ=255)
                                                    ttgttgGCCTGGAAGCGGTAAGCGACAGTTTTACCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAACGCGAGACGGCACAGGCGCTGGCAATACGCCTGGCa                                                                      >  2:142142/1‑140 (MQ=255)
                                                      gttgGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACg                                                                    <  1:503153/140‑1 (MQ=255)
                                                      gttgGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACg                                                                    <  2:598050/140‑1 (MQ=255)
                                                      gttgGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACg                                                                    <  1:368674/140‑1 (MQ=255)
                                                      gttgGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACg                                                                    <  1:25561/140‑1 (MQ=255)
                                                        tgGCCTGGAAGCGGTCAGCGACAGTCTTAACCCAATGCAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGcc                                                                  <  1:202423/140‑1 (MQ=255)
                                                        tgGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGcc                                                                  <  1:499047/140‑1 (MQ=255)
                                                              ggAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCggtgg                                                            >  1:456482/1‑140 (MQ=255)
                                                                   cGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTg                                                       <  2:434781/140‑1 (MQ=255)
                                                                              aGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGTATTACGACTGGCAAGCTCGGTGGGAGTGATCGATAAAAt                                            >  1:226320/1‑140 (MQ=255)
                                                                                        cAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATAGGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAg                                  >  1:411926/1‑140 (MQ=255)
                                                                                        cAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAg                                  <  1:586862/140‑1 (MQ=255)
                                                                                          aTGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAggg                                >  1:654595/1‑139 (MQ=255)
                                                                                           tGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAggtg                               <  1:632162/140‑1 (MQ=255)
                                                                                           tGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAggtg                               <  1:501254/140‑1 (MQ=255)
                                                                                                    aaaagaaaagTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGTGGTGACCAtt                      >  1:102805/1‑140 (MQ=255)
                                                                                                           aagTGACGGTGTCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGAGGTGACCATtgctgctg               >  2:277209/1‑140 (MQ=255)
                                                                                                                   gtgtCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGTGGTGACCATTGCTGCTGCAGCGGTg       >  1:476548/1‑140 (MQ=255)
                                                                                                                   gtgtCACGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGTGGTGACCATTGCTGCTGCAGCGGTg       >  1:112456/1‑140 (MQ=255)
                                                                                                                        aCGGAAATTGACGATGCTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGTGGTGACCATTGCTGCTGCAGCGGGGAAccc  >  2:522151/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                           
GGCAAGAAAAGTGGGCTGGGCGTGTACGACTGGCGCGCGGAACGTGAGGCAGTTGTTGGCCTGGAAGCGGTAAGCGACAGTTTTAGCCCAATGAAAGTAGAAAAGAAAAGTGACGGTGTCACGGAAATTGACGATGTTTTATTGATTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCC  >  NC_000913/1459876‑1460135

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