Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,316,025 G→T 10.6% 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,0250GT10.6% 149.1 / 7.2 54A36A (GCG→GCT) PP_3788non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (24/24);  new base T (3/3);  total (27/27)
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.91e-01

AATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCCGTCG  >  NC_002947/4315882‑4316167
                                                                                                                                               |                                                                                                                                               
aaTCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCAtt                                                                                                                                            <  6:104768/149‑1 (MQ=255)
  tCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTcg                                                                                                                                          >  5:259161/1‑149 (MQ=255)
    aGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCTTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTcgac                                                                                                                                        <  1:44153/149‑1 (MQ=255)
     gTCTGGTACGCCACACATCAAGGAGCGGAGCATGTACCAGTCATTACTGACGATCACGGGAGGTCCTGACTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcgat                                                                                                                                        <  8:122626/148‑2 (MQ=255)
      tctcGTACGCCACACAGCAATGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCgg                                                                                                                                                <  4:310944/139‑1 (MQ=255)
      tctcGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCgg                                                                                                                                                >  3:310944/1‑139 (MQ=255)
        tcGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGATTACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAgg                                                                                                                                    >  8:75671/1‑149 (MQ=255)
          gTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAgg                                                                                                                                    <  5:166688/147‑1 (MQ=255)
             cGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCa                                                                                                                                <  6:122455/148‑1 (MQ=255)
              gCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAgc                                                                                                                              <  6:142705/149‑1 (MQ=255)
                cacacaGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcgat                                                                                                                                        >  7:223657/1‑136 (MQ=255)
                cacacaGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcgat                                                                                                                                        <  8:223657/137‑2 (MQ=255)
                     aGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGAc                                                                                                                       >  3:336982/1‑149 (MQ=255)
                         aGGAGCGCAGCATGAACCAGGCATTTCTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCa                                                                                                                                <  6:352436/136‑1 (MQ=255)
                         aGGAGCGCAGCATGAACCAGGCATTTCTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCa                                                                                                                                >  5:352436/1‑136 (MQ=255)
                           gAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAg                                                                                                                 <  3:235309/149‑1 (MQ=255)
                                    aTGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAg                                                                                                        >  7:23373/1‑149 (MQ=255)
                                         ccAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATc                                                                                                    >  7:280137/1‑148 (MQ=255)
                                          cAGGCATTACTGACGATCACGGGAGGGCATGCCTACCCCGTGCACTTCTCGCTTGTCGAATACATCGAGCGTATCATCCACCTCTACCCATCGATTGCAGCGTCATTCTACGAGGACCAGCGCTTGACCTACAGCGAACTaa                                                                                                         <  3:28959/142‑1 (MQ=255)
                                          cAGGCATTACTGACGATCACGGGAGGGACTGCCTACACCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCAACCACCGCTACACATCGAATGCAGCGGCATTCAACGAGGACCAGCGCTTGACCTACAGCGAACTAAgacatag                                                                                                  >  4:28959/1‑143 (MQ=255)
                                               aTTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGc                                                                                            <  4:336982/149‑1 (MQ=255)
                                                ttACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTAc                                                                                                                   <  6:321624/126‑1 (MQ=255)
                                                      aCGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACt                                                                                     <  6:259161/149‑1 (MQ=255)
                                                            aCGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGc                                                                               >  5:199834/1‑149 (MQ=255)
                                                               ggAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCg                                                                            >  6:115751/1‑149 (MQ=255)
                                                                    gcctgcctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGc                                                                       <  6:199834/149‑1 (MQ=255)
                                                                         cctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGaa                                                                  >  5:130799/1‑149 (MQ=255)
                                                                           tACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACt                                                                >  8:229101/1‑149 (MQ=255)
                                                                             ccccGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCg                                                              <  6:130799/149‑1 (MQ=255)
                                                                               ccGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGc                                                            >  1:12804/1‑149 (MQ=255)
                                                                               ccGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGc                                                            <  2:12804/149‑1 (MQ=255)
                                                                                cGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCAt                                                                     <  1:213776/139‑1 (MQ=255)
                                                                                cGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCAt                                                                     >  2:213776/1‑139 (MQ=255)
                                                                                     aCTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGaa                                                                  >  8:210123/1‑137 (MQ=255)
                                                                                            tcTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACCCGGCGTGCGCAGCggg                                                >  8:3825/2‑148 (MQ=255)
                                                                                            gCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGAc                                                                                                                       <  4:27572/78‑1 (MQ=255)
                                                                                            gCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGAc                                                                                                                       >  3:27572/1‑78 (MQ=255)
                                                                                            gCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCg                                                  <  5:115751/146‑1 (MQ=255)
                                                                                                  cGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGAt                                             >  7:299920/1‑145 (MQ=255)
                                                                                                  cGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGAt                                             <  8:299920/145‑1 (MQ=255)
                                                                                                     aaaCATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TCATCACACTGCACGCTCGGCTGCATGAACTCCCCGTCCGCAGCGGGGATGTCGTTg                                      >  2:26259/3‑149 (MQ=255)
                                                                                                         aTCGAGCATATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGt                                  >  8:15710/1‑149 (MQ=255)
                                                                                                               cGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGc                            <  7:229101/149‑1 (MQ=255)
                                                                                                                   tcatcCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACCCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACa                        <  7:3825/149‑1 (MQ=255)
                                                                                                                       ccACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTg                     >  2:151748/1‑148 (MQ=255)
                                                                                                                         aCCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCgg                                                 >  7:248374/1‑118 (MQ=255)
                                                                                                                         aCCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCgg                                                 <  8:248374/118‑1 (MQ=255)
                                                                                                                                 ccATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACa                                                                                                                  <  4:286607/46‑1 (MQ=255)
                                                                                                                                 ccATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACa                                                                                                                  >  3:286607/1‑46 (MQ=255)
                                                                                                                                 ccATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACa                        >  4:61626/1‑135 (MQ=255)
                                                                                                                                 ccATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACa                        <  3:61626/135‑1 (MQ=255)
                                                                                                                                  cgtCGAATGCCGCTGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATg                             <  8:178924/127‑1 (MQ=255)
                                                                                                                                  cgtCGAATGCCGCTGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATg                             >  7:178924/3‑129 (MQ=255)
                                                                                                                                  cATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTAt         >  7:71884/1‑149 (MQ=255)
                                                                                                                                    tCGAATGCCGCTGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACGAGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATc        >  1:274409/1‑149 (MQ=255)
                                                                                                                                         tGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCcgtcg  <  8:71884/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
AATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCCGTCG  >  NC_002947/4315882‑4316167

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