Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 358,583 C→A 7.5% P16H (CCC→CAC)  gbdR → AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947358,5830CA 7.5% 170.1 / 2.4 53P16H (CCC→CAC) gbdRAraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (27/22);  new base A (2/2);  total (29/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.95e-01

GTTGAGATTGGCAAGCCTACCCCCGGAAAATGACCGAGGGGTCATAACAAGAAGTAAGTGCCCGCCTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGTGTGGGCCA  >  NC_002947/358434‑358732
                                                                                                                                                     |                                                                                                                                                     
gTTGAGATTGGCAAGCCTACCCCCGGAAAATGACCGAGGGGTCATAACAAGAAGTAAGTGCCCGCCTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACAcc                                                                                                                                                       >  1:939061/1‑150 (MQ=255)
 ttGAGATTGGCAAGCCTACCCCCGGAAAATGACCGAGGGGTCATAACAAGAAGTAAGTGCCCGCCTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACAccc                                                                                                                                                      <  2:445649/150‑1 (MQ=255)
          gCAAGCCTACCCCCGGAAAATGACCGAGGGGTCATAACAAGAAGTAAGTGCCCGCCTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATc                                                                                                                                             >  1:646125/1‑150 (MQ=255)
                 tACCCCCGGAAAATGACCGAGGGGTCATAACAAGAAGTAAGTGCCCGCCTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGACCCCCACACCACCCCAGTCCATctgtgttc                                                                                                                                      <  2:565693/150‑8 (MQ=255)
                              tGACCGAGGGGTCATAACAAGAAGTAAGTGCCCGCCTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACaa                                                                                                                         <  2:1149268/150‑1 (MQ=255)
                                   gAGGGGTCATAACAAGAAGTAAGTGCCCGCCTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCa                                                                                                                    <  2:643341/150‑1 (MQ=255)
                                         tCATAACAAGAAGTAAGTGCCCGCCTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCActttga                                                                                                              >  2:1035935/1‑145 (MQ=255)
                                            tAACAAGAAGTAAGTGCCCGCCTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGa                                                                                                                            >  2:721269/1‑133 (MQ=255)
                                            tAACAAGAAGTAAGTGCCCGCCTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGa                                                                                                                            <  1:721269/133‑1 (MQ=255)
                                              aCAAGAAGTAAGTGCCCGCCTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTcc                                                                                                         >  2:1147943/1‑150 (MQ=255)
                                                aagaagTAAGTGCCCGCCTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCt                                                                                                       >  1:714095/1‑150 (MQ=255)
                                                 agaagTAAGTGCCCGCCTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTg                                                                                                      >  1:722716/1‑150 (MQ=255)
                                                 agaagGAAGTGCCCGCCTGCGGCCCAGCCTGAATTTGTGCGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTg                                                                                                      <  2:1103131/150‑1 (MQ=255)
                                                            cccGCCTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGt                                                                                           <  2:646125/150‑1 (MQ=255)
                                                             ccGCCTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCgg                                                                                                                                           <  2:623751/101‑1 (MQ=255)
                                                             ccGCCTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCgg                                                                                                                                           >  1:623751/1‑101 (MQ=255)
                                                             ccGCCTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTc                                                                                          >  2:796892/1‑150 (MQ=255)
                                                             ccGCCTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTc                                                                                          <  2:939061/150‑1 (MQ=255)
                                                                ccTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCAcc                                                                                                                  >  2:1070698/1‑123 (MQ=255)
                                                                ccTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCAcc                                                                                                                  <  1:1070698/123‑1 (MQ=255)
                                                                 cTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTg                                                                                                      <  2:170930/134‑1 (MQ=255)
                                                                 cTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTg                                                                                                      >  1:170930/1‑134 (MQ=255)
                                                                        acacCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTgcg                                                                               >  2:1171741/1‑150 (MQ=255)
                                                                          acCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCgg                                                                                                                                           <  2:537254/88‑1 (MQ=255)
                                                                          acCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCgg                                                                                                                                           >  1:537254/1‑88 (MQ=255)
                                                                          acCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGa                                                                                                                            >  2:1084782/1‑103 (MQ=255)
                                                                          acCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGa                                                                                                                            <  1:1084782/103‑1 (MQ=255)
                                                                            ccTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATg                                                                           <  1:418669/150‑1 (MQ=255)
                                                                            ccTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTACTGGACAATTTCCCCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATg                                                                           <  2:266329/150‑1 (MQ=255)
                                                                            ccTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACACCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATg                                                                           >  2:785373/1‑150 (MQ=255)
                                                                                      tgtgAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGc                                                                 >  1:720996/1‑150 (MQ=255)
                                                                                            gagaAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTg                                                           >  2:461706/1‑150 (MQ=255)
                                                                                            aagaAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACACCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTg                                                           >  2:828179/2‑150 (MQ=255)
                                                                                                     cgATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGc                                                  >  1:664111/1‑150 (MQ=255)
                                                                                                            gtcgtACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCg                                           >  1:417588/1‑150 (MQ=255)
                                                                                                                tACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTgg                                       <  1:1171741/150‑1 (MQ=255)
                                                                                                                tACACCTCCGGGAATCCAACCCAGAACCGCACAACACACCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTgg                                       <  1:785373/150‑1 (MQ=255)
                                                                                                                     cTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCacac                                  >  2:1207490/1‑150 (MQ=255)
                                                                                                                      tCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACg                                 >  2:652490/1‑150 (MQ=255)
                                                                                                                           gAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGAc                            >  2:974572/1‑150 (MQ=255)
                                                                                                                             aTCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGAcga                          <  1:461706/150‑1 (MQ=255)
                                                                                                                              tCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGAcgat                         <  2:714095/150‑1 (MQ=255)
                                                                                                                                cAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGAcgatcg                       >  1:330695/1‑150 (MQ=255)
                                                                                                                                cAACCCAGAACCGCACAACACCCCAGTAAATTGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGAcgatcg                       >  1:822393/1‑150 (MQ=255)
                                                                                                                                 aaCCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGcc                                                        >  1:330809/1‑116 (MQ=255)
                                                                                                                                 aaCCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGcc                                                        <  2:330809/116‑1 (MQ=255)
                                                                                                                                 aaCCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGa                             <  2:874372/143‑1 (MQ=255)
                                                                                                                                 aaCCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGa                             >  1:874372/1‑143 (MQ=255)
                                                                                                                                    ccAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATgg                   <  2:722716/150‑1 (MQ=255)
                                                                                                                                      aGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATggcg                 >  1:971065/1‑150 (MQ=255)
                                                                                                                                         aCCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATc                                            >  1:51181/1‑120 (MQ=255)
                                                                                                                                         aCCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATc                                            <  2:51181/120‑1 (MQ=255)
                                                                                                                                          ccGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCa             >  2:821455/1‑150 (MQ=255)
                                                                                                                                            gCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAgg           <  1:652490/150‑1 (MQ=255)
                                                                                                                                               caacaCACCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGgtgt        <  1:828179/150‑1 (MQ=255)
                                                                                                                                                     cccAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGTGTGGGCCa  >  2:811679/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTGAGATTGGCAAGCCTACCCCCGGAAAATGACCGAGGGGTCATAACAAGAAGTAAGTGCCCGCCTGAGGCACACCCTGCATTTGTGTGAGGAGAAATCGCGATGACGTCGTACACCTCCGGGAATCCAACCCAGAACCGCACAACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGTGTGGGCCA  >  NC_002947/358434‑358732

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: ‑