Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,269,946 G→T 8.0% L64M (CTG→ATG)  PP_2868 ← Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,269,9460GT 8.0% 153.6 / 3.5 50L64M (CTG→ATG) PP_2868Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (24/22);  new base T (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.06e-01

GCAGTTGATCGACGATGCCATCTTCGTAATGCACGCTCAGGCGGCTGTTGCGCCGCACCACATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCTCCGCCAAGGTCATGCCGCGGCGCTTGCGCAGGCCGCGGATGCGGGTGCCGAGGAACTGCGGCGCGGC  >  NC_002947/3269799‑3270093
                                                                                                                                                    |                                                                                                                                                   
gCAGTTGATCGACGATGCCATCTTCGTAATGCACGCTCAGGCGGCTGTTGCGCCGCACCACATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAgcg                                                                                                                                                   >  1:122592/1‑150 (MQ=255)
 cAGTTGATCGACGATGCCATCTTCGTAATGCACGCTCAGGCGGCTGTTGCGCCGCACCACATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAgcgc                                                                                                                                                  >  1:430197/1‑150 (MQ=255)
    ttGATCGACGATGCCATCTTCGTAATGCACGCTCAGGCGGCTGTTGCGCCGCACCACATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCgg                                                                                                                                               >  2:1025290/1‑150 (MQ=255)
       aTCGACGATGCCATCTTCGTAATGCACGCTCAGGCGGCTGTTGCGCCGCACCACATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGAt                                                                                                                                            >  1:1060835/1‑150 (MQ=255)
            cgaTGCCATCTTCGTAATGCACGCTCAGGCGGCTGTTGCGCCGCACCACATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACg                                                                                                                                       >  2:886194/1‑150 (MQ=255)
                  cATCTTCGTAATGCACGCTCAGGCGGCTGTTGCGCCGCACCACATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAgg                                                                                                                                 <  1:169474/150‑1 (MQ=255)
                  cATCTTCGTAATGCACGCTCAGGCGGCTGTTGCGCCGCACCACATAGCCTAGGTCTTCGGGCGAGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAgg                                                                                                                                 >  2:618200/1‑150 (MQ=255)
                     cTTCGTAATGCACGCTCAGGCGGCTGTTGCGCCGCACCACATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGgaga                                                                                                                              <  2:41453/150‑1 (MQ=255)
                               tgggtgcAGGCGGCTGTTGCGCCGCACCACAGGGCCTGGGTCTTCGGGCGCGGTGGTGGGCTCTCTGGCGAAGAACCACTGGTTGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTc                                                                                                                    <  2:215156/144‑1 (MQ=255)
                               cACGCTCAGGCGGCTGTTGCGCCGCACCACATAGCCTTGGTCTTCGGGCGCGGCGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTc                                                                                                                    <  2:707567/150‑1 (MQ=255)
                                aCGCTCAGGCGGCTGTTGCGCCGCACCACATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTcc                                                                                                                   >  1:358499/1‑150 (MQ=255)
                                 cGCTCAGGCGGCTGTTGCGCCGCACCACATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCa                                                                                                                  <  2:122592/150‑1 (MQ=255)
                                 cGCTCAGGCGGCTGTTGCGCCGCACCACATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCa                                                                                                                  <  1:529210/150‑1 (MQ=255)
                                    tCAGGCGGCTGTTGCGCCGCACCACATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAgct                                                                                                               >  1:663371/1‑150 (MQ=255)
                                    tCAGGCGGCTGTTGCGCCGCACCACATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAgct                                                                                                               >  1:663904/1‑150 (MQ=255)
                                        gcggcTGTTGCGCCGCACCACATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGATGCGGGCGAGTGGTAACAGCGCCGGGATCCACGGGTAGGAGAGGTTGCGCTCCAgctggct                                                                                                           >  2:578279/1‑150 (MQ=255)
                                              gTTGCGCCGCACCACATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTa                                                                                                     <  2:663371/150‑1 (MQ=255)
                                              gTTGCGCCGCACCACATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTa                                                                                                     <  2:663904/150‑1 (MQ=255)
                                               ttGCGCCGCACCACATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAg                                                                                                    >  2:660090/1‑150 (MQ=255)
                                                   gccgcACCACATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATTTAGCCTg                                                                                                >  2:240232/1‑150 (MQ=255)
                                                        accacATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCgg                                                                                                                                               >  1:847042/1‑98 (MQ=255)
                                                        accacATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCgg                                                                                                                                               <  2:847042/98‑1 (MQ=255)
                                                           acaTAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAg                                                                                         >  1:372417/1‑149 (MQ=255)
                                                           acaTAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAg                                                                                         <  2:372417/149‑1 (MQ=255)
                                                             aTAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTc                                                                                      >  2:442547/1‑150 (MQ=255)
                                                               aGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGCTGTTCAACATCGcc                                                                                                                                                 <  2:878798/89‑1 (MQ=255)
                                                               aGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGCTGTTCAACATCGcc                                                                                                                                                 >  1:878798/1‑89 (MQ=255)
                                                                         ttCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGAt                                                                                                        <  1:28050/120‑1 (MQ=255)
                                                                         ttCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGAt                                                                                                        >  2:28050/1‑120 (MQ=255)
                                                                         ttCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGAGGGTCCCGTCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTCACTGTGTGCGGc                                                                          <  2:1051866/150‑1 (MQ=255)
                                                                              gcgcGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCt                                                                     >  1:729867/1‑150 (MQ=255)
                                                                                 cgGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGATGCGAGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCTCCg                                                                  >  2:271095/1‑150 (MQ=255)
                                                                                         gCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCTCCGCCAAGGTc                                                          <  2:610902/150‑1 (MQ=255)
                                                                                           cTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATTTAGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCTCCGCCAAGGTCAt                                                        <  1:240232/150‑1 (MQ=255)
                                                                                           cTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCTCCGCCAAGGTCAt                                                        >  2:37321/1‑150 (MQ=255)
                                                                                                      agaaCCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCTCCGCCAAGGTCATGCCGCGGCGCt                                             <  2:729867/150‑1 (MQ=255)
                                                                                                          ccACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCTCCGCCAAGGTCATGCCGCGGCGCTTgcg                                         >  2:437328/1‑150 (MQ=255)
                                                                                                           cACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCTCCGCCAAGGTCATGCCGCGGCGCTTgcgc                                        >  1:220395/1‑150 (MQ=255)
                                                                                                              tGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTGGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCTCCGCCAAGGTCATGCCGCGGCGCTTGCGCAgg                                     >  2:590800/1‑150 (MQ=255)
                                                                                                                      cACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCTCCGCCAAGGTCATGCCGCGGCGCTTGCGCAGGCCGCGGAt                             <  1:437328/150‑1 (MQ=255)
                                                                                                                        cGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCTCCGCCAAGGTCATGCCGCGGCGCTTGCGCAGGCCGCGGATGc                           <  2:1060835/150‑1 (MQ=255)
                                                                                                                        cGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCTCCGCCAAGGTCATGCCGCGGCGCTTGCGCAGGCCGCGGATGc                           <  2:430197/150‑1 (MQ=255)
                                                                                                                               gCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCTCCGCCAAGGTCATGCCGCGGCGCTTGCGCAGGCCGCGGATGCGGGTGcc                    >  2:43898/1‑150 (MQ=255)
                                                                                                                                     ggCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCTCCAGCAGGGTCATGCCGCGGCGCTTGCGCAGGCCGCGGATGCGGGTGCCGAGGaa              >  1:1092612/1‑150 (MQ=255)
                                                                                                                                          tGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCTCCGCCAAGGTCATGCCGCGGCGCTTGCGCAGGCCGCGGATGCGGGTGCCGAGGAACTgcg         >  2:958782/1‑150 (MQ=255)
                                                                                                                                            ttGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCTCCGCCAAGGTCATGCCGCGGCGCTTGCGCAGGCCGCGGATGCGGGTGCCGAGGAACTgcggc       <  2:152586/150‑1 (MQ=255)
                                                                                                                                              tcaaCATCGCCCGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCa                                                                                          <  2:159114/64‑1 (MQ=255)
                                                                                                                                              tcaaCATCGCCCGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCa                                                                                          >  1:159114/3‑66 (MQ=255)
                                                                                                                                                  cAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCa                                                                                          <  1:977213/62‑1 (MQ=255)
                                                                                                                                                  cAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCa                                                                                          >  2:977213/1‑62 (MQ=255)
                                                                                                                                                  cAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCTCCGCCAAGGTCATGCCGCGGCGCTTGCGCAGGCCGCGGATGCGGGTGCCGAGGAACTgcggcgcggc  <  1:886194/150‑1 (MQ=255)
                                                                                                                                                  cAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCTCCGCCAAGGTCATGCCGCGGCGCTTGCGCAGGCCGCGGATGCGGGTGCCGAGGAACTgcggcgcggc  <  1:43898/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCAGTTGATCGACGATGCCATCTTCGTAATGCACGCTCAGGCGGCTGTTGCGCCGCACCACATAGCCTTGGTCTTCGGGCGCGGTGGTGGCCTCGCTGGCGAAGAACCACTGGATGGTCACGCCCAGGCTGCGGGCGATGTTGAACAGCGCCGGGATCGACGGGTAGGAGAGGTTGCGCTCCAGCTGGCTGATGTAGCCTGCGGTCAGTTCACTGTGTGCGGCCAGCTCCGCCAAGGTCATGCCGCGGCGCTTGCGCAGGCCGCGGATGCGGGTGCCGAGGAACTGCGGCGCGGC  >  NC_002947/3269799‑3270093

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: