Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,134,498 G→A V385V (GTC→GTT)  wzc ← protein‑tyrosine kinase Wzc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,134,4980GA96.8% 94.6 / ‑4.3 31V385V (GTC→GTT) wzcprotein‑tyrosine kinase Wzc
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (19/11);  total (19/12)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

AATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTG  >  NC_000913/2134261‑2134736
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
aaTACCGTGTTCTGCCAGGACTTGCGGGCGTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCATGATAGAAAGCATCAGGCCAAGGATAATCGCTCCGATGATAATCAGCCCTGTCTGCGGTTTCAGCACGCCATGCTGAGTGATTGCCTGGTCAACAATGCTCACATCGCCGACGGTGCTGGCCTCGGGGATGTTCAGCTCCTGCACTTTATTCGGCAGTTTGATATAGACCTGGTGAcc                                                                                                                                                                                                                                     <  2:413963/249‑1 (MQ=255)
        gtTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGTTTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGAc                                                                                                                                                                                                                            >  2:140818/1‑250 (MQ=255)
                              ttCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAAtct                                                                                                                                                                                                     >  1:108633/1‑251 (MQ=255)
                                 aaTGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCct                                                                                                                                                                                                  >  1:245239/1‑251 (MQ=255)
                                 aaTGCAGCGATTACACAACGCGCGCCGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGCCGGACAATCTCct                                                                                                                                                                                                  >  2:454064/1‑251 (MQ=255)
                                   tGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGACAGCATCCGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCACCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCACTTTATTCAGCAGTTGCATATAAACCTGCTGAACAGACTCGACATCGCGGGTCAGACGGACAATATCctgc                                                                                                                                                                                                >  2:41153/1‑251 (MQ=255)
                                          aTTAAAAAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGttt                                                                                                                                                                                         >  1:109948/1‑251 (MQ=255)
                                              aacaacGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACCATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCgg                                                                                                                                                                                     >  2:113728/1‑251 (MQ=255)
                                                     aGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGcc                                                                                                                                                                              <  1:140818/251‑1 (MQ=255)
                                                         cagcagcaCACCCACGATGGAAAGCCTCAGGCCAAGGATAATCGCGCCGCGGATAATCAGCCCTTTCTTCGGTTTAAGCACGCAAGGATGAGTGATTGACGGGTCAACAATGCGCACATCGACGACGGTGATGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGAATATCAACCTTCTGACCAGACTCGACATAGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGtt                                                                                                                                                                           >  2:419864/1‑250 (MQ=255)
                                                         cagcagcaCACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGCCTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGtt                                                                                                                                                                           >  2:344631/1‑250 (MQ=255)
                                                                acacCCACGATAGAAAGCATAAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGAGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGTTTTTCGGCATCGCCGTTACGCGAcc                                                                                                                                                                   >  2:335915/1‑251 (MQ=255)
                                                                     cACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTa                                                                                                                                                               >  2:104445/1‑250 (MQ=255)
                                                                         aTAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCTACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGttt                                                                                                                                                          <  1:263178/251‑1 (MQ=255)
                                                                                            ggATAATCGCGCCGAGGATAAGCAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAgc                                                                                                                                       <  2:172349/251‑1 (MQ=255)
                                                                                                   cgcgCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGAc                                                                                                                                >  1:445633/1‑251 (MQ=255)
                                                                                                            gATAATCAGCCCTTTCTTCGGTGTCAGCACGCCAGGCTGATTGATTGCCGGGTCAACAATGCGCTCATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGGGCCTGACGTTTCTCca                                                                                                                       <  1:120817/251‑1 (MQ=255)
                                                                                                                  cAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCGCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTAAGACGGACAATCCCCTGCTGGGTTTTCGGCATCGCCGTTACGCGATCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCg                                                                                                                 <  1:74/251‑1 (MQ=255)
                                                                                                                   aGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGATTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGt                                                                                                                >  1:279688/1‑251 (MQ=255)
                                                                                                                          tcttcGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTGTCTCCAGCAGCGTGCGGTa                                                                                                          <  1:335915/250‑1 (MQ=255)
                                                                                                                                                      gATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTgg                                                                             <  1:454064/251‑1 (MQ=255)
                                                                                                                                                         tGCCGGGTCAAGAATGCGCACATCGCCCACGGTGCTTGCCTCGGTGATTTTCAGCTCCGTCTCTTTATTCAGCAGTGTCATATAAGCCGTCGGACCAGTCGCGACATCGCGGGTCAGACGGTCAATCGCCTACTGGTTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTGGGCTTTTTCGTCTTCCAGTTCCTGTCGTCTCTCTAGCAGCGGGCTGTACGCCGGGTGAACtgtgttgttgagcgtggag                                                                           <  2:414098/250‑19 (MQ=255)
                                                                                                                                                                                 gCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgtt                                                  >  1:175773/1‑251 (MQ=255)
                                                                                                                                                                                             ggCCTCGGTGATTTTCAGCTACTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTAAGTTCGTTCAACTgcg                                      >  2:450202/1‑251 (MQ=255)
                                                                                                                                                                                                cTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCAt                                   >  1:345657/1‑251 (MQ=255)
                                                                                                                                                                                                  cggagaatttCAGCACCTGCTCTTTATTCAGTAGTTGCATCCCACCCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATAGCCGTTACGCGAACGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTCTTGTACAGCTTGGAGATTTCTGCCTCTTTAACGGTCAGTTCGTTCAACTGCGCATCg                                 >  2:294513/8‑251 (MQ=255)
                                                                                                                                                                                                      gATTTACATCTCCTTCTCTTCATTCAGCAGTTGCATATAAACCTGCTGCCCCGACCCGACATCGCGGGTCAGACGGACAATCTCCGGTTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGTCGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATg                              <  2:252194/250‑1 (MQ=255)
                                                                                                                                                                                                              gCTCCTGCTCTTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCGTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCatc                     <  1:344631/251‑1 (MQ=255)
                                                                                                                                                                                                                    gCTCTTTATTCATCAGTTGTATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGTGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTTCGCATCGATGTTCACCatcgaatcg               <  2:434144/251‑1 (MQ=255)
                                                                                                                                                                                                                      tctTTATTCAGCAGTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCatcgaatcg               >  2:389098/1‑249 (MQ=255)
                                                                                                                                                                                                                                 cagTTGCATATAAACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTg  >  2:365946/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
AATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTG  >  NC_000913/2134261‑2134736

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