Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,316,025 G→T 9.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,0250GT 9.3% 159.0 / 6.7 58A36A (GCG→GCT) PP_3788non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (27/25);  new base T (2/3);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.95e-01

AGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCCGTCG  >  NC_002947/4315878‑4316167
                                                                                                                                                   |                                                                                                                                               
aGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCgg                                                                                                                                                >  7:258550/1‑149 (MQ=255)
        aGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTcga                                                                                                                                         <  4:176903/148‑1 (MQ=255)
        aGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACACCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTcgac                                                                                                                                        >  7:150170/1‑149 (MQ=255)
         gTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcgatg                                                                                                                                       <  7:259969/149‑3 (MQ=255)
         gTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcgatg                                                                                                                                       <  2:5907/149‑3 (MQ=255)
           ctcGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTcgacga                                                                                                                                      <  2:12046/148‑1 (MQ=255)
              gTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGAc                                                                                                                                  >  2:6890/1‑149 (MQ=255)
              gTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAAGACATCGAGCGTATCCTCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGAc                                                                                                                                  >  4:276410/1‑149 (MQ=255)
               tACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTAGCTGGTCGAATACATCGAGCGTATCATCCACCGCTAACCATCGAATGCAGCGGCCTTCGACGAGGAcc                                                                                                                                 >  8:145689/1‑149 (MQ=255)
                 cGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCa                                                                                                                                <  8:150170/148‑1 (MQ=255)
                     acacaGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCt                                                                                                                           >  6:123721/1‑149 (MQ=255)
                       acaGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCtt                                                                                                                          >  1:56092/1‑148 (MQ=255)
                        caGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGa                                                                                                                        <  8:258550/149‑1 (MQ=255)
                        caGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGa                                                                                                                        >  4:305828/1‑149 (MQ=255)
                         aGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGAc                                                                                                                       <  1:6890/149‑1 (MQ=255)
                          gCAAGGAGCGCAGCATGAACCAGGCATTACTGACGCTCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGAc                                                                                                                       >  5:202030/1‑148 (MQ=255)
                                      gcaTGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTAc                                                                                                                   >  4:201812/1‑140 (MQ=255)
                                      gcaTGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTAc                                                                                                                   <  3:201812/140‑1 (MQ=255)
                                             ccAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTTGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATc                                                                                                    >  5:39493/1‑148 (MQ=255)
                                              cAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCg                                                                                                   >  4:305133/1‑148 (MQ=255)
                                                  cATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTg                                                                                              <  2:93452/149‑1 (MQ=255)
                                                                aCGGGAGGGCCTGCCTACCCCGTGCACTTCTTGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGc                                                                                <  6:39493/149‑1 (MQ=255)
                                                                      gggCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGaa                                                                                                             >  4:170085/1‑114 (MQ=255)
                                                                      gggCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGaa                                                                                                             <  3:170085/114‑1 (MQ=255)
                                                                       ggCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCt                                                                                                                     >  2:284944/1‑105 (MQ=255)
                                                                       ggCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCt                                                                                                                     <  1:284944/105‑1 (MQ=255)
                                                                        gcctgcctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGc                                                                        <  2:56092/149‑1 (MQ=255)
                                                                        gcctgcctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGc                                                                        <  5:123721/149‑1 (MQ=255)
                                                                          ctgcctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCAt                                                                      >  5:21362/1‑149 (MQ=255)
                                                                             cctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGaa                                                                   >  6:265470/1‑149 (MQ=255)
                                                                                 ccccGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATGGAATGCAGCGGCATTCGACGAGTACCAGCGGTTGAGCTTCAGCGAACTAAGGATCTAACTGCTCACACTGCACGCTCGGCTGCATTAACTCg                                                               <  1:234092/149‑1 (MQ=255)
                                                                                  cccGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACTCTCCGCTGCATGAACTCgg                                                              >  6:203362/1‑149 (MQ=255)
                                                                                  cccGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATAGAATGCAGCGGCATTAGACGAGGACCAGCGCTTGACCTACAGCGAAATAAGGATCGAACTGATCACACTGCAAGCTACGCTGCATGAACTCg                                                               >  3:256926/1‑148 (MQ=255)
                                                                                           ttCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAg                                                    >  4:53803/1‑149 (MQ=255)
                                                                                           ttCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAg                                                    >  2:163340/1‑149 (MQ=255)
                                                                                             ctcGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTg                                                                         <  8:247608/127‑1 (MQ=255)
                                                                                             ctcGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTg                                                                         >  7:247608/1‑127 (MQ=255)
                                                                                                 cTGGTCGAATACATCGACCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGACCTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGa                                              <  1:129341/149‑1 (MQ=255)
                                                                                                     tCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGAtgtg                                          <  6:21362/149‑1 (MQ=255)
                                                                                                     tCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGAtgtg                                          >  1:124003/1‑149 (MQ=255)
                                                                                                       gAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTgg                                         <  5:265470/148‑1 (MQ=255)
                                                                                                             aTCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGt                                  <  3:305133/149‑1 (MQ=255)
                                                                                                             aTCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGt                                  >  7:299971/1‑149 (MQ=255)
                                                                                                                gAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGaa                               <  3:305828/149‑1 (MQ=255)
                                                                                                                 atcGTATCATCCACCGCTAACCATCGAATGCAGCGGCATTCGAAGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGTCGCAGCGGGGATGTGGTTGCGGTGaa                               <  1:163340/147‑1 (MQ=255)
                                                                                                                 aGCGTATCATCCACCGCTACCCATCGAATGCAGCAGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGAAt                              <  3:143503/149‑1 (MQ=255)
                                                                                                                          tcCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCg                                    >  4:102375/1‑134 (MQ=255)
                                                                                                                          tcCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCg                                    <  3:102375/134‑1 (MQ=255)
                                                                                                                            cACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGACTGCGACAGTGc                    >  7:291980/1‑148 (MQ=255)
                                                                                                                                       gtCGAATGCCGCTGCATTCGATGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGt                                  <  1:115814/122‑1 (MQ=255)
                                                                                                                                       gtCGAATGCCGCTGCATTCGATGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGt                                  >  2:115814/2‑123 (MQ=255)
                                                                                                                                       aTCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACGTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCTGCGTG‑CGCAGCGGGGATGTGGTTGCGGt                                  <  2:85386/123‑1 (MQ=255)
                                                                                                                                       aTCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATc        <  8:299971/149‑1 (MQ=255)
                                                                                                                                       aTCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATc        <  4:256926/149‑1 (MQ=255)
                                                                                                                                       aTCGAATGCAGCGGCATTCGACGAGGAACAGCGCTTGACCTACAGCCAACTAAGGATCGAACTCCTCACACTGCACGCTCCGCTGCATGAACTCGGCCTG‑CGCAGCGGGGATGTGGTTGCGGt                                  >  1:85386/1‑123 (MQ=255)
                                                                                                                                        tCGAATGCCGCTGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATcc       >  4:34289/1‑149 (MQ=255)
                                                                                                                                          gAATGCAGCGGCATTCGACGAGGACCAGGGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCcg     >  1:177099/1‑149 (MQ=255)
                                                                                                                                                aGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCcgtcg  <  8:102613/146‑1 (MQ=255)
                                                                                                                                                aGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCcgtcg  >  7:102613/1‑146 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
AGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCCGTCG  >  NC_002947/4315878‑4316167

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