Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,278,931 C→G 9.5% H29Q (CAC→CAG)  PP_3750 → GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,278,9310CG 9.5% 156.7 / 8.4 53H29Q (CAC→CAG) PP_3750GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (23/24);  new base G (3/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

CTGAGTCGCAGCCAATCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCGGCTGACCTACGATAACCTGCTCAGCGGCGTGGTCGGCCGTGGCACATTCGT  >  NC_002947/4278792‑4279071
                                                                                                                                           |                                                                                                                                            
cTGAGTCGCAGCCAATCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGccgccg                                                                                                                                    >  2:1186496/1‑150 (MQ=255)
 tcactCGCAGCCAATCCTCTATCCTGAGCAATCCTCCCCTCCGCACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGAAAAGACCTGGACTATCAGCTTTACATGTAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCAgg                                                                                                                                   <  2:42068/146‑1 (MQ=255)
 tGAGTCGCAGCCAATCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCgg                                                                                                                                   >  2:811578/1‑150 (MQ=255)
  gAGTCGCAGCCAATCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGc                                                                                                                                  >  1:575899/1‑150 (MQ=255)
       gCAGCCAATCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCa                                                                                                                             >  1:925717/1‑150 (MQ=255)
        cAGCCAATCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCaa                                                                                                                            >  2:975738/1‑150 (MQ=255)
         aGCCAATCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCAc                                                                                                                                              >  2:771089/1‑131 (MQ=255)
         aGCCAATCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCAc                                                                                                                                              <  1:771089/131‑1 (MQ=255)
             aaTCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTgcc                                                                                                                       >  1:555732/1‑150 (MQ=255)
                  gCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTa                                                                                                                  <  2:925717/150‑1 (MQ=255)
                        cTGAGCAATCCTCCGCTCCGGACACGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTc                                                                                                            >  2:198674/1‑150 (MQ=255)
                         tGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCAGGTGgccgctgt                                                                                                                                   >  1:644322/1‑123 (MQ=255)
                         tGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCAGGTGgccgctgt                                                                                                                                   <  2:644322/126‑4 (MQ=255)
                               aTCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTgctggc                                                                                                     <  2:575899/150‑1 (MQ=255)
                                tcctccGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCAGGTGgccgctgtcga                                                                                                                                <  1:282575/122‑7 (MQ=255)
                                tcctccGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCAGGTGgccgctgtcga                                                                                                                                >  2:282575/1‑116 (MQ=255)
                                tcctccGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTgctggct                                                                                                    >  1:419828/1‑150 (MQ=255)
                                tcctccGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTgctggct                                                                                                    >  2:276041/1‑150 (MQ=255)
                                tcctccGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTgctgg                                                                                                      >  2:266338/1‑148 (MQ=255)
                                tcctccGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTgctgg                                                                                                      <  1:266338/148‑1 (MQ=255)
                                   tccgctccgGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGTGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGaa                                                                                                 >  1:894432/1‑150 (MQ=255)
                                   tccgctccgGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGaa                                                                                                 <  1:276041/150‑1 (MQ=255)
                                      gctccgGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGaacaa                                                                                              >  2:474903/1‑150 (MQ=255)
                                          cgGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTgg                                                                                          <  1:905058/150‑1 (MQ=255)
                                          cgGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTgg                                                                                          <  2:1011101/150‑1 (MQ=255)
                                           ggACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCCCGGGCTCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTggg                                                                                         <  1:811578/150‑1 (MQ=255)
                                                 gCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTctgctg                                                                                                       <  2:976198/130‑1 (MQ=255)
                                                 gCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTctgctg                                                                                                       >  1:976198/1‑130 (MQ=255)
                                                 gCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATc                                                                                   >  2:151288/1‑150 (MQ=255)
                                                  cccATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATcg                                                                                  <  2:555732/150‑1 (MQ=255)
                                                     aTGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCATACCCGTCTGCTGGCTGAACAACTGGGCGTATCgcgc                                                                               <  1:975738/150‑1 (MQ=255)
                                                       ggAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTgg                                                                                                                                          <  1:977709/89‑1 (MQ=255)
                                                       ggAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTgg                                                                                                                                          >  2:977709/1‑89 (MQ=255)
                                                         aaCTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCCCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGa                                                                           <  1:752849/150‑1 (MQ=255)
                                                            tGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCg                                                                        <  1:629514/150‑1 (MQ=255)
                                                            tGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCg                                                                        >  1:1158225/1‑150 (MQ=255)
                                                             gCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGaac                                                                                                >  1:1159614/1‑125 (MQ=255)
                                                             gCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGaac                                                                                                <  2:1159614/125‑1 (MQ=255)
                                                               acaTTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGTGCGATGGTATCGACAGCGGCCAACTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTgg                                                                     >  1:811198/1‑150 (MQ=255)
                                                                 aTTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGcc                                                                   >  1:939837/1‑150 (MQ=255)
                                                                    cGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGTGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAg                                                                <  2:811198/150‑1 (MQ=255)
                                                                                  cAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTgg                                                                                          >  1:488186/1‑110 (MQ=255)
                                                                                  cAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTgg                                                                                          <  2:488186/110‑1 (MQ=255)
                                                                                    aaGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCGGCTGa                                                <  2:1158225/150‑1 (MQ=255)
                                                                                    aaGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCGGCTGa                                                <  2:229743/150‑1 (MQ=255)
                                                                                              cGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGAGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCGGCTGACCTACGATaa                                      <  1:198674/150‑1 (MQ=255)
                                                                                               gATCAGCTTTACCAGCAATTGCTCGATGGTCTCGAACGCGGCCACCTGGCGGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTCAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCGGCTGACCTACGATGAc                                     <  2:83823/150‑1 (MQ=255)
                                                                                                aTCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCGGCTGACCTACGATAAcc                                    <  2:939837/150‑1 (MQ=255)
                                                                                                 tCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTgg                                                                                          >  1:1099720/1‑95 (MQ=255)
                                                                                                 tCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTgg                                                                                          <  2:1099720/95‑1 (MQ=255)
                                                                                                     cTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCGGCTGACCTACGATAACCTGCTc                               >  2:373987/1‑150 (MQ=255)
                                                                                                        tACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCGGCTGACCTACGATAACCTGCTCAgc                            >  1:212209/1‑150 (MQ=255)
                                                                                                                              tCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCGGCTGACCTACGATAACCTGCTCAGCGGCGTGGTCGGCCGTGGCACAt      >  1:150757/1‑150 (MQ=255)
                                                                                                                               agccggCGGCCAGGTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCGGCTGACCTACGATAACCTGCTCAGCGGCGTGGTCGGCCGTGGCACAtt     >  1:82898/6‑150 (MQ=255)
                                                                                                                                  cAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCGGCTGACCTACGATAACCTGCTCAGCGGCGTGGTCGGCCGTGGCACAGTCGt  <  1:1186496/150‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                            
CTGAGTCGCAGCCAATCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCGGCTGACCTACGATAACCTGCTCAGCGGCGTGGTCGGCCGTGGCACATTCGT  >  NC_002947/4278792‑4279071

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