Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,601,818 Δ1 bp 7.8% coding (1664/5418 nt) PP_4924 ← subtilase family serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,601,8180G. 7.8% 160.1 / 7.6 51coding (1664/5418 nt)PP_4924subtilase family serine protease
Reads supporting (aligned to +/- strand):  ref base G (22/25);  new base . (2/2);  total (24/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 = 9.48e-01

GCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGC  >  NC_002947/5601673‑5601957
                                                                                                                                                  |                                                                                                                                             
gCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTgcg                                                                                                                                            <  5:58820/148‑1 (MQ=255)
 cTGTCACTGATCACGAATACTTCGTCGCCCGCGCCGCCGTTTTAGGTGTTTCTGACACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGTTCAAGACATCTTCGCCCGCGCTGCCGTATATCGTCTTTACGCCGAGCGCTgcggc                                                                                                                                          <  8:161829/149‑1 (MQ=255)
 cTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTgcggc                                                                                                                                          <  7:300774/149‑1 (MQ=255)
  tGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTgcggcg                                                                                                                                         >  2:357624/1‑149 (MQ=255)
      aCTGATAACGAATACATCGTCGACAGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACAAGCCAACTGCCCTCTGGTTCCGCGGTCAAGACATCATCGCCCGCGCTCCCGTAAAGCGTCTGTACGACGAGCGCTGCGGCGAcc                                                                                                                                      >  4:123908/1‑148 (MQ=255)
          aTCACGACTACATCGTCGCCCGCCCCCCCGTGATAGGTGGGGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGGCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGt                                                                                                                                 <  6:19683/149‑1 (MQ=255)
           tCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGtg                                                                                                                                >  1:26883/1‑149 (MQ=255)
             aCGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGTCATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGc                                                                                                                                   >  1:140846/1‑144 (MQ=255)
             aCGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGtgct                                                                                                                              >  2:269443/1‑149 (MQ=255)
                  tACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCtt                                                                                                                        <  4:316657/149‑1 (MQ=255)
                  tACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCtt                                                                                                                        <  4:103494/149‑1 (MQ=255)
                     aTCGTCGCCCGCCCCGCCGTGATAGGTGTTGCAGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCAGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCCGCGTCCAGCGTGCTGCCTTGAg                                                                                                                     >  5:244719/1‑149 (MQ=255)
                        gtcgCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCa                                                                                                                   <  1:319749/148‑1 (MQ=255)
                           gcccgcccCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCActct                                                                                                               <  1:269443/149‑1 (MQ=255)
                           gcccgcccCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCActc                                                                                                                <  1:218490/148‑1 (MQ=255)
                              cgccccgccGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGc                                                                                                            >  2:89741/1‑149 (MQ=255)
                                ccccgccGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTg                                                                                                          <  2:85501/149‑1 (MQ=255)
                                   cgccgTGATAGGGGGTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCAGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATg                                                                                                       <  6:244719/149‑1 (MQ=255)
                                     ccgTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTa                                                                                                     <  1:58359/149‑1 (MQ=255)
                                              gtgtTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTgcg                                                                                                                                            <  7:138184/102‑1 (MQ=255)
                                              gtgtTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTgcg                                                                                                                                            >  8:138184/1‑102 (MQ=255)
                                                 ttGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCggg                                                                                          <  1:89741/149‑1 (MQ=255)
                                                 ttGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCggg                                                                                          >  2:70983/1‑148 (MQ=255)
                                                     agccaccgccaccAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTg                                                                                     >  2:283376/2‑149 (MQ=255)
                                                       ccaccgccaccAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTgg                                                                                    <  2:320295/148‑1 (MQ=255)
                                                         accgccaccAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGCAGCGCTCGGCGtac                                                                                                                                      <  6:291042/97‑3 (MQ=255)
                                                         accgccaccAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGCAGCGCTCGGCGtac                                                                                                                                      >  5:291042/1‑95 (MQ=255)
                                                         accgccaccAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCa                                                                                                                   <  1:95068/115‑1 (MQ=255)
                                                         accgccaccAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCa                                                                                                                   >  2:95068/1‑115 (MQ=255)
                                                                ccaaccaGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGa                                                                          <  1:357624/149‑1 (MQ=255)
                                                                  aaccaGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGAc                                                                         >  2:7328/1‑148 (MQ=255)
                                                                     caGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGATCACTCTGGCTGATGTAGGCCCTTCGTCAAGCGGGCCCACTTCCCGAAg                                                                     <  3:123908/149‑1 (MQ=255)
                                                                        ccAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGGCGGCTCTGCCGTAAAGCGTCTTTACGCGGAGCGCTGCTGCGTCCAGCGTTCTGCCTTCAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCa                                                                  <  7:129627/149‑1 (MQ=255)
                                                                        ccAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCa                                                                  <  1:7328/149‑1 (MQ=255)
                                                                           aCTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGGACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGATGTCAg                                                                 >  8:270916/1‑147 (MQ=255)
                                                                           aCTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTa                                                               <  1:283376/149‑1 (MQ=255)
                                                                                       ggCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTa                                                                                                     >  5:33787/1‑99 (MQ=255)
                                                                                       ggCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACcct                                                   >  1:145381/1‑149 (MQ=255)
                                                                                       ggCCGCGGGCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTa                                                                                                     <  6:33787/99‑1 (MQ=255)
                                                                                        gCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGa                                                                          <  1:250362/125‑1 (MQ=255)
                                                                                        gCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGa                                                                          >  2:250362/1‑125 (MQ=255)
                                                                                         ccGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTTATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGaa                                                                      <  5:230475/128‑1 (MQ=255)
                                                                                         ccGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGaa                                                                      >  6:230475/1‑128 (MQ=255)
                                                                                           gcgGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTg                                               >  1:296881/1‑149 (MQ=255)
                                                                                              gTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGc                                            <  2:26883/149‑1 (MQ=255)
                                                                                                              ccGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTAcc                             >  1:89478/1‑149 (MQ=255)
                                                                                                              ccGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTg                                                                                                          >  2:295482/1‑71 (MQ=255)
                                                                                                              ccGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTg                                                                                                          <  1:295482/71‑1 (MQ=255)
                                                                                                                      ccGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCa                                                                  >  7:249715/1‑103 (MQ=255)
                                                                                                                      ccGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCa                                                                  <  8:249715/103‑1 (MQ=255)
                                                                                                                           aaGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCga               <  4:288214/149‑1 (MQ=255)
                                                                                                                                tCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAAc          <  2:89478/149‑1 (MQ=255)
                                                                                                                                tCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAAc          >  3:346906/1‑149 (MQ=255)
                                                                                                                                  tggaCGCCGCAGCGCTCGGCGTACAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCga               >  6:328379/4‑142 (MQ=255)
                                                                                                                                  tggaCGCCGCAGCGCTCGGCGTACAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCga               <  5:328379/139‑1 (MQ=255)
                                                                                                                                         cGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGc  >  8:18673/1‑148 (MQ=255)
                                                                                                                                          gAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGc  >  3:241597/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGAGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACCGC  >  NC_002947/5601673‑5601957

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: