Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,726,657 G→A 7.3% P122L (CCG→CTG)  opgH ← glucans biosynthesis glucosyltransferase H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,726,6570GA 7.3% 154.2 / 2.9 55P122L (CCG→CTG) opgHglucans biosynthesis glucosyltransferase H
Reads supporting (aligned to +/- strand):  ref base G (26/25);  new base A (2/2);  total (28/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 = 5.72e-01

ACCAGCCCTGGTACGGCAGGATGCCTTTCATGTACCAGCCAGCGACGATGGTCTGGCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGA  >  NC_002947/5726520‑5726798
                                                                                                                                          |                                                                                                                                             
aCCAGCCCTGGTACGGCAGGATGCCTTTCATGTACCAGCCAGCGACGATGGTCTGGCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCAt                                                                                                                                   >  1:401876/1‑150 (MQ=255)
        tGGTACGGCAGGATGCCTTTCATGTACCAGCCAGCGACGATGGTCTGGCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCgggg                                                                                                                                      >  2:858399/1‑139 (MQ=255)
        tGGTACGGCAGGATGCCTTTCATGTACCAGCCAGCGACGATGGTCTGGCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCgggg                                                                                                                                      <  1:858399/139‑1 (MQ=255)
               gCAGGATGCCTTTCATGTACCAGCCAGCGACGATGGTCTGGCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGggcg                                                                                                                                               >  2:916818/1‑123 (MQ=255)
               gCAGGATGCCTTTCATGTACCAGCCAGCGACGATGGTCTGGCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGggcg                                                                                                                                               <  1:916818/123‑1 (MQ=255)
               gCAGGATGCCTTTCATGTACCAGCCAGCGACGATGGTCTGGCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCa                                                                                                                    >  1:671306/1‑150 (MQ=255)
                         tttCATGTACCGGCCAGCGACGATGGTCTGGCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCca                                                                                                          <  2:56954/150‑1 (MQ=255)
                         tttCATGTACCAGCCAGCGACGATGGTCTGGCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCca                                                                                                          >  1:149311/1‑150 (MQ=255)
                         tttCATGTACCAGCCAGCGACGATGGTCTGGCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTCTGCTAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCca                                                                                                          >  1:912358/1‑150 (MQ=255)
                              tGTACCAGCCAGCGACGATGGTCTGGCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGc                                                                                                     <  1:590054/150‑1 (MQ=255)
                                tACCAGCCAGCGACGATGGTCTGGCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCgg                                                                                                   >  2:822597/1‑150 (MQ=255)
                                               aTGGTCTGGCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTg                                                                                    >  1:210402/1‑150 (MQ=255)
                                               aTGGTCTGGCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTg                                                                                    <  2:37666/150‑1 (MQ=255)
                                                     tGGCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGt                                                                                                               <  1:44303/117‑1 (MQ=255)
                                                     tGGCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGt                                                                                                               >  2:44303/1‑117 (MQ=255)
                                                     tGGCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTg                                                                                       <  2:312876/141‑1 (MQ=255)
                                                     tGGCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTg                                                                                       >  1:312876/1‑141 (MQ=255)
                                                     gggCCGATCATCAGCGACAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCAt                                                                              <  1:109469/149‑1 (MQ=255)
                                                       gCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATATAAAGTACTGGGTTCTATCGGGCCGGCGGCAGCCCCCGGTTGGGGCGGATGGTTTCAGTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGAACATTGTTTCGCCAttg                                                                            >  1:523907/1‑148 (MQ=255)
                                                          gATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTc                                                                         >  2:225721/1‑150 (MQ=255)
                                                            tcatcaGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCgg                                                                       <  2:452773/150‑1 (MQ=255)
                                                            tcatcaGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCgg                                                                       <  1:813774/150‑1 (MQ=255)
                                                            tcatcaGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCgg                                                                       <  2:557900/150‑1 (MQ=255)
                                                              atcaGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGc                                                                     >  1:846528/1‑150 (MQ=255)
                                                                caGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCAc                                                                   >  1:871542/1‑150 (MQ=255)
                                                                      cagcagAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCtt                                                             >  1:238714/1‑150 (MQ=255)
                                                                      cagcagAATGTAACGGCGGATAGAGCCGAAAGTGCGCTTGTGTTTCGTTTAATTGTCTTGCTTGGGCGGCTCTAGTGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCtt                                                             >  1:193278/1‑150 (MQ=255)
                                                                           gAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGt                                                          >  1:498576/1‑148 (MQ=255)
                                                                           gAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGt                                                          <  2:498576/148‑1 (MQ=255)
                                                                           gAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGc                                                        >  2:963974/1‑150 (MQ=255)
                                                                           gAATGTAACGGCGGATAGAGCCGACCGTGCGCAAGAATGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGc                                                        >  1:366647/1‑150 (MQ=255)
                                                                            aaTGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCg                                                       >  2:523367/1‑150 (MQ=255)
                                                                             aTGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCGt                                                      <  1:404820/150‑1 (MQ=255)
                                                                               gTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCGttt                                                    <  1:822597/150‑1 (MQ=255)
                                                                                  aCGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCAGCTGCgg                                                                                                                                        >  1:1120466/1‑63 (MQ=255)
                                                                                  aCGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCAGCTGCgg                                                                                                                                        <  2:538492/63‑1 (MQ=255)
                                                                                  aCGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCAGCTGCgg                                                                                                                                        <  2:1120466/63‑1 (MQ=255)
                                                                                  aCGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCAGCTGCgg                                                                                                                                        >  1:538492/1‑63 (MQ=255)
                                                                                       ggATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATcggcg                                            <  1:209873/150‑1 (MQ=255)
                                                                                        gATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATcggcgg                                           >  2:1006058/1‑150 (MQ=255)
                                                                                         aTAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCggg                                          >  1:902328/1‑150 (MQ=255)
                                                                                          tAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCgggg                                         <  1:963974/150‑1 (MQ=255)
                                                                                           agagCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGc                                        <  2:322196/150‑1 (MQ=255)
                                                                                               ccGNCCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTTCGCCATGGCTCGGGCACGACCTTAGTGCGTTTGATCGGCGGGGCGATc                                    >  2:694865/1‑150 (MQ=255)
                                                                                               ccGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATc                                    <  2:149311/150‑1 (MQ=255)
                                                                                               ccGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATc                                    >  1:1071862/1‑150 (MQ=255)
                                                                                                    cGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAg                               <  2:1071862/150‑1 (MQ=255)
                                                                                                          ccAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAg                         <  1:1006058/150‑1 (MQ=255)
                                                                                                           cAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAgg                        <  1:923370/150‑1 (MQ=255)
                                                                                                               ccgctggcggcAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAggcggc                    >  1:854971/6‑150 (MQ=255)
                                                                                                                  ccGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGc                 >  2:165145/1‑150 (MQ=255)
                                                                                                                         cAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGCGCGGCCTGTCAGGCGTCGCTACATGCGGATCAGCACATTGGTGCTCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGCCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTc          <  1:144805/150‑1 (MQ=255)
                                                                                                                           gCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGa        <  2:172396/150‑1 (MQ=255)
                                                                                                                               gcccgtcgggCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGc     <  2:907180/143‑1 (MQ=255)
                                                                                                                                cgcTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGAt                                                                                                 <  2:1064645/56‑1 (MQ=255)
                                                                                                                                cgcTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGAt                                                                                                 >  1:1064645/1‑56 (MQ=255)
                                                                                                                                 gcTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGa  <  2:217133/150‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                             
ACCAGCCCTGGTACGGCAGGATGCCTTTCATGTACCAGCCAGCGACGATGGTCTGGCCGATCATCAGCGCCAGCAGAATGTAACGGCGGATAGAGCCGACCGTGCGCCAGCGGGCCGGCGGCAGCTCGCGCTTGGGCGGCTGCGGGGCATTGGGGCGGCCGGTCATGCGTCGCCACATGCGGATCAGCACATTGGTGCGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGA  >  NC_002947/5726520‑5726798

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: