Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,641,817 T→A 9.4% Q463L (CAG→CTG)  PP_3209 ← isoxanthopterin deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,641,8170TA 9.4% 149.2 / 4.7 53Q463L (CAG→CTG) PP_3209isoxanthopterin deaminase
Reads supporting (aligned to +/- strand):  ref base T (24/24);  new base A (3/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

CCAGTCGAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGACCGAAG  >  NC_002947/3641673‑3641963
                                                                                                                                                |                                                                                                                                                  
ccAGTCGAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCg                                                                                                                                                >  6:202081/1‑149 (MQ=255)
ccAGTCGAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCAGGca                                                                                                                                                >  7:223621/1‑148 (MQ=255)
      gAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTg                                                                                                                                          >  3:234556/1‑149 (MQ=255)
       aaGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTgg                                                                                                                                         <  5:202081/149‑1 (MQ=255)
        aGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGc                                                                                                                                        <  6:312587/149‑1 (MQ=255)
         gggCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCAGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCAGGCACGCCAGGCt                                                                                                                                       <  8:223621/149‑1 (MQ=255)
                   gCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCg                                                                                                                             >  8:98946/1‑149 (MQ=255)
                   gCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCg                                                                                                                             <  3:148500/149‑1 (MQ=255)
                    cAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCgg                                                                                                                            <  7:295177/149‑1 (MQ=255)
                       gCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGc                                                                                                                           >  3:191186/1‑147 (MQ=255)
                               ggggCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAg                                                                                                                 <  8:185028/149‑1 (MQ=255)
                                gggCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGATGTACAGCCAGGCACGCCAGGCTgta                                                                                                                                    <  1:65971/129‑3 (MQ=255)
                                gggCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGATGTACAGCCAGGCACGCCAGGCTgta                                                                                                                                    >  2:65971/1‑127 (MQ=255)
                                  gCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCa                                                                                                              <  5:154216/149‑1 (MQ=255)
                                         ccACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATc                                                                                                        <  3:149681/148‑1 (MQ=255)
                                          cacaCCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCg                                                                                                       >  7:57604/1‑148 (MQ=255)
                                               ccAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGt                                                                                                 >  5:179173/1‑149 (MQ=255)
                                                cAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGtt                                                                                                >  3:283094/1‑149 (MQ=255)
                                                 aGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGttt                                                                                               >  1:45833/1‑149 (MQ=255)
                                                  gCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGtttt                                                                                              <  2:45833/149‑1 (MQ=255)
                                                    aCGTCGTCTGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTcc                                                                                            <  3:293995/149‑1 (MQ=255)
                                                    aCGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTcc                                                                                            <  8:315170/149‑1 (MQ=255)
                                                      gtcgtcGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCAc                                                                                          <  4:234556/149‑1 (MQ=255)
                                                               ttCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAggcg                                                                                 >  2:180877/1‑149 (MQ=255)
                                                                        cAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAg                                                                                    >  5:70227/1‑137 (MQ=255)
                                                                        cAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAg                                                                                    <  6:70227/137‑1 (MQ=255)
                                                                            gCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAg                                                                     >  2:264113/1‑148 (MQ=255)
                                                                             cccGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGc                                                                    >  8:156884/1‑148 (MQ=255)
                                                                                gCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGCACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAg                                                                                                                        >  3:279271/1‑93 (MQ=255)
                                                                                gCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGCACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAg                                                                                                                        <  4:279271/93‑1 (MQ=255)
                                                                                  aCCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTACAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGcc                                                              >  3:160079/1‑149 (MQ=255)
                                                                                    cTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTcc                                                                                            <  4:255760/117‑1 (MQ=255)
                                                                                    cTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTcc                                                                                            >  3:255760/1‑117 (MQ=255)
                                                                                    cTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCg                                                            >  6:300287/1‑149 (MQ=255)
                                                                                      gCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCg                                                            >  3:115575/1‑147 (MQ=255)
                                                                                        agagGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGc                                                        >  4:285585/1‑149 (MQ=255)
                                                                                        agagGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGc                                                        <  7:98946/149‑1 (MQ=255)
                                                                                             cGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCgg                                                   <  8:314728/149‑1 (MQ=255)
                                                                                                           gcgTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCaa                                     <  7:209680/149‑1 (MQ=255)
                                                                                                             gTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAAc                                    >  5:225773/1‑148 (MQ=255)
                                                                                                                    cgcgGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTgggg                            <  5:300287/149‑1 (MQ=255)
                                                                                                                           tGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCg                          >  3:284750/1‑144 (MQ=255)
                                                                                                                           tGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCg                          <  4:284750/144‑1 (MQ=255)
                                                                                                                                 aGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCa                   >  5:149788/1‑145 (MQ=255)
                                                                                                                                 aGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCa                   <  6:149788/145‑1 (MQ=255)
                                                                                                                                  ggCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCCCGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCa              >  6:79855/1‑149 (MQ=255)
                                                                                                                                  ggCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCa              <  4:115575/149‑1 (MQ=255)
                                                                                                                                  ggCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCa              <  3:285585/149‑1 (MQ=255)
                                                                                                                                   gCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCAt             <  4:283094/149‑1 (MQ=255)
                                                                                                                                        gggcAGCCAGGCACGCCAGGCTGCCCAGTTCGGCGAGGTCCAAGACAGGAATCGCCACGTGTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGACGTACGGACACCGCCAACCAATTGGGCCATGCCCAGATCAAGCAGa        >  7:285398/4‑149 (MQ=255)
                                                                                                                                         tACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGAc       <  6:225773/149‑1 (MQ=255)
                                                                                                                                         tACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGAc       >  4:147590/1‑149 (MQ=255)
                                                                                                                                              ccTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGACCGAAg  >  8:126977/1‑149 (MQ=255)
                                                                                                                                                tGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGACCGAAg  >  5:318749/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CCAGTCGAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGACCGAAG  >  NC_002947/3641673‑3641963

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