Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,267,497 T→A 9.8% L28Q (CTG→CAG)  PP_4643 → xanthine/uracil permease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,267,4970TA 9.8% 145.8 / 3.8 52L28Q (CTG→CAG) PP_4643xanthine/uracil permease family protein
Reads supporting (aligned to +/- strand):  ref base T (24/22);  new base A (3/2);  total (27/25)
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

TAGGAAGTGTCTACAGTAGCGGCACAACAATAAACAGAGAGCCTGCTCCATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCGGCATGCTGATCGCCGCCTCGCTGTTCGCCGGTGGCCTTGCCACCCTGCTGCAGACCCTTGGC  >  NC_002947/5267366‑5267645
                                                                                                                                   |                                                                                                                                                     
tAGGAAGTGTCTACAGTAGCGGCACAACAATAAACAGAGAGCCTGCTCCATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCa                                                                                                                                     >  3:171881/1‑149 (MQ=255)
   gAAGTGTCTACCGTAGCGGCACAACCATAAACAGAGAGCCTGCTCCATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGt                                                                                                                                  <  6:176332/149‑1 (MQ=255)
       tgtCTACAGTAGCGGCACAACAATAAACAGAGAGCCTGCTCCATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTAtgg                                                                                                                              <  3:108871/149‑1 (MQ=255)
         tCTACAGTAGCGGCACAACAATAAACAGAGAGCCTGCTCCATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAACACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTAtggtg                                                                                                                            <  2:257913/149‑1 (MQ=255)
         tCTACAGTAGCGGCACAACAATAAACAGAGAGCCTGCTCCATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTAtggt                                                                                                                             <  3:2000/148‑1 (MQ=255)
           tACAGTAGCGGCACAACAATAAACAGAGAGCCTGCTCCATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGt                                                                                                                                  >  6:77319/1‑141 (MQ=255)
           tACAGTAGCGGCACAACAATAAACAGAGAGCCTGCTCCATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGt                                                                                                                                  <  5:77319/141‑1 (MQ=255)
           tACAGTAGCGGCACAACAATAAACAGAGAGCCTGCTCCATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTAtggtgg                                                                                                                           >  7:222762/1‑148 (MQ=255)
            aCAGTAGCGGCACAACAATAAACAGAGAGCCTGCTCCATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGa                                                                                                                         >  3:263581/1‑149 (MQ=255)
                       acaacaATAAACAGAGAGCCTGCTCCATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTa                                                                                                              >  8:63317/1‑149 (MQ=255)
                         aacaaTAAACAGAGAGCCTGCTCCATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTAcc                                                                                                            >  4:41074/1‑149 (MQ=255)
                         aacaaTAAACAGAGAGCCTGCTCCATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTAcc                                                                                                            >  4:160607/1‑149 (MQ=255)
                           caaTAAACAGAGAGCCTGCTCCATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTAcccc                                                                                                          <  4:263581/149‑1 (MQ=255)
                             aTAAACAGAGAGCCTGCTCCATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCt                                                                                                         >  5:38327/1‑148 (MQ=255)
                                 aCAGAGAGCCTACTCCATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATcc                                                                                                    >  6:237267/1‑149 (MQ=255)
                                   agagagCCTGCTCCATGACTACCTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTg                                                                                                  <  8:222762/149‑1 (MQ=255)
                                             cTCCATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCcg                                                                                        >  1:172839/1‑149 (MQ=255)
                                                cATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCAGCACATGCTGCAGACCa                                                                                                                                     >  1:124625/1‑101 (MQ=255)
                                                cATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCAGCACATGCTGCAGACCa                                                                                                                                     <  2:124625/101‑1 (MQ=255)
                                                      tACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTg                                                                                >  4:60673/1‑148 (MQ=255)
                                                        cGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATTGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAAt                                                                             <  3:60673/149‑1 (MQ=255)
                                                        cGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGaa                                                                              <  6:38327/148‑1 (MQ=255)
                                                         gTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCAGCACATGCTGCAGACCAtag                                                                                                                                  >  5:63637/1‑93 (MQ=38)
                                                         gTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCAGCACATGCTGCAGACCAtag                                                                                                                                  <  6:63637/95‑3 (MQ=38)
                                                          tCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATgg                                                                           >  7:114505/1‑149 (MQ=255)
                                                           ccccTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGt                                                                          <  7:63317/149‑1 (MQ=255)
                                                             ccTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCAGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGc                                                                        >  4:298311/1‑149 (MQ=255)
                                                              cTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGt                                                                          <  6:111846/146‑1 (MQ=255)
                                                              cTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGt                                                                          >  5:111846/1‑146 (MQ=255)
                                                                  cACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCACCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGaa                                                                    >  3:27799/1‑148 (MQ=255)
                                                                     gTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCg                                                                >  3:55084/1‑149 (MQ=255)
                                                                     gTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCg                                                                <  5:279028/149‑1 (MQ=255)
                                                                        tcCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCGGCa                                                             >  2:250365/1‑149 (MQ=255)
                                                                         cccGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCGGCAt                                                            >  4:336425/1‑149 (MQ=255)
                                                                          ccGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCGGc                                                              <  3:298311/146‑1 (MQ=255)
                                                                          ccGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCGGCATg                                                           >  3:121513/1‑149 (MQ=255)
                                                                             ccAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCGGCATGCTg                                                        <  3:160607/149‑1 (MQ=255)
                                                                                      ccGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTg                                                                         >  6:315136/1‑123 (MQ=255)
                                                                                         aGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCGGCATGCTGATCGCCGCCTCg                                            <  1:250365/149‑1 (MQ=255)
                                                                                               aaaaCCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCGGCATGCTGATCGCCGCCTCGCTGTTc                                      >  3:1370/1‑149 (MQ=255)
                                                                                                aaaCCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCGGCATGCTGATCGCAGCCTCGCTGTTCg                                     >  3:173510/1‑149 (MQ=255)
                                                                                                  aCCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCGGCATGCTGATCGCCGCCTCGCTGTTCGcc                                   <  3:41074/149‑1 (MQ=255)
                                                                                                       ggCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCGGCATGCTGATCGCCGCCTCGCTGTTCGCCGGTgg                              >  1:295927/1‑149 (MQ=255)
                                                                                                         ccTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCGGCATGCTGATCGCCGCCTCGCTGTTCGCCGGTGGcc                            <  8:114505/149‑1 (MQ=255)
                                                                                                                  cAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGAtcgtc                                                                                                                    >  1:304040/1‑52 (MQ=255)
                                                                                                                  cAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGAtcgtc                                                                                                                    <  2:304040/52‑1 (MQ=255)
                                                                                                                   aaCCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCGGCATGCTGATCGCCGCCTCGCTGTTCGCCGGTGGCCTTGCCACCct                  <  4:173510/149‑1 (MQ=255)
                                                                                                                         gCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCGGCATGCTGATCGCCGCCTCGCTGTTCGCCGGTGGCCTTGCCACCctgctg              <  4:121513/147‑1 (MQ=255)
                                                                                                                           cTATGGTCCGCACCATGTGCTGCCCATGTATGGAGCGATCGTCCCAGTACCCCTGATCCTTTGGAATGCCTCCGGCCTGAATGGTGCCGAACACTGCATGCTTCTCGCCGCCTCTCTTTTCGCCGGTGGCCTTGCCACCCTGCTGCAGa          <  6:294437/149‑1 (MQ=255)
                                                                                                                                gTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCGGCATGCTGATCGCCGCCTCGCTGTTCGCCGGTGGCCTTGCCACCCTGCTGCAGACCCtt     <  2:172839/149‑1 (MQ=255)
                                                                                                                                 tCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCGGCATGCTGATCGCCGCCTCGCTGTTCGCCGGTGGCCTTGCCACCCTGCTGCAGACCCTTg    <  5:237267/149‑1 (MQ=255)
                                                                                                                                 tCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCGGCATGCTGATCGCCGCCTCGCTGTTCGCCGGTGGCCTTGCCACCCTGCTGCAGACACTTg    >  6:352082/1‑149 (MQ=255)
                                                                                                                                   tGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCGGCATGCTGATCGCCGCCTCGCTGTTCGCCGGTGGCCTTGCCACCCGGCTGCAGACCCTTGGc  >  4:272139/1‑149 (MQ=255)
                                                                                                                                   |                                                                                                                                                     
TAGGAAGTGTCTACAGTAGCGGCACAACAATAAACAGAGAGCCTGCTCCATGACTACGTCCCCTAGCACGTCTCCCGCCAAGCGACCCGAGGATGAAAACCTCGGCCTTGGTGCCAACCTGGCCTATGGTCTGCAGCATGTGCTGACCATGTATGGTGGGATCGTCGCAGTACCCCTGATCCTGGGGCAGGCCGCGGGCCTGAATGGTGCCGAAATCGGCATGCTGATCGCCGCCTCGCTGTTCGCCGGTGGCCTTGCCACCCTGCTGCAGACCCTTGGC  >  NC_002947/5267366‑5267645

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