Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,047,960 C→G 13.2% P224A (CCC→GCC) ‡ PP_2660 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,047,9600CG13.2% 146.6 / 14.1 53P224A (CCC→GCC) ‡PP_2660hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (19/27);  new base G (3/4);  total (22/31)
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.50e-01

AGGGTTACTTCAACCCGTCCTCGGTCATACGCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGCCAGGTACGCATTTCGTTCGATATCGACGACACGCTGGCCTGCCTGCCCGAA  >  NC_002947/3047811‑3048102
                                                                                                                                                     |                                                                                                                                              
aGGGTTACTTCAACCCGTCCTCGGTCATACGCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGc                                                                                                                                                <  1:1195458/150‑1 (MQ=255)
  ggTTACTTCAACCCGTCCTCGGTCATACGCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGccc                                                                                                                                              >  1:431645/1‑150 (MQ=255)
   gTTACTTCAACCCGTCCTCGGTCATACGCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCg                                                                                                                                             >  2:126845/1‑150 (MQ=255)
      aCTTCAACCCGTCCTCGGTCATACGCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGtgc                                                                                                                                          <  2:511244/150‑1 (MQ=255)
        ttCAACCCGTCCTCGGTCATACGCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGGGCGtgccg                                                                                                                                        >  1:371668/1‑149 (MQ=255)
        ttCAACCCGTCCTCGGTCATACGCAACTCCGCAAAACACGGTGTCAAGGCTTTCTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGACAGCAGCTGAGCGCGGCAAGTTCTTCCTGTACGGCACCCCCGtccct                                                                                                                                        >  2:619239/1‑150 (MQ=255)
          cAACCCGTCCTCGGTCATACGCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGtgcctgc                                                                                                                                      <  2:1071244/150‑1 (MQ=255)
            aCCCGTCCTCGGTCATACGCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGGGCGtgccgttca                                                                                                                                    <  2:371668/150‑6 (MQ=255)
             cccGTCCTCGGTCATACGCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGGGCGtgccgttcag                                                                                                                                   <  2:912829/150‑7 (MQ=255)
                gTCCTCGGTCATACGCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGa                                                                                                                                <  2:399516/150‑1 (MQ=255)
                        tCATACGCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGcc                                                                                                                        <  1:126845/150‑1 (MQ=255)
                        tCATACGCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGcc                                                                                                                        <  2:1137025/150‑1 (MQ=255)
                          aTACGCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCTTGCCTGCGAGGGGAGCCGc                                                                                                                      <  1:203544/150‑1 (MQ=255)
                          aTACGCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGc                                                                                                                      <  2:1389179/150‑1 (MQ=255)
                           tACGCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCt                                                                                                                     >  2:1229615/1‑150 (MQ=255)
                              gCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGAt                                                                                                                  >  1:155779/1‑150 (MQ=255)
                              gCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGAt                                                                                                                  <  2:502045/150‑1 (MQ=255)
                              gCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGAt                                                                                                                  <  1:1111207/150‑1 (MQ=255)
                              gCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTGTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGAt                                                                                                                  <  1:526259/150‑1 (MQ=255)
                                       ccAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTc                                                                                                         >  2:47227/1‑150 (MQ=255)
                                           acacGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCACGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATcc                                                                                                     <  2:1087625/150‑1 (MQ=255)
                                             acGATGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCt                                                                                                   >  2:1279238/1‑150 (MQ=255)
                                               ggTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGt                                                                                                 <  2:431645/150‑1 (MQ=255)
                                                     aaGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGGGCGtgccgttca                                                                                                                                    >  2:1294464/1‑104 (MQ=255)
                                                     aaGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGGGCGtgccgttca                                                                                                                                    <  1:1294464/109‑6 (MQ=255)
                                                      aGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGc                                                                                                                      >  2:425295/1‑122 (MQ=255)
                                                      aGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGc                                                                                                                      <  1:425295/122‑1 (MQ=255)
                                                       ggCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGa                                                                                         <  1:47227/150‑1 (MQ=255)
                                                          tttGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGa                                                                                                                                >  2:703913/1‑108 (MQ=255)
                                                          tttGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGa                                                                                                                                <  1:703913/108‑1 (MQ=255)
                                                                cGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCaa                                                                                <  1:868698/150‑1 (MQ=255)
                                                                         cAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGaa                                                                       >  2:139301/1‑150 (MQ=255)
                                                                             aTCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCa                                                                   >  2:435082/1‑150 (MQ=255)
                                                                                        gTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGt                                                        <  1:619239/150‑1 (MQ=255)
                                                                                            gccgcTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGcc                                                    >  1:504745/1‑150 (MQ=255)
                                                                                                cTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGCCAGGt                                                <  1:1229615/150‑1 (MQ=255)
                                                                                                       ggAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGCCAGGTACGCAtt                                         >  1:539503/1‑150 (MQ=255)
                                                                                                       ggAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGCCAGGTACGCAtt                                         <  2:539503/150‑1 (MQ=255)
                                                                                                        gAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGCCAGGTACGCAttt                                        >  1:470476/1‑150 (MQ=255)
                                                                                                         aaGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGCCAGGTACGCATttc                                       <  1:139301/150‑1 (MQ=255)
                                                                                                           gCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGTCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGCCAGGTACGCATttcgt                                     >  1:392380/1‑150 (MQ=255)
                                                                                                             gAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGCCAGGTACGCATttcgttc                                   <  1:948933/150‑1 (MQ=255)
                                                                                                                     cTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATAGAACACAAGCTCAAGCGTGGCCAGGTACGCATTTCGTTAGATATcga                           >  2:610888/1‑150 (MQ=255)
                                                                                                                           ccGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGCCAGGTACGCATTTCGTTCGATATCGACGacac                     >  2:237570/1‑150 (MQ=255)
                                                                                                                            cGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGCCAGGTACGCATTTCGTTCGATATCGACGACACg                    >  2:538771/1‑150 (MQ=255)
                                                                                                                                 aGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGCCAGGTACGCATTTCGTTCGATATCGACGACACGCTGgc               <  1:1279238/150‑1 (MQ=255)
                                                                                                                                      ttcCTGAATGGCACGCACGTACATGCCTGCCAGGTAAGCCGGTGAGGGAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGCCAGGTACGCATTTCGTTCGATATCGACGACACGCTGgcttgcc          <  2:155779/150‑1 (MQ=255)
                                                                                                                                      ttcCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGCCAGGTACGCATTTCGTTCGATATCGACGACACGCTGgcctgcc          <  1:435082/150‑1 (MQ=255)
                                                                                                                                          tGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGCCAGGTACGCATTTCGTTCGATATCGACGACACGCTGgc               <  2:1126255/141‑1 (MQ=255)
                                                                                                                                          tGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGCCAGGTACGCATTTCGTTCGATATCGACGACACGCTGgc               >  1:1126255/1‑141 (MQ=255)
                                                                                                                                              cGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGCCAGGTACGCATTTCGTTCGATATCGACGACACGCTGGCCTGCCTGCCCGaa  <  2:175534/150‑1 (MQ=255)
                                                                                                                                              aggCACGGGCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGCCAGGTACGCATttcgttc                                   >  1:1141464/2‑117 (MQ=255)
                                                                                                                                              aggCACGGGCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGCCAGGTACGCATttcgttc                                   <  2:1141464/116‑1 (MQ=255)
                                                                                                                                                    gCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATc                                                                          <  2:668602/72‑1 (MQ=255)
                                                                                                                                                    gCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATc                                                                          >  1:668602/1‑72 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
AGGGTTACTTCAACCCGTCCTCGGTCATACGCAACTCCGCCAAACACGGTGTCAAGGCTTTGTACGAACTCCCCAAGATCAAGAAAAAGTACGCCGCTTTCCTGGAAGCGAAAGCACCTGAGCCCGGCAAGTTCTTCCTGAACGGCACGCCCGTGCCTGCCTGCGAGGGGAGCCGCTGATGCAAGGTTCATCCCTGTTCGCACGAGGGCGGCAAGCCATCGAACGCAAGCTCAAGCGTGGCCAGGTACGCATTTCGTTCGATATCGACGACACGCTGGCCTGCCTGCCCGAA  >  NC_002947/3047811‑3048102

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