Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,316,022 A→C 10.2% A35A (GCA→GCC)  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,0220AC10.2% 161.6 / 6.3 57A35A (GCA→GCC) PP_3788non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (26/25);  new base C (3/3);  total (29/28)
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

TCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCCGTCG  >  NC_002947/4315876‑4316167
                                                                                                                                                  |                                                                                                                                                  
tcaGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGc                                                                                                                                                  >  5:122455/1‑149 (MQ=255)
      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                                                                                                                                        <  8:223657/137‑2 (MQ=255)
                      cacacaGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcgat                                                                                                                                        >  7:223657/1‑136 (MQ=255)
                           aGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGAc                                                                                                                       >  3:336982/1‑149 (MQ=255)
                               aGGAGCGCAGCATGAACCAGGCATTTCTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCa                                                                                                                                >  5:352436/1‑136 (MQ=255)
                               aGGAGCGCAGCATGAACCAGGCATTTCTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCa                                                                                                                                <  6:352436/136‑1 (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                                                                                                                       >  3:27572/1‑78 (MQ=255)
                                                                                                  gCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGAc                                                                                                                       <  4:27572/78‑1 (MQ=255)
                                                                                                  gCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCg                                                  <  5:115751/146‑1 (MQ=255)
                                                                                                        cGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGAt                                             <  8:299920/145‑1 (MQ=255)
                                                                                                        cGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGAt                                             >  7:299920/1‑145 (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                                                                                                                  >  3:286607/1‑46 (MQ=255)
                                                                                                                                       ccATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACa                                                                                                                  <  4:286607/46‑1 (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)
                                                                                                                                                  |                                                                                                                                                  
TCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAAC‑TGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGCGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCCGTCG  >  NC_002947/4315876‑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: ‑