Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,041,186 G→A 7.2% A269V (GCA→GTA)  PP_2653 ← Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,041,1860GA 7.2% 157.6 / 2.4 56A269V (GCA→GTA) PP_2653Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (22/29);  new base A (2/2);  total (25/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 = 7.94e-01

GCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAGGTGGTGAACGTGGACAGGAACCGCAGCTCACCATACG  >  NC_002947/3041040‑3041333
                                                                                                                                                  |                                                                                                                                                    
gCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGcc                                                                                                                                                    <  3:205488/149‑1 (MQ=255)
           cGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTCGATTTGCATGAAATAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTc                                                                                                                                         >  8:307304/1‑149 (MQ=255)
                gTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACg                                                                                                                                    >  7:4216/1‑149 (MQ=255)
                     cAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGctgcc                                                                                                                               <  4:304294/149‑1 (MQ=255)
                      aGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGctgcct                                                                                                                              <  5:45625/149‑1 (MQ=255)
                              aaGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCtt                                                                                                                      <  1:321106/149‑1 (MQ=255)
                              aaGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCtt                                                                                                                      <  2:260812/149‑1 (MQ=255)
                                    ggggATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACGCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACTCTGCCTGCATTCTTTGCCgg                                                                                                                <  2:334011/149‑1 (MQ=255)
                                    ggggATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCgg                                                                                                                >  1:104771/1‑149 (MQ=255)
                                           tCGTGACCACCGGGGCCATGTGCCAGCCGTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAg                                                                                                         >  7:18900/1‑149 (MQ=255)
                                            cGTGACCACCGGGGCCATGTGCCAGCAGTTCGATATGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCATGGCCAGACCCTGACGTTTACGCTTTAACTCCCGCCCAATGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAgc                                                                                                        >  3:37423/1‑148 (MQ=255)
                                             gTGACCACCGGGGCCATGTGCCAGCCGTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAggg                                                                                                       >  5:364940/1‑149 (MQ=255)
                                              tGACCACCGGGGCCATGTGCCAGCCGTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGct                                                                                                                                  <  5:9545/121‑1 (MQ=255)
                                              tGACCACCGGGGCCATGTGCCAGCCGTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGct                                                                                                                                  >  6:9545/1‑121 (MQ=255)
                                               gACCACCGGGGCCATGTGCCAGCCGTTCGATTTGCATCAAACAGGCAACCAGGTCAGCGCCAGGATCCTAGTCCAGACCCTGACGTTTACGCTGTAACTGCCGAGCATTGCTCGAACGCTGCATGCATTCTTAGCCGGGTCTCAGGGTc                                                                                                    >  3:213427/1‑149 (MQ=255)
                                                accaccGGGGCCATGTGCCAGCCGTGCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTc                                                                                                    <  7:298835/148‑1 (MQ=255)
                                                     cGGGGCCATGTGCCAGCCGTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCggg                                                                                              <  2:14038/149‑1 (MQ=255)
                                                      ggggCCATGTGCCAGCCGTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCgggg                                                                                             <  3:24633/149‑1 (MQ=255)
                                                           cATGTGCCAGCCGTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAg                                                                                        >  5:49820/1‑149 (MQ=255)
                                                           cATGTGCCAGCCGTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAg                                                                                        <  2:104771/149‑1 (MQ=255)
                                                                    gCCGTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGAt                                                                               <  8:18900/149‑1 (MQ=255)
                                                                       gTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATgcg                                                                            >  8:11136/1‑149 (MQ=255)
                                                                        ttCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATgcgc                                                                           <  6:364940/149‑1 (MQ=255)
                                                                           gATTTGCATGAAATAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAg                                                                         <  7:307304/148‑1 (MQ=255)
                                                                             tttGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTcggcgg                                                                                               <  3:30753/124‑1 (MQ=255)
                                                                             tttGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTcggcgg                                                                                               >  4:30753/1‑124 (MQ=255)
                                                                             tttGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGGGCAGTGa                                                                      <  6:173131/149‑1 (MQ=255)
                                                                                 cATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTg                                                                   >  2:111387/1‑148 (MQ=255)
                                                                                  aTGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTACCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCa                                                                 >  4:27814/1‑149 (MQ=255)
                                                                                  aTGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTACCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCa                                                                 >  8:27264/1‑149 (MQ=255)
                                                                                        cAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGa                                                           <  1:176472/149‑1 (MQ=255)
                                                                                         aGGCCACCAGGTAAGCGCCAGGTTCCTTGGTCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGAc                                                          <  7:128368/149‑1 (MQ=255)
                                                                                             caccaGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATcc                                                      >  1:268045/1‑149 (MQ=255)
                                                                                              accaGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCt                                                     >  2:52060/1‑149 (MQ=255)
                                                                                                     cAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCAt                                              >  5:30955/1‑149 (MQ=255)
                                                                                                     cAGCGCCAGGTTCCTTGGCCAGACCCTGACGGTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGGCTGAGGGTCGGCGGGGATCAGGGGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCAt                                              <  4:37423/149‑1 (MQ=255)
                                                                                                          ccAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCaa                                         <  7:23935/149‑1 (MQ=255)
                                                                                                          ccAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCa                                          >  1:751/1‑148 (MQ=255)
                                                                                                            aGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAg                                       <  1:52060/149‑1 (MQ=255)
                                                                                                               ttCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAggtg                                    <  6:30955/149‑1 (MQ=255)
                                                                                                                   ttGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAGgtggtg                                 >  5:239746/1‑148 (MQ=255)
                                                                                                                     ggCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGCTAGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAGgtggtg                                 <  1:111387/147‑1 (MQ=255)
                                                                                                                           aCCCTGACGTTTACGCTGTAACTGCCGCCCATTGCTCGAACGCTGCCTGGATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAGGTGGTGAACGTGGAc                        >  2:344152/1‑149 (MQ=255)
                                                                                                                           aCCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGGGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAGGTGGTGAACGTGGAc                        <  8:4216/149‑1 (MQ=255)
                                                                                                                             ccTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAGGTGGTGAACGTGGACAg                      <  6:49820/149‑1 (MQ=255)
                                                                                                                                 aCGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAGGTGGTGAACGTGGACAGGAAc                  <  2:268045/149‑1 (MQ=255)
                                                                                                                                   gTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAGGTGGTGAACGTGGACAGGAACCg                >  5:347582/1‑149 (MQ=255)
                                                                                                                                          cTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAGGTGGTGAACGTGGACAGGAACCGCAGCTCa         <  2:751/149‑1 (MQ=255)
                                                                                                                                           tGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAggtgg                                   >  7:152713/1‑122 (MQ=255)
                                                                                                                                           tGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAggtgg                                   <  8:152713/122‑1 (MQ=255)
                                                                                                                                             tAACTACCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAGGTGGTGAACGTGGACAGGAACCGCAGCTCAc        <  7:27264/147‑1 (MQ=255)
                                                                                                                                             tAACTACCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAGGTGGTGAACGTGGACAGGAACCGCAGCTCAc        <  3:27814/147‑1 (MQ=255)
                                                                                                                                               aCTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAGGTGGTGAACGTGGAc                        <  6:209965/129‑1 (MQ=255)
                                                                                                                                               aCTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAGGTGGTGAACGTGGAc                        >  5:209965/1‑129 (MQ=255)
                                                                                                                                               aCTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAGGTGGTGAACGTGGACAGGAACCGCAGCTCACCATa    >  8:277079/1‑149 (MQ=255)
                                                                                                                                                 tGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAGGTGGTGAACGTGGACAGGAACCGCAGCTCACCATACg  >  7:3411/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCACACGCTGCCGAGCGTGGCCAGACCTTGAAGGCCGGGGATGTCGTGACCACCGGGGCCATGTGCCAGCCGTTCGATTTGCATGAAACAGGCCACCAGGTCAGCGCCAGGTTCCTTGGCCAGACCCTGACGTTTACGCTGTAACTGCCGACCATTGCTCGAACGCTGCCTGCATTCTTTGCCGGGTCTGAGGGTCGGCGGGGATCAGGTGTTCGATGCGCAGTGACTGCACGGTGACATCCTGCGGCATGCCAAAGGTGGTGAACGTGGACAGGAACCGCAGCTCACCATACG  >  NC_002947/3041040‑3041333

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: