Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 904,937 Δ1 bp 8.0% coding (1193/2160 nt) PP_0788 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947904,9370C. 8.0% 154.8 / 8.2 50coding (1193/2160 nt)PP_0788hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/22);  new base . (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 = 8.82e-01

GCGTTCATCCACAACGACGACAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACTTCC  >  NC_002947/904789‑905077
                                                                                                                                                    |                                                                                                                                             
gcgTTCATCCACAACGACGACAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGcc                                                                                                                                              >  1:5362/1‑150 (MQ=255)
             acgacgACAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTg                                                                                                                                <  1:106980/150‑1 (MQ=255)
              cgacgaCAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTgg                                                                                                                               <  1:446418/150‑1 (MQ=255)
                  gacAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTgg                                                                                                                           >  1:269370/1‑150 (MQ=255)
                     aaCCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAg                                                                                                                        >  2:787481/1‑150 (MQ=255)
                        cTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGACTGGGAGTGc                                                                                                                     <  2:408362/150‑1 (MQ=255)
                              gCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGc                                                                                                               >  2:215094/1‑150 (MQ=255)
                                 agacTGCGGGCGCTGATTACGCGCTTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAg                                                                                                            >  2:683795/4‑150 (MQ=255)
                                  gCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGt                                                                                                           <  1:60852/150‑1 (MQ=255)
                                        gggCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGc                                                                                                     <  2:957186/150‑1 (MQ=255)
                                           cgcTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCgg                                                                                                                   >  2:817272/1‑133 (MQ=255)
                                           cgcTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCgg                                                                                                                   <  1:817272/133‑1 (MQ=255)
                                           cgcTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTg                                                                                                  <  2:411242/150‑1 (MQ=255)
                                            gcTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCGGCGCCGGCGAtttc                                                                                                                               >  2:403211/1‑118 (MQ=255)
                                            gcTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCGGCGCCGGCGAtttc                                                                                                                               <  1:403211/120‑3 (MQ=255)
                                                    cgcgCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGc                                                                                         <  1:180634/150‑1 (MQ=255)
                                                     gcgcCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGca                                                                                        >  1:1092782/1‑150 (MQ=255)
                                                         cTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATg                                                                                    >  1:838674/1‑150 (MQ=255)
                                                                     cTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGCGCCGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGt                                                                        >  1:906083/1‑150 (MQ=255)
                                                                         tCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGCGCGGCCGAGAGGGCCAGCGGTAGATCGTGCACATGCGGGTAGAGCGTATCg                                                                    >  2:1029675/1‑150 (MQ=255)
                                                                             ggTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCa                                                                <  2:5362/150‑1 (MQ=255)
                                                                                   aTTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCtt                                                          >  2:537592/1‑150 (MQ=255)
                                                                                         ggACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCaa                                                    <  1:1028979/150‑1 (MQ=255)
                                                                                               ggggCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAgcgc                                                            <  1:525839/136‑1 (MQ=255)
                                                                                               ggggCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAgcgc                                                            >  2:525839/1‑136 (MQ=255)
                                                                                               ggggCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGAttt                                              >  1:1053313/1‑150 (MQ=255)
                                                                                                  gCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGc                                           >  2:451832/1‑150 (MQ=255)
                                                                                                   cTGTTCCGCCGGAGTAGCCCCCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACAGGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGcc                                          <  1:1029675/150‑1 (MQ=255)
                                                                                                        cAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCcgacga                                     >  1:286920/1‑150 (MQ=255)
                                                                                                          gCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGc                                                                                         >  2:827836/1‑96 (MQ=255)
                                                                                                          gCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGc                                                                                         <  1:827836/96‑1 (MQ=255)
                                                                                                             ggAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCa                                <  2:535787/150‑1 (MQ=255)
                                                                                                              gAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGTGTGCGGCAGCGAGTGCCAGCGTGACAGCtt                                                                                           >  1:959462/1‑90 (MQ=255)
                                                                                                              gAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGTGTGCGGCAGCGAGTGCCAGCGTGACAGCtt                                                                                           <  2:959462/90‑1 (MQ=255)
                                                                                                                 tAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCt                            <  2:286920/150‑1 (MQ=255)
                                                                                                                   gCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGt                          >  2:837098/1‑150 (MQ=255)
                                                                                                                      aCCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCg                       >  2:109682/1‑150 (MQ=255)
                                                                                                                            cgcgcCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAgc                 <  1:215094/150‑1 (MQ=255)
                                                                                                                            cgcgcCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCATGATAATTGATCTGTCCGACGATGCCAGCCTGTTCGCCGAgc                 >  2:820509/1‑150 (MQ=255)
                                                                                                                               gcCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAgctgc              <  2:838674/150‑1 (MQ=255)
                                                                                                                               gcCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAgctgc              <  1:109682/150‑1 (MQ=255)
                                                                                                                                      gCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTcgacga       <  1:451832/150‑1 (MQ=255)
                                                                                                                                        tCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCGTGCTCAAGTAGGTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACt     >  1:249282/1‑150 (MQ=255)
                                                                                                                                         cAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACtt    >  1:616903/1‑150 (MQ=255)
                                                                                                                                         cAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACtt    <  2:616903/150‑1 (MQ=255)
                                                                                                                                          cagccaaTCGCGGCGCCGGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGAtttt                                             <  1:857108/103‑1 (MQ=255)
                                                                                                                                          cagccaaTCGCGGCGCCGGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGAtttt                                             >  2:857108/6‑108 (MQ=255)
                                                                                                                                          aaCGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGTGTTACAGCTTGCACATGCGGGGAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACTTc   >  1:577433/1‑150 (MQ=255)
                                                                                                                                          aaCGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTcggcga       >  1:200067/1‑146 (MQ=255)
                                                                                                                                          aaCGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTcgacga       <  2:200067/146‑1 (MQ=255)
                                                                                                                                           aCGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACTTcc  >  1:28231/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GCGTTCATCCACAACGACGACAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCCGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACTTCC  >  NC_002947/904789‑905077

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: