Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 904,944:1 +G 8.0% coding (1200/2160 nt) PP_0788 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947904,9441.G 8.0% 188.0 / 10.9 50coding (1200/2160 nt)PP_0788hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (24/22);  new base G (2/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

ACGACGACAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACTTCCTGGCATTCAATCAGG  >  NC_002947/904802‑905092
                                                                                                                                               |                                                                                                                                                    
acgacgACAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTg                                                                                                                                               <  1:106980/150‑1 (MQ=255)
 cgacgaCAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTgg                                                                                                                                              <  1:446418/150‑1 (MQ=255)
     gacAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTgg                                                                                                                                          >  1:269370/1‑150 (MQ=255)
        aaCCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAg                                                                                                                                       >  2:787481/1‑150 (MQ=255)
           cTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGACTGGGAGTGc                                                                                                                                    <  2:408362/150‑1 (MQ=255)
                 gCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGc                                                                                                                              >  2:215094/1‑150 (MQ=255)
                    agacTGCGGGCGCTGATTACGCGCTTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAg                                                                                                                           >  2:683795/4‑150 (MQ=255)
                     gCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGt                                                                                                                          <  1:60852/150‑1 (MQ=255)
                           gggCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGc                                                                                                                    <  2:957186/150‑1 (MQ=255)
                              cgcTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCgg                                                                                                                                  <  1:817272/133‑1 (MQ=255)
                              cgcTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCgg                                                                                                                                  >  2:817272/1‑133 (MQ=255)
                              cgcTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTg                                                                                                                 <  2:411242/150‑1 (MQ=255)
                               gcTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCG‑CGGCGCCGGCGAtttc                                                                                                                                              >  2:403211/1‑118 (MQ=255)
                               gcTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCG‑CGGCGCCGGCGAtttc                                                                                                                                              <  1:403211/120‑3 (MQ=255)
                                       cgcgCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGc                                                                                                        <  1:180634/150‑1 (MQ=255)
                                        gcgcCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGca                                                                                                       >  1:1092782/1‑150 (MQ=255)
                                            cTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATg                                                                                                   >  1:838674/1‑150 (MQ=255)
                                                        cTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGCGCCGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGt                                                                                       >  1:906083/1‑150 (MQ=255)
                                                            tCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGCGCGGCCGAGAGGGCCAGCGGTAGATCGTGCACATGCGGGTAGAGCGTATCg                                                                                   >  2:1029675/1‑150 (MQ=255)
                                                                ggTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCa                                                                               <  2:5362/150‑1 (MQ=255)
                                                                      aTTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCtt                                                                         >  2:537592/1‑150 (MQ=255)
                                                                            ggACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCaa                                                                   <  1:1028979/150‑1 (MQ=255)
                                                                                  ggggCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAgcgc                                                                           <  1:525839/136‑1 (MQ=255)
                                                                                  ggggCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAgcgc                                                                           >  2:525839/1‑136 (MQ=255)
                                                                                  ggggCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGAttt                                                             >  1:1053313/1‑150 (MQ=255)
                                                                                     gCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGc                                                          >  2:451832/1‑150 (MQ=255)
                                                                                      cTGTTCCGCCGGAGTAGCCCCCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACAGGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGcc                                                         <  1:1029675/150‑1 (MQ=255)
                                                                                           cAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCcgacga                                                    >  1:286920/1‑150 (MQ=255)
                                                                                             gCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGc                                                                                                        <  1:827836/96‑1 (MQ=255)
                                                                                             gCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGc                                                                                                        >  2:827836/1‑96 (MQ=255)
                                                                                                ggAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCa                                               <  2:535787/150‑1 (MQ=255)
                                                                                                 gAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGTGTGCGGCAGCGAGTGCCAGCGTGACAGCtt                                                                                                          <  2:959462/90‑1 (MQ=255)
                                                                                                 gAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGTGTGCGGCAGCGAGTGCCAGCGTGACAGCtt                                                                                                          >  1:959462/1‑90 (MQ=255)
                                                                                                    tAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCt                                           <  2:286920/150‑1 (MQ=255)
                                                                                                      gCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGt                                         >  2:837098/1‑150 (MQ=255)
                                                                                                         aCCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCg                                      >  2:109682/1‑150 (MQ=255)
                                                                                                               cgcgcCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAgc                                <  1:215094/150‑1 (MQ=255)
                                                                                                               cgcgcCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCATGATAATTGATCTGTCCGACGATGCCAGCCTGTTCGCCGAgc                                >  2:820509/1‑150 (MQ=255)
                                                                                                                  gcCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAgctgc                             <  1:109682/150‑1 (MQ=255)
                                                                                                                  gcCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAgctgc                             <  2:838674/150‑1 (MQ=255)
                                                                                                                         gCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTcgacga                      <  1:451832/150‑1 (MQ=255)
                                                                                                                           tCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCGTGCTCAAGTAGGTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACt                    >  1:249282/1‑150 (MQ=255)
                                                                                                                            cAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACtt                   <  2:616903/150‑1 (MQ=255)
                                                                                                                            cAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACtt                   >  1:616903/1‑150 (MQ=255)
                                                                                                                             cagccaaTCG‑CGGCGCCGGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGAtttt                                                            <  1:857108/103‑1 (MQ=255)
                                                                                                                             cagccaaTCG‑CGGCGCCGGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGAtttt                                                            >  2:857108/6‑108 (MQ=255)
                                                                                                                             aaCGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGTGTTACAGCTTGCACATGCGGGGAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACTTc                  >  1:577433/1‑150 (MQ=255)
                                                                                                                             aaCGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTcggcga                      >  1:200067/1‑146 (MQ=255)
                                                                                                                             aaCGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTcgacga                      <  2:200067/146‑1 (MQ=255)
                                                                                                                              aCGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACTTcc                 >  1:28231/1‑150 (MQ=255)
                                                                                                                                            gcc‑gcGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACTTCCTGGCATTCAATCAg   <  2:269370/150‑1 (MQ=255)
                                                                                                                                             cc‑gcGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACTTCCTGGCATTCAATCAgg  <  1:537592/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
ACGACGACAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACTTCCTGGCATTCAATCAGG  >  NC_002947/904802‑905092

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑