Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,579,050 G→A A199V (GCC→GTC)  yhhW ← quercetin 2,3‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,579,0500GA93.1% 77.4 / ‑5.9 29A199V (GCC→GTC) yhhWquercetin 2,3‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (15/12);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CATCGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGGCTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCG  >  NC_000913/3578840‑3579263
                                                                                                                                                                                                                  |                                                                                                                                                                                                                     
cATCGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTg                                                                                                                                                                               <  1:102451/251‑1 (MQ=255)
   cGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGCTGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGCGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGAt                                                                                                                                                                            >  2:83240/1‑251 (MQ=255)
         gtgCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGAc                                                                                                                                                                      >  1:8202/1‑251 (MQ=255)
                    gcgGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTc                                                                                                                                                             >  1:64240/1‑249 (MQ=255)
                           aaCCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAAt                                                                                                                                                    >  2:248681/1‑251 (MQ=255)
                                                 aCGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCt                                                                                                                              >  2:64070/1‑251 (MQ=255)
                                                  cgcgACCGGAGCGAGCGAGGAAGCCGGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTAGCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGGCTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCtt                                                                                                                             <  2:236455/251‑1 (MQ=255)
                                                                      aaGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTAc                                                                                                         <  1:87464/251‑1 (MQ=255)
                                                                            gCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCGCATCTTTCAGCAACGCCCAGCGGTACAGTTcc                                                                                                   <  1:111571/251‑1 (MQ=255)
                                                                              aaaaaGGCAGGCTTAACATATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACCGTGACACTTCGCTGTCGCTATCCGCAGGGATTGAGATTTCCGGCGAATCCCCGATTGCCAGACAATCGCTGGTCGAGACTGTCAGGCCGTTAATGGAGACCTTGCCTTTCACCGCCAAGATCCAGACGCGCCGTTCGGCGGCAATCTGATGCACCGCCTGCTCCTCTTTCAGCAACGACCGTCGGTACAGTTCCat                                                                                                 >  2:38049/1‑251 (MQ=255)
                                                                                 aaGGCAGGCTTGAAGGATTCCGTTGATTTTTAACCCGGCGGCTGATCGAACAGTAACTCGTCGCTGTCGCAATCGGCATGGAGGGATATTTCCTTCTCATCCGAGATTGCCAGACCATTACGGGTCGAGGCTTTCACGCCGTTAATGGTGACATTGCTTGTCACCACCTGGATCTAGAGGCGTCGTTCAGCGGCAGTCTGATGCACCGACGGCTCGTGTTTCAGCAACGCCCAGCGGTACAGTTCCATATc                                                                                              <  2:145460/251‑1 (MQ=255)
                                                                                           tCAAAGATGCCGTTGAGTTGTAAACCGGCGGCAGATCGAACAGTACCACTTCGCTGTCGCTATCCGCATGGAGGGAGATTGCCTGCTCATCCCAGATGGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGGTCAGCGGCAATCTGATGCACCTACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCAc                                                                                     <  2:144276/250‑1 (MQ=255)
                                                                                                   gCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGCAACACTTCGCTGTCGCTCTCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATg                                                                             >  1:128112/1‑250 (MQ=255)
                                                                                                      gTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCg                                                                                                                                         <  2:84432/187‑1 (MQ=255)
                                                                                                      gTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCg                                                                                                                                         >  1:84432/1‑187 (MQ=255)
                                                                                                       ttGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAAcc                                                                         >  2:231201/1‑250 (MQ=255)
                                                                                                        tGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTAATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTAATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGCTGCACTTTCAATGAACCAt                                                                       >  1:278558/1‑251 (MQ=255)
                                                                                                          aGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCAtct                                                                     <  1:34930/251‑1 (MQ=255)
                                                                                                                               cGAACAGTCACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCAcc                                                >  2:267479/1‑251 (MQ=255)
                                                                                                                                        acacTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTg                                       >  2:63663/1‑251 (MQ=255)
                                                                                                                                        acacTTCGCTGTCGCTACCCGAATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTg                                       <  2:229362/251‑1 (MQ=255)
                                                                                                                                                 tGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACg                              >  1:240382/1‑251 (MQ=255)
                                                                                                                                                 ggtggCTATCCGCATGGATGGAGATTGCCGGCTCATCCCAGATTGCCGGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTGTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACg                              <  2:44160/247‑1 (MQ=255)
                                                                                                                                                  gTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACgg                             >  1:8979/1‑251 (MQ=255)
                                                                                                                                                                  atggaGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTg             <  1:248681/251‑1 (MQ=255)
                                                                                                                                                                         ttGCCTGCTCACCCCAGATTGCCAAACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCTCCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATa       <  2:472/250‑1 (MQ=255)
                                                                                                                                                                           gCCTGCTCCTCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTGCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGATCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCGTAAcg    <  2:181272/251‑1 (MQ=255)
                                                                                                                                                                             cTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGACGCGACGCTGTTCATAAcgcg  >  2:194140/1‑251 (MQ=255)
                                                                                                                                                                               gCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGACTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCGCTGTACGGCATCGAAGCGACGCTGTTCATAAcgcg  <  2:19351/249‑1 (MQ=255)
                                                                                                                                                                                                                  |                                                                                                                                                                                                                     
CATCGCATTGTGCTAACGACGCGGCAAAACCCCGGAAGCGTACATCAGCACGCGACCGGAGCGAGCGAGGAAGCCTGCAAAAAGGCAGGCTTCAAAGATGCCGTTGAGTTTTAAACCGGCGGCAGATCGAACAGTAACACTTCGCTGTCGCTATCCGCATGGATGGAGATTGCCTGCTCATCCCAGATTGCCAGACCATCGCTGGTCGAGGCTTTCACGCCGTTAATGGTGACATTGCCTTTCACCACCTGGATCCAGACGCGGCGTTCAGCGGCAATCTGATGCACCGACTGCTCATCTTTCAGCAACGCCCAGCGGTACAGTTCCATATCCTGATGCACTTTCAATGAACCATCTCGCGCATCCGGCGAGAGCACCAGCTGTTTGCCCTGTACGGCATCGAAGCGACGCTGTTCATAACGCG  >  NC_000913/3578840‑3579263

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