Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,316,025 G→T 13.3% A36A (GCG→GCT)  PP_3788 → non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,316,0250GT13.3% 128.3 / 9.9 51A36A (GCG→GCT) PP_3788non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (29/15);  new base T (2/5);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 9.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

AGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCCGTCGT  >  NC_002947/4315878‑4316168
                                                                                                                                                   |                                                                                                                                               
aGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCg                                                                                                                                                 <  6:280843/148‑1 (MQ=255)
        aGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcgat                                                                                                                                        <  6:202834/149‑2 (MQ=255)
        aGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcgat                                                                                                                                        <  5:167807/149‑2 (MQ=255)
        aGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTcgac                                                                                                                                        >  3:325660/1‑149 (MQ=255)
        aGTCTCGTACGCCACACACCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGGCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTcgac                                                                                                                                        <  4:82608/149‑1 (MQ=255)
         gTCTCGTACGCCACACTGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGATCGAATACATCGAGCGTATCATCCACCGCTACCCGTCGAATGCCGCTGCATTcgatg                                                                                                                                       <  5:198287/149‑3 (MQ=255)
         gTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTACCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTcgacg                                                                                                                                       <  4:275383/149‑1 (MQ=255)
            tcGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAg                                                                                                                                     >  4:178428/1‑148 (MQ=255)
                      cacaGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCtt                                                                                                                          <  3:378137/149‑1 (MQ=255)
                       acaGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCtt                                                                                                                          >  1:138436/1‑148 (MQ=255)
                         aGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGAc                                                                                                                       <  2:138436/149‑1 (MQ=255)
                            aaGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGAc                                                                                                                       <  1:160671/146‑1 (MQ=255)
                            aaGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGAc                                                                                                                       >  2:160671/1‑146 (MQ=255)
                            aaGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTa                                                                                                                    >  8:286413/1‑149 (MQ=255)
                                    cagcaTGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACAACGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAAc                                                                                                            >  4:114291/1‑149 (MQ=255)
                                                  cATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTg                                                                                              >  1:121315/1‑149 (MQ=255)
                                                        tGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTcac                                                                                         >  3:129607/1‑148 (MQ=255)
                                                                  gggAGGCCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACCAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTc                                                                              >  4:113597/1‑149 (MQ=255)
                                                                       ggCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcga                                                                                                                                         >  1:368722/1‑85 (MQ=255)
                                                                       ggCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcga                                                                                                                                         <  2:368722/85‑1 (MQ=255)
                                                                        gcctgcctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTg                                                                         >  8:2030/1‑148 (MQ=255)
                                                                          ctgcctACGCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCAt                                                                      <  7:45990/149‑1 (MQ=255)
                                                                           tgcctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGAACGAACTGCTCACACTGCACGCTCAGCTGCAt                                                                      >  3:255151/1‑148 (MQ=255)
                                                                               tACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACt                                                                 >  7:177187/1‑149 (MQ=255)
                                                                                aCCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTc                                                                >  1:149659/1‑149 (MQ=255)
                                                                                  cccGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCgg                                                              >  6:337253/1‑149 (MQ=255)
                                                                                   ccGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGc                                                             <  6:37551/149‑1 (MQ=255)
                                                                                      tGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCTGCGTg                                                          >  7:70223/1‑149 (MQ=255)
                                                                                           ttCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAg                                                     >  2:122013/1‑149 (MQ=255)
                                                                                           ttCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAg                                                     >  2:52397/1‑149 (MQ=255)
                                                                                                     tCGAATACATCGAGCGTATCATCTACCGCTACCCATAGAATGACGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCAGGACCACGGCGTGCGCAGCGGGGAtgtg                                           >  2:247925/1‑149 (MQ=255)
                                                                                                        aaTACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCggg                                                 >  7:118512/1‑140 (MQ=255)
                                                                                                        aaTACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCggg                                                 <  8:118512/140‑1 (MQ=255)
                                                                                                                     tATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAAt                               >  1:303341/1‑145 (MQ=255)
                                                                                                                     tATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAAt                               <  2:303341/145‑1 (MQ=255)
                                                                                                                          tcCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTg                      >  4:314954/1‑149 (MQ=255)
                                                                                                                          tcCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTg                      >  4:136992/1‑149 (MQ=255)
                                                                                                                            cACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGc                     >  3:345919/1‑148 (MQ=255)
                                                                                                                            cACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGc                     <  2:149659/148‑1 (MQ=255)
                                                                                                                               cGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACt                 <  3:314954/149‑1 (MQ=255)
                                                                                                                               cGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACt                 >  5:328555/1‑149 (MQ=255)
                                                                                                                               cGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACt                 <  3:178428/149‑1 (MQ=255)
                                                                                                                                 cTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTgc               >  2:308863/1‑149 (MQ=255)
                                                                                                                                    cccATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTgcgc             <  5:337253/148‑1 (MQ=255)
                                                                                                                                     ccATCGAATTCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGAGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGGGGTTGCGGTGAATGCGACAGTGCAACTGCGCTa           <  1:247925/149‑1 (MQ=255)
                                                                                                                                      cATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTCACCTACAGCGAACTAAGGATCGAAATGATCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACAGCGCTAt          >  1:116637/1‑149 (MQ=255)
                                                                                                                                        tCGAATGCCGCTGCATTCGATGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTgc               <  2:275170/142‑1 (MQ=255)
                                                                                                                                        tCGAATGCCGCTGCATTCGATGAGGACCAGCGCTTGACCTACAGCCAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTgc               >  1:275170/1‑142 (MQ=255)
                                                                                                                                          gAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCcg      >  1:301005/1‑149 (MQ=255)
                                                                                                                                            aTGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCcgt     >  7:108246/1‑148 (MQ=255)
                                                                                                                                            aTGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCcgt     <  8:108246/148‑1 (MQ=255)
                                                                                                                                              gCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCcgtcgt  >  6:302026/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
AGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCCGTCGT  >  NC_002947/4315878‑4316168

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: ‑